Impact of the choice of surface mass balance models and their calibration on large-scale glacier change projections

Author(s):  
Lilian Schuster ◽  
David Rounce ◽  
Fabien Maussion

<p>A recent large model intercomparison study (GlacierMIP) showed that differences between the glacier models is a dominant source of uncertainty for future glacier change projections, in particular in the first half of the century.  Each glacier model has their own unique set of process representations and climate forcing methodology, which makes it impossible to determine the model components that contribute most to the projection uncertainty. This study aims to improve our understanding of the sources of large scale glacier model uncertainty using the Open Global Glacier Model (OGGM), focussing on the surface mass balance (SMB) in a first step. We calibrate and run a set of interchangeable SMB model parameterizations (e.g. monthly vs. daily, constant vs. variable lapse rates, albedo, snowpack evolution and refreezing) under controlled boundary conditions. Based on ensemble approaches, we explore the influence of (i) the parameter calibration strategy and (ii) SMB model complexity on regional to global glacier change. These uncertainties are then put in relation to a qualitative selection of other model design choices, such as the forcing climate dataset and ice dynamics model parameters. </p>

2018 ◽  
Author(s):  
Harry Zekollari ◽  
Matthias Huss ◽  
Daniel Farinotti

Abstract. Glaciers in the European Alps play an important role in the hydrological cycle, act as a source for hydroelectricity and have a large touristic importance. The future evolution of these glaciers is driven by surface mass balance and ice flow processes, which the latter is to date not included in regional glacier projections for the Alps. Here, we model the future evolution of glaciers in the European Alps with GloGEMflow, an extended version of the Global Glacier Evolution Model (GloGEM), in which both surface mass balance and ice flow are explicitly accounted for. The mass balance model is calibrated with glacier-specific geodetic mass balances, and forced with high-resolution regional climate model (RCM) simulations from the EURO-CORDEX ensemble. The evolution of the total glacier volume in the coming decades is relatively similar under the various representative concentrations pathways (RCP2.6, 4.5 and 8.5), with volume losses of about 47–52 % in 2050 with respect to 2017. We find that under RCP2.6, the ice loss in the second part of the 21st century is relatively limited and that about one-third (36.8 % ± 11.1 %) of the present-day (2017) ice volume will still present in 2100. Under a strong warming (RCP8.5) the future evolution of the glaciers is dictated by a substantial increase in surface melt, and glaciers are projected to largely disappear by 2100 (94.4 ± 4.4 % volume loss vs. 2017). For a given RCP, differences in future changes are mainly determined by the driving global climate model, rather than by the RCM that is coupled to it, and these differences are larger than those arising from various model parameters. We find that under a limited warming, the inclusion of ice dynamics reduces the projected mass loss and that this effect increases with the glacier elevation range, implying that the inclusion of ice dynamics is likely to be important for global glacier evolution projections.


2019 ◽  
Vol 13 (4) ◽  
pp. 1125-1146 ◽  
Author(s):  
Harry Zekollari ◽  
Matthias Huss ◽  
Daniel Farinotti

Abstract. Glaciers in the European Alps play an important role in the hydrological cycle, act as a source for hydroelectricity and have a large touristic importance. The future evolution of these glaciers is driven by surface mass balance and ice flow processes, of which the latter is to date not included explicitly in regional glacier projections for the Alps. Here, we model the future evolution of glaciers in the European Alps with GloGEMflow, an extended version of the Global Glacier Evolution Model (GloGEM), in which both surface mass balance and ice flow are explicitly accounted for. The mass balance model is calibrated with glacier-specific geodetic mass balances and forced with high-resolution regional climate model (RCM) simulations from the EURO-CORDEX ensemble. The evolution of the total glacier volume in the coming decades is relatively similar under the various representative concentrations pathways (RCP2.6, 4.5 and 8.5), with volume losses of about 47 %–52 % in 2050 with respect to 2017. We find that under RCP2.6, the ice loss in the second part of the 21st century is relatively limited and that about one-third (36.8 % ± 11.1 %, multi-model mean ±1σ) of the present-day (2017) ice volume will still be present in 2100. Under a strong warming (RCP8.5) the future evolution of the glaciers is dictated by a substantial increase in surface melt, and glaciers are projected to largely disappear by 2100 (94.4±4.4 % volume loss vs. 2017). For a given RCP, differences in future changes are mainly determined by the driving global climate model (GCM), rather than by the RCM, and these differences are larger than those arising from various model parameters (e.g. flow parameters and cross-section parameterisation). We find that under a limited warming, the inclusion of ice dynamics reduces the projected mass loss and that this effect increases with the glacier elevation range, implying that the inclusion of ice dynamics is likely to be important for global glacier evolution projections.


