Ice Sheet Model Research Articles

Modeling the timing of Patagonian Ice Sheet retreat in the Chilean Lake District from 22–10 ka

Abstract. Studying the retreat of the Patagonian Ice Sheet (PIS) during the last deglaciation represents an important opportunity to understand how ice sheets outside the polar regions have responded to deglacial changes in temperature and large-scale atmospheric circulation. At the northernmost extension of the PIS during the Last Glacial Maximum (LGM), the Chilean Lake District (CLD) was influenced by the southern westerly winds (SWW), which strongly modulated the hydrologic and heat budgets of the region. Despite progress in constraining the nature and timing of deglacial ice retreat across this area, considerable uncertainty in the glacial history still exists due to a lack of geologic constraints on past ice margin change. Where the glacial chronology is lacking, ice sheet models can provide important insight into our understanding of the characteristics and drivers of deglacial ice retreat. Here we apply the Ice Sheet and Sea-level System Model (ISSM) to simulate the LGM and last deglacial ice history of the PIS across the CLD at high spatial resolution (450 m). We present a transient simulation of ice margin change across the last deglaciation using climate inputs from the National Center for Atmospheric Research Community Climate System Model (CCSM3) Trace-21ka experiment. At the LGM, the simulated ice extent across the CLD agrees well with the most comprehensive reconstruction of PIS ice history (PATICE). Coincident with deglacial warming, ice retreat ensues after 19 ka, with large-scale ice retreat occurring across the CLD between 18 and 16.5 ka. By 17 ka, the northern portion of the CLD becomes ice free, and by 15 ka, ice only persists at high elevations as mountain glaciers and small ice caps. Our simulated ice history agrees well with PATICE for early deglacial ice retreat but diverges at and after 15 ka, where the geologic reconstruction suggests the persistence of an ice cap across the southern CLD until 10 ka. However, given the high uncertainty in the geologic reconstruction of the PIS across the CLD during the later deglaciation, this work emphasizes a need for improved geologic constraints on past ice margin change. While deglacial warming drove the ice retreat across this region, sensitivity tests reveal that modest variations in wintertime precipitation (∼10 %) can modulate the pacing of ice retreat by up to 2 ka, which has implications when comparing simulated outputs of ice margin change to geologic reconstructions. While we find that TraCE-21ka simulates large-scale changes in the SWW across the CLD that are consistent with regional paleoclimate reconstructions, the magnitude of the simulated precipitation changes is smaller than what is found in proxy records. From our sensitivity analysis, we can deduce that larger anomalies in precipitation, as found in paleoclimate proxies, may have had a large impact on modulating the magnitude and timing of deglacial ice retreat. This fact highlights an additional need for better constraints on the deglacial change in strength, position, and extent of the SWW as it relates to understanding the drivers of deglacial PIS behavior.

Glacial inception through rapid ice area increase driven by albedo and vegetation feedbacks

Abstract. We present transient simulations of the last glacial inception using the Earth system model CLIMBER-X with dynamic vegetation, interactive ice sheets, and visco-elastic solid Earth responses. The simulations are initialized at the middle of the Eemian interglacial (125 kiloyears before present, ka) and run until 100 ka, driven by prescribed changes in Earth's orbital parameters and greenhouse gas concentrations from ice core data. CLIMBER-X simulates a rapid increase in Northern Hemisphere ice sheet area through MIS5d, with ice sheets expanding over northern North America and Scandinavia, in broad agreement with proxy reconstructions. While most of the increase in ice sheet area occurs over a relatively short period between 119 and 117 ka, the larger part of the increase in ice volume occurs afterwards with an almost constant ice sheet extent. We show that the vegetation feedback plays a fundamental role in controlling the ice sheet expansion during the last glacial inception. In particular, with prescribed present-day vegetation the model simulates a global sea level drop of only ∼ 20 m, compared with the ∼ 35 m decrease in sea level with dynamic vegetation response. The ice sheet and carbon cycle feedbacks play only a minor role during the ice sheet expansion phase prior to ∼ 115 ka but are important in limiting the deglaciation during the following phase characterized by increasing summer insolation. The model results are sensitive to climate model biases and to the parameterization of snow albedo, while they show only a weak dependence on changes in the ice sheet model resolution and the acceleration factor used to speed up the climate component. Overall, our simulations confirm and refine previous results showing that climate–vegetation–cryosphere feedbacks play a fundamental role in the transition from interglacial to glacial states characterizing Quaternary glacial cycles.

