Dwindling relevance of large volcanic eruptions for global glacier changes in the anthropocene

Large volcanic eruptions impact climate through the injection of ash and sulfur-containing gases into the atmosphere. While the ash particles fall out rapidly, the gases are converted to sulfate aerosols that reflect solar radiation in the stratosphere and cause a lowering of the global mean surface temperature. Earlier studies have suggested that major volcanic eruptions resulted in positive mass balances and advances of glaciers. Here, we perform a multivariate analysis to decompose global glacier mass changes from 1961 to 2005 into components associated with anthropogenic influences, volcanic and solar activities, and the El Niño-Southern Oscillation. We find that the global glacier mass loss was mainly driven by the anthropogenic forcing, interrupted by a few years of intermittent mass gains following large volcanic eruptions. The relative impact of volcanic eruptions has dwindled due to strongly increasing greenhouse gas concentrations since the mid-20th century

Reconstruction of Past Glacier Changes with an Ice-Flow Glacier Model: Proof of Concept and Validation

Estimations of global glacier mass changes over the course of the 20th century require automated initialization methods, allowing the reconstruction of past glacier states from limited information. In a previous study, we developed a method to initialize the Open Global Glacier Model (OGGM) from past climate information and present-day geometry alone. Tested in an idealized framework, this method aimed to quantify how much information present-day glacier geometry carries about past glacier states. The method was not applied to real-world cases, and therefore, the results were not comparable with observations. This study closes the gap to real-world cases by introducing a glacier-specific calibration of the mass balance model. This procedure ensures that the modeled present-day geometry matches the observed area and that the past glacier evolution is consistent with bias-corrected past climate time series. We apply the method to 517 glaciers, spread globally, for which either mass balance observations or length records are available, and compare the observations to the modeled reconstructed glacier changes. For the validation of the initialization method, we use multiple measures of reconstruction skill (e.g., MBE, RMSE, and correlation). We find that the modeled mass balances and glacier lengths are in good agreement with the observations, especially for glaciers with many observation years. These results open the door to a future global application

The land ice contribution to sea level during the satellite era

Since 1992, there has been a revolution in our ability to quantify the land ice contribution to SLR using a variety of satellite missions and technologies. Each mission has provided unique, but sometimes conflicting, insights into the mass trends of land ice. Over the last decade, over fifty estimates of land ice trends have been published, providing a confusing and often inconsistent picture. The IPCC Fifth Assessment Report (AR5) attempted to synthesise estimates published up to early 2013. Since then, considerable advances have been made in understanding the origin of the inconsistencies, reducing uncertainties in estimates and extending time series. We assess and synthesise results published, primarily, since the AR5, to produce a consistent estimate of land ice mass trends during the satellite era (1992 to 2016). We combine observations from multiple missions and approaches including sea level budget analyses. Our resulting synthesis is both consistent and rigorous, drawing on i) the published literature, ii) expert assessment of that literature, and iii) a new analysis of Arctic glacier and ice cap trends combined with statistical modelling. 
 We present annual and pentad (five-year mean) time series for the East, West Antarctic and Greenland Ice Sheets and glaciers separately and combined. When averaged over pentads, covering the entire period considered, we obtain a monotonic trend in mass contribution to the oceans, increasing from 0.31±0.35 mm of sea level equivalent for 1992-1996 to 1.85±0.13 for 2012-2016. Our integrated land ice trend is lower than many estimates of GRACE-derived ocean mass change for the same periods. This is due, in part, to a smaller estimate for glacier and ice cap mass trends compared to previous assessments. We discuss this, and other likely reasons, for the difference between GRACE ocean mass and land ice trends

Impact of frontal ablation on the ice thickness estimation of marine-terminating glaciers in Alaska

Frontal ablation is a major component of the mass budget of calving glaciers, strongly affecting their dynamics. Most global-scale ice volume estimates to date still suffer from considerable uncertainties related to (i) the implemented frontal ablation parameterization or (ii) not accounting for frontal ablation at all in the glacier model. To improve estimates of the ice thickness distribution of glaciers, it is thus important to identify and test low-cost and robust parameterizations of this process. By implementing such parameterization into the ice thickness estimation module of the Open Global Glacier Model (OGGM v1.1.2), we conduct a first assessment of the impact of accounting for frontal ablation on the estimate of ice stored in glaciers in Alaska. We find that inversion methods based on mass conservation systematically underestimate the mass turnover and, therefore, the thickness of tidewater glaciers when neglecting frontal ablation. This underestimation can amount to up to 19 % on a regional scale and up to 30 % for individual glaciers. The effect is independent of the size of the glacier. Additionally, we perform different sensitivity experiments to study the influence of (i) a constant of proportionality (k) used in the frontal ablation parameterization, (ii) Glen’s temperature-dependent creep parameter (A) and (iii) a sliding velocity parameter (fs) on the regional dynamics of Alaska tidewater glaciers. OGGM is able to reproduce previous regional frontal ablation estimates, applying a number of combinations of values for k, Glen’s A and fs. Our sensitivity studies also show that differences in thickness between accounting for and not accounting for frontal ablation occur mainly at the lower parts of the glacier, both above and below sea level. This indicates that not accounting for frontal ablation will have an impact on the estimate of the glaciers’ potential contribution to sea-level rise. Introducing frontal ablation increases the volume estimate of Alaska marine-terminating glaciers from 9.18±0.62 to 10.61±0.75 mm s.l.e, of which 1.52±0.31 mm s.l.e (0.59±0.08 mm s.l.e when ignoring frontal ablation) are found to be below sea level

