The near-surface electron radiation environment of Saturn's moon Mimas

Saturn's inner mid-size moons are exposed to a number of external weathering processes, including charged particle bombardment and UV photolysis, as well as deposition of E-ring grains and interplanetary dust. While optical remote sensing observations by several instruments onboard the Cassini spacecraft have revealed a number of weathering patterns across the surfaces of these moons, it is not entirely clear which external process is responsible for which observed weathering pattern. Here we focus on Saturn's moon Mimas and model the effect of energetic electron bombardment across its surface. By using a combination of a guiding center, bounce-averaged charged particle tracing approach and a particle physics code, we investigate how the radiation dose due to energetic electrons is deposited with depth at different locations. We predict a lens-shaped electron energy deposition pattern that extends down to ∼cm depths at low latitudes centered around the apex of the leading hemisphere (90° W). These results are consistent with previous remote sensing observations of a lens-shaped color anomaly observed by the Imaging Science Subsystem (ISS) instrument as well as a thermal inertia anomaly observed by the Visual and Infrared Mapping Spectrometer (VIMS) and the Composite Infrared Spectrometer (CIRS). Our results confirm that these features are produced by MeV electrons that have a penetration depth into the surface comparable to the effective sampling depths of these instruments. On the trailing hemisphere we predict a similar lens-shaped electron energy deposition pattern, whose effects have to date not been observed by the Cassini remote sensing instruments. We suggest that no corresponding lens-shaped weathering pattern has been observed on the trailing hemisphere because of the comparatively short range of lower energy (<1 MeV) electrons into surface ice, as well as competing effects from cold plasma, neutral, and dust bombardment

Constraining the surface properties of Helene

We analyze two sets of observations of Dione's co-orbital satellite Helene taken by Cassini's Composite Infrared Spectrometer (CIRS). The first observation was a CIRS FP3 (600 to 1100 cm−1, 9.1 to 16.7 μm) stare of Helene's trailing hemisphere, where two of the ten FP3 pixels were filled. The daytime surface temperatures derived from these observations were 83.3 ± 0.9 K and 88.8 ± 0.8 K at local times 223° to 288° and 180° to 238° respectively. When these temperatures were compared to a 1-D thermophysical model only albedos between 0.25 and 0.70 were able to fit the data, with a mean and standard deviation of 0.43 ± 0.12. All thermal inertias tested between 1 and 2000 J m−2 K−1 s-1/2 could fit the data (i.e. thermal inertia was not constrained). The second observation analyzed was a FP3 and FP4 (1100 to 1400 cm−1, 7.1 to 9.1 μm) scan of Helene's leading hemisphere. Temperatures between 77 and 89 K were observed with FP3, with a typical error between 5 and 10 K. The surface temperatures derived from FP4 were higher, between 98 and 106 K, but with much larger errors (between 10 and 30 K) and thus the FP3- and FP4-derived temperature largely agree within their uncertainty. Dione's disk-integrated bolometric Bond albedos have been found to be between 0.63 ± 0.15 (Howett et al. 2010) and 0.44 ± 0.13 (Howett et al. 2014). Thus Helene may be darker than Dione, which is the opposite of the trend found at shorter wavelengths (c.f. Hedman et al. 2020; Royer et al., 2021). However few conclusions can be drawn since the albedos of Dione and Helene agree within their uncertainty

Bolometric bond albedo and thermal inertia maps of Mimas

In 2011 a thermally anomalous region was discovered on Mimas, Saturn's innermost major icy satellite (Howett et al., 2011). The anomalous region is a lens-like shape located at low latitudes on Mimas' leading hemisphere. It manifests as a region with warmer nighttime temperatures, and cooler daytime ones than its surroundings. The thermally anomalous region is spatially correlated with a darkening in Mimas' IR/UV surface color (Schenk et al., 2011) and the region preferentially bombarded by high-energy electrons (Paranicas et al., 2012, Paranicas et al., 2014; Nordheim et al., 2017). We use data from Cassini's Composite Infrared Spectrometer (CIRS) to map Mimas' surface temperatures and its thermophysical properties. This provides a dramatic improvement on the work in Howett et al. (2011), where the values were determined at only two regions on Mimas (one inside, and another outside of the anomalous region). We use all spatially-resolved scans made by CIRS' focal plane 3 (FP3, 600 to 1100 cm−1) of Mimas' surface, which are largely daytime observations but do include one nighttime one. The resulting temperature maps confirm the presence and location of Mimas' previously discovered thermally anomalous region. No other thermally anomalous regions were discovered, although we note that the surface coverage is incomplete on Mimas' leading and anti-Saturn hemisphere. The thermal inertia map confirms that the anomalous region has a notably higher thermal inertia than its surroundings: 98 ± 42 J m−2 K−1 s-1/2 inside of the anomaly, compared to 34 ± 32 J m−2 K−1 s-1/2 outside. The albedo inside and outside of the anomalous region agrees within their uncertainty: 0.45 ± 0.08 inside compared to 0.41 ± 0.07 outside the anomaly. Interestingly the albedo appears brighter inside the anomaly region, which may not be surprising given this region does appear brighter at some UV wavelengths (0.338 μm, see Schenk et al., 2011). However, this result should be treated with caution because, as previously stated, statistically the albedo of these two regions is the same when their uncertainties are considered. These thermal inertia and albedo values determined here are consistent with those found by Howett et al. (2011), who determined the thermal inertia inside the anomaly to be 66 ± 23 J m−2 K−1 s-1/2 and <16 J m−2 K−1 s-1/2 outside, with albedos that varied from 0.49 to 0.70

