scholarly journals Greenland Ice Sheet runoff reduced by meltwater refreezing in bare ice

Author(s):  
Matthew Cooper ◽  
Laurence Smith ◽  
Åsa Rennermalm ◽  
Kang Yang ◽  
Glen Liston ◽  
...  

Abstract The Greenland Ice Sheet’s contribution to global sea-level rise is accelerating1 due to increased melting of its bare-ice ablation zone2–6, but there is growing evidence that climate models overestimate runoff from this critical area of the ice sheet7–12. Current climate models assume all bare ice runoff escapes to the ocean, unlike snow covered areas where some fraction of runoff is retained and/or refrozen in porous firn13–15. Here we use in situ measurements and numerical modeling to reveal extensive retention and refreezing of liquid meltwater in bare glacial ice, explaining chronic runoff overestimation by climate models. From 2009–2018, refreezing of liquid meltwater in bare, porous glacial ice reduced meltwater runoff by 11–23 Gt a-1 in southwest Greenland alone, equivalent to 10–20% of annual meltwater production. This mass retention is commensurate with current estimates of climate model ice sheet meltwater runoff uncertainty, and may represent an overlooked buffer on projected runoff increases for the coming century16. Inclusion of bare-ice retention and refreezing processes in climate models therefore has immediate potential to improve forecasts of ice sheet runoff and its contribution to global sea-level rise.

2021 ◽  
Author(s):  
Max Brils ◽  
Peter Kuipers Munneke ◽  
Willem Jan van de Berg ◽  
Achim Heilig ◽  
Baptiste Vandercrux ◽  
...  

<p>Recent studies indicate that a declining surface mass balance will dominate the Greenland Ice Sheet’s (GrIS) contribution to 21<sup>st</sup> century sea level rise. It is therefore crucial to understand the liquid water balance of the ice sheet and its response to increasing temperatures and surface melt if we want to accurately predict future sea level rise. The ice sheet firn layer covers ~90% of the GrIS and provides pore space for storage and refreezing of meltwater. Because of this, the firn layer can retain up to ~45% of the surface meltwater and thus act as an efficient buffer to ice sheet mass loss. However, in a warming climate this buffer capacity of the firn layer is expected to decrease, amplifying meltwater runoff and sea-level rise. Dedicated firn models are used to understand how firn layers evolve and affect runoff. Additionally, firn models are used to estimate the changing thickness of the firn layer, which is necessary in altimetry to convert surface height change into ice sheet mass loss.</p><p>Here, we present the latest version of our firn model IMAU-FDM. With respect to the previous version, changes have been made to the handling of the freshly fallen snow, the densification rate of the firn and the conduction of heat. These changes lead to an improved representation of firn density and temperature. The results have been thoroughly validated using an extensive dataset of density and temperature measurements that we have compiled covering 126 different locations on the GrIS. Meltwater behaviour in the model is validated with upward-looking GPR measurements at Dye-2. Lastly, we present an in-depth look at the evolution firn characteristics at some typical locations in Greenland.</p><p>Dedicated, stand-alone firn models offer various benefits to using a regional climate model with an embedded firn model. Firstly, the vertical resolution for buried snow and ice layers can be larger, improving accuracy. Secondly, a stand-alone firn model allows for spinning up the model to a more accurate equilibrium state. And thirdly, a stand-alone model is more cost- and time-effective to use. Firn models are increasingly capable of simulating the firn layer, but areas with large amounts of melt still pose the greatest challenge.</p>


2020 ◽  
Author(s):  
Michiel van den Broeke ◽  
Brice Noël ◽  
Leo van Kampenhout ◽  
Willem-Jan van de Berg