Future sea level rise constrained by observations and long-term commitment

Sea level has been steadily rising over the past century, predominantly due to anthropogenic climate change. The rate of sea level rise will keep increasing with continued global warming, and, even if temperatures are stabilized through the phasing out of greenhouse gas emissions, sea level is still expected to rise for centuries. This will affect coastal areas worldwide, and robust projections are needed to assess mitigation options and guide adaptation measures. Here we combine the equilibrium response of the main sea level rise contributions with their last century’s observed contribution to constrain projections of future sea level rise. Our model is calibrated to a set of observations for each contribution, and the observational and climate uncertainties are combined to produce uncertainty ranges for 21st century sea level rise. We project anthropogenic sea level rise of 28–56 cm, 37–77 cm, and 57–131 cm in 2100 for the greenhouse gas concentration scenarios RCP26, RCP45, and RCP85, respectively. Our uncertainty ranges for total sea level rise overlap with the process-based estimates of the Intergovernmental Panel on Climate Change. The “constrained extrapolation” approach generalizes earlier global semiempirical models and may therefore lead to a better understanding of the discrepancies with processbased projections

Global sea-level budget and ocean-mass budget, with a focus on advanced data products and uncertainty characterisation

Studies of the global sea-level budget (SLB) and the global ocean-mass budget (OMB) are essential to assess the reliability of our knowledge of sea-level change and its contributors. Here we present datasets for times series of the SLB and OMB elements developed in the framework of ESA's Climate Change Initiative. We use these datasets to assess the SLB and the OMB simultaneously, utilising a consistent framework of uncertainty characterisation. The time series, given at monthly sampling and available at https://doi.org/10.5285/17c2ce31784048de93996275ee976fff (Horwath et al., 2021), include global mean sea-level (GMSL) anomalies from satellite altimetry, the global mean steric component from Argo drifter data with incorporation of sea surface temperature data, the ocean-mass component from Gravity Recovery and Climate Experiment (GRACE) satellite gravimetry, the contribution from global glacier mass changes assessed by a global glacier model, the contribution from Greenland Ice Sheet and Antarctic Ice Sheet mass changes assessed by satellite radar altimetry and by GRACE, and the contribution from land water storage anomalies assessed by the global hydrological model WaterGAP (Water Global Assessment and Prognosis). Over the period January 1993–December 2016 (P1, covered by the satellite altimetry records), the mean rate (linear trend) of GMSL is 3.05 ± 0.24 mm yr−1. The steric component is 1.15 ± 0.12 mm yr−1 (38 % of the GMSL trend), and the mass component is 1.75 ± 0.12 mm yr−1 (57 %). The mass component includes 0.64  ± 0.03 mm yr−1 (21 % of the GMSL trend) from glaciers outside Greenland and Antarctica, 0.60 ± 0.04 mm yr−1 (20 %) from Greenland, 0.19 ± 0.04 mm yr−1 (6 %) from Antarctica, and 0.32 ± 0.10 mm yr−1 (10 %) from changes of land water storage. In the period January 2003–August 2016 (P2, covered by GRACE and the Argo drifter system), GMSL rise is higher than in P1 at 3.64 ± 0.26 mm yr−1. This is due to an increase of the mass contributions, now about 2.40 ± 0.13 mm yr−1 (66 % of the GMSL trend), with the largest increase contributed from Greenland, while the steric contribution remained similar at 1.19 ± 0.17 mm yr−1 (now 33 %). The SLB of linear trends is closed for P1 and P2; that is, the GMSL trend agrees with the sum of the steric and mass components within their combined uncertainties. The OMB, which can be evaluated only for P2, shows that our preferred GRACE-based estimate of the ocean-mass trend agrees with the sum of mass contributions within 1.5 times or 0.8 times the combined 1σ uncertainties, depending on the way of assessing the mass contributions. Combined uncertainties (1σ) of the elements involved in the budgets are between 0.29 and 0.42 mm yr−1, on the order of 10 % of GMSL rise. Interannual variations that overlie the long-term trends are coherently represented by the elements of the SLB and the OMB. Even at the level of monthly anomalies the budgets are closed within uncertainties, while also indicating possible origins of remaining misclosures

Exploring the impact of a frontal ablation parameterization on projected 21st-century mass change for Northern Hemisphere glaciers

