Regional and global projections of twenty-first century glacier mass changes in response to climate scenarios from global climate models

A large component of present-day sea-level rise is due to the melt of glaciers other than the ice sheets. Recent projections of their contribution to global sea-level rise for the twenty-first century range between 70 and 180 mm, but bear significant uncertainty due to poor glacier inventory and lack of hypsometric data. Here, we aim to update the projections and improve quantification of their uncertainties by using a recently released global inventory containing outlines of almost every glacier in the world. We model volume change for each glacier in response to transient spatially-differentiated temperature and precipitation projections from 14 global climate models with two emission scenarios (RCP4.5 and RCP8.5) prepared for the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. The multi-model mean suggests sea-level rise of 155 ± 41 mm (RCP4.5) and 216 ± 44 mm (RCP8.5) over the period 2006–2100, reducing the current global glacier volume by 29 or 41 %. The largest contributors to projected global volume loss are the glaciers in the Canadian and Russian Arctic, Alaska, and glaciers peripheral to the Antarctic and Greenland ice sheets. Although small contributors to global volume loss, glaciers in Central Europe, low-latitude South America, Caucasus, North Asia, and Western Canada and US are projected to lose more than 80 % of their volume by 2100. However, large uncertainties in the projections remain due to the choice of global climate model and emission scenario. With a series of sensitivity tests we quantify additional uncertainties due to the calibration of our model with sparsely observed glacier mass changes. This gives an upper bound for the uncertainty range of ±84 mm sea-level rise by 2100 for each projection.     

Contribution of glacier melt to sea-level rise since AD 1865: a regionally differentiated calculation

 The contribution of glacier melt, including the Greenland ice-sheet, to sea-level change since AD 1865 is estimated on the basis of modelled sensitivity of glacier mass balance to climate change and historical temperature data. Calculations are done in a regionally differentiated manner to overcome the inhomogeneity of the global distribution of glaciers. A distinction is made between changes in summer temperature and in temperature over the rest of the year. Our best estimate of the ice melt in the period 1865–1990 in terms of sea-level change equivalent is 5.7 cm (2.7 cm for glaciers and 3.0 cm for the Greenland ice-sheet). Additional calculations show that simpler methods, like using annual or even global mean temperature anomaly give estimates that differ by up to 55%. Consequently, a regionally differentiating approach is advised for making projections of glacier melt with GCM output. Received: 6 December 1996/Accepted: 30 May 1997     

Projecting twenty-first century regional sea-level changes

We present regional sea-level projections and associated uncertainty estimates for the end of the 21 ^st century. We show regional projections of sea-level change resulting from changing ocean circulation, increased heat uptake and atmospheric pressure in CMIP5 climate models. These are combined with model- and observation-based regional contributions of land ice, groundwater depletion and glacial isostatic adjustment, including gravitational effects due to mass redistribution. A moderate and a warmer climate change scenario are considered, yielding a global mean sea-level rise of 0.54 ±0.19 m and 0.71 ±0.28 m respectively (mean ±1σ). Regionally however, changes reach up to 30 % higher in coastal regions along the North Atlantic Ocean and along the Antarctic Circumpolar Current, and up to 20 % higher in the subtropical and equatorial regions, confirming patterns found in previous studies. Only 50 % of the global mean value is projected for the subpolar North Atlantic Ocean, the Arctic Ocean and off the western Antarctic coast. Uncertainty estimates for each component demonstrate that the land ice contribution dominates the total uncertainty.     

Modelling the response of glaciers to climate change by applying volume-area scaling in combination with a high resolution GCM

 A seasonally and regionally differentiated glacier model is used to estimate the contribution that glaciers are likely to make to global sea level rise over a period of 70 years. A high resolution general circulation model (ECHAM4 T106) is used to estimate temperature and precipitation changes for a doubled CO₂ climate and serves as input for the glacier model. Volume-area relations are used to take into account the reduction of glacier area resulting from greenhouse warming. Each glacieriated region has a specified glacier size distribution, defined by the number of glaciers in a size class and a mean area. Changes in glacier volume are calculated by a precipitation dependent mass balance sensitivity. The model predicts a global sea level rise of 57 mm over a period of 70 years. This corresponds to a sensitivity of 0.86 mm yr⁻¹K⁻¹. Assuming a constant glacier area as done in earlier work leads to an overestimation of 19% for the contribution to sea level rise. Received: 16 August 2000 / Accepted: 21 May 2001     

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6?±?1 Gta^?1, which decreases about sixfold to -35?±?2 Gta^?1 by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12?±?2 Gta^?1 by the 2080s. We also find that 10.6 and 27 % of the glaciers could face ‘eventual disappearance’ by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.     