A framework for estimating the anthropogenic part of Antarctica’s sea level contribution in a synthetic setting

The relative contributions of anthropogenic climate change and internal variability in sea level rise from the West Antarctic Ice Sheet are yet to be determined. Even the way to address this question is not yet clear, since these two are linked through ice-ocean feedbacks and probed using ice sheet models with substantial uncertainty. Here we demonstrate how their relative contributions can be assessed by simulating the retreat of a synthetic ice sheet setup using an ice sheet model. Using a Bayesian approach, we construct distributions of sea level rise associated with this retreat. We demonstrate that it is necessary to account for both uncertainties arising from both a poorly-constrained model parameter and stochastic variations in climatic forcing, and our distributions of sea level rise include these two. These sources of uncertainty have only previously been considered in isolation. We identify characteristic effects of climate change on sea level rise distributions in this setup, most notably that climate change increases both the median and the weight in tails of distributions. From these findings, we construct metrics quantifying the role of climate change on both past and future sea level rise, suggesting that its attribution is possible even for unstable marine ice sheets.

Multimethod dating of ice-rafted dropstones reveals hidden localized glacial erosion in Wilkes Subglacial Basin, Antarctica

Abstract The Antarctic ice sheet blankets >99% of the continent and limits our ability to study how subglacial geology and topography have evolved through time. Ice-rafted dropstones derived from the Antarctic subglacial continental interior at different times during the late Cenozoic provide valuable thermal history proxies to understand this geologic history. We applied multiple thermochronometers covering a range of closure temperatures (60–800 °C) to 10 dropstones collected during Integrated Ocean Drilling Program (IODP) Expedition 318 in order to explore the subglacial geology and thermal and exhumation history of the Wilkes Subglacial Basin. The Wilkes Subglacial Basin is a key target for study because ice-sheet models show it was an area of ice-sheet retreat that significantly contributed to sea-level rise during past warm periods. Depositional ages of dropstones range from early Oligocene to late Pleistocene and have zircon U-Pb or 40Ar/39Ar ages indicating sources from the Mertz shear zone, Adélie craton, Ferrar large igneous province, and Millen schist belt. Dropstones from the Mertz shear zone and Adélie craton experienced three cooling periods (1700–1500 Ma; 500–280 Ma; 34–0 Ma) and two periods of extremely slow cooling rates (1500–500 Ma; 280–34 Ma). Low-temperature thermochronometers from seven of the dropstones record cooling during the Paleozoic, potentially recording the Ross or Pan-African orogenies, and during the Mesozoic, potentially recording late Paleozoic to Mesozoic rifting. These dropstones then resided within ~500 m of the surface since the late Paleozoic and early Mesozoic. In contrast, two dropstones deposited during the mid-Pliocene, one from the Mertz shear zone and one from Adélie craton, show evidence for localized post-Eocene glacial erosion of ≥2 km.

Local forcing mechanisms challenge parameterizations of ocean thermal forcing for Greenland tidewater glaciers

Abstract. Frontal ablation has caused 32 %–66 % of Greenland Ice Sheet mass loss since 1972, and despite its importance in driving terminus change, ocean thermal forcing remains crudely incorporated into large-scale ice sheet models. In Greenland, local fjord-scale processes modify the magnitude of thermal forcing at the ice–ocean boundary but are too small scale to be resolved in current global climate models. For example, simulations used in the Ice Sheet Intercomparison Project for CMIP6 (ISMIP6) to predict future ice sheet change rely on the extrapolation of regional ocean water properties into fjords to drive terminus ablation. However, the accuracy of this approach has not previously been tested due to the scarcity of observations in Greenland fjords, as well as the inability of fjord-scale models to realistically incorporate icebergs. By employing the recently developed IceBerg package within the Massachusetts Institute of Technology general circulation model (MITgcm), we here evaluate the ability of ocean thermal forcing parameterizations to predict thermal forcing at tidewater glacier termini. This is accomplished through sensitivity experiments using a set of idealized Greenland fjords, each forced with equivalent ocean boundary conditions but with varying tidal amplitudes, subglacial discharge, iceberg coverage, and bathymetry. Our results indicate that the bathymetric obstruction of external water is the primary control on near-glacier thermal forcing, followed by iceberg submarine melting. Despite identical ocean boundary conditions, we find that the simulated fjord processes can modify grounding line thermal forcing by as much as 3 °C, the magnitude of which is largely controlled by the relative depth of bathymetric sills to the Polar Water–Atlantic Water thermocline. However, using a common adjustment for fjord bathymetry we can still predict grounding line thermal forcing within 0.2 °C in our simulations. Finally, we introduce new parameterizations that additionally account for iceberg-driven cooling that can accurately predict interior fjord thermal forcing profiles both in iceberg-laden simulations and in observations from Kangiata Sullua (Ilulissat Icefjord).