Marine-terminating glaciers cover more than one-fourth of the total glacierized area in the Northern Hemisphere outside the Greenland ice sheet. It is therefore crucial to ensure an adequate representation of these glaciers when projecting large-scale glacier mass changes. We investigate how the introduction of marine frontal processes in the modeling chain influences the results of mass change projections, compared to projections neglecting such processes. We find that including frontal processes reduces the projected glacier mass loss, since incorporating frontal ablation in the model's mass-balance calibration results in a decrease in marine-terminating glaciers’ sensitivity to atmospheric temperatures. We also find that retrograde bed slopes lead to increased frontal ablation as the atmosphere warms, while frontal ablation decreases if bed slopes are prograde. These opposing effects have the potential to partly cancel each other when considering large glacier ensembles. Although we do not account for potential future changes in oceanic climate yet, any effect of these would be moderated by around half of today's marine-terminating glaciers becoming land-terminating in the course of the 21st century. While we find a significant influence of ice flow parameters on our results, boundary conditions remain the largest source of uncertainty in our projections

Spatial correlation structures of wind speed and irradiance in Europe as modelled in regional climate models and the ERA5 reanalysis

Objective & Background For comprehensive energy systems analysis considering high shares of weather-dependent renewables, the spatiotemporal characteristics of meteorological input data (in particular, wind speed and irradiance) are of considerable interest. Studies aiming to evaluate the potential effects of future climate change rely on model-based projections of these variables. In order to reach as high a spatiotemporal resolution as possible, global climate models (GCM) can be used as drivers for regional climate models (RCM). In doing so, global climate effects are considered on large scales, while the higher resolution computations are limited to the domain of interest. Here, we analyse such regional climate model outputs as well as reanalysis data with regards to their suitability in energy systems analysis, using spatial correlation structures of the wind and solar resource anomalies. Method Based on 10 years of historical experiments of EURO-CORDEX regional climate model output and ERA5 reanalyses, we compute climate anomalies of both wind speed and irradiance for each point in space and time by subtracting the multi-year mean of each point from the instantaneous values at that point. In doing so, we largely remove effects of seasonal and diurnal cycles, especially in the case of irradiance. From each of these wind speed and irradiance anomaly fields, we sample distinct two-dimensional correlation structures (i.e., maps of the correlation coefficient between each pixel’s time series and a reference pixel’s series). Increasing the number of reference pixels hence increases the number of spatial correlation structures available in the analysis. We finally estimate the degree of similarity between the spatial correlation structures of each EURO-CORDEX member and its corresponding ERA5 counterpart using the coefficient of pattern correlation (i.e., the correlation coefficient of the mapped spatial correlation structures in vector form). Principal Findings Kernel density estimates and box-plot statistics of all available pattern correlations show that the wind speed anomaly spatial correlation structures tend to be better represented than their irradiance counterparts. Between different driving global models, the performance is relatively similar for wind speed anomaly correlation structures, while solar irradiance anomaly structures can be associated with a few unusually low pattern correlations for some global models. In terms of regional climate models, these cases of low pattern correlation values are largely associated with the SMHI-RCA4 regional model during the summer. Discussion This initial comparison of CORDEX model output and ERA5 reanalysis based on correlation structures in wind speed and solar irradiance anomalies allows to help in the selection of global and regional models when compiling an ensemble for use in energy systems analysis. In this context, the spatial correlation structures of wind speed need less attention when selecting model outputs, while the irradiance anomaly projections of the regional climate model SMHI-RCA4 should be handled with caution. Increasing the number of reference points used in the spatial correlation structure calculations may lead to further insights regarding potential regional differences in CORDEX model performances

GlacierMIP – A model intercomparison of global-scale glacier mass-balance models and projections

Global-scale 21st-century glacier mass change projections from six published global glacier models are systematically compared as part of the Glacier Model Intercomparison Project. In total 214 projections of annual glacier mass and area forced by 25 General Circulation Models (GCMs) and four Representative Concentration Pathways (RCP) emission scenarios and aggregated into 19 glacier regions are considered. Global mass loss of all glaciers (outside the Antarctic and Greenland ice sheets) by 2100 relative to 2015 averaged over all model runs varies from 18 ± 7% (RCP2.6) to 36 ± 11% (RCP8.5) corresponding to 94 ± 25 and 200 ± 44 mm sea-level equivalent (SLE), respectively. Regional relative mass changes by 2100 correlate linearly with relative area changes. For RCP8.5 three models project global rates of mass loss (multi-GCM means) of >3 mm SLE per year towards the end of the century. Projections vary considerably between regions, and also among the glacier models. Global glacier mass changes per degree global air temperature rise tend to increase with more pronounced warming indicating that mass-balance sensitivities to temperature change are not constant. Differences in glacier mass projections among the models are attributed to differences in model physics, calibration and downscaling procedures, initial ice volumes and varying ensembles of forcing GCMs