Subsurface Thermophysical Properties of Europa’s Leading and Trailing Hemispheres as Revealed by ALMA

Abstract

We present best-fit values of porosity – and the corresponding effective thermal inertiae – determined from three different depths in Europa’s near-subsurface ($\sim 1-20$ cm). The porosity of the upper $\sim 20$ cm of Europa’s subsurface varies between 75–50$\%$ ($\Gamma_{eff}\approx 50-140$ J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$) on the leading hemisphere and 50–40$\%$ ($\Gamma_{eff}\approx 140-180$ J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$) on the trailing hemisphere. Residual maps produced by comparison with these models reveal thermally anomalous features that cannot be reproduced by globally homogeneous porosity models. These regions are compared to Europa’s surface terrain and known compositional variations. We find that some instances of warm thermal anomalies are co-located with known geographical or compositional features on both the leading and trailing hemisphere; cool temperature anomalies are well correlated with surfaces previously observed to contain pure, crystalline water ice and the expansive rays of Pwyll crater. Anomalous regions correspond to locations with subsurface properties different from those of our best-fit models, such as potentially elevated thermal inertia, decreased emissivity, or more porous regolith. We also find that ALMA observations at $\sim 3$ mm sound below the thermal skin depth of Europa ($\sim 10-15$ cm) for a range of porosity values, and thus do not exhibit features indicative of diurnal variability or residuals similar to other frequency bands. Future observations of Europa at higher angular resolution may reveal additional locations of variable subsurface thermophysical properties, while those at other wavelengths will inform our understanding of the regolith compaction length and the effects of external processes on the shallow subsurface.

1 Introduction

The surface of Europa, the smallest of the Galilean Satellites, is notable for its varied terrain units and hemispheric asymmetries – particularly when compared to its sister moons that are resurfaced by active volcanoes (Io), heavily cratered and relatively dark (Callisto), or somewhat intermediary (Ganymede) (see the reviews in McEwen et al. 2004, Greeley et al. 2004, Moore et al. 2004, Pappalardo et al. 2004, de Pater et al. 2021b and references therein). The presence of ridged plains, chaotic terrain (comprised of small, incoherent ice latticework), hydrated salts and sulfuric compounds across Europa’s icy surface indicates the crust above its subsurface ocean may be relatively young and tectonically active, while also being exogenically weathered (Smith et al., 1979; McCord et al., 1998; Zahnle et al., 1998; Cooper et al., 2001; Paranicas et al., 2001; Zahnle et al., 2003; Schenk & Pappalardo, 2004; Bierhaus et al., 2009; Doggett et al., 2009; Carlson et al., 2009). Tidally locked and orbiting slower than Jupiter’s magnetic field, which is tied to the planet’s rapid rotation ($\sim 10$ hr), the trailing hemisphere (centered at 270${}^{\circ}$ W) is constantly bombarded by charged particles and heavy ions (e.g. S${}^{+}$, O${}^{+}$) sourced from Io’s plasma torus that are entrained in Jupiter’s magnetic field (Paranicas et al., 2009). The leading hemisphere (centered at 90${}^{\circ}$ W) is exposed to the highest energy particles from the Jovian magnetosphere (Nordheim et al., 2022) and micrometeorite gardening (Zahnle et al., 1998), and presents a brighter, less altered surface. However, the influence of exogenic erosion of the upper layers of Europa’s regolith and the slow resurfacing from beneath are not easily related to latitude or longitude, resulting in a complex surface whose composition and structure are influenced from both the subsurface ocean and exogenic environment (Anderson et al., 1998; Carr et al., 1998; Pappalardo et al., 1999; Kivelson et al., 2000). Hydrated minerals and salts have been detected across the varied surface terrain, possibly originating in the subsurface while radiolysis provides the formation of sulfur-bearing species, hydrogen peroxide, carbon dioxide, among others (McCord et al., 1998; Carlson et al., 2005, 2009; Brown & Hand, 2013; Trumbo et al., 2019a, b, 2022; Trumbo & Brown, 2023; Villanueva et al., 2023).