Numerical simulation of heterogeneous ice sheet-structure interaction based on cohesive element method

The coupling process between sea ice and offshore structures is complex, which involves the influences of physical and mechanical properties of sea ice, structural dynamics, and fluid–structure interaction. In this work, a numerical model was established to simulate the ice-structure interactions, which is based on the cohesive element method combined with random field. The complicated spatial heterogeneity of ice sheet materials was considered, and the fracture energy, which is one of the material properties of the ice sheet model, was assumed to be a Gaussian random field. Numerical simulations of the ice sheet-vertical structures interactions were conducted with two mesh patterns, respectively. A more suitable element was adopted to build the heterogeneous ice field to investigate the effect of heterogeneity on the interactions. Moreover, the parametric sensitivity analysis of Gaussian random fields was performed. It has been found that the heterogeneity of ice sheet led to significant uncertainties in the numerical simulations. The location of ice sheet failure and the shape of ice fragments are associated with the distribution of the properties, and the cracks emerge and propagate in those weaker regions. Furthermore, the mean of ice load and the collision frequency show a decreasing trend with an increase in correlation length.

Timing of Cordilleran-Laurentide ice-sheet separation: Implications for sea-level rise

During the last deglaciation, collapse of the saddle between the North American Cordilleran and Laurentide ice sheets led to rapid ice-sheet mass loss and separation, with meltwater discharge contributing to deglacial sea-level rise. We directly date ice-sheet separation at the end of the saddle collapse using 64 10Be exposure ages along an ∼1200-km transect of the ice-sheet suture zone. Collapse began in the south by 15.4 ± 0.4 ka and ended by 13.8 ± 0.1 ka at ∼56°N. Ice-sheet model simulations consistent with the 10Be ages find that the saddle collapse contributed 6.2–7.2 m to global mean sea-level rise from ∼15.5 ka to ∼14.0 ka, or approximately one third of global mean sea-level rise over this period. We determine 3.1–3.6 m of the saddle collapse meltwater was released during Meltwater Pulse 1A ∼14.6-14.3 ka, constituting 20–40% of this meltwater pulse's volume. Because the separation of the Cordilleran and Laurentide ice sheets occurred over 1–2 millennia, the associated release of meltwater during the saddle collapse supplied a smaller contribution to the magnitude of Meltwater Pulse 1A than has been recently proposed.

Holocene southwest Greenland ice sheet behavior constrained by sea-level modeling

The melting of the Greenland ice sheet (GrIS) is a major contributor to past and future global sea-level rise. Understanding the response of the GrIS to times in the past when temperatures were as warm or warmer than today offers insights into its current and future response to climate change. In the southwest sector, the GrIS retreated inland beyond its current margin during the (at least regionally) warmer-than-present mid-Holocene, before it readvanced to the historical maximum position during the Little Ice Age. This was then followed by a slight retreat to its current position. To investigate the timing and magnitude of southwest GrIS retreat and readvance in response to Holocene warmth, we model the response of the solid Earth and local relative sea level (RSL) to past ice sheet change. We test a suite of eleven ice sheet scenarios that are based on ICE-6G_C but are modified in the timing and magnitude of ice retreat and readvance and pair them with four different viscoelastic Earth structures. We compare model predictions to observations of paleo sea level, present-day sea-level change, and present-day vertical land motion (VLM) around Nuuk, Greenland. We find that the modeled timing and magnitude of the Holocene retreat and readvance have a significant impact on modern sea-level change and VLM in Nuuk. Models that assume a readvance approaching the southwest GrIS’ historical maximum between 2 and 1 ka are most consistent with observations. The RSL response, however, is less sensitive to the timing of the minimum GrIS extent. Nonetheless, better data-model fits are generally obtained when the minimum ice sheet extent is reached between 5 and 3 ka, within the tested range of 6-3 ka. Comparing this timing to local and regional records of temperature and ice-sheet change suggest that the evolution of the southwestern GrIS presented here was in-phase with the likely evolution of southwestern GrIS mass balance through the Holocene. Our results have implications for future ice sheet modeling studies targeting southwestern Greenland by providing additional constraints and strengthening existing ones. Moreover, this work provides a deeper understanding of the interactions between the climate and the cryosphere and thus of future ice sheet change.