2016 ◽  
Vol 62 (236) ◽  
pp. 1083-1092 ◽  
Author(s):  
SHUN TSUTAKI ◽  
SHIN SUGIYAMA ◽  
DAIKI SAKAKIBARA ◽  
TAKANOBU SAWAGAKI

ABSTRACTTo quantify recent thinning of marine-terminating outlet glaciers in northwestern Greenland, we carried out field and satellite observations near the terminus of Bowdoin Glacier. These data were used to compute the change in surface elevation from 2007 to 2013 and this rate of thinning was then compared with that of the adjacent land-terminating Tugto Glacier. Comparing DEMs of 2007 and 2010 shows that Bowdoin Glacier is thinning more rapidly (4.1 ± 0.3 m a−1) than Tugto Glacier (2.8 ± 0.3 m a−1). The observed negative surface mass-balance accounts for <40% of the elevation change of Bowdoin Glacier, meaning that the thinning of Bowdoin Glacier cannot be attributable to surface melting alone. The ice speed of Bowdoin Glacier increases down-glacier, reaching 457 m a−1 near the calving front. This flow regime causes longitudinal stretching and vertical compression at a rate of −0.04 a−1. It is likely that this dynamically-controlled thinning has been enhanced by the acceleration of the glacier since 2000. Our measurements indicate that ice dynamics indeed play a predominant role in the rapid thinning of Bowdoin Glacier.


2018 ◽  
Author(s):  
Andreas Plach ◽  
Kerim H. Nisancioglu ◽  
Sébastien Le clec’h ◽  
Andreas Born ◽  
Petra M. Langebroek ◽  
...  

Abstract. Understanding the behavior of the Greenland ice sheet in a warmer climate, and particularly its surface mass balance (SMB), is important for assessing Greenland’s potential contribution to future sea level rise. The Eemian interglacial, the most recent warmer-than-present period in Earth’s history approximately 125 000 years ago, provides an analogue for a warm summer climate over Greenland. The Eemian is characterized by a positive Northern Hemisphere summer insolation anomaly, which introduces uncertainties in Eemian SMB when using positive degree day estimates. In this study, we use Eemian global and regional climate simulations in combination with three types of SMB models – a simple positive degree day, an intermediate complexity, and a full surface energy balance model – to evaluate the importance of regional climate and model complexity for estimates of Greenland SMB. We find that all SMB models perform well under the relatively cool pre-industrial and late Eemian. For the relatively warm early Eemian, the differences between SMB models are large which is associated with the representation of insolation in the respective models. For all simulated time slices there is a systematic difference between globally and regionally forced SMB models, due to the different representation of the regional climate over Greenland. We conclude that both the resolution of the simulated climate as well as the method used to estimate the SMB, are important for an accurate simulation of Greenland’s SMB. Whether model resolution or SMB method is most important depends on the climate state and in particular the prevailing insolation pattern. We suggest that future Eemian climate model inter-comparison studies are combined with different SMB models to quantify Eemian SMB uncertainty estimates.