As far back as the early 20${}^{th}$ century questions regarding the hemispheric dichotomy of Europa’s surface brightness and properties began to arise, initially from ground-based observations (Stebbins, 1927; Stebbins & Jacobsen, 1928). These questions persisted into the 1970’s based on further ground-based observations and data from the Pioneer flybys of Jupiter (Fimmel et al., 1974). Subsequently, the coloration, non-icy material composition, mixtures of amorphous or crystalline ice, and weathering by Jovian magnetospheric ions have been investigated in-depth using near-infrared through ultraviolet wavelength instruments onboard the spacecraft venturing near and into the Jovian system. Initial observations from the Voyager spacecraft determined differences in color, albedo, and water ice distributions between the leading and trailing hemispheres (Pilcher et al., 1972; Lucchitta & Soderblom, 1982; McEwen, 1986; Spencer, 1987), while subsequent Galileo measurements revealed compositional and thermal variations using the Ultraviolet Spectrometer (UVS; Hendrix et al., 1998), Near Infrared Mass Spectrometer (NIMS; Carlson et al., 1996; McCord et al., 1998; Hansen & McCord, 2004), Photopolarimeter-radiometer (PPR; Spencer et al., 1999; Rathbun et al., 2010; Rathbun & Spencer, 2020), and Solid-state Imaging (SSI; Fanale et al., 2000; Leonard et al., 2018) instruments. Flybys with the Cassini and New Horizons spacecraft allowed for further study of the water and hydrated mineral composition of the surface ice via the Visual and Infrared Mapping Spectrometer (VIMS), Linear Etalon Imaging Spectral Array (LEISA), and LOng-Range Reconnaissance Imager (LORRI) instruments (Brown et al., 2003; McCord et al., 2004; Grundy et al., 2007). Recently, high spatial resolution observations of the surface with the Jovian InfraRed Auroral Mapper (JIRAM) onboard the Juno spacecraft allowed for constraints on the ice grain size, while in situ magnetometer measurements helped to better characterize the charged particle environment at Europa (Filacchione et al., 2019; Mishra et al., 2021; Addison et al., 2023).

Observations of the Galilean Satellites in support of these missions (and in-between) have been conducted with ground- and space-based assets, improving our understanding of the distinct coloration, albedo differences, and hydrate absorption features across the surface. These include compositional and thermal studies utilizing the International Ultraviolet Explorer (IUE; Lane et al. 1981; Domingue & Lane 1998), Hubble Space Telescope (HST; Noll et al. 1995; Brown & Hand 2013; Trumbo et al. 2020, 2022), the airborne SOFIA observatory (de Pater et al., 2021a), and various ground-based facilities such as the Very Large Telescope (VLT), Infrared Telescope Facility (IRTF), and Keck (Hansen, 1973; de Pater et al., 1989; Spencer & Calvin, 2002; Spencer et al., 2006; Fischer et al., 2015; Ligier et al., 2016; Fischer et al., 2017; Trumbo et al., 2017b; King et al., 2022). Recent results from the James Webb Space Telescope (JWST) also show evidence for the endogenous origin of previously detected surface CO${}_{2}$, potentially sourced from the subsurface ocean (Trumbo & Brown, 2023; Villanueva et al., 2023). Though the variability in surface terrain and reddened trailing hemisphere have now been well characterized, the endo- and exogenic processes that have influenced Europa’s surface composition and evolution are currently poorly understood, and will likely remain so until the arrival of the JUpier ICy moons Explorer (JUICE) and Europa Clipper spacecraft in the future (which will undoubtedly provide many additional questions of their own).