Coupling MAR (Modèle Atmosphérique Régional) with PISM (Parallel Ice Sheet Model) mitigates the positive melt–elevation feedback

Abstract. The Greenland Ice Sheet is a key contributor to sea level rise. By melting, the ice sheet thins, inducing higher surface melt due to lower surface elevations, accelerating the melt coming from global warming. This process is called the melt–elevation feedback and can be considered by using two types of models: either (1) atmospheric models, which can represent the surface mass balance (SMB), or SMB estimates resulting from simpler models such as positive degree day models or (2) ice sheet models representing the surface elevation evolution. The latter ones do not represent the surface mass balance explicitly as well as polar-oriented climate models. A new coupling between the MAR (Modèle Atmosphérique Régional) regional climate model and the PISM (Parallel Ice Sheet Model) ice sheet model is presented here following the CESM2 (Community Earth System Model; SSP5-8.5, Shared Socioeconomic Pathway) scenario until 2100 at the MAR lateral boundaries. The coupling is extended to 2200 with a stabilised climate (+7 ∘C compared to 1961–1990) by randomly sampling the last 10 years of CESM2 to force MAR and reaches a sea level rise contribution of 64 cm. The fully coupled simulation is compared to a one-way experiment where surface topography remains fixed in MAR. However, the surface mass balance is corrected for the melt–elevation feedback when interpolated on the PISM grid by using surface mass balance vertical gradients as a function of local elevation variations (offline correction). This method is often used to represent the melt–elevation feedback and prevents a coupling which is too expensive in computation time. In the fully coupled MAR simulation, the ice sheet morphology evolution (changing slope and reducing the orographic barrier) induces changes in local atmospheric patterns. More specifically, wind regimes are modified, as well as temperature lapse rates, influencing the melt rate through modification of sensible heat fluxes at the ice sheet margins. We highlight mitigation of the melt lapse rate on the margins by modifying the surface morphology. The lapse rates considered by the offline correction are no longer valid at the ice sheet margins. If used (one-way simulation), this correction implies an overestimation of the sea level rise contribution of 2.5 %. The mitigation of the melt lapse rate on the margins can only be corrected by using a full coupling between an ice sheet model and an atmospheric model.

Disentangling the drivers of future Antarctic ice loss with a historically calibrated ice-sheet model

Abstract. We use an observationally calibrated ice-sheet model to investigate the future trajectory of the Antarctic ice sheet related to uncertainties in the future balance between sub-shelf melting and ice discharge, on the one hand, and the surface mass balance, on the other. Our ensemble of simulations, forced by a panel of climate models from the sixth phase of the Coupled Model Intercomparison Project (CMIP6), suggests that the ocean will be the primary driver of short-term Antarctic mass loss, initiating ice loss in West Antarctica already during this century. The atmosphere initially plays a mitigating role through increased snowfall, leading to an Antarctic contribution to global mean sea-level rise by 2100 of 6 (−8 to 15) cm under a low-emission scenario and 5.5 (−10 to 16) cm under a very high-emission scenario. However, under the very high-emission pathway, the influence of the atmosphere shifts beyond the end of the century, becoming an amplifying driver of mass loss as the ice sheet's surface mass balance decreases. We show that this transition occurs when Antarctic near-surface warming exceeds a critical threshold of +7.5 ∘C, at which the increase in surface runoff outweighs the increase in snow accumulation, a signal that is amplified by the melt–elevation feedback. Therefore, under the very high-emission scenario, oceanic and atmospheric drivers are projected to result in a complete collapse of the West Antarctic ice sheet along with significant grounding-line retreat in the marine basins of the East Antarctic ice sheet, leading to a median global mean sea-level rise of 2.75 (6.95) m by 2300 (3000). Under a more sustainable socio-economic pathway, we find that the Antarctic ice sheet may still contribute to a median global mean sea-level rise of 0.62 (1.85) m by 2300 (3000). However, the rate of sea-level rise is significantly reduced as mass loss is likely to remain confined to the Amundsen Sea Embayment, where present-day climate conditions seem sufficient to commit to a continuous retreat of Thwaites Glacier.