2005 ◽  
Vol 41 ◽  
pp. 97-104 ◽  
Author(s):  
W.J. Van De Berg ◽  
M.R. Van Den Broeke ◽  
C.H. Reijmer ◽  
E. Van Meijgaard

AbstractTemporal and spatial characteristics of the Antarctic specific surface mass balance (SSMB) are presented, including its components solid precipitation, sublimation/deposition and melt. For this purpose, we use the output of a regional atmospheric climate model (RACMO2/ANT, horizontal resolution of ~55 km) for the period 1958–2002. RACMO2/ANT uses European Centre for Medium-Range Weather Forecasts (ECMWF) 40 year re-analysis (ERA-40) fields as forcing at the lateral boundaries. RACMO2/ANT underestimates SSMB in the high interior of East and West Antarctica and overestimates SSMB on the steep coastal slopes. Otherwise, the modeled spatial pattern of SSMB is in good qualitative agreement with recent compilations of in situ observations. Large-scale patterns, like the precipitation shadow effect of the Antarctic Peninsula, are well reproduced, and mesoscale SSMB patterns, such as the strong precipitation gradients on Law Dome, are well represented in the model. The integrated SSMB over the grounded ice sheet is 153mmw.e. a–1 for the period 1958–2002, which agrees within 5% with the latest measurement compilations. Sublimation and melt remove 7% and <1% respectively of the solid precipitation. We found significant seasonality of solid precipitation, with a maximum in autumn and a minimum in summer. No meaningful trend was identified for the SSMB, because the time series of solid precipitation and SSMB are affected by an inhomogeneity in 1980 within the ERA-40 fields that drive RACMO2/ANT. Sublimation, melt and liquid precipitation increase in time, which is related to a modeled increase in 2m temperature.


2016 ◽  
Vol 62 (233) ◽  
pp. 525-534 ◽  
Author(s):  
DENIS CALLENS ◽  
REINHARD DREWS ◽  
EMMANUEL WITRANT ◽  
MORGANE PHILIPPE ◽  
FRANK PATTYN

ABSTRACTIce rises are locally grounded parts of Antarctic ice shelves that play an important role in regulating ice flow from the continent towards the ocean. Because they protrude out of the otherwise horizontal ice shelves, ice rises induce an orographic uplift of the atmospheric flow, resulting in an asymmetric distribution of the surface mass balance (SMB). Here, we combine younger and older internal reflection horizons (IRHs) from radar to quantify this distribution in time and space across Derwael Ice Rise (DIR), Dronning Maud Land, Antarctica. We employ two methods depending on the age of the IRHs, i.e. the shallow layer approximation for the younger IRHs near the surface and an optimization technique based on an ice flow model for the older IRHs. We identify an SMB ratio of 2.5 between the flanks and the ice divide with the SMB ranging between 300 and 750 kg m−2 a−1. The SMB maximum is located on the upwind side, ~4 km offset to today's topographic divide. The large-scale asymmetry is consistently observed in time until 1966. The SMB from older IRHs is less-well constrained, but the asymmetry has likely persisted for >ka, indicating that DIR has been a stable features over long time spans.


2002 ◽  
Vol 35 ◽  
pp. 67-72 ◽  
Author(s):  
Edward Hanna ◽  
Philippe Huybrechts ◽  
Thomas L. Mote

AbstractWe used surface climate fields from high-resolution (~0.5660.56˚) European Centre for Medium-RangeWeather Forecasts (ECMWF) operational analyses (1992–98), together with meteorological and glaciological models of snow accumulation and surface meltwater runoff/retention, to produce novel maps of Greenland ice sheet (GIS) net accumulation, net runoff and surface mass balance (SMB). We compared our runoff maps with similar-scaled runoff (melt minus refreezing) maps based on passive-microwave satellite data. Our gross spatial/temporal patterns of runoff compared well with those from the satellite data, although amounts of modelled runoff are likely too low. Mean accumulation was 0.287 (0.307)ma–1, and mean runoff was 0.128 (0.151)ma–1, averaged across the W. Abdalati (T. L. Mote) GIS mask. Corresponding mean SMB was 0.159 (0.156)ma–1, with considerable interannual variability (standard deviation ~0.11ma–1) primarily due to variations in runoff. Considering best estimates of current iceberg calving, overall the GIS is probably currently losing mass. Our study shows great promise for meaningfully modelling SMB based on forthcoming ``second-generation’’ ECMWF re-analysis (ERA-40) data, and comparing the results with ongoing laser/radarmeasurements of surface elevation. This should help elucidate to what extent surface elevation changes are caused by short-term SMB variations or other factors (e.g. ice dynamics).