k-means aperture optimization applied to Kepler K2 time series photometry of Titan

Motivated by the Kepler K2 time series of Titan, we present an aperture optimization technique for extracting photometry of saturated moving targets with high temporally and spatially varying backgrounds. Our approach uses k-means clustering to identify interleaved families of images with similar point-spread function and saturation properties, optimizes apertures for each family independently, then merges the time series through a normalization procedure. By applying k-means aperture optimization to the K2 Titan data, we achieve ≤0.33% photometric scatter in spite of background levels varying from 15% to 60% of the target's flux. We find no compelling evidence for signals attributable to atmospheric variation on the timescales sampled by these observations. We explore other potential applications of the k-means aperture optimization technique, including testing its performance on a saturated K2 eclipsing binary star. We conclude with a discussion of the potential for future continuous high-precision photometry campaigns for revealing the dynamical properties of Titan's atmosphere

Optical constants of ammonium hydrosulfide ice and ammonia ice

Thin-film transmission spectra of ammonium hydrosulfide (NH4SH) ice and ammonia (NH3) ice between 1300 and 12,000 cm-1 were used to determine the ice's optical constants. The films were grown on a sapphire substrate, and a Fourier-transform spectrometer and a grating spectrometer were used together to record the spectra. Lambert's law was used to directly determine the imaginary component of the complex refractive indices; from this, the real component was derived using the Kramers-Kronig algorithm. It is shown that, contrary to what is expected, the optical constants determined for NH3 ice at 80 K are in good agreement with those in the cubic phase, rather than the metastable one. The phase of the NH4SH ice was observed to change from amorphous to polycrystalline as the film was annealed to 160 K. © 2006 Optical Society of America

Maps of Tethys’ thermophysical properties

On 11th April 2015 Cassini's Composite Infrared Spectrometer (CIRS) made a series of observations of Tethys&#x2019; daytime anti-Saturn hemisphere over a nine-hour time period. During this time the sub-spacecraft position was remarkably stable (0.3&#xB0; S to 3.9&#xB0; S; 153.2&#xB0; W to 221.8&#xB0; W), and so these observations provide unprecedented coverage of diurnal temperature variations on Tethys&#x2019; anti-Saturn hemisphere. In 2012 a thermal anomaly was discovered at low latitudes on Tethys&#x2019; leading hemisphere; it appears cooler during the day and warmer at night than its surroundings (Howett et al., 2012) and is spatially correlated with a decrease in the IR3/UV3 visible color ratio (Schenk et al., 2011). The cause of this anomaly is believed to be surface alteration by high-energy electrons, which preferentially bombard low-latitudes of Tethys&#x2019; leading hemisphere (Schenk et al., 2011; Howett et al., 2012; Paranicas et al. 2014; Schaible et al., 2017). The thermal anomaly was quickly dubbed &#x201C;Pac-Man&#x201D; due to its resemblance to the 1980s video game icon. We use these daytime 2015 CIRS data, along with two sets of nighttime CIRS observations of Tethys (from 27 June 2007 and 17 August 2015) to make maps of bolometric Bond albedo and thermal inertia variations across the anti-Saturn hemisphere of Tethys (including the edge of its Pac-Man region). These maps confirm the presence of the Pac-Man thermal anomaly and show that while Tethys&#x2019; bolometric Bond albedo varies negligibly outside and inside the anomaly (0.69 &#xB1; 0.02 inside, compared to 0.71 &#xB1; 0.04 outside) the thermal inertia varies dramatically (29 &#xB1; 10 J m&#x2212;2 K&#x2212;1 s&#x2212;1/2 inside, compared to 9 &#xB1; 4 J m&#x2212;2 K&#x2212;1 s&#x2212;1/2 outside). These thermal inertias are in keeping with previously published values: 25 &#xB1; 3 J m&#x2212;2 K&#x2212;1 s&#x2212;1/2 inside, and 5 &#xB1; 1 J m&#x2212;2 K&#x2212;1 s&#x2212;1/2 outside the anomaly (Howett et al., 2012). A detailed analysis shows that on smaller spatial-scales the bolometric Bond albedo does vary: increasing to a peak value at 180&#xB0; W. For longitudes between &#x223C;100&#xB0; W and &#x223C;160&#xB0; W the thermal inertia increases from northern to southern latitudes, while the reverse is true for bolometric Bond albedo. The thermal inertia on Tethys generally increases towards the center of its leading hemisphere but also displays other notable small-scale variations. These thermal inertia and bolometric Bond albedo variations are perhaps due to differences in competing surface modification by E ring grains and high-energy electrons which both bombard Tethys&#x2019; leading hemisphere (but in different ways). A comparison between the observed temperatures and our best thermal model fits shows notable discrepancies in the morning warming curve, which may provide evidence of regional variations in surface roughness effects, perhaps again due to variations in surface alteration mechanisms