<p>The mass balance of the Greenland ice sheet (GrIS, units Gt per year) equals the surface mass balance (SMB) minus solid ice discharge across the grounding line. As the latter is definite positive, an important threshold for irreversible GrIS mass loss occurs when long-term average SMB becomes negative. For this to happen, runoff (mainly meltwater, some rain) must exceed mass accumulation (mainly snowfall minus sublimation). Even for a single year, this threshold has not been passed since at least 1958, the first year with reliable estimates of SMB components, although recent years with warm summers (e.g. 2012 and 2019) came close. Simply extrapolating the recent (1991-present) negative SMB trend into the future suggests that the SMB = 0 threshold could be reached before ~2040, but such predictions are extremely uncertain given the very large interannual SMB variability, the relative brevity of the time series and the uncertainty in future warming. In this study we use a cascade of models, extensively evaluated with in-situ and remotely sensed (GRACE) SMB observations, to better constrain the future regional warming threshold for the 5-year average GrIS SMB to become negative. To this end, a 1950-2100 climate change run with the global model CESM2 (app. 100 km resolution) was dynamically downscaled using the regional climate model RACMO2 (app. 11 km), which in turn was statistically downscaled to 1 km resolution. The result is a threshold regional Greenland warming of close to 4 degrees. We then use a range of CMIP5 and CMIP6 global climate models to translate the regional value into a global warming threshold for various warming scenarios, including its timing this century. We find substantial differences, ranging from stabilization before the threshold is reached in the RCP/SSP2.6 scenarios with a limited but still significant sea-level rise contribution (< 5 cm by 2100) to an imminent crossing of the warming threshold for the RCP/SSP8.5 scenarios with substantial and ever-growing contributions to sea level rise (> 10 cm by 2100). These results stress the need for strong mitigation to avoid irreversible GrIS mass loss. We finish by discussing the caveats and uncertainties of our approach.</p>


2013 ◽  
Vol 9 (1) ◽  
pp. 353-366 ◽  
Author(s):  
A. Quiquet ◽  
C. Ritz ◽  
H. J. Punge ◽  
D. Salas y Mélia

Abstract. As pointed out by the forth assessment report of the Intergovernmental Panel on Climate Change, IPCC-AR4 (Meehl et al., 2007), the contribution of the two major ice sheets, Antarctica and Greenland, to global sea level rise, is a subject of key importance for the scientific community. By the end of the next century, a 3–5 °C warming is expected in Greenland. Similar temperatures in this region were reached during the last interglacial (LIG) period, 130–115 ka BP, due to a change in orbital configuration rather than to an anthropogenic forcing. Ice core evidence suggests that the Greenland ice sheet (GIS) survived this warm period, but great uncertainties remain about the total Greenland ice reduction during the LIG. Here we perform long-term simulations of the GIS using an improved ice sheet model. Both the methodologies chosen to reconstruct palaeoclimate and to calibrate the model are strongly based on proxy data. We suggest a relatively low contribution to LIG sea level rise from Greenland melting, ranging from 0.7 to 1.5 m of sea level equivalent, contrasting with previous studies. Our results suggest an important contribution of the Antarctic ice sheet to the LIG highstand.


2013 ◽  
Vol 38 (1) ◽  
pp. 19-54 ◽  
Author(s):  
Vena W. Chu

Understanding Greenland ice sheet (GrIS) hydrology is essential for evaluating response of ice dynamics to a warming climate and future contributions to global sea level rise. Recently observed increases in temperature and melt extent over the GrIS have prompted numerous remote sensing, modeling, and field studies gauging the response of the ice sheet and outlet glaciers to increasing meltwater input, providing a quickly growing body of literature describing seasonal and annual development of the GrIS hydrologic system. This system is characterized by supraglacial streams and lakes that drain through moulins, providing an influx of meltwater into englacial and subglacial environments that increases basal sliding speeds of outlet glaciers in the short term. However, englacial and subglacial drainage systems may adjust to efficiently drain increased meltwater without significant changes to ice dynamics over seasonal and annual scales. Both proglacial rivers originating from land-terminating glaciers and subglacial conduits under marine-terminating glaciers represent direct meltwater outputs in the form of fjord sediment plumes, visible in remotely sensed imagery. This review provides the current state of knowledge on GrIS surface water hydrology, following ice sheet surface meltwater production and transport via supra-, en-, sub-, and proglacial processes to final meltwater export to the ocean. With continued efforts targeting both process-level and systems analysis of the hydrologic system, the larger picture of how future changes in Greenland hydrology will affect ice sheet glacier dynamics and ultimately global sea level rise can be advanced.