A stochastic parameterization of ice sheet surface mass balance for the Stochastic Ice-Sheet and Sea-Level System Model (StISSM v1.0)

Abstract. Many scientific and societal questions that draw on ice sheet modeling necessitate sampling a wide range of potential climatic changes and realizations of internal climate variability. For example, coastal planning literature demonstrates a demand for probabilistic sea level projections with quantified uncertainty. Further, robust attribution of past and future ice sheet change to specific processes or forcings requires a full understanding of the space of possible ice sheet behaviors. The wide sampling required to address such questions is computationally infeasible with sophisticated numerical climate models at the resolution required to accurately force ice sheet models. Stochastic generation of climate forcing of ice sheets offers a complementary alternative. Here, we describe a method to construct a stochastic generator for ice sheet surface mass balance varying in time and space. We demonstrate the method with an application to Greenland Ice Sheet surface mass balance for 1980–2012. We account for spatial correlations among glacier catchments using sparse covariance techniques, and we apply an elevation-dependent downscaling to recover gridded surface mass balance fields suitable for forcing an ice sheet model while including feedback from changing ice sheet surface elevation. The efficiency gained in the stochastic method supports large-ensemble simulations of ice sheet change in a new stochastic ice sheet model. We provide open source Python workflows to support use of our stochastic approach for a broad range of applications.

Abrupt Holocene ice loss due to thinning and ungrounding in the Weddell Sea Embayment

The extent of grounded ice and buttressing by the Ronne Ice Shelf, which provides resistance to the outflow of ice streams, moderate West Antarctic Ice Sheet stability. During the Last Glacial Maximum, the ice sheet advanced and was grounded near the Weddell Sea continental shelf break. The timing of subsequent ice sheet retreat and the relative roles of ice shelf buttressing and grounding line changes remain unresolved. Here we use an ice core record from grounded ice at Skytrain Ice Rise to constrain the timing and speed of early Holocene ice sheet retreat. Measured δ18O and total air content suggest that the surface elevation of Skytrain Ice Rise decreased by about 450 m between 8.2 and 8.0 kyr before 1950 ce (±0.13 kyr). We attribute this elevation change to dynamic thinning due to flow changes induced by the ungrounding of ice in the area. Ice core sodium concentrations suggest that the ice front of this ungrounded ice shelf then retreated about 270 km (±30 km) from 7.7 to 7.3 kyr before 1950 ce. These centennial-scale changes demonstrate how quickly ice mass can be lost from the West Antarctic Ice Sheet due to changes in grounded ice without extensive ice shelf calving. Our findings both support and temporally constrain ice sheet models that exhibit rapid ice loss in the Weddell Sea sector in the early Holocene.

Uncertainties originating from GCM downscaling and bias correction with application to the MIS-11c Greenland Ice Sheet

Abstract. The Marine Isotope Stage 11c (MIS-11c) interglacial is an enigmatic period characterized by a long duration of relatively weak insolation forcing, but it is thought to have been coincident with a large global sea-level rise of 6–13 m. The configuration of the Greenland Ice Sheet during the MIS-11c interglacial highstand is therefore of great interest. Given the constraints of limited data, model-based analysis may be of use but only if model uncertainties are adequately accounted for. A particularly under-addressed issue in coupled climate and ice-sheet modeling is the coupling of surface air temperatures to the ice model. Many studies apply a uniform “lapse rate” accounting for the temperature differences at different altitudes over the ice surface, but this uniformity neglects both regional and seasonal differences in near-surface temperature dependencies on altitude. Herein we provide the first such analysis for MIS-11c Greenland that addresses these uncertainties by comparing one-way coupled Community Earth System Model (CESM) and ice-sheet model results from several different downscaling methodologies. In our study, a spatially and temporally varying temperature downscaling method produced the greatest success rate in matching the constraints of limited paleodata, and it suggests a peak ice volume loss from Greenland during MIS-11c of approximately 50 % compared to present day (∼ 3.9 m contribution to sea-level rise). This result is on the lower bound of existing data- and model-based studies, partly as a consequence of the applied one-way coupling methodology that neglects some feedbacks. Additional uncertainties are examined by comparing two different present-day regional climate analyses for bias correction of temperatures and precipitation, a spread of initialization states and times, and different spatial configurations of precipitation bias corrections. No other factor exhibited greater influence over the simulated Greenland ice sheet than the choice of temperature downscaling scheme.