Limits on Dione's activity using Cassini/CIRS data

We use nighttime Cassini Composite Infrared Spectrometer (CIRS) data to look for discrete regions of elevated nighttime temperatures indicative of endogenic activity on Dione's surface. This is achieved by producing low latitude and midlatitude (less than 60°) maps of Dione's nighttime surface temperature, derived from 10 to 1,100-cm−1 CIRS data. The surface temperatures observed do not show evidence of any small discrete regions of elevated nighttime temperatures and are comparable to temperatures predicted by a passive thermophysical model of Dione's surface. Thus, we conclude that no evidence for activity exists on Dione at midlatitude to low latitude. Using the derived surface temperature maps, we set upper limits for the temperature at which a 50-, 100-, or 200-km2 hot spot would remain undetected by this study. We find the mean temperature of such a hot spot would be 117.1 ± 47.2 K (−249 F), 104.8 ± 27.7 K (−272 F), and 95.4 ± 19.5 K (−288 F) for a 50-, 100-, and 200-km2 hot spot, respectively, corresponding to endogenic emission of 1.07, 0.68, and 0.47 GW

Enceladus plume structure and time variability: comparison of Cassini observations

During three low-altitude (99, 66, 66 km) flybys through the Enceladus plume in 2010 and 2011, Cassini's ion neutral mass spectrometer (INMS) made its first high spatial resolution measurements of the plume's gas density and distribution, detecting in situ the individual gas jets within the broad plume. Since those flybys, more detailed Imaging Science Subsystem (ISS) imaging observations of the plume's icy component have been reported, which constrain the locations and orientations of the numerous gas/grain jets. In the present study, we used these ISS imaging results, together with ultraviolet imaging spectrograph stellar and solar occultation measurements and modeling of the three-dimensional structure of the vapor cloud, to constrain the magnitudes, velocities, and time variability of the plume gas sources from the INMS data. Our results confirm a mixture of both low and high Mach gas emission from Enceladus' surface tiger stripes, with gas accelerated as fast as Mach 10 before escaping the surface. The vapor source fluxes and jet intensities/densities vary dramatically and stochastically, up to a factor 10, both spatially along the tiger stripes and over time between flyby observations. This complex spatial variability and dynamics may result from time-variable tidal stress fields interacting with subsurface fissure geometry and tortuosity beyond detectability, including changing gas pathways to the surface, and fluid flow and boiling in response evolving lithostatic stress conditions. The total plume gas source has 30% uncertainty depending on the contributions assumed for adiabatic and nonadiabatic gas expansion/acceleration to the high Mach emission. The overall vapor plume source rate exhibits stochastic time variability up to a factor ∼5 between observations, reflecting that found in the individual gas sources/jets

Extreme exospheric dynamics at Charon: Implications for the red spot

Charon's exosphere may exhibit extreme seasonal dynamics, with centuries of quiescence punctuated by short lived (∼4 earth years) exospheric surges near the equinoxes, as spring sunrise bi-annually drives frozen methane off the polar night zones. Charon's pole-centric red spot has been proposed to be the product of Ly-α photolysis of frozen methane into refractory hydrocarbon “tholins”, but the role of exospheric dynamics in the red material's formation has not been investigated. We show with exospheric modeling that methane “polar-swap”, in which exospheric CH4 sublimated from the spring polar zone is rapidly re-frozen onto the autumn hemisphere, deposits ∼30 μm polar frosts too thick for Ly-α light to penetrate. Ethane, the primary methane photoproduct under these conditions, may unlike methane remain frozen decades after polar sunrise under solar wind exposure. Solar wind radiolysis of polar ethane frost synthesizes higher-order refractories that may contribute to the coloration of Charon's polar zones

Extreme exospheric dynamics at Charon: Implications for the red spot

Charon's exosphere may exhibit extreme seasonal dynamics, with centuries of quiescence punctuated by short lived (∼4 earth years) exospheric surges near the equinoxes, as spring sunrise bi-annually drives frozen methane off the polar night zones. Charon's pole-centric red spot has been proposed to be the product of Ly-α photolysis of frozen methane into refractory hydrocarbon “tholins”, but the role of exospheric dynamics in the red material's formation has not been investigated. We show with exospheric modeling that methane “polar-swap”, in which exospheric CH4 sublimated from the spring polar zone is rapidly re-frozen onto the autumn hemisphere, deposits ∼30 μm polar frosts too thick for Ly-α light to penetrate. Ethane, the primary methane photoproduct under these conditions, may unlike methane remain frozen decades after polar sunrise under solar wind exposure. Solar wind radiolysis of polar ethane frost synthesizes higher-order refractories that may contribute to the coloration of Charon's polar zones