Stability of polar frosts in spherical bowl-shaped craters on the Moon, Mercury, and Mars☆
Abstract
Following Svitek (Martian Water Frost: Control of Global Distribution by Small-Scale Processes, Ph.D Thesis, California Institute of Technology, 1992), analytic solutions are presented for the effective albedo, the effective emissivity, and the radiative equilibrium temperature in the shadowed portions of a spherical bowl-shaped crater. The model assumes that the surface is a Lambert scatterer with visual albedo and infrared emissivity each independent of wavelength across their respective spectral ranges. Absorption, emission, and multiple scattering from the walls of the crater are treated rigorously to all orders. For airless bodies whose surfaces are in radiative equilibrium, all shadowed portions of any individual crater have the same temperature, whose value depends on four quantities: the insolation (product of the solar constant and the sine of the solar elevation angle), the depth/diameter ratio of the crater, the visual albedo, and the infrared emissivity. As long as the crater is deep enough to have shadows, the lowest temperatures are for the shallowest crater—those with the smallest depth/diameter ratio. The model is applied first to the Moon and Mercury using a depth/diameter ratio of 0.2, which is typical of the lunar highlands according to Pike (Geophys. Res. Lett. 1, 291–294 (1974); in Impact and Explosion Cratering (Roddy et al., Eds.), pp. 489–509, Pergamon, New York, 1977). For Mercury and the Moon, temperatures in shadows in polar craters are below 102 K, so the sublimation rate of water ice calculated according to the model of Watson et al. (J. Geophys. Res. 66, 3033-3015 (1961)) is less than 1 cm per byr. The latitudinal extent of the cold zone on the Moon is greater than that on Mercury, although temperatures at the poles of the two planets are similar. The other application is to polar frosts on Mars. Illuminated water frosts in radiative equilibrium grow rougher, because the average temperature of a depression is greater than that of flat ground. Subliming CO2 frosts, which are always at the same temperature, grow rougher at low solar elevation angles because the heat flux absorbed by a depression is greater than that for a flat surface. At high insolation rates (high Sun near perihelion) the average heat flux to a depression is less than for a flat surface. The latter evaporates faster, which makes the average surface smoother and leads to a high average albedo. This behavior helps explain the fact that the south CO2 cap, which receives its greatest insolation near perihelion, has a higher effective albedo and therefore can survive the summer, whereas the north CO2 cap has a lower effective albedo and disappears each year around summer solstice.
References (19)
- J.E. Colwell et al.
Evolution of topography on comets. II. Icy craters and trenches
Icarus
(1990) - B. Hapke
Bidirectional reflectance spectroscopy
Icarus
(1984) - A.J. Verbiscer et al.
Scattering properties of natural snow and frost: Comparison with icy satellite photometry
Icarus
(1990) - A.S. Adorjan
Temperature distribution in shadowed lunar craters
J. Spacecraft
(1970) - C.W. Allen
Astrophysical Quantities
(1973) - C.E. Bryson et al.
Sublimation rates and vapor pressures of H2O, CO2, N2O, and Xe
J. Chem. Eng. Data
(1974) - B. Butler et al.
Mercury: Full disk radar images, and the detection and stability of ice at the North Pole
(1992) - J.K. Harmon et al.
Radar mapping of Mercury: Full-disk delay-doppler images
(1992) - R.R. Hodges
Lunar cold traps and their influence on argon-40
Cited by (132)
Multi-mission, multi-sensor study of the Shackleton Crater constrained for volatiles with emphasis on albedo distribution of the Lunar South Pole
2024, Advances in Space ResearchShackleton Crater on the Moon is named in honor of a pivotal figure of the Heroic Age of Antarctic Exploration Sir Ernest Henry Shackleton. The Crater was formed about 3.6 billion years ago at the Lunar South Pole. Due to the low (1.5°) tilt of the lunar axis of rotation, the interiors of the crater receive almost no direct sunlight and observations by passive hyperspectral sensors are only possible due to reflected light. Active sensors are however able to illuminate the interior surface of the crater and make observations. Here we present the results of study of the crater using data from multiple sensors onboard two missions of Indian Space Research Organization (ISRO) and National Aeronautics and Space Administration (NASA). Hyperspectral observations from Moon Mineralogy Mapper (M3) on-board Chandrayaan-1, temperature from Diviner Lunar Radiometer Experiment (DLRE), albedo from Lunar Orbiter Laser Altimeter (LOLA) and Off/On band ratios from Lyman-Alpha Mapping Project (LAMP) on board Lunar Reconnaissance Orbiter (LRO) have been used. Permanently Shaded Regions (PSRs) are identified and Shackleton has been found, as reported, to be a PSR. Emphasis has been laid on LOLA albedo analysis of lunar south pole (80°-90° South). It is observed that the Shackleton Crater exhibits significantly higher albedo compared to its surrounding regions. Some regions exhibit albedo even higher than the Shackleton Crater in the lunar south pole region. The findings of high albedo, coupled with high Off/On band ratio of LAMP, extremely low DLRE temperatures in the permanently shaded regions and detection of hydration features in the M3 hyperspectral data suggests higher probability of hydration in the Shackleton Crater.
Fast hierarchical low-rank view factor matrices for thermal irradiance on planetary surfaces
2023, Journal of Computational Physics: XWe present an algorithm for compressing the radiosity view factor model commonly used in radiation heat transfer and computer graphics. We use a format inspired by the hierarchical off-diagonal low rank format, where elements are recursively partitioned using a quadtree or octree and blocks are compressed using a sparse singular value decomposition—the hierarchical matrix is assembled using dynamic programming. The motivating application is time-dependent thermal modeling on vast planetary surfaces, with a focus on permanently shadowed craters which receive energy through indirect irradiance. In this setting, shape models are comprised of a large number of triangular facets which conform to a rough surface. At each time step, a quadratic number of triangle-to-triangle scattered fluxes must be summed; that is, as the sun moves through the sky, we must solve the same view factor system of equations for a potentially unlimited number of time-varying righthand sides. We first conduct numerical experiments with a synthetic spherical cap-shaped crater, where the equilibrium temperature is analytically available. We also test our implementation with triangle meshes of planetary surfaces derived from digital elevation models recovered by orbiting spacecraft. Our results indicate that the compressed view factor matrix can be assembled in quadratic time, which is comparable to the time it takes to assemble the full view matrix itself. Memory requirements during assembly are reduced by a large factor. Finally, for a range of compression tolerances, the size of the compressed view factor matrix and the speed of the resulting matrix vector product both scale linearly (as opposed to quadratically for the full matrix), resulting in orders of magnitude savings in processing time and memory space.
Modeling the thermal environment of Mercury's north pole using MLA. Implications for locations of water ice
2023, IcarusRadar observations made by the Arecibo Observatory have revealed radar-bright deposits in the north polar region of Mercury which are consistent with water ice. These deposits are generally associated with floors of deep craters in permanent shadow. In this paper, we revisit topographic data, illumination, and associated thermal conditions of the north polar region of Mercury to investigate these radar-bright deposits. We used Mercury Laser Altimeter (MLA) data from the MESSENGER mission to derive a dedicated, adjusted high-resolution Digital Terrain Model (DTM) of the north polar region at 250 m/pixel. Based on this DTM, illumination conditions and temperature were derived for the first time on the fully available MLA dataset. Generally, craters equatorward of 86 N are found to be too hot to contain surface water ice due to effective self-heating. For those warm craters that do contain radar-bright deposits, temperatures within the upper 10–20 cm of regolith are cold enough for water ice to be stable. Using our improved DTM, we confirm previously MLA-modeled illumination conditions in Prokofiev crater and argue that radar-bright patches outside the PSR in the area may be explained by icy deposits in cold and shallow (within 5 cm) sub-surface layers. In the radar observations we find a strong, unphysical correlation between radar-bright pixels and the number of radar observations per pixel. Consequently, several PSRs which were classified as free of radar-bright material (mainly at longitudes 120–340 E) are most likely subject to a bias in the averaging process in deriving the composite radar image and might contain water ice after all.
The origins of low-albedo material on Europa’s surface have been of interest since Voyager first returned close-up images of the icy moon. Material ejected from Io is known to contribute an exogenic flux of dark material to Europa’s trailing hemisphere, and hydrated salt compounds concentrated within chaos terrain, ridges, and pits may be endogenous to the subsurface ocean. Many of Europa’s impact craters also exhibit dark ejecta, the origins of which are unknown. Our study examines the ejecta of several large impact craters to determine possible origins for this dark material. We compared the dark material found in crater ejecta to other surface materials using Galileo Near-Infrared Mapping Spectrometer data to assess similarities in composition between ejecta material and other dark materials on Europa’s surface. Our analysis shows that dark material found in crater ejecta has similar composition to other dark features on Europa and may be the result of comparable sources or alteration processes. We also considered dark impactors as sources for the dark ejecta material. Using crater scaling laws, we estimated the impactor size for each crater and determined the impactor’s potential contribution of dark material to the ejecta. We then compared these quantities to those derived using a radiative transfer model and the measured reflectance values of each dark-ejecta crater. Our model results show that the lower albedo of the ejecta of these craters cannot be solely attributed to an intimate mixture of the impactor material with the target material. In contrast, modeled impact heating and vaporization suggest sufficient amounts of ice could be removed in order to explain the observed low-albedo patterns, if preexisting or impactor-derived dark material is just 0.1% more resistant to vaporization than the ice. Given the lack of spatial correlation, and the presence of similar-sized craters without dark ejecta, these results point to either localized differences in the concentration of dark non-ice materials in Europa’s shell, or variations in impact velocities and geometries leading to differences in the amount of vaporized ejecta.
Volatile interactions with the lunar surface
2022, GeochemistryThe Moon is generally depleted in volatile elements and this depletion extends to the surface where the most abundant mineral, anorthite, features <6 ppm H2O. Presumably the other nominally anhydrous minerals that dominate the mineral composition of the global surface—olivine and pyroxene—are similarly depleted in water and other volatiles. Thus the Moon is tabula rasa for the study of volatiles introduced in the wake of its origin. Since the formation of the last major basin (Orientale), volatiles from the solar wind, from impactors of all sizes, and from volatiles expelled from the interior during volcanic eruptions have all interacted with the lunar surface, leaving a volatile record that can be used to understand the processes that enable processing, transport, sequestration, and loss of volatiles from the lunar system. Recent discoveries have shown the lunar system to be complex, featuring emerging recognition of chemistry unanticipated from the Apollo era, confounding issues regarding transport of volatiles to the lunar poles, the role of the lunar regolith as a sink for volatiles, and the potential for active volatile dynamics in the polar cold traps. While much has been learned since the overturn of the “Moon is dry” paradigm by innovative sample and spacecraft measurements, the data point to a more complex lunar volatile environment than is currently perceived.
Sub-field of view surface thermal modeling of Cassini CIRS observations of Rhea during south polar winter
2022, IcarusRhea's exosphere is thought to originate from sources of carbon, water ice and other volatiles that arrived at Rhea by bombardment. Its seasonal variability is directly driven by polar surface temperatures allowing surface adsorption, and the persistence of source volatiles require that seasonal temperatures remain sufficiently cold to retain them. Cassini CIRS detected temperatures in Rhea's winter polar region of 23 K, amongst the coldest measured in the Solar System, but with a relatively large footprint we seek to add value to these observations by modeling sub field-of-view (FOV) temperature distribution and examining how the coldest scenes evolve on a seasonal basis.
A simple, rough surface 1-dimensional thermal model is developed using a digital elevation map as input to a thermal model. We compared averaged rough and flat modeled FOVs to CIRS temperatures for a set of case study observations in the south polar region in winter darkness and found they both agree within expected CIRS error in all cases. We develop an asymmetric estimate of CIRS FOV temperature uncertainty, which is particularly important for very cold spectra to accurately represent the sensitivity of the instrument. This approach provides a conservative upper temperature limit to spectral fits whilst highlighting there is often little information to constrain the lower bound of temperature uncertainty in these cases.
The distribution of modeled sub-FOV facet temperatures is explored for the full range of azimuth angles and slope gradients. More extensive and cooler temperatures were found in the rough model, complemented by fewer but warmer areas than the flat model. We find temperature contrasts of individual model facets of up to 15 K warmer and 11 K colder within some CIRS FOV when the rough model was compared flat, in a case study of scenes in winter darkness. This is due to the persistence of the seasonal thermal wave beneath the surface. We tested model sensitivity to thermal input parameters (thermal inertia and bolometric Bond albedo). These values are challenging to constrain due to limited observations and measurement noise and are expected to vary with subtle changes in surface characteristics. We found that within 20° of the pole, temperatures do not rise significantly above 80 K all year, implying that a variety of simple organic molecules, linear amides and other carbon-containing compounds would remain stable on a quasi-permanent basis, potentially until photolytic or radiolytic mechanisms liberate by-products which provide a source for exospheric constituents. The simple rough model indicates that some facets experience very different summer/winter season lengths than a flat model can represent, which is important in terms of the exospheric sequestration process.
We use threshold temperatures of 55 K for CO2 and 30 K for O2 ice as proxies for their relative stability at the surface. The surface coverage of these temperatures with time was compared to modeled exospheric abundance by Teolis and Waite (2016). They are anti-correlated which is expected, but with slight differences in timings of minima and the equinoctial maxima. We find that the rough model suggests larger areas able to adsorb these species than the flat model, and that the cumulative influence of topography would have a direct relationship on the timing and abundance of the exosphere.
We place Rhea's polar environment in context with other icy moons in the Saturnian system using the flat model. The satellites Tethys, Dione and Enceladus experience similar polar thermal regimes to Rhea (with exceptionally cold winters) and would potentially benefit from the consideration of topography in relation to exospheric modeling.
- ☆
Contribution number 5186 from the Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California 91125.