An official website of the United States government
Abstract. The effect of the North Atlantic Ocean on the Greenland Ice Sheet through submarine melting of Greenland's tidewater glacier calving fronts is thought to be a key driver of widespread glacier retreat, dynamic mass loss and sea level contribution from the ice sheet. Despite its critical importance, problems of process complexity and scale hinder efforts to represent the influence of submarine melting in ice-sheet-scale models. Here we propose parameterizing tidewater glacier terminus position as a simple linear function of submarine melting, with submarine melting in turn estimated as a function of subglacial discharge and ocean temperature. The relationship is tested, calibrated and validated using datasets of terminus position, subglacial discharge and ocean temperature covering the full ice sheet and surrounding ocean from the period 1960–2018. We demonstrate a statistically significant link between multi-decadal tidewater glacier terminus position change and submarine melting and show that the proposed parameterization has predictive power when considering a population of glaciers. An illustrative 21st century projection is considered, suggesting that tidewater glaciers in Greenland will undergo little further retreat in a low-emission RCP2.6 scenario. In contrast, a high-emission RCP8.5 scenario results in a median retreat of 4.2 km, with a quarter of tidewater glaciers experiencing retreat exceeding 10 km. Our study provides a long-term and ice-sheet-wide assessment of the sensitivity of tidewater glaciers to submarine melting and proposes a practical and empirically validated means of incorporating ocean forcing into models of the Greenland ice sheet.
more »
« less
Submarine melting of glaciers in Greenland amplified by atmospheric warming
Abstract Rapid ice loss from the Greenland ice sheet since 1992 is due in equal parts to increased surface melting and accelerated ice flow. The latter is conventionally attributed to ocean warming, which has enhanced submarine melting of the fronts of Greenland’s marine-terminating glaciers. Yet, through the release of ice sheet surface meltwater into the ocean, which excites near-glacier ocean circulation and in turn the transfer of heat from ocean to ice, a warming atmosphere can increase submarine melting even in the absence of ocean warming. The relative importance of atmospheric and oceanic warming in driving increased submarine melting has, however, not been quantified. Here, we reconstruct the rate of submarine melting at Greenland’s marine-terminating glaciers from 1979 to 2018 and estimate the resulting dynamic mass loss. We show that in south Greenland, variability in submarine melting was indeed governed by the ocean, but, in contrast, the atmosphere dominated in the northwest. At the ice sheet scale, the atmosphere plays a first-order role in controlling submarine melting and the subsequent dynamic mass loss. Our results challenge the attribution of dynamic mass loss to ocean warming alone and show that a warming atmosphere has amplified the impact of the ocean on the Greenland ice sheet.
more »
« less
- PAR ID:
- 10376019
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Geoscience
- Volume:
- 15
- Issue:
- 10
- ISSN:
- 1752-0894
- Page Range / eLocation ID:
- p. 794-799
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Submarine melting at Greenland's marine terminating glaciers is a crucial, yet poorly constrained process in the coupled ice‐ocean system. Application of Antarctic melt rate representations, derived for floating glacier tongues, to non‐floating marine terminating glaciers commonly found in Greenland, results in a dramatic underestimation of submarine melting. Here, we revisit the physical theory underlying melt rate parameterizations and leverage recently published observational data to derive a novel melt rate parameterization. This is the first parameterization that (a) consistently comprises both convective‐ and shear‐dominated melt regimes, (b) includes coefficients quantitatively constrained using observational data, and (c) is applicable to any vertical glacier front. We show that, compared to the current state‐of‐the‐art approach, the scheme provides an improved fit to observed melt rates on the scale of the terminating front, offering an opportunity to incorporate this critical missing forcing into ocean circulation models.more » « less
-
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
-
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
-
null (Ed.)Abstract Understanding marine-terminating ice sheet response to past climate transitions provides valuable long-term context for observations of modern ice sheet change. Here, we reconstruct the last deglaciation of marine-terminating Cordilleran Ice Sheet (CIS) margins in Southeast Alaska and explore potential forcings of western CIS retreat. We combine 27 new cosmogenic 10 Be exposure ages, 13 recently published 10 Be ages, and 25 new 14 C ages from raised marine sediments to constrain CIS recession. Retreat from the outer coast was underway by 17 ka, and the inner fjords and sounds were ice-free by 15 ka. After 15 ka, the western margin of the CIS became primarily land-terminating and alpine glaciers disappeared from the outer coast. Isolated alpine glaciers may have persisted in high inland peaks until the early Holocene. Our results suggest that the most rapid phase of CIS retreat along the Pacific coast occurred between ~17 and 15 ka. This retreat was likely driven by processes operating at the ice-ocean interface, including sea level rise and ocean warming. CIS recession after ~15 ka occurred during a time of climatic amelioration in this region, when both ocean and air temperatures increased. These data highlight the sensitivity of marine-terminating CIS regions to deglacial climate change.more » « less