Validating ensemble historical simulations of Upernavik Isstrøm (1985–2019) using observations of surface velocity and elevation

Abstract The future of tidewater glaciers in response to climate warming is one of the largest sources of uncertainty in the contribution of the Greenland ice sheet to global sea-level rise. In this study, we investigate the ability of an ice-sheet model to reproduce the past evolution of the velocity and surface elevation of a tidewater glacier, Upernavik Isstrøm, by prescribing front positions. To achieve this, we run two ensembles of simulations with a Weertman and a regularised-Coulomb friction law. We show that the ice-flow model has to include a reduction in friction in the first 15 km upstream of the ice front in fast-flowing regions to capture the trends observed during the 1985–2019 period. Without this process, the ensemble model overestimates the ice flow before the retreat of the front in 2005 and does not fully reproduce its acceleration during the retreat. This results in an overestimation of the total mass loss between 1985 and 2019 of 50% (300 vs 200 Gt). Using a variance-based sensitivity analysis, we show that uncertainties in the friction law and the ice-flow law have a greater impact on the model results than surface mass balance and initial surface elevation.

Miocene Antarctic Ice Sheet area adapts significantly faster than volume to CO2-induced climate change

Abstract. The strongly varying benthic δ18O levels of the early and mid-Miocene (23 to 14 Myr ago) are primarily caused by a combination of changes in Antarctic Ice Sheet (AIS) volume and deep-ocean temperatures. These factors are coupled since AIS changes affect deep-ocean temperatures. It has recently been argued that this is due to changes in ice sheet area rather than volume because area changes affect the surface albedo. This finding would be important when the transient AIS grows relatively faster in extent than in thickness, which we test here. We analyse simulations of Miocene AIS variability carried out using the three-dimensional ice sheet model IMAU-ICE forced by warm (high CO2, no ice) and cold (low CO2, large East AIS) climate snapshots. These simulations comprise equilibrium and idealized quasi-orbital transient runs with strongly varying CO2 levels (280 to 840 ppm). Our simulations show a limited direct effect of East AIS changes on Miocene orbital-timescale benthic δ18O variability because of the slow build-up of volume. However, we find that relative to the equilibrium ice sheet size, the AIS area adapts significantly faster and more strongly than volume to the applied forcing variability. Consequently, during certain intervals the ice sheet is receding at the margins, while ice is still building up in the interior. That means the AIS does not adapt to a changing equilibrium size at the same rate or with the same sign everywhere. Our results indicate that the Miocene Antarctic Ice Sheet affects deep-ocean temperatures more than its volume suggests.

Evaluating different geothermal heat-flow maps as basal boundary conditions during spin-up of the Greenland ice sheet

Abstract. There is currently poor scientific agreement on whether the ice–bed interface is frozen or thawed beneath approximately one third of the Greenland ice sheet. This disagreement in basal thermal state results, at least partly, from differences in the subglacial geothermal heat-flow basal boundary condition used in different ice-flow models. Here, we employ seven widely used Greenland geothermal heat-flow maps in 10 000-year spin-ups of the Community Ice Sheet Model (CISM). We perform two spin-ups: one nudged toward thickness observations and the other unconstrained. Across the seven heat-flow maps, and regardless of unconstrained or nudged spin-up, the spread in basal ice temperatures exceeds 10 ∘C over large areas of the ice–bed interface. For a given heat-flow map, the thawed-bed ice-sheet area is consistently larger under unconstrained spin-ups than nudged spin-ups. Under the unconstrained spin-up, thawed-bed area ranges from 33.5 % to 60.0 % across the seven heat-flow maps. Perhaps counterintuitively, the highest iceberg calving fluxes are associated with the lowest heat flows (and vice versa) for both unconstrained and nudged spin-ups. These results highlight the direct, and non-trivial, influence of the heat-flow boundary condition on the simulated equilibrium thermal state of the ice sheet. We suggest that future ice-flow model intercomparisons should employ a range of basal heat-flow maps, and limit direct intercomparisons with simulations using a common heat-flow map.

Accelerating Subglacial Hydrology for Ice Sheet Models With Deep Learning Methods

AbstractSubglacial 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.95) 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.

Relative importance of the mechanisms triggering the Eurasian ice sheet deglaciation in the GRISLI2.0 ice sheet model