Complementary to the aforementioned observations at shorter wavelengths are those in the radio and (sub)millimeter regime, which probe the near-surface crust down to $\sim 10$s of cm to m depths; beyond, the deeper layers of the crust may be probed by microwave and radar observations down to $\sim 10$ km (Ostro, 1982; Ostro et al., 1992; Chyba et al., 1998; Bruzzone et al., 2013), including recent in situ remote sensing with the Juno MicroWave Radiometer (MWR; Janssen et al., 2017; Zhang et al., 2023), and future thermal imaging, submillimeter, and radar observations from the Europa Clipper and JUICE spacecraft (Hartogh et al., 2013; Phillips & Pappalardo, 2014; Pappalardo et al., 2017). Millimeter wave observations at different wavelengths permit the measurement of thermal radiation as a function of subsurface depth, the modification of which is governed by the thermophysical properties of the surface. These include the millimeter emissivity, subsurface thermal inertia, porosity, dust fraction, and grain size, all of which inform our understanding of how the various endo- and exogenic processes have altered the surface, and to what extent they change the subsurface structure and composition. Initial characterization of the subsurface properties and thermal emission of the Galilean Satellites were made with a number of long-wavelength facilities throughout the last half-century, including single dish facilities such as the 2.24-m telescope on Maunakea, the 12-m dish at Kitt Peak, the Effelsberg 100-m telescope, and the Institut de Radioastronomie Millimétrique (IRAM) 30-m telescope (Morrison et al., 1972; Morrison & Cruikshank, 1973; Ulich & Conklin, 1976; Pauliny-Toth et al., 1977; Ulich et al., 1984; Altenhoff et al., 1988); the Owens Valley Radio Observatory (OVRO) 3-element array (Berge & Muhleman, 1975; Muhleman & Berge, 1991); the the SubMillimeter Array (SMA) and Very Large Array (VLA) interferometers (de Pater et al., 1982, 1984; Muhleman et al., 1986); and the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space-based telescope (Müller et al., 2016). Often, the Galilean Satellites were also used for flux calibration observations for (sub)millimeter facilities, along with Saturn’s largest moon, Titan (Ulich, 1981; Moreno, 2007; Butler, 2012).

Radio and (sub)millimeter interferometric observations from modern telescopes can spatially resolve small Solar System bodies, such as Europa, and thus enable the measurement of thermophysical properties as a function of location on the body, by modeling the thermal radiation from the subsurface (cm-m depths). Utilizing the Atacama Large Millimeter/submillimeter Array (ALMA), Trumbo et al. (2018) mapped the thermal inertia of Europa’s surface using 1.3 mm (233 GHz) observations; they also investigated the correlation of thermal anomalies observed with ALMA with potential plume locations (Trumbo et al., 2017a). These studies revealed that a global thermal inertia of 95 J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$ and emissivity of 0.75 provided good fits to the ALMA observations. They found that anomalously cold locations in the ALMA observations, such as around Pwyll crater ($\sim 271^{\circ}$W, $25^{\circ}$S) and a region on the leading hemisphere ($90^{\circ}$W, 23${}^{\circ}$N), were indicative of localized, high thermal inertia regions or low emissivity; thermal inertia values ranging from 40–300 J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$ or emissivities from 0.67–0.84 were found to characterize outlying regions in the residual maps, though thermal anomalies were not correlated with geological or morphological features (excepting Pwyll). The retrieved thermal inertias are comparable to those derived for the surface from Galileo/PPR observations, which provided constraints on Europa’s thermal inertia from 40–150 J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$, with elevated measurements in similarly anomalous regions such as near Pwyll (Spencer et al., 1999; Rathbun et al., 2010; Rathbun & Spencer, 2020). Recent analyses have also been carried out for Ganymede (de Kleer et al., 2021a) and Callisto (Camarca et al., 2023) using ALMA to investigate the change in porosity or thermal inertia as a function of depth and correlate brightness temperatures to geographically distinct surface regions. On Ganymede, de Kleer et al. (2021a) found that a porosity gradient between 10–40$\%$ provided good fits to ALMA observations sounding the upper $\sim$0.5 m of the subsurface. From ALMA Band 7 data, Camarca et al. (2023) derived a mixture of high (1200–2000 J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$) and low (15–50 J m${}^{-2}$ K${}^{-1}$ s${}^{-1/2}$) thermal inertia components to correctly model the thermal emission from Callisto’s leading hemisphere. Both studies found cold thermal anomalies co-located with the locations of crater basins or complexes. Generally, these studies revealed higher thermal inertias on the near subsurface of Ganymede and Callisto than Europa.

Here, we present the analysis of ALMA observations of Europa at three wavelengths (0.88, 1.25, and 3.05 mm) that probe distinct depths in Europa’s subsurface, which allows us to investigate the change in thermophysical properties with depth and latitude, and ascertain their potential link to exogenic sources and the evolution of Europa’s ice shell. These observations complement the recent studies of Ganymede and Callisto with ALMA, and provide context for Juno observations of Europa with infrared and microwave instruments. In Section 2, we detail the ALMA observations, reduction and imaging procedures, followed by the radiative transfer modeling in Section 3. A discussion of the modeling results is presented in Sections 4, followed by our conclusions in Section 5.

