Two outstanding questions for future Greenland predictions are (1) how enhanced meltwater draining beneath the ice sheet will impact the behavior of large tidewater glaciers, and (2) to what extent tidewater glacier velocity is driven by changes at the terminus versus changes in sliding velocity due to meltwater input. We present a two-way coupled framework to simulate the nonlinear feedbacks of evolving subglacial hydrology and ice dynamics using the Subglacial Hydrology And Kinetic, Transient Interactions (SHAKTI) model within the Ice-sheet and Sea-level System Model (ISSM). Through coupled simulations of Helheim Glacier, we find that terminus effects dominate the seasonal velocity pattern up to 15 km from the terminus, while hydrology primarily drives the velocity response upstream. With increased melt, the hydrology influence yields seasonal acceleration of several hundred meters per year in the interior, suggesting that hydrologic forcing will play an important role in future mass balance of tidewater glaciers.

The transport of meltwater through porous snow is a fundamental process in hydrology that remains poorly understood but essential for more robust prediction of how the cryosphere will respond under climate change. Here we propose a continuum model that resolves the nonlinear coupling of preferential melt flow and the nonequilibrium thermodynamics of ice-melt phase change at the Darcy scale. We assume that the commonly observed unstable melt infiltration is due to the gravity fingering instabililty, and capture it using the modified Richards equation that is extended with a higher-order term in saturation. Our model accounts for changes in porosity and the thermal budget of the snowpack caused by melt refreezing at the continuum scale, based on a mechanistic estimate of the ice-water phase change kinetics formulated at the pore scale. We validate the model in 1D against field data and laboratory experiments of infiltration in snow and find generally good agreement. Compared to existing theory of stable melt infiltration, our 2D simulation results show that preferential infiltration delivers melt faster to deeper depths, and as a result, changes in porosity and temperature can occur at deeper parts of the snow. The simulations also capture the formation of vertical low porosity annulus known as ice pipes, which have been observed in the field but lack mechanistic understanding to date. Our results demonstrate how melt refreezing and unstable infiltration reshape the porosity structure of snow and impacts thermal and mass transport in highly nonlinear ways, which are not captured by simpler models.

Ocean worlds have been identified as high-priority astrobiology targets due to the link between life and liquid water. Young surface terrain on many icy bodies indicates they support active geophysical cycles that may facilitate ocean-surface transport that could provide observables for upcoming missions. Accurately interpreting spacecraft observations requires constraining the relationship between ice shell characteristics and interior dynamics. On Earth, the composition, physical characteristics, and bioburden of ocean-derived ices are related to their formation history and parent fluid composition. In such systems the ice-ocean interface, which exists as a multiphase mushy layer, dictates the overlying ice’s properties and evolution. Inclusion of the physics governing these boundaries is a novel strategy in modeling planetary ices, and thus far has been limited to 1D approaches. Here we present results from 2D simulations of an archetypal ice-ocean world. We track the evolution of temperature, salinity, porosity, and brine velocity within a thickening ice shell enabling us to place improved constraints on ice-ocean world properties, including: the composition of planetary ice shells, the thickness and hydraulic connectivity of ice-ocean interfaces, and heterogeneous dynamics/structures in the interfacial mushy layer. We show that stable eutectic horizons are likely a common feature of ice-ocean worlds and that ocean composition plays an important role in governing the structure and dynamics of the interface, including the formation of chemical gradient-rich regions within the mushy layer. We discuss the geophysical and astrobiological implications of our results and highlight how they can be validated by instrument specific measurements.

Although the nature of the early Martian climate is a matter of considerable debate, the presence of valley networks (VN) provides unambiguous evidence for the presence of liquid water on Mars’ surface. A subaerial fluvial origin of VN is at odds with the expected phase instability of near-surface water in the cold, dry Late Noachian climate. Furthermore, many observed geomorphometric properties of VN are inconsistent with surface water flow. Conversely, subglacial channels exhibit many of these characteristics and could have persisted beneath ice sheets even in a cold climate. Here we model basal melting beneath a Late Noachian Icy Highlands ice sheet and map subglacial hydrological flow paths to investigate the distribution and geomorphometry of subglacial channels. We show that subglacial processes produce enough melt water to carve Mars’ VN; that predicted channel distribution is consistent with observations; and corroborate geomorphometric measurements of VN consistent with subglacial formation mechanisms. We suggest that subglacial hydrology may have played a key role in the surface modification of Mars.

