scholarly journals Simulating the growth of supra-glacial lakes at the western margin of the Greenland ice sheet

2012 ◽  
Vol 6 (2) ◽  
pp. 1307-1336
Author(s):  
A. A. Leeson ◽  
A. Shepherd ◽  
S. Palmer ◽  
A. Sundal ◽  
X. Fettweis

Abstract. We present a method of modelling the growth of supra-glacial lakes at the western margin of the Greenland ice sheet, based on routeing runoff estimated by a Regional Climate Model (RCM) across a digital elevation model (DEM) of the ice sheet surface. Using data acquired during the 2003 melt season, we demonstrate that the model is 18 times more likely to correctly predict the presence or absence of lakes identified in MODIS satellite imagery within an elevation range of 1000 to 1600 metres above sea level (m.a.s.l.) than it is to make incorrect predictions. Our model does not, however, simulate processes leading to lake stagnation or decay, such as refreezing or drainage – a process which affects approximately 17% of lakes in our study area (Selmes et al., 2011). This likely explains much of why our model over-predicts cumulative area by 32% although other factors including uncertainty in the DEM and in the MODIS derived observations used for validation contribute to this error. Simulated lake filling tends to lead observations by approximately 5 days which could be related to a filling period required to saturate cracks, crevasses and other porous space within the ice. We find that the maximum modelled lake covered ice sheet area is 6% and suggest that this is a topographic limitation for this sector. We can take this as an upper bound; given the absence of drainage in the model. In 2003, the difference between RCM estimates of runoff and the maximum volume of water simulated to be stored in lakes was 12.49 km3. This can be taken as a measure of potential water available for lubrication and is calculated to be 1.86 m3 per square metre of ice. This study has proved a good first step towards capturing the variability of supra-glacial lake evolution with a numerical model; we are optimistic that the model will develop further into a useful tool for use in analysing the behaviour of supra-glacial lakes on the Greenland ice sheet in the present day and beyond.

2012 ◽  
Vol 6 (5) ◽  
pp. 1077-1086 ◽  
Author(s):  
A. A. Leeson ◽  
A. Shepherd ◽  
S. Palmer ◽  
A. Sundal ◽  
X. Fettweis

Abstract. We present a new method of modelling the growth of supraglacial lakes at the western margin of the Greenland ice sheet, based on routing runoff estimated by a regional climate model across a digital elevation model (DEM) of the ice sheet surface. Using data acquired during the 2003 melt season, we demonstrate that the model is 19 times more likely to correctly predict the presence (or absence) of lakes than it is to make incorrect predictions, within an elevation range of 1100 to 1700 metres above sea level (m a.s.l.), when compared with MODIS satellite imagery. Of the 66% of observed lake locations which the model correctly reproduces, the simulated lake onset day is found to be correlated with that observed with a Pearson correlation coefficient of 0.76. Our model accurately simulates maximum cumulative lake area with only a 1.5% overestimate. However, because our model does not simulate processes leading to lake stagnation or decay, such as refreezing or drainage, at present we do not simulate absolute daily lake area. We find that the maximum potential lake-covered ice sheet area is limited by topography to 6.4%. We estimate that this corresponds to a volume of 1.49 km3, 12% of the runoff produced in 2003. This can be taken as an upper bound given uncertainty in the DEM. This study has proved a good first step towards capturing the variability of supraglacial lake evolution with a numerical model. These initial results are promising and suggest that the model is a useful tool for use in analysing the behaviour of supraglacial lakes on the Greenland ice sheet in the present day and potentially beyond.


2016 ◽  
Vol 10 (5) ◽  
pp. 2361-2377 ◽  
Author(s):  
Brice Noël ◽  
Willem Jan van de Berg ◽  
Horst Machguth ◽  
Stef Lhermitte ◽  
Ian Howat ◽  
...  

Abstract. This study presents a data set of daily, 1 km resolution Greenland ice sheet (GrIS) surface mass balance (SMB) covering the period 1958–2015. Applying corrections for elevation, bare ice albedo and accumulation bias, the high-resolution product is statistically downscaled from the native daily output of the polar regional climate model RACMO2.3 at 11 km. The data set includes all individual SMB components projected to a down-sampled version of the Greenland Ice Mapping Project (GIMP) digital elevation model and ice mask. The 1 km mask better resolves narrow ablation zones, valley glaciers, fjords and disconnected ice caps. Relative to the 11 km product, the more detailed representation of isolated glaciated areas leads to increased precipitation over the southeastern GrIS. In addition, the downscaled product shows a significant increase in runoff owing to better resolved low-lying marginal glaciated regions. The combined corrections for elevation and bare ice albedo markedly improve model agreement with a newly compiled data set of ablation measurements.


2016 ◽  
Author(s):  
Brice Noël ◽  
Willem Jan van de Berg ◽  
Horst Machguth ◽  
Stef Lhermitte ◽  
Ian Howat ◽  
...  