2 Observations

The ALMA Main Array is an interferometer consisting of up to 50 12-m antennas located in the Atacama Desert, Chile. Every pair of antennas acts as a two-element interferometer, measuring a single complex component (often called a “visibility”) of the Fourier transform of the sky brightness. Together, the collection of visibilities allows for the reconstruction of the full sky brightness in both dimensions via image deconvolution techniques (see Thompson et al., 2001, and references therein). As part of ALMA Project Code 2016.1.00691.S, the leading and trailing hemispheres of each of the Galilean Satellites were observed in three distinct frequency bands that probe different subsurface depths: ALMA Band 3 (97.5 GHz; 3.05 mm), Band 6 (233 GHz; 1.25 mm), and Band 7 (343.5 GHz, 0.88 mm). Europa was observed 8 times between 2016 and 2017. As the angular resolution of interferometric observations depends on the distances between antennas in the array, these observations were executed using different antenna configurations so as to obtain relatively consistent resolution across all three frequency bands. A configuration with maximum antenna separation of 6.4 km was used for Band 3 observations to achieve comparable resolution to data from higher frequency bands, while a configuration with a shorter maximum antenna separation of 1.3 km was used for Bands 6 and 7. Separate observations in each frequency band were executed to target both the leading and trailing hemispheres of Europa, with typical integration times of $\sim 120-300$ s; as such, longitudinal smearing over this time period was well below the size of a resolution element. All observations were carried out using between 40 and 45 antennas. In some cases, multiple executions (i.e. observing integrations) were acquired for each hemisphere in a single band, allowing for additional longitude coverage and higher constraints on thermophysical properties. An additional execution in ALMA Band 7 was performed, but was set to incorrect sky coordinates, and as such was not analyzed here. The observation parameters for each integration are detailed in Table 2.

Table 1: Observational Parameters

Obs. Date	Freq.${}^{a}$	$\lambda$	Tag$\#^{b}$	Ang. Diam.	Spatial Res.	Pos. Ang.${}^{c}$	Lat.	W Lon.	Corr.${}^{d}$
(UTC)	(GHz)	(mm)		(arcsec)	(arcsec)	(${}^{\circ}$)	(${}^{\circ}$N)	(${}^{\circ}$W)	Factor
2017 Sep 19 17:15	97.5	3.05	3L	0.684	$0.110\times 0.086$	52.23	-3.08	103.3	1.0
2017 Sep 28 14:43${}^{\dagger}$			3T	0.677	$0.148\times 0.085$	54.87	-3.12	283.2	1.0
2017 Aug 07 19:15	233	1.25	6L	0.737	$0.107\times 0.082$	87.24	-2.92	80.63	0.985
2017 Jul 09 01:03			6T0${}^{e}$	0.797	$0.280\times 0.107$	-73.13	-2.89	307.7	1.068
2017 Jul 30 00:05			6T1	0.754	$0.163\times 0.078$	-70.25	-2.91	270.1	1.028
2017 Aug 16 22:55			6T2	0.722	$0.127\times 0.079$	-70.64	-2.95	287.0	0.937
2017 Jul 06 23:56${}^{\dagger}$	343.5	0.88	7L	0.803	$0.139\times 0.068$	-68.88	-2.90	100.3	0.956
2016 Oct 25 12:28			7T${}^{e}$	0.678	$0.196\times 0.151$	66.85	-2.50	226.7	0.938

Note. — ${}^{a}$Averaged frequency of all continuum windows. Frequencies correspond to ALMA Band 3 (97.5 GHz), Band 6 (233 GHz), and Band 7 (343.5 GHz). ${}^{b}$Tag denoting the ALMA frequency band, targeted hemisphere (L = leading, T = trailing), and observation number; exact longitudes vary slightly for each execution. Hemispheres with multiple integrations are denoted with separate labels for each individual execution. ${}^{c}$The position angle of the synthesized ALMA beam, denoted in degrees counter-clockwise from the positive vertical. ${}^{d}$Correction factor derived from variability of quasars used for flux density calibrations. ${}^{e}$6T0: Though data from this execution were reduced and modeled, the beam dimensions prevents the data from yielding meaningful longitudinal information regarding Europa’s surface properties. 7T: A second execution for the trailing hemisphere in Band 7 was not used. ${}^{{\dagger}}$Denotes observations where interloping satellites were present in the ALMA field - see Appendix A.

