Abstract. Many marine-terminating outlet glaciers have retreated rapidly in recent decades, but these changes have not been formally attributed to anthropogenic climate change. A key challenge for such an attribution assessment is that if glacier termini are sufficiently perturbed from bathymetric highs, ice-dynamic feedbacks can cause rapid retreat even without further climate forcing. In the presence of internal climate variability, attribution thus depends on understanding whether (or how frequently) these rapid retreats could be triggered by climatic noise alone. Our simulations with idealized glaciers show that in a noisy climate, rapid retreat is a stochastic phenomenon. We therefore propose a probabilistic approach to attribution and present a framework for analysis that uses ensembles of many simulations with independent realizations of random climate variability. Synthetic experiments show that century-scale climate trends substantially increase the likelihood of rapid glacier retreat. This effect depends on the timescales over which ice dynamics integrate forcing. For a population of synthetic glaciers with different topographies, we find that external trends increase the number of large retreats triggered within the population, offering a metric for regional attribution. Our analyses suggest that formal attribution studies are tractable and should be further pursued to clarify the human role in recent ice-sheet change. We emphasize that early-industrial-era constraints on glacier and climate state are likely to be crucial for such studies.
Abstract. We introduce the first version of the Stochastic Ice-sheet and Sea-level System Model (StISSM v1.0), which adds stochastic parameterizations within a state-of-the-art large-scale ice sheet model. In StISSM v1.0, stochastic parameterizations target climatic fields with internal variability, as well as glaciological processes exhibiting variability that cannot be resolved at the spatiotemporal resolution of ice sheet models: calving and subglacial hydrology. Because both climate and unresolved glaciological processes include internal variability, stochastic parameterizations allow StISSM v1.0 to account for the impacts of their high-frequency variability on ice dynamics and on the long-term evolution of modeled glaciers and ice sheets. StISSM v1.0 additionally includes statistical models to represent surface mass balance and oceanic forcing as autoregressive processes. Such models, once appropriately calibrated, allow users to sample irreducible uncertainty in climate prediction without the need for computationally expensive ensembles from climate models. When combined together, these novel features of StISSM v1.0 enable quantification of irreducible uncertainty in ice sheet model simulations and of ice sheet sensitivity to noisy forcings. We detail the implementation strategy of StISSM v1.0, evaluate its capabilities in idealized model experiments, demonstrate its applicability at the scale of a Greenland ice sheet simulation, and highlight priorities for future developments. Results from our test experiments demonstrate the complexity of ice sheet response to variability, such as asymmetric and/or non-zero mean responses to symmetric, zero-mean imposed variability. They also show differing levels of projection uncertainty for stochastic variability in different processes. These features are in line with results from stochastic experiments in climate and ocean models, as well as with the theoretical expected behavior of noise-forced non-linear systems.
Abstract. Increasing melt of ice sheets at their floating or vertical interfaces with the ocean is a major driver of marine ice sheet retreat and sea level rise. However, the extent to which warm, salty seawater may drive melting under the grounded portions of ice sheets is still not well understood. Previous work has explored the possibility that dense seawater intrudes beneath relatively light subglacial freshwater discharge, similar to the âsalt wedgeâ observed in many estuarine systems. In this study, we develop a generalized theory of layered seawater intrusion under grounded ice, including where subglacial hydrology occurs as a macroporous water sheet over impermeable beds or as microporous Darcy flow through permeable till. Using predictions from this theory, we show that seawater intrusion over flat or reverse-sloping impermeable beds may feasibly occur up to tens of kilometers upstream of a glacier terminus or grounding line. On the other hand, seawater is unlikely to intrude more than tens of meters through permeable till. Simulations using the Ice-sheet and Sea-level System Model (ISSM) show that even just a few hundred meters of basal melt caused by seawater intrusion upstream of marine ice sheet grounding lines can cause projections of marine ice sheet volume loss to be 10â%â50â% higher. Kilometers of intrusion-induced basal melt can cause projected ice sheet volume loss to more than double. These results suggest that further observational, experimental and numerical investigations are needed to determine the conditions under which seawater intrusion occurs and whether it will indeed drive rapid marine ice sheet retreat and sea level rise in the future.
Ice sheet models use observations to infer basal shear stress, but the variety of methods and datasets available has resulted in a wide range of estimates. Radar-based metrics such as reflectivity and specularity content have been used to characterize subglacial hydrologic conditions that are linked to spatial variations in basal shear stress. We explore whether radar metrics can be used to inform models about basal shear stress. At Thwaites Glacier, West Antarctica, we sample basal shear stress inversions across a wide range of ice sheet models to see how the basal shear stress distribution changes in regions of varying relative reflectivity and specularity content. Our results reveal three key findings: (1) Regions of high specularity content exhibit lower mean basal shear stresses (2) Wet and bumpy regions, as characterized by high relative reflectivity and low specularity content, exhibit higher mean basal shear stresses (3) Models disagree about what basal shear stress should be at the onset of rapid ice flow and high basal melt where relative reflectivity and specularity content are low.
We develop a statistical method to generate ocean forcing boundary conditions for Greenland ice sheet model simulations.• The method bias-corrects and extrapolates global climate model output using reanalysis products and high-resolution model results.• Stochastic time series models reproduce the spatiotemporal variability of ocean conditions at negligible computational expense.
Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on prior developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.
Subglacial drainage networks regulate the response of ice sheet flow to surface meltwater input to the subglacial environment. Simulating subglacial hydrology evolution is critical to projecting ice sheet sensitivity to climate, and contribution to sea-level change. However, current numerical subglacial hydrology models are computationally expensive, and, consequently, evolving subglacial hydrology is neglected in large-scale ice sheet simulations. We present a deep learning emulator of a state-of-the-art subglacial hydrology model, trained at multiple Greenland glaciers. Our emulator performs strongly in both temporal (R2>0.99) and spatial (R2>0.96) generalization, offers high computational savings, and can be used to force numerical ice sheet models. This will enable century- and large-scale ice sheet model simulations, including interactions between ice flow and increased meltwater input to the subglacial environment. Generally, our work demonstrates that machine learning can further improve ice sheet models, reduce computational bottlenecks, and exploit information from high-fidelity models and novel observational platforms.
Abstract. Increasing melt of ice sheets at their floating or vertical interfaces with the ocean is a major driver of marine ice sheet retreat and sea level rise. However, the extent to which warm, salty seawater may drive melting under the grounded portions of ice sheets is still not well understood. Previous work has explored the possibility that dense seawater intrudes beneath relatively light subglacial freshwater discharge, similar to the âsalt wedgeâ observed in many estuarine systems. In this study, we develop a generalized theory of layered seawater intrusion under grounded ice, including where subglacial hydrology occurs as a macroporous water sheet over impermeable beds or as microporous Darcy flow through permeable till. Using predictions from this theory, we show that seawater intrusion over flat or reverse-sloping impermeable beds may feasibly occur up to tens of kilometers upstream of a glacier terminus or grounding line. On the other hand, seawater is unlikely to intrude more than tens of meters through permeable till. Simulations using the Ice-sheet and Sea-level System Model (ISSM) show that even just a few hundred meters of basal melt caused by seawater intrusion upstream of marine ice sheet grounding lines can cause projections of marine ice sheet volume loss to be 10â%â50â% higher. Kilometers of intrusion-induced basal melt can cause projected ice sheet volume loss to more than double. These results suggest that further observational, experimental and numerical investigations are needed to determine the conditions under which seawater intrusion occurs and whether it will indeed drive rapid marine ice sheet retreat and sea level rise in the future.