Abstract. The last deglaciation (21 to 8 ka) of the Eurasian ice sheet (EIS) is thought to have been responsible for a sea level rise of about 20 m. While many studies have examined the timing and rate of the EIS retreat during this period, many questions remain about the key processes that triggered the EIS deglaciation 21 kyr ago. Due to its large marine-based parts in the Barents–Kara (BKIS) and British Isles sectors, the BKIS is often considered to be a potential analogue of the current West Antarctic ice sheet (WAIS). Identifying the mechanisms that drove the EIS evolution might provide a better understanding of the processes at play in the West Antarctic destabilization. To investigate the relative impact of key drivers on the EIS destabilization, we used the three-dimensional ice sheet model GRISLI (GRenoble Ice Shelf and Land Ice) (version 2.0) forced by climatic fields from five Paleoclimate Modelling Intercomparison Project phases 3 and 4 (PMIP3, PMIP4) Last Glacial Maximum (LGM) simulations. In this study, we performed sensitivity experiments to test the response of the simulated Eurasian ice sheets to surface climate, oceanic temperatures (and thus basal melting under floating ice tongues), and sea level perturbations. Our results highlight that the EIS retreat simulated with the GRISLI model is primarily triggered by atmospheric warming. Increased atmospheric temperatures further amplify the sensitivity of the ice sheets to sub-shelf melting. These results contradict those of previous modelling studies mentioning the central role of basal melting on the deglaciation of the marine-based Barents–Kara ice sheet. However, we argue that the differences with previous works are mainly related to differences in the methodology followed to generate the initial LGM ice sheet. Due to the strong sensitivity of EIS to the atmospheric forcing highlighted with the GRISLI model and the limited extent of the confined ice shelves during the LGM, we conclude by questioning the analogy between EIS and the current WAIS. However, because of the expected rise in atmospheric temperatures, the risk of hydrofracturing is increasing and could ultimately put the WAIS in a configuration similar to the past Eurasian ice sheet.

Net evolution of subglacial sediment transport in the Quebec–Labrador sector of the Laurentide Ice Sheet

The Laurentide Ice Sheet's (LIS) interior had a dynamic polythermal base, but the spatiotemporal variations of subglacial processes related to ice divide migration and other transient changes remain largely unknown, limiting our understanding of regional glacial dynamics. Previous studies focused on the regional glacial landform record, while ice sheet models lacked detailed parameterization within these regions, leading to an overestimation of cold-based subglacial conditions' extent and duration. In this study, glacial sediment dispersal patterns as identified by heavy minerals, clasts, and multivariate statistics of till matrix geochemistry were used to assess ice sheet dynamics within the Quebec–Labrador sector of the LIS. The earliest ice-flow phase produced and transported till across the study area (>175 km). However, major oxide data from till matrix geochemistry show a correlation with underlying bedrock, and this relationship is relatively common in areas of thin till cover and resistant bedrock lithologies. These results suggest a switch from an early phase of widespread erosion and long, sustained sediment transport to one of more limited erosion, perhaps abrasion dominant and shorter transport. Till compositional data and related dispersal patterns add supporting evidence to earlier ice sheet reconstructions based on ice-flow indicators and 10Be data together suggesting a transition from widespread uniform warm-based conditions during the earliest ice flow, followed by the development of an ice divide, its migration, and more sporadic warm-based conditions. Consequently, a thorough understanding of ice-flow history is essential for ice sheet modelling and future mineral exploration programs in inner ice sheet regions of the LIS.

Tidally Modulated Glacial Slip and Tremor at Helheim Glacier, Greenland

AbstractNumerical modeling of ice sheet motion and hence projections of global sea level rise require information about the evolving subglacial environment, which unfortunately remains largely unknown due to its difficulty of access. Here we advance such subglacial observations by reporting multi‐year observations of seismic tremor likely associated with glacier sliding at Helheim Glacier. This association is confirmed by correlation analysis between tremor power and multiple environmental forcings on different timescales. Variations of the observed tremor power indicate that different factors affect glacial sliding on different timescales. Effective pressure may control glacial sliding on long (seasonal/annual) timescales, while tidal forcing modulates the sliding rate and tremor power on short (hourly/daily) timescales. Polarization results suggest that the tremor source comes from an upstream subglacial ridge. This observation provides insights on how different factors should be included in ice sheet modeling and how their timescales of variability play an essential role.