Ice sintering is a form of metamorphism that drives the microstructural evolution of an aggregate of grains through surface and volume diffusion. This leads to an increase in the grain-to-grain contact area ("neck") and density of the aggregate over time, resulting in the evolution of its strength, porosity, thermal conductivity, and other properties. This process plays an important role in the evolution of icy planetary surfaces, though its rate and nature are not well constrained. In this study, we explore the model of Swinkels and Ashby (1981), and assess the extent to which it can be used to quantify sintering timescales for water ice. We compare predicted neck growth rates to new and historical observations of ice sintering, and find agreement to some studies at the order of magnitude level. First-order estimates of neck growth timescales on planetary surfaces show that ice may undergo significant modification over geologic timescales, even in the outer solar system. Densification occurs over much longer timescales, suggesting some surfaces may develop cohesive, but porous, crusts. Sintering rates are extremely sensitive to temperature and grain size, occurring faster in warmer aggregates of smaller grains. This suggests that the microstructural evolution of ices may vary not only throughout the solar system, but also spatially across the surface and in the near-surface of a given body. Our experimental observations of complex grain growth and mass redistribution in ice aggregates point to components of the model that may benefit from improvement, and areas where additional laboratory studies are needed.

This work investigates the macroscopic thermomechanical behavior of lunar boulders by modeling their response to diurnal thermal forcing. Our results reveal a bimodal, spatiotemporally-complex stress response. During sunrise, stresses occur in the boulders' interiors that are associated with large-scale temperature gradients developed due to overnight cooling. During sunset, stresses occur at the boulders' exteriors due to the cooling and contraction of the surface. Both kinds of stresses are on the order of 10 MPa in 1 m boulders and decrease for smaller diameters, suggesting that larger boulders break down more quickly. Boulders <30 cm exhibit a weak response to thermal forcing, suggesting a threshold below which crack propagation may not occur. Boulders of any size buried by regolith are shielded from thermal breakdown. As boulders increase in size (>1 m), stresses increase to several 10s of MPa as the behavior of their surfaces approaches that of an infinite halfspace. As the thermal wave loses contact with the boulder interior, stresses become limited to the near-surface. This suggests that the survival time of a boulder is not only controlled by the amplitude of induced stress, but also by its diameter as compared to the diurnal skin depth. While stresses on the order of 10 MPa are enough to drive crack propagation in terrestrial environments, crack propagation rates in vacuum are not well constrained. We explore the relationship between boulder size, stress, and the direction of crack propagation, and discuss the implications for the relative breakdown rates and estimated lifetimes of boulders on airless body surfaces.

Thermomechanical processes such as fatigue and shock have been suggested to cause and contribute to rock breakdown on Earth, and on other planetary bodies, particularly airless bodies in the inner solar system. In this study, we modeled grain-scale stresses induced by diurnal temperature variations on simple microstructures made of pyroxene and plagioclase on various solar system bodies. We found that a heterogeneous microstructure on the Moon experiences peak tensile stresses on the order of 100 MPa. The stresses induced are controlled by the coefficient of thermal expansion and Young's modulus of the mineral constituents, and the average stress within the microstructure is determined by relative volume of each mineral. Amplification of stresses occurs at surface-parallel boundaries between adjacent mineral grains and at the tips of pore spaces. We also found that microscopic spatial and temporal surface temperature gradients do not correlate with high stresses, making them inappropriate proxies for investigating microcrack propagation. Although these results provide very strong evidence for the significance of thermomechanical processes on airless bodies, more work is needed to quantify crack propagation and rock breakdown rates.