A precipitation-dominated,mid-latitude glacier system: Mount Shasta,California

Temperature is often seen as the dominant control on inter-decadal glacier volume changes. However, despite regional warming over the past half-century, the glaciers of Mount Shasta have continued to expand following a contraction during a prolonged drought in the early twentieth century, indicating a greater sensitivity to precipitation than temperature. We use the 110 year record of fluctuations in Mount Shasta’s glaciers and climate to calibrate numerical glacier models of the two largest glaciers. The reconstructed balance and volume histories show a much greater correlation to precipitation than temperature and significant correlation to oscillatory modes of Pacific Ocean climate. An approximately 20% increase in precipitation is needed for every 1°C increase in temperature to maintain stability. Under continued historical trends, oscillations in climate modes and random variability will dominate inter-decadal variability in ice volume. Under the strong warming trend predicted by a regional climate model, the temperature trend will be the dominant forcing resulting in near total loss of Mount Shasta’s glaciers by the end of the twenty-first century.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Projected change in East Asian summer monsoon precipitation under RCP scenario

Future changes in East Asian summer monsoon precipitation climatology, frequency, and intensity are analyzed using historical climate simulations and future climate simulations under the RCP4.5 scenario using the World Climate Research Programme’s (WCRP) Coupled Model Intercomparison Project 5 (CMIP5) multi-model dataset. The model reproducibility is evaluated, and well performance in the present-day climate simulation can be obtained by most of the studied models. However, underestimation is obvious over the East Asian region for precipitation climatology and precipitation intensity, and overestimation is observed for precipitation frequency. The overestimation of precipitation frequency is mainly due to the large positive bias of the light precipitation (precipitation <10 mm/day) days, and the underestimation of precipitation intensity is mainly caused by the negative bias of the intense precipitation (precipitation >10 mm/day) intensity. For the future climate simulations, simple multi-model ensemble (MME) averages using all of the models show increases in precipitation and its intensity over almost all of East Asia, while the precipitation frequency is projected to decrease over eastern China and around Japan and increase in other regions. When the weighted MME is considered, no large difference can be observed compared with the simple MME. For the MME using the six best models that have good performance in simulating the present-day climate, the future climate changes over East Asia are very similar to those predicted using all of the models. Further analysis shows that the frequency and intensity of intense precipitation events are also projected to significantly increase over East Asia. Increases in precipitation frequency and intensity are the main contributors to increases in precipitation, and the contribution of frequency increases (contributed by 40.8 % in the near future and by 58.9 % by the end of the twenty-first century) is much larger than that of intensity increases (contributed by 29.9 % in the near future and by 30.1 % by the end of the twenty-first century). This finding also implies an increased risk of intense precipitation events over the East Asian region under global warming scenario. These results regarding future climate simulations show much greater reliability than those using CMIP3 simulations.     

Climate scenarios of sea level rise for the northeast Atlantic Ocean: a study including the effects of ocean dynamics and gravity changes induced by ice melt

Here we present a set of regional climate scenarios of sea level rise for the northeast Atlantic Ocean. In this study, the latest observations and results obtained with state-of-the-art climate models are combined. In addition, regional effects due to ocean dynamics and changes in the Earth’s gravity field induced by melting of land-based ice masses have been taken into account. The climate scenarios are constructed for the target years 2050 and 2100, for both a moderate and a large rise in global mean atmospheric temperature (2 °C and 4 °C in 2100 respectively). The climate scenarios contain contributions from changes in ocean density (global thermal expansion and local steric changes related to changing ocean dynamics) and changes in ocean mass (melting of mountain glaciers and ice caps, changes in the Greenland and Antarctic ice sheets, and (minor) terrestrial water-storage contributions). All major components depend on the global temperature rise achieved in the target periods considered. The resulting set of climate scenarios represents our best estimate of twenty-first century sea level rise in the northeast Atlantic Ocean, given the current understanding of the various contributions. For 2100, they yield a local rise of 30 to 55 cm and 40 to 80 cm for the moderate and large rise in global mean atmospheric temperature, respectively.     

Recent precipitation trends over the polar ice sheets