2021 ◽  
Author(s):  
Nicolaj Hansen ◽  
Peter L. Langen ◽  
Fredrik Boberg ◽  
Rene Forsberg ◽  
Sebastian B. Simonsen ◽  
...  

Abstract. Antarctic surface mass balance (SMB) is largely determined by precipitation over the continent and subject to regional climate variability related to the Southern Annular Mode (SAM) and other climatic drivers at the large scale. Locally however, firn and snow pack processes are important in determining SMB and the total mass balance of Antarctica and global sea level. Here, we examine factors that influence Antarctic SMB and attempt to reconcile the outcome with estimates for total mass balance determined from the GRACE satellites. This is done by having the regional climate model HIRHAM5 forcing two versions of an offline subsurface model, to estimate Antarctic ice sheet (AIS) SMB from 1980 to 2017. The Lagrangian subsurface model estimates AIS SMB of 2473.5 ± 114.4 Gt per year, while the Eulerian subsurface model variant results in slightly higher modelled SMB of 2564.8 ± 113.7 Gt per year. The majority of this difference in modelled SMB is due to melt and refreezing over ice shelves and demonstrates the importance of firn modelling in areas with substantial melt. Both the Eulerian and the Lagrangian SMB estimates are within uncertainty ranges of each other and within the range of other SMB studies. However, the Lagrangian version has better statistics when modelling the densities. There is a mean bias in modelled density of −24.0 ± 18.4 kg m−3 and −8.2 ± 15.3 kg m−3 for layers less than 550 kg m−3 for the Eulerian and Lagrangian framework, respectively. For layers with a density above 550 kg m−3 the bias is −31.7 ± 23.4 kg m−3 and −35.0 ± 23.7 kg m−3 for the Eulerian and Lagrangian framework, respectively. The mean firn 10 m temperature bias is 0.42–0.52 °C. Further, analysis of the relationship between SMB in individual drainage basins and the SAM, is carried out using a bootstrapping approach. This shows a robust relationship between SAM and SMB in half of the basins (13 out of 27). In general, when SAM is positive there is a lower SMB over the Plateau and a higher SMB on the westerly side of the Antarctic Peninsula, and vice versa when the SAM is negative. Finally, we compare the modelled SMB to GRACE data by subtracting the solid ice discharge, and find that there is a good agreement in East Antarctica, but large disagreements over the Antarctic Peninsula.There is a large difference between published estimates of discharge that make it challenging to use mass reconciliation in evaluating SMB models on the basin scale.


2018 ◽  
Author(s):  
Cécile Agosta ◽  
Charles Amory ◽  
Christoph Kittel ◽  
Anais Orsi ◽  
Vincent Favier ◽  
...  