Abstract. This study presents a dataset of daily, 1-km resolution Greenland ice sheet (GrIS) surface mass balance (SMB) covering the period 1958–2015. Using elevation dependence, the high-resolution product is statistically downscaled from the native daily output of the polar regional climate model RACMO2.3 at 11-km. The dataset includes all individual SMB components projected to a down-sampled version of the Greenland Ice Mapping Project (GIMP) digital elevation model and ice mask. The 1-km mask better resolves narrow ablation zones, valley glaciers, fjords and disconnected ice caps. Relative to the 11-km product, the more detailed representation of confined glaciated areas leads to increased precipitation over the southeastern GrIS. In addition, the downscaled product shows a significant increase in runoff owing to better resolved low-lying marginal glaciated regions. The combined corrections for elevation and bare ice albedo markedly improve model agreement with a newly compiled dataset of ablation measurements.


2013 ◽  
Vol 59 (216) ◽  
pp. 733-749 ◽  
Author(s):  
H. Goelzer ◽  
P. Huybrechts ◽  
J.J. Fürst ◽  
F.M. Nick ◽  
M.L. Andersen ◽  
...  

AbstractPhysically based projections of the Greenland ice sheet contribution to future sea-level change are subject to uncertainties of the atmospheric and oceanic climatic forcing and to the formulations within the ice flow model itself. Here a higher-order, three-dimensional thermomechanical ice flow model is used, initialized to the present-day geometry. The forcing comes from a high-resolution regional climate model and from a flowline model applied to four individual marine-terminated glaciers, and results are subsequently extended to the entire ice sheet. The experiments span the next 200 years and consider climate scenario SRES A1B. The surface mass-balance (SMB) scheme is taken either from a regional climate model or from a positive-degree-day (PDD) model using temperature and precipitation anomalies from the underlying climate models. Our model results show that outlet glacier dynamics only account for 6–18% of the sea-level contribution after 200 years, confirming earlier findings that stress the dominant effect of SMB changes. Furthermore, interaction between SMB and ice discharge limits the importance of outlet glacier dynamics with increasing atmospheric forcing. Forcing from the regional climate model produces a 14–31 % higher sea-level contribution compared to a PDD model run with the same parameters as for IPCC AR4.


2020 ◽  
Author(s):  
Christiaan T. van Dalum ◽  
Willem Jan van de Berg ◽  
Michiel R. van den Broeke

Abstract. This study evaluates the impact of a new snow and ice albedo and radiative transfer scheme on the surface mass and energy budget for the Greenland ice sheet in the latest version of the regional climate model RACMO2, version 2.3p3. We also evaluate the modeled (sub)surface temperature and snow melt, as subsurface heating by radiation penetration now occurs. The results are compared to the previous model version and are evaluated against stake measurements and automatic weather station data of the K-transect and PROMICE projects. In addition, subsurface snow temperature profiles are compared at the K-transect, Summit and southeast Greenland. The surface mass balance is in good agreement with observations, and only changes considerably with respect to the previous RACMO2 version around the ice margins and in the percolation zone. Snow melt and refreezing, on the other hand, are changed more substantially in various regions due to the changed albedo representation, subsurface energy absorption and melt water percolation. Internal heating leads to considerably higher snow temperatures in summer, in agreement with observations, and introduces a shallow layer of subsurface melt.


2015 ◽  
Vol 9 (5) ◽  
pp. 1831-1844 ◽  
Author(s):  
B. Noël ◽  
W. J. van de Berg ◽  
E. van Meijgaard ◽  
P. Kuipers Munneke ◽  
R. S. W. van de Wal ◽  
...  

Abstract. We discuss Greenland Ice Sheet (GrIS) surface mass balance (SMB) differences between the updated polar version of the RACMO climate model (RACMO2.3) and the previous version (RACMO2.1). Among other revisions, the updated model includes an adjusted rainfall-to-snowfall conversion that produces exclusively snowfall under freezing conditions; this especially favours snowfall in summer. Summer snowfall in the ablation zone of the GrIS has a pronounced effect on melt rates, affecting modelled GrIS SMB in two ways. By covering relatively dark ice with highly reflective fresh snow, these summer snowfalls have the potential to locally reduce melt rates in the ablation zone of the GrIS through the snow-albedo-melt feedback. At larger scales, SMB changes are driven by differences in orographic precipitation following a shift in large-scale circulation, in combination with enhanced moisture to precipitation conversion for warm to moderately cold conditions. A detailed comparison of model output with observations from automatic weather stations, ice cores and ablation stakes shows that the model update generally improves the simulated SMB-elevation gradient as well as the representation of the surface energy balance, although significant biases remain.


2012 ◽  
Vol 6 (2) ◽  
pp. 255-272 ◽  
Author(s):  
M. M. Helsen ◽  
R. S. W. van de Wal ◽  
M. R. van den Broeke ◽  
W. J. van de Berg ◽  
J. Oerlemans