Summary Meteorological and glaciological analyses are integrated to examine the precipitation trends during the last three decades over the ice sheets covering Antarctica and Greenland. For Antarctica, the best data source is provided by glaciologically-measured trends of snow accumulation, and for limited sectors of East Antarctica consistency with precipitation amounts calculated from the atmospheric water balance equation is obtained. For Greenland, precipitation rates parameterized from atmospheric analyses yield the only comprehensive depiction. The precipitation rate over Antarctica appears to have increased by about 5% over a time period spanning the accumulation means for the 1955–65 to 1965–75 periods, while over Greenland it has decreased by about 15% since 1983 with a secondary increase over the southern part of the ice sheet starting in 1977. At the end of the 10-year overlapping period, the global sea-level impact of the precipitation changes over Antarctica dominates that for Greenland and yields a net ice-sheet precipitation contribution of roughly 0.02 mm yr^–1. These changes are likely due to marked variations in the cyclonic forcing affecting the ice sheets, but are only weakly reflected in the temperature regime, consistent with the episodic nature of cyclonic precipitation. These conclusions are not founded on high quality data bases. The importance of such changes for understanding global sea-level variations argues for a modest research effort to collect simultaneous meteorological and glaciological observations in order to describe and understand the current precipitation variations over both ice sheets. Some suggestions are offered for steps that could be taken.With 8 Figures     

High-resolution precipitation and temperature downscaling for glacier models

The spatial resolution gap between global or regional climate models and the requirements for local impact studies motivates the need for climate downscaling. For impact studies that involve glacier modelling, the sparsity or complete absence of climate monitoring activities within the regions of interest presents a substantial additional challenge. Downscaling methods for this application must be independent of climate observations and cannot rely on tuning to station data. We present new, computationally-efficient methods for downscaling precipitation and temperature to the high spatial resolutions required to force mountain glacier models. Our precipitation downscaling is based on an existing linear theory for orographic precipitation, which we modify for large study regions by including moist air tracking. Temperature is downscaled using an interpolation scheme that reconstructs the vertical temperature structure to estimate surface temperatures from upper air data. Both methods are able to produce output on km to sub-km spatial resolution, yet do not require tuning to station measurements. By comparing our downscaled precipitation (1 km resolution) and temperature (200 m resolution) fields to station measurements in southern British Columbia, we evaluate their performance regionally and through the annual cycle. Precipitation is improved by as much as 30% (median relative error) over the input reanalysis data and temperature is reconstructed with a mean bias of 0.5°C at locations with high vertical relief. Both methods perform best in mountainous terrain, where glaciers tend to be concentrated.     

Uncertainty in climate change projections: the role of internal variability

Uncertainty in future climate change presents a key challenge for adaptation planning. In this study, uncertainty arising from internal climate variability is investigated using a new 40-member ensemble conducted with the National Center for Atmospheric Research Community Climate System Model Version 3 (CCSM3) under the SRES A1B greenhouse gas and ozone recovery forcing scenarios during 2000–2060. The contribution of intrinsic atmospheric variability to the total uncertainty is further examined using a 10,000-year control integration of the atmospheric model component of CCSM3 under fixed boundary conditions. The global climate response is characterized in terms of air temperature, precipitation, and sea level pressure during winter and summer. The dominant source of uncertainty in the simulated climate response at middle and high latitudes is internal atmospheric variability associated with the annular modes of circulation variability. Coupled ocean-atmosphere variability plays a dominant role in the tropics, with attendant effects at higher latitudes via atmospheric teleconnections. Uncertainties in the forced response are generally larger for sea level pressure than precipitation, and smallest for air temperature. Accordingly, forced changes in air temperature can be detected earlier and with fewer ensemble members than those in atmospheric circulation and precipitation. Implications of the results for detection and attribution of observed climate change and for multi-model climate assessments are discussed. Internal variability is estimated to account for at least half of the inter-model spread in projected climate trends during 2005–2060 in the CMIP3 multi-model ensemble.     

近30a青藏高原气候与冰川变化中的两种特殊现象

近30 a全球强烈变暖,水循环加快,冰川也加剧退缩。青藏高原以其特殊的地理位置与下垫面,既对全球变暖有正常的反应,也出现了异常特殊现象。这种特殊现象已发现两处：1) 青藏高原北部偏西冰芯记录降温0.6℃,相应的冰川退缩微弱,融水径流降低;2) 青藏高原东南部以岗日嘎布山区为代表,出现较多的冰川前进,可能指示降水量有较大的增加。上述事实指示气候变化与冰川响应的复杂性。     