Abstract. The Antarctic ice sheet mass balance is a major component of the sea level budget and results from the difference of two fluxes of a similar magnitude: ice flow discharging in the ocean and net snow accumulation on the ice sheet surface, i.e. the surface mass balance (SMB). Separately modelling ice dynamics and surface mass balance is the only way to project future trends. In addition, mass balance studies frequently use regional climate models (RCMs) outputs as an alternative to observed fields because SMB observations are particularly scarce on the ice sheet. Here we evaluate new simulations of the polar RCM MAR forced by three reanalyses, ERA-Interim, JRA-55 and MERRA2, for the period 1979–2015, and we compare our results to the last outputs of the RCM RACMO2 forced by ERA-Interim. We show that MAR and RACMO2 perform similarly well in simulating coast to plateau SMB gradients, and we find no significant differences in their simulated SMB when integrated over the ice sheet or its major basins. More importantly, we outline and quantify missing processes in both RCMs. Along stake transects, we show that both models accumulate too much snow on crests, and not enough snow in valleys, as a result of erosion-deposition processes not included in MAR, where the drifting snow module has been switched off, and probably underestimated in RACMO2 by a factor of three. As a consequence, the amount of drifting snow sublimating in the atmospheric boundary layer remains a potentially large mass sink needed to be better constrained. Moreover, MAR generally simulates larger SMB and snowfall amounts than RACMO2 inland, whereas snowfall rates are significantly lower in MAR than in RACMO2 at the ice sheet margins. This divergent behaviour at the margins results from differences in model parameterisations, as MAR explicitly advects precipitating particles through the atmospheric layers and sublimates snowflakes in the undersaturated katabatic layer, whereas in RACMO2 precipitation is added to the surface without advection through the atmosphere. Consequently, we corroborate a recent study concluding that sublimation of precipitation in the low-level atmospheric layers is a significant mass sink for the Antarctic SMB, as it may represent ∼ 240 ± 25 Gt yr-1 of difference in snowfall between RACMO2 and MAR for the period 1979–2015, which is 10 % of the simulated snowfall loaded on the ice sheet and more than twice the surface snow sublimation as currently simulated by MAR.


2020 ◽  
Author(s):  
Marie G. P. Cavitte ◽  
Quentin Dalaiden ◽  
Hugues Goosse ◽  
Jan T. M. Lenaerts ◽  
Elizabeth R. Thomas

Abstract. Ice cores are an important record of the past surface mass balance (SMB) of ice sheets, with SMB mitigating the ice sheets’ sea level impact over the recent decades. For the Antarctic Ice Sheet (AIS), SMB is dominated by large-scale atmospheric circulation, which collects warm moist air from further north and releases it in the form of snow as widespread accumulation or focused atmospheric rivers on the continent. This implies that the snow deposited at the surface of the AIS should record strongly coupled SMB and surface air temperature (SAT) variations. Ice cores use δ18O as a proxy for SAT as they do not record SAT directly. Here, using isotope-enabled global climate models and the RACMO2.3 regional climate model, we calculate positive SMB-SAT and δ18O-SMB correlations over ∼90 % of the AIS. The high spatial resolution of the RACMO2.3 model allows us to highlight a number of areas where SMB and SAT are not correlated, and show that wind-driven processes acting locally, such as Foehn and katabatic effects, can overwhelm the large-scale atmospheric input in SMB and SAT responsible for the positive SMB-SAT correlations. We focus in particular on Dronning Maud Land, East Antarctica, where the ice promontories clearly show these wind-induced effects. However, using the PAGES2k ice core compilations of SMB and δ18O of Thomas et al. (2017) and Stenni et al. (2017), we obtain a weak correlation, on the order of 0.1, between SMB and δ18O over the past ~150 years. We obtain an equivalently weak correlation between ice core SMB and the SAT reconstruction of Nicolas and Bromwich (2014) over the past ~50 years, although the ice core sites are not spatially co-located with the areas displaying a low SMB-SAT correlation in the models. To resolve the discrepancy between the measured and modeled signals, we show that averaging the ice core records in close spatial proximity increases their SMB-SAT correlation. This increase shows that the weak measured correlation likely results from random noise present in the ice core records, but is not large enough to match the correlation calculated in the models. Our results indicate thus a positive correlation between SAT and SMB in models and ice core reconstructions but with a weaker value in observations that may be due to missing processes in models or some systematic biases in ice core data that are not removed by a simple average.


Sign in / Sign up

Export Citation Format

Share Document