Data from each integration were reduced using the Common Astronomy Software Applications (CASA) package ver. 4.7 (Jaeger, 2008) and the provided ALMA pipeline scripts. Continuum images were produced by flagging channels with telluric contamination and then averaging to channel bins of 125 or 256 MHz to reduce data volume. The resulting data were then combined using multi-frequency synthesis imaging methods to produce a single, high signal-to-noise ratio (SNR) broadband image of the thermal continuum emission. Phase self-calibration was performed on each observation to compensate for tropospheric phase fluctuations, which improves image coherence and SNR for each observation (see the discussion in Cornwell & Fomalont, 1999; Butler & Bastian, 1999; Brogan et al., 2018, and ALMA Memo 620 by Richards et al.). Similar procedures were applied to the accompanying observations of Ganymede and Callisto (de Kleer et al., 2021a; Camarca et al., 2023).

Final image deconvolution, which removed interferometric artifacts induced by the lack of complete antenna coverage on the sky, was performed using the CASA tclean task with image sizes of $1000\times 1000$ pixels of $0.01^{\prime\prime}$ size (note that this is not the effective resolution, which is shown in Table 2, but simply the pixel size). Briggs weighting was applied with a “robust” factor of 0, which slightly increases the weight of data from larger antenna separations (Briggs, 1995). The removal of interferometric artifacts – and thus the improvement of the final image quality – for two of the ALMA observations was facilitated by accounting for the emission from nearby Galilean Satellites (Ganymede, Callisto) that intervened on the relatively large ALMA Field-of-View (FOV), introducing additional signal in the sidelobes. These procedures and the improvements in the images are detailed in Appendix A.

For each ALMA integration, the disk-averaged flux density of Europa was determined by fitting a disk model to the calibrated visibility data, often excluding data from larger antenna spacings (e.g. >100–200 m), which are sensitive to smaller scale thermal structure (such as surface variations) and not the total flux density. A correction to this value was made based on the variability of measured quasar brightnesses for each quasar used for each ALMA observation, as was done for previous ALMA observations (Trumbo et al., 2018; de Kleer et al., 2021a). The flux density for each quasar was interpolated based on the measurements from the nearest dates in the cases of Band 3 and 7 observations, where quasars were commonly monitored. For Band 6 observations, quasar flux density curves were derived based on the functional form detailed in Ennis et al. (1982), using contemporaneous quasar observations in both Bands 3 and 7 to determine the variability of flux density with frequency. No corrections were needed for Band 3 data because the quasar flux densities were determined on the same date as the observations. For Bands 6 and 7, we found correction factors from $1.5-6.8\%$ were needed (Table 2). As found previously, the dependence of ALMA on quasar observations can result in higher flux density scale calibration uncertainties (Francis et al., 2020); as a result, our quoted uncertainties on the disk-averaged flux densities, temperatures, and emissivities are no less than $5\%$, which are often larger than the statistical uncertainties derived from the model fit for the flux density.

After converting from flux density units (Jy) to brightness temperature (K, the expected thermal temperature the surface would emit if it was solely parameterized by the Planck function; see also de Kleer et al. 2021a and Camarca et al. 2023), the final emission maps were compared to radiative transfer models generated using a variety of thermophysical properties and global porosity or thermal inertia conditions. The Europa continuum image maps are shown in Figure 1. The measured flux densities and brightness temperatures are listed in Table 2.

Figure 1: Brightness temperature maps of Europa’s leading hemisphere ($\sim 90^{\circ}$ W longitude; top row) and trailing hemisphere ($\sim 270^{\circ}$ W longitude; bottom row) from Band 3 (3.05 mm; A, E), Band 6 (1.25 mm; B, F-H), and Band 7 (0.88 mm; C, G). The ALMA beam (the full-width at half-maximum of the ALMA point spread function) is shown as a hashed ellipse in the bottom left corner. All images are aligned with Europa’s north pole along the vertical axis.

Table 2: Derived Properties and Results

Note. — ${}^{a}$Flux densities and brightness temperatures listed here are derived as hemispheric averages. ${}^{*}$Porosity and effective thermal inertia values derived for ALMA Band 3 observations were inferred through bounds as discussed in Appendix B.