Present and future climates of the Greenland ice sheet according to the IPCC AR4 models

The atmosphere?Cocean general circulation models (AOGCMs) used for the IPCC 4th Assessment Report (IPCC AR4) are evaluated for the Greenland ice sheet (GrIS) current climate modelling. The most suited AOGCMs for Greenland climate simulation are then selected on the basis of comparison between the 1970?C1999 outputs of the Climate of the twentieth Century experiment (20C3M) and reanalyses (ECMWF, NCEP/NCAR). This comparison indicates that the representation quality of surface parameters such as temperature and precipitation are highly correlated to the atmospheric circulation (500?hPa geopotential height) and its interannual variability (North Atlantic oscillation). The outputs of the three most suitable AOGCMs for present-day climate simulation are then used to assess the changes estimated by three IPCC greenhouse gas emissions scenarios (SRES) over the GrIS for the 2070?C2099 period. Future atmospheric circulation changes are projected to dampen the zonal flow, enhance the meridional fluxes and therefore provide additional heat and moisture to the GrIS, increasing temperature over the whole ice sheet and precipitation over its northeastern area. We also show that the GrIS surface mass balance anomalies from the SRES A1B scenario amount to ?300?km³/year with respect to the 1970?C1999 period, leading to a global sea-level rise of 5?cm by the end of the 21st century. This work can help to select the boundaries conditions for AOGCMs-based downscaled future projections.     

针对冰盖的定向地球工程研究

针对冰盖的定向地球工程研究旨在增强冰盖稳定性和减缓冰盖物质流失,从源头上减少冰盖对海平面上升的贡献,有望为应对气候变化和保护海岸线争取几百年的时间。冰盖地球工程主要作用在冰底和冰架-海洋界面上,主要途径如下：(1)排干或冻结冰盖底部水来干燥冰床,增强冰盖底部摩擦力;(2)在海洋中建造人造岛来支撑漂浮的冰架;(3)在冰架前端建造水下隔离墙,阻止温暖的海水到达冰川底部以减缓其融化。冰盖地球工程包括数值模拟、方案设计、工程试验和政治法律等诸多方面的研究。国际上的研究团队正在开展数值模拟和方案设计方面的研究,工程试验和政治法律等方面的研究尚未起步。预计工程试验的难度阶梯很可能是从实验室试验开始,到小尺度的野外试验,接着到格陵兰冰盖的入海冰川,最后到南极冰盖的入海冰川。针对冰盖的定向地球工程研究很有可能成为21世纪全球变化领域新兴的研究方向。     

Hydrological response to climate change in a glacierized catchment in the Himalayas

The analysis of climate change impact on the hydrology of high altitude glacierized catchments in the Himalayas is complex due to the high variability in climate, lack of data, large uncertainties in climate change projection and uncertainty about the response of glaciers. Therefore a high resolution combined cryospheric hydrological model was developed and calibrated that explicitly simulates glacier evolution and all major hydrological processes. The model was used to assess the future development of the glaciers and the runoff using an ensemble of downscaled climate model data in the Langtang catchment in Nepal. The analysis shows that both temperature and precipitation are projected to increase which results in a steady decline of the glacier area. The river flow is projected to increase significantly due to the increased precipitation and ice melt and the transition towards a rain river. Rain runoff and base flow will increase at the expense of glacier runoff. However, as the melt water peak coincides with the monsoon peak, no shifts in the hydrograph are expected.     

Greenland��s contribution to global sea-level rise by the end of the 21st century

The Greenland ice sheet holds enough water to raise the global sea level with ??7 m. Over the last few decades, observations manifest a substantial increase of the mass loss of this ice sheet. Both enhanced melting and increase of the dynamical discharge, associated with calving at the outlet-glacier fronts, are contributing to the mass imbalance. Using a dynamical and thermodynamical ice-sheet model, and taking into account speed up of outlet glaciers, we estimate Greenland??s contribution to the 21st-century global sea-level rise and the uncertainty of this estimate. Boundary fields of temperature and precipitation extracted from coupled climate-model projections used for the IPCC Fourth Assessment Report, are applied to the ice-sheet model. We implement a simple parameterization for increased flow of outlet glaciers, which decreases the bias of the modeled present-day surface height. It also allows for taking into account the observed recent increase in dynamical discharge, and it can be used for future projections associated with outlet-glacier speed up. Greenland contributes 0?C17?cm to global sea-level rise by the end of the 21st century. This range includes the uncertainties in climate-model projections, the uncertainty associated with scenarios of greenhouse-gas emissions, as well as the uncertainties in future outlet-glacier discharge. In addition, the range takes into account the uncertainty of the ice-sheet model and its boundary fields.     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