Abstract. It is notoriously difficult to couple surface mass balance (SMB) results from climate models to the changing geometry of an ice sheet model. This problem is traditionally avoided by using only accumulation from a climate model, and parameterizing the meltwater run-off as a function of temperature, which is often related to surface elevation (Hs). In this study, we propose a new strategy to calculate SMB, to allow a direct adjustment of SMB to a change in ice sheet topography and/or a change in climate forcing. This method is based on elevational gradients in the SMB field as computed by a regional climate model. Separate linear relations are derived for ablation and accumulation, using pairs of Hs and SMB within a minimum search radius. The continuously adjusting SMB forcing is consistent with climate model forcing fields, also for initially non-glaciated areas in the peripheral areas of an ice sheet. When applied to an asynchronous coupled ice sheet – climate model setup, this method circumvents traditional temperature lapse rate assumptions. Here we apply it to the Greenland Ice Sheet (GrIS). Experiments using both steady-state forcing and glacial-interglacial forcing result in realistic ice sheet reconstructions.


2020 ◽  
Author(s):  
Christiaan van Dalum ◽  
Willem Jan van de Berg ◽  
Stef Lhermitte ◽  
Michiel van den Broeke

<p>Snow and ice albedo schemes in present day climate models often lack a sophisticated radiation penetration scheme and are limited to a broadband albedo. In this study, we evaluate a new snow albedo scheme in the regional climate model RACMO2 that uses the two-stream radiative transfer in snow model TARTES and the spectral-to-narrowband albedo module SNOWBAL for the Greenland ice sheet. Additionally, the bare ice albedo parameterization has been updated. The snow and ice albedo output of the updated version of RACMO2, referred to as RACMO2.3p3, is evaluated using PROMICE and K-transect in-situ data and MODIS remote-sensing observations. Generally, the RACMO2.3p3 albedo is in very good agreement with satellite observations, leading to a domain-averaged bias of only -0.012. Some discrepancies are, however, observed for regions close to the ice margin. Compared to the previous iteration RACMO2.3p2, the albedo of RACMO2.3p3 is considerably higher in the bare ice zone during the ablation season, as atmospheric conditions now alter the bare ice albedo. For most other regions, however, the albedo of RACMO2.3p3 is lower due to spectral effects, radiation penetration, snow metamorphism or a delayed firn-ice transition. Furthermore, a white-out effect during cloudy conditions is captured and the snow albedo shows a low sensitivity to low soot concentrations. The surface mass balance of RACMO2.3p3 compares well with observations. Subsurface heating, however, now leads to increased melt and refreezing in south Greenland, changing the snow structure.</p>


2014 ◽  
Vol 8 (1) ◽  
pp. 1057-1093
Author(s):  
R. T. W. L. Hurkmans ◽  
J. L. Bamber ◽  
C. H. Davis ◽  
I. R. Joughin ◽  
K. S. Khvorostovsky ◽  
...  

Abstract. Mass changes of the Greenland ice sheet may be estimated by the Input Output Method (IOM), satellite gravimetry, or via surface elevation change rates (dH / dt). Whereas the first two have been shown to agree well in reconstructing mass changes over the last decade, there are few decadal estimates from satellite altimetry and none that provide a time evolving trend that can be readily compared with the other methods. Here, we interpolate radar and laser altimetry data between 1995 and 2009 in both space and time to reconstruct the evolving volume changes. A firn densification model forced by the output of a regional climate model is used to convert volume to mass. We consider and investigate the potential sources of error in our reconstruction of mass trends, including geophysical biases in the altimetry, and the resulting mass change rates are compared to other published estimates. We find that mass changes are dominated by SMB until about 2001, when mass loss rapidly accelerates. The onset of this acceleration is somewhat later, and less gradual, compared to the IOM. Our time averaged mass changes agree well with recently published estimates based on gravimetry, IOM, laser altimetry, and with radar altimetry when merged with airborne data over outlet glaciers. We demonstrate, that with appropriate treatment, satellite radar altimetry can provide reliable estimates of mass trends for the Greenland ice sheet. With the inclusion of data from CryoSat II, this provides the possibility of producing a continuous time series of regional mass trends from 1992 onward.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Thomas Slater ◽  
Andrew Shepherd ◽  
Malcolm McMillan ◽  
Amber Leeson ◽  
Lin Gilbert ◽  
...  

AbstractRunoff from the Greenland Ice Sheet has increased over recent decades affecting global sea level, regional ocean circulation, and coastal marine ecosystems, and it now accounts for most of the contemporary mass imbalance. Estimates of runoff are typically derived from regional climate models because satellite records have been limited to assessments of melting extent. Here, we use CryoSat-2 satellite altimetry to produce direct measurements of Greenland’s runoff variability, based on seasonal changes in the ice sheet’s surface elevation. Between 2011 and 2020, Greenland’s ablation zone thinned on average by 1.4 ± 0.4 m each summer and thickened by 0.9 ± 0.4 m each winter. By adjusting for the steady-state divergence of ice, we estimate that runoff was 357 ± 58 Gt/yr on average – in close agreement with regional climate model simulations (root mean square difference of 47 to 60 Gt/yr). As well as being 21 % higher between 2011 and 2020 than over the preceding three decades, runoff is now also 60 % more variable from year-to-year as a consequence of large-scale fluctuations in atmospheric circulation. Because this variability is not captured in global climate model simulations, our satellite record of runoff should help to refine them and improve confidence in their projections.


Sign in / Sign up

Export Citation Format

Share Document