2016 ◽  
Vol 10 (3) ◽  
pp. 1147-1159 ◽  
Author(s):  
Andreas Bech Mikkelsen ◽  
Alun Hubbard ◽  
Mike MacFerrin ◽  
Jason Eric Box ◽  
Sam H. Doyle ◽  
...  

Abstract. It has been argued that the infiltration and retention of meltwater within firn across the percolation zone of the Greenland ice sheet has the potential to buffer up to  ∼  3.6 mm of global sea-level rise (Harper et al., 2012). Despite evidence confirming active refreezing processes above the equilibrium line, their impact on runoff and proglacial discharge has yet to be assessed. Here, we compare meteorological, melt, firn stratigraphy and discharge data from the extreme 2010 and 2012 summers to determine the relationship between atmospheric forcing and melt runoff at the land-terminating Kangerlussuaq sector of the Greenland ice sheet, which drains into the Watson River. The 6.8 km3 bulk discharge in 2012 exceeded that in 2010 by 28 %, despite only a 3 % difference in net incoming melt energy between the two years. This large disparity can be explained by a 10 % contribution of runoff originating from above the long-term equilibrium line in 2012 caused by diminished firn retention. The amplified 2012 response was compounded by catchment hypsometry; the disproportionate increase in area contributing to runoff as the melt-level rose high into the accumulation area.Satellite imagery and aerial photographs reveal an extensive supraglacial network extending 140 km from the ice margin that confirms active meltwater runoff originating well above the equilibrium line. This runoff culminated in three days with record discharge of 3100 m3 s−1 (0.27 Gt d−1) that peaked on 11 July and washed out the Watson River Bridge. Our findings corroborate melt infiltration processes in the percolation zone, though the resulting patterns of refreezing are complex and can lead to spatially extensive, perched superimposed ice layers within the firn. In 2012, such layers extended to an elevation of at least 1840 m and provided a semi-impermeable barrier to further meltwater storage, thereby promoting widespread runoff from the accumulation area of the Greenland ice sheet that contributed directly to proglacial discharge and global sea-level rise.


2015 ◽  
Vol 11 (3) ◽  
pp. 403-424 ◽  
Author(s):  
A. M. Dolan ◽  
S. J. Hunter ◽  
D. J. Hill ◽  
A. M. Haywood ◽  
S. J. Koenig ◽  
...  