The regional nature of global challenges: a need and strategy for integrated regional modeling

Glaciers of the conterminous United States have been receding for the past century. Since 1900 the recession has varied from a 24 % loss in area (Mt. Rainier, Washington) to a 66 % loss in the Lewis Range of Montana. The rates of retreat are generally similar with a rapid loss in the early decades of the 20th century, slowing in the 1950s–1970s, and a resumption of rapid retreat starting in the 1990s. Decadal estimates of changes in glacier area for a subset of 31 glaciers from 1900 to 2000 are used to test a snow water equivalent model that is subsequently employed to examine the effects of temperature and precipitation variability on annual glacier area changes for these glaciers. Model results indicate that both winter precipitation and winter temperature have been important climatic factors affecting the variability of glacier variability during the 20th Century. Most of the glaciers analyzed appear to be more sensitive to temperature variability than to precipitation variability. However, precipitation variability is important, especially for high elevation glaciers. Additionally, glaciers with areas greater than 1 km² are highly sensitive to variability in temperature.     

Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble

The recent increase in the rate of the Greenland ice sheet melting has raised with urgency the question of the impact of such a melting on the climate. As former model projections, based on a coarse representation of the melting, show very different sensitivity to this melting, it seems necessary to consider a multi-model ensemble to tackle this question. Here we use five coupled climate models and one ocean-only model to evaluate the impact of 0.1 Sv (1 Sv = 10⁶ m³/s) of freshwater equally distributed around the coast of Greenland during the historical era 1965–2004. The ocean-only model helps to discriminate between oceanic and coupled responses. In this idealized framework, we find similar fingerprints in the fourth decade of hosing among the models, with a general weakening of the Atlantic Meridional Overturning Circulation (AMOC). Initially, the additional freshwater spreads along the main currents of the subpolar gyre. Part of the anomaly crosses the Atlantic eastward and enters into the Canary Current constituting a freshwater leakage tapping the subpolar gyre system. As a consequence, we show that the AMOC weakening is smaller if the leakage is larger. We argue that the magnitude of the freshwater leakage is related to the asymmetry between the subpolar-subtropical gyres in the control simulations, which may ultimately be a primary cause for the diversity of AMOC responses to the hosing in the multi-model ensemble. Another important fingerprint concerns a warming in the Nordic Seas in response to the re-emergence of Atlantic subsurface waters capped by the freshwater in the subpolar gyre. This subsurface heat anomaly reaches the Arctic where it emerges and induces a positive upper ocean salinity anomaly by introducing more Atlantic waters. We found similar climatic impacts in all the coupled ocean–atmosphere models with an atmospheric cooling of the North Atlantic except in the region around the Nordic Seas and a slight warming south of the equator in the Atlantic. This meridional gradient of temperature is associated with a southward shift of the tropical rains. The free surface models also show similar sea-level fingerprints notably with a comma-shape of high sea-level rise following the Canary Current.     

Estimates of twenty-first century sea-level changes for Norway

In this work we establish a framework for estimating future regional sea-level changes for Norway. Following recently published works, we consider how different physical processes drive non-uniform sea-level changes by accounting for spatial variations in (1) ocean density and circulation (2) ice and ocean mass changes and associated gravitational effects on sea level and (3) vertical land motion arising from past surface loading change and associated gravitational effects on sea level. An important component of past and present sea-level change in Norway is glacial isostatic adjustment. Central to our study, therefore, is a reassessment of vertical land motion using a far larger set of new observations from a permanent GNSS network. Our twenty-first century sea-level estimates are split into two parts. Firstly, we show regional projections largely based on findings from the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4) and dependent on the emission scenarios A2, A1B and B1. These indicate that twenty-first century relative sea-level changes in Norway will vary between ?0.2 to 0.3 m (1-sigma ± 0.13 m). Secondly, we explore a high-end scenario, in which a global atmospheric temperature rise of up to 6 °C and emerging collapse for some areas of the Antarctic ice sheets are assumed. Using this approach twenty-first century relative sea-level changes in Norway are found to vary between 0.25 and 0.85 m (min/max ± 0.45 m). We attach no likelihood to any of our projections owing to the lack of understanding of some of the processes that cause sea-level change.