An official website of the United States government
Abstract. The northern sector of the Greenland Ice Sheet is considered to beparticularly susceptible to ice mass loss arising from increased glacierdischarge in the coming decades. However, the past extent and dynamics ofoutlet glaciers in this region, and hence their vulnerability to climatechange, are poorly documented. In the summer of 2019, the Swedish icebreakerOden entered the previously unchartered waters of Sherard Osborn Fjord, whereRyder Glacier drains approximately 2 % of Greenland's ice sheet into theLincoln Sea. Here we reconstruct the Holocene dynamics of Ryder Glacier andits ice tongue by combining radiocarbon dating with sedimentary faciesanalyses along a 45 km transect of marine sediment cores collected betweenthe modern ice tongue margin and the mouth of the fjord. The resultsillustrate that Ryder Glacier retreated from a grounded position at thefjord mouth during the Early Holocene (> 10.7±0.4 ka cal BP) and receded more than 120 km to the end of Sherard Osborn Fjord by theMiddle Holocene (6.3±0.3 ka cal BP), likely becoming completelyland-based. A re-advance of Ryder Glacier occurred in the Late Holocene,becoming marine-based around 3.9±0.4 ka cal BP. An ice tongue,similar in extent to its current position was established in the LateHolocene (between 3.6±0.4 and 2.9±0.4 ka cal BP) andextended to its maximum historical position near the fjord mouth around 0.9±0.3 ka cal BP. Laminated, clast-poor sediments were deposited duringthe entire retreat and regrowth phases, suggesting the persistence of an icetongue that only collapsed when the glacier retreated behind a prominenttopographic high at the landward end of the fjord. Sherard Osborn Fjordnarrows inland, is constrained by steep-sided cliffs, contains a number ofbathymetric pinning points that also shield the modern ice tongue andgrounding zone from warm Atlantic waters, and has a shallowing inlandsub-ice topography. These features are conducive to glacier stability andcan explain the persistence of Ryder's ice tongue while the glacier remainedmarine-based. However, the physiography of the fjord did not halt thedramatic retreat of Ryder Glacier under the relatively mild changes inclimate forcing during the Holocene. Presently, Ryder Glacier is groundedmore than 40 km seaward of its inferred position during the Middle Holocene,highlighting the potential for substantial retreat in response to ongoingclimate change.
more »
« less
This content will become publicly available on March 20, 2026
Simulating the Holocene evolution of Ryder Glacier, North Greenland
Abstract. The Greenland Ice Sheet's negative mass balance is driven by a sensitivity to both a warming atmosphere and ocean. The fidelity of ice-sheet models in accounting for ice-ocean interaction is inherently uncertain and often constrained against recent fluctuations in the ice-sheet margin from the previous decades. The geological record can be utilised to contextualise ice-sheet mass loss and understand the drivers of changes at the marine margin across climatic shifts and previous extended warm periods, aiding our understanding of future ice-sheet behaviour. Here, we use the Ice-sheet and Sea-level System Model (ISSM) to explore the Holocene evolution of Ryder Glacier draining into Sherard Osborn Fjord, Northern Greenland. Our modelling results are constrained with terrestrial reconstructions of the paleo-ice sheet margin and an extensive marine sediment record from Sherard Osborn Fjord that details ice dynamics over the past 12.5 ka years. By employing a consistent mesh resolution of <1 km at the ice-ocean boundary, we assess the importance of atmospheric and oceanic changes to Ryder Glacier's Holocene behaviour. Our simulations show that the initial retreat of the ice margin after the Younger Dryas cold period was driven by a warming climate and the resulting fluctuations in Surface Mass Balance. Changing atmospheric conditions remain the first order control in the timing of ice retreat during the Holocene. We find ice-ocean interactions become increasingly fundamental to Ryder's retreat in the mid-Holocene; with higher than contemporary melt rates required to force grounding line retreat and capture the collapse of the ice tongue during the Holocene Thermal Maximum. Regrowth of the tongue during the neo-glacial cooling of the late Holocene is necessary to advance both the terrestrial and marine margins of the glacier. Our results stress the importance of accurately resolving the ice-ocean interface in modelling efforts over centennial and millennial time scales, in particular the role of floating ice tongues and submarine melt, and provide vital analogous for the future evolution of Ryder in a warming climate.
more »
« less
- Award ID(s):
- 2105960
- PAR ID:
- 10628549
- Publisher / Repository:
- The cryosphere
- Date Published:
- Format(s):
- Medium: X
- Institution:
- UCLA Joint institute for regional earth system science and engineering
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Numerical simulations of the Greenland Ice Sheet (GrIS) over geologictimescales can greatly improve our knowledge of the critical factors drivingGrIS demise during climatically warm periods, which has clear relevance forbetter predicting GrIS behavior over the upcoming centuries. To assess thefidelity of these modeling efforts, however, observational constraints ofpast ice sheet change are needed. Across southwestern Greenland, geologicrecords detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorouslybenchmark model simulations against geologic reconstructions of ice sheetchange. Here, we present regional ice sheet modeling results using theIce-sheet and Sea-level System Model (ISSM) of Holocene ice sheet historyacross an extensive fjord region in southwestern Greenland covering thelandscape around the Kangiata Nunaata Sermia (KNS) glacier and extendingoutward along the 200 km Nuup Kangerula (Godthåbsfjord). Oursimulations, forced by reconstructions of Holocene climate and recentlyimplemented calving laws, assess the sensitivity of ice retreat across theKNS region to atmospheric and oceanic forcing. Our simulations reveal thatthe geologically reconstructed ice retreat across the terrestrial landscapein the study area was likely driven by fluctuations in surface mass balancein response to Early Holocene warming – and was likely not influencedsignificantly by the response of adjacent outlet glaciers to calving andocean-induced melting. The impact of ice calving within fjords, however,plays a significant role by enhancing ice discharge at the terminus, leadingto interior thinning up to the ice divide that is consistent withreconstructed magnitudes of Early Holocene ice thinning. Our results,benchmarked against geologic constraints of past ice-margin change, suggestthat while calving did not strongly influence Holocene ice-margin migrationacross terrestrial portions of the KNS forefield, it strongly impactedregional mass loss. While these results imply that the implementation andresolution of ice calving in paleo-ice-flow models is important towardsmaking more robust estimations of past ice mass change, they also illustratethe importance these processes have on contemporary and future long-term icemass change across similar fjord-dominated regions of the GrIS.more » « less
-
Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.more » « less
-
Ocean temperatures have warmed in the fjords surrounding the Greenland Ice Sheet, causing increased melt along their ice fronts and rapid glacier retreat and contributing to rising global sea levels. However, there are many physical mechanisms that can mediate the glacier response to ocean warming and variability. Warm ocean waters can directly cause melt at horizontal and vertical ice interfaces or promote iceberg calving by weakening proglacial melange or undercutting the glacier front. Sermeq Kujalleq (also known as Jakobshavn Isbræ) is the largest and fastest glacier in Greenland and has undergone substantial retreat, which started in the late 1990s. In this study, we use an ensemble modeling approach to disentangle the dominant mechanisms that drive the retreat of Sermeq Kujalleq. Within this ensemble, we vary the sensitivity of three different glaciological parameters to ocean temperature: frontal melt, subshelf melt, and a calving-stress threshold. Comparing results to the observed retreat behavior from 1985 to 2018, we select a best-fitting simulation which reproduces the observed retreat well. In this simulation, the arrival of warm water at the front of Sermeq Kujalleq in the late 1990s led to enhanced rates of subshelf melt, triggering the disintegration of the floating ice tongue over a decade. The recession of the calving front into a substantially deeper bed trough around 2010 accelerated the calving-driven retreat, which continued nearly unabated despite local ocean cooling in 2016. An extended ensemble of simulations with varying calving thresholds shows evidence of hysteresis in the calving rate, which can only be inhibited by a substantial increase in the calving-stress threshold beyond the values suggested for the historical period. Our findings indicate that accurate simulation of rapid calving-driven glacier retreats requires more sophisticated models of iceberg mélange and calving evolution coupled to ice flow models.more » « less
-
Abstract Variability in oceanic conditions directly impacts ice loss from marine outlet glaciers in Greenland, influencing the ice sheet mass balance. Oceanic conditions are available from Atmosphere‐Ocean Global Climate Model (AOGCM) output, but these models require extensive computational resources and lack the fine resolution needed to simulate ocean dynamics on the Greenland continental shelf and close to glacier marine termini. Here, we develop a statistical approach to generate ocean forcing for ice sheet model simulations, which incorporates natural spatiotemporal variability and anthropogenic changes. Starting from raw AOGCM ocean heat content, we apply: (a) a bias‐correction using ocean reanalysis, (b) an extrapolation accounting for on‐shelf ocean dynamics, and (c) stochastic time series models to generate realizations of natural variability. The bias‐correction reduces model errors by ∼25% when compared to independent in‐situ measurements. The bias‐corrected time series are subsequently extrapolated to fjord mouth locations using relations constrained from available high‐resolution regional ocean model results. The stochastic time series models reproduce the spatial correlation, characteristic timescales, and the amplitude of natural variability of bias‐corrected AOGCMs, but at negligible computational expense. We demonstrate the efficiency of this method by generating >6,000 time series of ocean forcing for >200 Greenland marine‐terminating glacier locations until 2100. As our method is computationally efficient and adaptable to any ocean model output and reanalysis product, it provides flexibility in exploring sensitivity to ocean conditions in Greenland ice sheet model simulations. We provide the output and workflow in an open‐source repository, and discuss advantages and future developments for our method.more » « less