Abstract. During an interval of the Late Pliocene, referred to here as the mid-Pliocene Warm Period (mPWP; 3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the end of this century, and atmospheric carbon dioxide concentrations were higher than pre-industrial levels. Sea level was also higher than today, implying a significant reduction in the extent of the ice sheets. Thus, the mPWP provides a natural laboratory in which to investigate the long-term response of the Earth's ice sheets and sea level in a warmer-than-present-day world. At present, our understanding of the Greenland ice sheet during the mPWP is generally based upon predictions using single climate and ice sheet models. Therefore, it is essential that the model dependency of these results is assessed. The Pliocene Model Intercomparison Project (PlioMIP) has brought together nine international modelling groups to simulate the warm climate of the Pliocene. Here we use the climatological fields derived from the results of the 15 PlioMIP climate models to force an offline ice sheet model. We show that mPWP ice sheet reconstructions are highly dependent upon the forcing climatology used, with Greenland reconstructions ranging from an ice-free state to a near-modern ice sheet. An analysis of the surface albedo variability between the climate models over Greenland offers insights into the drivers of inter-model differences. As we demonstrate that the climate model dependency of our results is high, we highlight the necessity of data-based constraints of ice extent in developing our understanding of the mPWP Greenland ice sheet.


2020 ◽  
Author(s):  
Andrew Shepherd ◽  

<p>In recent decades, the Antarctic and Greenland Ice Sheets have been major contributors to global sea-level rise and are expected to be so in the future. Although increases in glacier flow and surface melting have been driven by oceanic and atmospheric warming, the degree and trajectory of today’s imbalance remain uncertain. Here we compare and combine 26 individual satellite records of changes in polar ice sheet volume, flow and gravitational potential to produce a reconciled estimate of their mass balance. <strong>Since the early 1990’s, ice losses from Antarctica and Greenland have caused global sea-levels to rise by 18.4 millimetres, on average, and there has been a sixfold increase in the volume of ice loss over time. Of this total, 41 % (7.6 millimetres) originates from Antarctica and 59 % (10.8 millimetres) is from Greenland. In this presentation, we compare our reconciled estimates of Antarctic and Greenland ice sheet mass change to IPCC projection of sea level rise to assess the model skill in predicting changes in ice dynamics and surface mass balance.  </strong>Cumulative ice losses from both ice sheets have been close to the IPCC’s predicted rates for their high-end climate warming scenario, which forecast an additional 170 millimetres of global sea-level rise by 2100 when compared to their central estimate.</p>


2020 ◽  
Author(s):  
Samuel Helsen ◽  
Sam Vanden Broucke ◽  
Alexandra Gossart ◽  
Niels Souverijns ◽  
Nicole van Lipzig

<p>The Totten glacier is a highly dynamic outlet glacier, situated in E-Antarctica, that contains a potential sea level rise of about 3.5 meters. During recent years, this area has been influenced by sub-shelf intrusion of warm ocean currents, contributing to higher basal melt rates. Moreover, most of the ice over this area is grounded below sea level, which makes the ice shelf potentially vulnerable to the marine ice sheet instability mechanism. It is expected that, as a result of climate change, the latter mechanisms may contribute to significant ice losses in this region within the next decades, thereby contributing to future sea level rise. Up to now, most studies have been focusing on sub-shelf melt rates and the influence of the ocean, with much less attention for atmospheric processes (often ignored), which also play a key-role in determining the climatic conditions over this region. For example: surface melt is important because it contributes to hydrofracturing, a process that may lead to ice cliff instabilities. Also precipitation is an important atmospheric process, since it determines the input of mass to the ice sheet and contributes directly to the surface mass balance. In order to perform detailed studies on these processes, we need a well-evaluated climate model that represents all these processes well. Recently, the COSMO-CLM<sup>2</sup> (CCLM<sup>2</sup>) model was adapted to the climatological conditions over Antarctica. The model was evaluated by comparing a 30 year Antarctic-wide hindcast run (1986-2016) at 25 km resolution with meteorological observational products (Souverijns et al., 2019). It was shown that the model performance is comparable to other state-of-the-art regional climate models over the Antarctic region. We now applied the CCLM<sup>2</sup> model in a regional configuration over the Totten glacier area (E-Antarctica) at 5 km resolution and evaluated its performance over this region by comparing it to climatological observations from different stations. We show that the performance for temperature in the high resolution run is comparable to the performance of the Antarctic-wide run. Precipitation is, however, overestimated in the high-resolution run, especially over dome structures (Law-Dome). Therefore, we applied an orographic smoothening, which clearly improves the precipitation pattern with respect to observations. Wind speed is overestimated in some places, which is solved by increasing the surface roughness. This research frames in the context of the PARAMOUR project. Within PARAMOUR, CCLM<sup>2 </sup>is currently being coupled to an ocean model (NEMO) and an ice sheet model (f.ETISh/BISICLES) in order to understand decadal predictability over this region.</p>


2020 ◽  
Author(s):  
Xavier Fettweis ◽  

<p>The Greenland Ice Sheet (GrIS) mass loss has been accelerating at a rate of about 20 +/- 10 Gt/yr<sup>2</sup> since the end of the 1990's, with around 60% of this mass loss directly attributed to enhanced surface meltwater runoff. However, in the climate and glaciology communities, different approaches exist on how to model the different surface mass balance (SMB) components using: (1) complex physically-based climate models which are computationally expensive; (2) intermediate complexity energy balance models; (3) simple and fast positive degree day models which base their inferences on statistical principles and are computationally highly efficient. Additionally, many of these models compute the SMB components based on different spatial and temporal resolutions, with different forcing fields as well as different ice sheet topographies and extents, making inter-comparison difficult. In the GrIS SMB model intercomparison project (GrSMBMIP) we address these issues by forcing each model with the same data (i.e., the ERA-Interim reanalysis) except for two global models for which this forcing is limited to the oceanic conditions, and at the same time by interpolating all modelled results onto a common ice sheet mask at 1 km horizontal resolution for the common period 1980-2012. The SMB outputs from 13 models are then compared over the GrIS to (1) SMB estimates using a combination of gravimetric remote sensing data from GRACE and measured ice discharge, (2) ice cores, snow pits, in-situ SMB observations, and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Our results reveal that the mean GrIS SMB of all 13 models has been positive between 1980 and 2012 with an average of 340 +/- 112 Gt/yr, but has decreased at an average rate of -7.3 Gt/yr<sup>2</sup> (with a significance of 96%), mainly driven by an increase of 8.0 Gt/yr<sup>2</sup> (with a significance of 98%) in meltwater runoff. Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting the need for accurate representation of the GrIS ablation zone extent and processes driving the surface melt. In addition, a higher density of in-situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 mWE/yr due to large discrepancies in modelled snowfall accumulation. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of same order than RCMs with observations and remain then useful tools for long-term simulations. It is also interesting to note that the ensemble mean of the 13 models produces the best estimate of the present day SMB relative to observations, suggesting that biases are not systematic among models. Finally, results from MAR forced by ERA5 will be added in this intercomparison to evaluate the added value of using this new reanalysis as forcing vs the former ERA-Interim reanalysis (used in SMBMIP). </p>


2020 ◽  
Vol 11 (4) ◽  
pp. 953-976
Author(s):  
Jonas Van Breedam ◽  
Heiko Goelzer ◽  
Philippe Huybrechts

Abstract. The emphasis for informing policy makers on future sea-level rise has been on projections by the end of the 21st century. However, due to the long lifetime of atmospheric CO2, the thermal inertia of the climate system and the slow equilibration of the ice sheets, global sea level will continue to rise on a multi-millennial timescale even when anthropogenic CO2 emissions cease completely during the coming decades to centuries. Here we present global sea-level change projections due to the melting of land ice combined with steric sea effects during the next 10 000 years calculated in a fully interactive way with the Earth system model of intermediate complexity LOVECLIMv1.3. The greenhouse forcing is based on the Extended Concentration Pathways defined until 2300 CE with no carbon dioxide emissions thereafter, equivalent to a cumulative CO2 release of between 460 and 5300 GtC. We performed one additional experiment for the highest-forcing scenario with the inclusion of a methane emission feedback where methane is slowly released due to a strong increase in surface and oceanic temperatures. After 10 000 years, the sea-level change rate drops below 0.05 m per century and a semi-equilibrated state is reached. The Greenland ice sheet is found to nearly disappear for all forcing scenarios. The Antarctic ice sheet contributes only about 1.6 m to sea level for the lowest forcing scenario with a limited retreat of the grounding line in West Antarctica. For the higher-forcing scenarios, the marine basins of the East Antarctic Ice Sheet also become ice free, resulting in a sea-level rise of up to 27 m. The global mean sea-level change after 10 000 years ranges from 9.2 to more than 37 m. For the highest-forcing scenario, the model uncertainty does not exclude the complete melting of the Antarctic ice sheet during the next 10 000 years.


Sign in / Sign up

Export Citation Format

Share Document