The cause of Heinrich events and their relationship with Dansgaard-Oeschger (DO) events are not fully understood. Previous modeling studies have argued that Heinrich events result from either internal oscillations generated within ice sheets or ocean warming occurring during DO events. In this study, we present a coupled model of ice stream and ocean dynamics to evaluate the behavior of the coupled system with few degrees of freedom and minimal parameterizations. Both components of the model may oscillate independently, with stagnant versus active phases for the ice stream model and strong versus weak Atlantic Meridional Overturning Circulation (AMOC) phases for the ocean model. The ice sheet and ocean interact through submarine melt at the ice stream grounding line and freshwater flux into the ocean from ice sheet discharge. We show that these two oscillators have a strong tendency to synchronize, even when their interaction is weak, due to the amplification of small perturbations typical in nonlinear oscillators. In syn- chronized regimes with ocean-induced melt at the ice stream grounding line, Heinrich events always follow DO events by a constant time lag. We also introduce noise into the ocean system and find that ice-ocean interactions not only maintain a narrow distribu- tion of timing between Heinrich and DO events, but also regulate DO event periodic- ity against noise in the climate system. This synchronization persists across a broad range of parameters, indicating that it is a robust explanation for Heinrich events and their timing despite the significant uncertainty associated with past ice sheet conditions.

We investigate the structure and evolution of multiphase ice-ocean interfaces (‘mushy layers’) and the implications for the geophysics and habitability of ice-ocean worlds. Understanding the potential diversity of these multiphase layers across solar system bodies provides insight into the potential rates and mechanisms of heat and solute transport between their respective oceans and ice shells - which remain largely unconstrained. Additionally, variations in mushy layer properties may drive diverse geophysical processes unique to individual bodies or that may vary regionally on an individual icy world. We explore mushy layer evolution by analytically solving for the thickness of a simplified ice-ocean mushy layer system. We investigate two dynamic regimes, one driven by molecular diffusion and one driven by convection of brine within the mushy layer. We analyze the impact of gravity, thermal gradient, and ocean composition on the thickness of mushy layers. Additionally, a perturbation analysis is carried out to investigate the existence of mushy layer steady states. We show that stable mushy layers exist when ice shells are thickening, suggesting that mushy layers are likely persistent and common features of growing ice shells and accretionary regions of ice-ocean worlds.

Most of the mass loss from the Antarctic Ice Sheet occurs through glaciers and ice streams, where fast-flow is partially controlled by rapid ice deformation in the margins. Deformation drives thermomechanical and recrystallization processes that influence further deformation, a feedback which may destabilize glaciers. However, few models account for the feedback between deformation and recrystallization. We derive an idealized model for ice temperature and grain-size that partitions deformational work into dissipated heat and changes in strain and surface energy, all of which drive dynamic recrystallization. Under conditions common in glacier shear margins, we show that a large portion of deformational work is stored as elastic energy, with the remainder dissipated as heat. This result revises our current picture of the amount of heat generated in glacier shear margins and suggests that changes in internal strain through dynamic recrystallization of ice likely play an important role in facilitating fast-flowing glacial ice.

Non-ice impurities within the ice shells of ocean worlds (e.g., Europa, Enceladus, Titan) are believed to play a fundamental role in their geophysics and habitability and may become a surface expression of subsurface ocean properties. Heterogeneous entrainment and distribution of impurities within planetary ice shells have been proposed as mechanisms that can drive ice shell overturn, generate diverse geological features, and facilitate ocean-surface material transport critical for maintaining a habitable subsurface ocean. However, current models of ice shell composition suggest that impurity rejection at the ice-ocean interface of thick contemporary ice shells will be exceptionally efficient, resulting in relatively pure, homogeneous ice. As such, additional mechanisms capable of facilitating enhanced and heterogeneous impurity entrainment are needed to reconcile the observed physicochemical diversity of planetary ice shells. Here we investigate the potential for hydrologic features within planetary ice shells (sills and basal fractures), and the unique freezing geometries they promote, to provide such a mechanism. By simulating the two-dimensional thermal and physicochemical evolution of these hydrological features as they solidify, we demonstrate that bottom-up solidification at sill floors and horizontal solidification at fracture walls generate distinct ice compositions and provide mechanisms for both enhanced and heterogeneous impurity entrainment. We compare our results with magmatic and metallurgic analogs that exhibit similar micro- and macroscale chemical zonation patterns during solidification. Our results suggest variations in ice-ocean/brine interface geometry could play a fundamental role in introducing compositional heterogeneities into planetary ice shells and cryoconcentrating impurities in (re)frozen hydrologic features.