Glaciers in the Earth’s Hydrological Cycle: Assessments of Glacier Mass and Runoff Changes on Global and Regional Scales

Changes in mass contained by mountain glaciers and ice caps can modify the Earth’s hydrological cycle on multiple scales. On a global scale, the mass loss from glaciers contributes to sea-level rise. On regional and local scales, glacier meltwater is an important contributor to and modulator of river flow. In light of strongly accelerated worldwide glacier retreat, the associated glacier mass losses raise concerns over the sustainability of water supplies in many parts of the world. Here, we review recent attempts to quantify glacier mass changes and their effect on river runoff on regional and global scales. We find that glacier runoff is defined ambiguously in the literature, hampering direct comparison of findings on the importance of glacier contribution to runoff. Despite consensus on the hydrological implications to be expected from projected future warming, there is a pressing need for quantifying the associated regional-scale changes in glacier runoff and responses in different climate regimes.     

New fjords,new coasts,new landscapes: The geomorphology of paraglacial coasts formed after recent glacier retreat in Brepollen (Hornsund,southern Svalbard)

Changes in the properties and dynamics of tidewater glacier systems are key indicators of the state of Arctic climate and environment. Calving of tidewater glacier fronts is currently the dominant form of ice mass loss and a major contributor to global sea-level rise. An important yet under-studied aspect of this process is transformation of Arctic landscapes, where new lands and coastal systems are revealed due to the recession of marine-terminating ice masses. The evolution of those freshly exposed paraglacial coastal environments is controlled by nearshore marine, coastal and terrestrial geomorphic processes, which rework glacial-derived sediments to create new coastal paraglacial landforms and landscapes. Here, we present the first study of the paraglacial coasts of Brepollen, one of the youngest bays of Svalbard revealed by ice retreat. We describe and classify coastal systems and the variety of landforms (deltas, cliffs, tidal flats, beaches) developed along the shores of Brepollen during the last 100 years. We further discuss the main modes of sediment supply to the coast in different parts of the new bay, highlighting the fast rate of coastal transformation as a paraglacial response to rapid deglaciation in the Arctic. This study provides an exemplar of likely coastal responses to be anticipated in similar tidewater settings under future climate change. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Sea-level rise: A review of recent past and near-future trends

Global mean sea level is a potentially sensitive indicator of climate change. Global warming will contribute to worldwide sea-level rise (SLR) from thermal expansion of ocean water, melting of mountain glaciers and polar ice sheets. A number of studies, mostly using tide-gauge data from the Permanent Service for Mean Sea Level, Bidston Observatory, England, have obtained rates of global SLR within the last 100 years that range between 0·3 and 3 mm yr^?1, with most values concentrated between 1 and 2 mm yr^?1. However, the reliability of these results has been questioned because of problems with data quality and physical processes that introduce a high level of spatial and temporal variability. Sources of uncertainty in the sea-level data include variations in winds, ocean currents, river runoff, vertical earth movements, and geographically uneven distribution of long-term records. Crustal motions introduce a major source of error. To a large extent, these can be filtered by employing palaeo-sea-level proxies, and geophysical modelling to remove glacio-isostatic changes. Ultimately, satellite geodesy will help resolve the inherent ambiguity between the land and ocean level changes recorded by tide gauges. Future sea level is expected to rise by ～ 1 m, with a ‘best-guess’ value of 48 cm by the year 2100. Such rates represent an acceleration of four to seven times over present rates. Local land subsidence could substantially increase the apparent SLR. For example, Louisiana is currently experiencing SLR trends nearly 10 times the global mean rate. These recently reduced SLR estimates are based on climate models that predict a zero to negative contribution to SLR from Antarctica. Most global climate models (GCMs) indicate an ice accumulation over Antarctica, because in a warmer world, precipitation will exceed ablation/snow-melt. However, the impacts of attritional processes, such as thinning of the ice shelves, have been downplayed according to some experts. Furthermore, not all climate models are in agreement. Opposite conclusions may be drawn from the results of other GCMs. In addition, the West Antarctic Ice Sheet is potentially subject to dynamic and volcanic instabilities that are difficult to predict. Because of the great uncertainty in SLR projections, careful monitoring of future sea-level trends by upgraded tide-gauge networks and satellite geodesy will become essential. Finally, because of the high spatial variability in crustal subsidence rates, wave climates and tidal regimes, it will be the set of local conditions (especially the relative sea-level rise), rather than a single global mean sea-level trend, that will determine each locality's vulnerability to future SLR.     

Present State and Prospects of Ice Sheet and Glacier Modelling

Since the late 1970s, numerical modelling has become established as an important technique for the understanding of ice sheet and glacier dynamics, and several models have been developed over the years. Ice sheet models are particularly relevant for predicting the possible response of ice sheets to climate change. Recent observations suggest that ice dynamics could play a crucial role for the contribution of ice sheets to future sea level rise under global warming conditions, and the need for further research into the matter was explicitly stated in the Fourth Assessment Report (AR4) of the United Nations Intergovernmental Panel on Climate Change (IPCC). In this paper, we review the state of the art and current problems of ice sheet and glacier modelling. An outline of the underlying theory is given, and crucial processes (basal sliding, calving, interaction with the solid Earth) are discussed. We summarise recent progress in the development of ice sheet and glacier system models and their coupling to climate models, and point out directions for future work.     

Estimating the Glacier Contribution to Sea-Level Rise for the Period 1800�C2005

In this study, a new estimate of the contribution of glaciers and ice caps to the sea-level rise over the period 1800?C2005 is presented. We exploit the available information on changes in glacier length. Length records form the only direct evidence of glacier change that has potential global coverage before 1950. We calculate a globally representative signal from 349 glacier length records. By means of scaling, we deduce a global glacier volume signal, that is calibrated on the mass-balance and geodetic observations of the period 1950?C2005. We find that the glacier contribution to sea-level rise was 8.4 ± 2.1 cm for the period 1800?C2005 and 9.1 ± 2.3 cm for the period 1850?C2005.     

Impacts of climate change in three hydrologic regimes in British Columbia,Canada

Hydrologic modelling has been applied to assess the impacts of projected climate change within three study areas in the Peace, Campbell and Columbia River watersheds of British Columbia, Canada. These study areas include interior nival (two sites) and coastal hybrid nival–pluvial (one site) hydro‐climatic regimes. Projections were based on a suite of eight global climate models driven by three emission scenarios to project potential climate responses for the 2050s period (2041–2070). Climate projections were statistically downscaled and used to drive a macro‐scale hydrology model at high spatial resolution. This methodology covers a large range of potential future climates for British Columbia and explicitly addresses both emissions and global climate model uncertainty in the final hydrologic projections. Snow water equivalent is projected to decline throughout the Peace and Campbell and at low elevations within the Columbia. At high elevations within the Columbia, snow water equivalent is projected to increase with increased winter precipitation. Streamflow projections indicate timing shifts in all three watersheds, predominantly because of changes in the dynamics of snow accumulation and melt. The coastal hybrid site shows the largest sensitivity, shifting to more rainfall‐dominated system by mid‐century. The two interior sites are projected to retain the characteristics of a nival regime by mid‐century, although streamflow‐timing shifts result from increased mid‐winter rainfall and snowmelt, and earlier freshet onset. Copyright © 2012 John Wiley & Sons, Ltd.     

Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change

Glaciers have strongly contributed to sea-level rise during the past century and will continue to be an important part of the sea-level budget during the twenty-first century. Here, we review the progress in estimating global glacier mass change from in situ measurements of mass and length changes, remote sensing methods, and mass balance modeling driven by climate observations. For the period before the onset of satellite observations, different strategies to overcome the uncertainty associated with monitoring only a small sample of the world’s glaciers have been developed. These methods now yield estimates generally reconcilable with each other within their respective uncertainty margins. Whereas this is also the case for the recent decades, the greatly increased number of estimates obtained from remote sensing reveals that gravimetry-based methods typically arrive at lower mass loss estimates than the other methods. We suggest that strategies for better interconnecting the different methods are needed to ensure progress and to increase the temporal and spatial detail of reliable glacier mass change estimates.     

Estimation of uncertainty sources in the projections of Lithuanian river runoff

Particular attention is given to the reliability of hydrological modelling results. The accuracy of river runoff projection depends on the selected set of hydrological model parameters, emission scenario and global climate model. The aim of this article is to estimate the uncertainty of hydrological model parameters, to perform sensitivity analysis of the runoff projections, as well as the contribution analysis of uncertainty sources (model parameters, emission scenarios and global climate models) in forecasting Lithuanian river runoff. The impact of model parameters on the runoff modelling results was estimated using a sensitivity analysis for the selected hydrological periods (spring flood, winter and autumn flash floods, and low water). During spring flood the results of runoff modelling depended on the calibration parameters that describe snowmelt and soil moisture storage, while during the low water period—the parameter that determines river underground feeding was the most important. The estimation of climate change impact on hydrological processes in the Merkys and Neris river basins was accomplished through the combination of results from A1B, A2 and B1 emission scenarios and global climate models (ECHAM5 and HadCM3). The runoff projections of the thirty-year periods (2011–2040, 2041–2070, 2071–2100) were conducted applying the HBV software. The uncertainties introduced by hydrological model parameters, emission scenarios and global climate models were presented according to the magnitude of the expected changes in Lithuanian rivers runoff. The emission scenarios had much greater influence on the runoff projection than the global climate models. The hydrological model parameters had less impact on the reliability of the modelling results.     

Projected 21st century climate change on snow conditions over Shasta Dam watershed by means of dynamical downscaling

Snow is an important component of the Earth's climate system and is particularly vulnerable to global warming. It has been suggested that warmer temperatures may cause significant declines in snow water content and snow cover duration. In this study, snowfall and snowmelt were projected by means of a regional climate model that was coupled to a physically based snow model over Shasta Dam watershed to assess changes in snow water content and snow cover duration during the 21st century. This physically based snow model requires both physical data and future climate projections. These physical data include topography, soils, vegetation, and land use/land cover, which were collected from associated organizations. The future climate projections were dynamically downscaled by means of the regional climate model under 4 emission scenarios simulated by 2 general circulation models (fifth‐generation of the ECHAM general circulation model and the third‐generation atmospheric general circulation model). The downscaled future projections were bias corrected before projecting snowfall and snowmelt processes over Shasta Dam watershed during 2010–2099. This study's results agree with those of previous studies that projected snow water equivalent is decreasing by 50–80% whereas the fraction of precipitation falling as snowfall is decreasing by 15% to 20%. The obtained projection results show that future snow water content will change in both time and space. Furthermore, the results confirm that physical data such as topography, land cover, and atmospheric–hydrologic data are instrumental in the studies on the impact of climate change on the water resources of a region.     

Prospective scenarios of the saltwater intrusion in an estuary under climate change context using Bayesian neural networks

In recent decades, saltwater intrusion over some low-lying coastal regions was deteriorated by rising sea-level and decreasing streamflow in the context of climate change. Though physically-based hydrodynamic models are the most detailed means to simulate salinity processes, they are commonly restricted by data insufficiency issues both in spatial resolution and temporal lasting. This motivates us to build a statistical model enable simulation and scenario analysis for coastal salinity change with limited observations. A Bayesian neural network (BNN) model is built hereby to simulate salinity. It offers more precise estimation compared with the conventional artificial neural network. Meanwhile, the model gives the uncertainty behaviors of the final salinity simulation which is not available for other methods. Future scenarios of salinity change are constructed and analyzed in different time periods on the basis of the validated BNN model. Results indicate that the water quality over lower Pearl River is degrading along with more significant uncertainties. Further analysis suggests that streamflow alteration has a more direct impact on salinity variations than the sea-level change does. The method allows a profound analysis of the potential influence on water quality degradation in coastal and low-lying regions in support of water management and adaptation toward global climate change.     

Projection of glacier runoff in Yarkant River basin and Beida River basin,Western China

Potential changes in glacier area, mass balance and runoff in the Yarkant River Basin (YRB) and Beida River Basin (BRB) are projected for the period from 2011 to 2050 employing the modified monthly degree‐day model forced by climate change projection. Future monthly air temperature and precipitation were derived from the simple average of 17, 16 and 17 General Circulation Model (GCM) projections following the A1B, A2 and B1 scenarios, respectively. These data were downscaled to each station employing the Delta method, which computes differences between current and future GCM simulations and adds these changes to observed time series. Model parameters calibrated with observations or results published in the literature between 1961 and 2006 were kept unchanged. Annual glacier runoff in YRB is projected to increase until 2050, and the total runoff over glacier area in 1970 is projected to increase by about 13%–35% during 2011–2050 relative to the average during 1961–2006. Annual glacier runoff and the total runoff over glacier area in 1970 in BRB is projected to increase initially and then to reach a tipping point during 2011–2030. There are prominent increases in summer, but only small increase in May and October of glacier runoff in YRB, and significant increases during late spring and early summer and significant decreases in July and late summer of glacier runoff in BRB. This study highlights the great differences among basins in their response to future climate warming. The specific runoff from areas exposed after glacier retreat relative to 1970 is projected to general increasing, which must be considered when evaluating the potential change of glacier runoff. Copyright © 2011 John Wiley & Sons, Ltd.     

基于偏移追踪的青藏高原卓琼冰川流速监测及成因分析

冰川作为地球系统中重要的组成部分,是全球气候变化的敏感指示器和调节器.冰川运动的遥感监测也是进行冰川研究的重要内容.本文主要采用偏移追踪的方法,利用2003-2010年期间7对雷达单视复数数据监测了青藏高原珠穆朗玛峰地区的卓琼冰川运动情况.经验证,本研究获得的冰川流速结果可靠.卓琼冰川的流向自西向东,流速自冰川末端向上游积累区逐渐增加,流速大小主要来自距离向的贡献.研究结果显示,卓琼冰川在2003-2005年期间流速较大,最大流速达到45 m·a~(-1),而在2005-2010年期间流速有所降低,最大流速在35~40 m·a~(-1)范围内波动;卓琼冰川年际间流速变化基本一致,冰川中部存在流速突变情况,突变幅度为3~7 m·a~(-1),但该尺度的突变并不会影响卓琼冰川的总体运动趋势.此外,本文还分析了气象以及地理位置因素对卓琼冰川运动的影响.     

Analysis of uncertainties in regional climate model outputs over the Czech Republic

In this study we present results of uncertainty analysis in eight regional climate model (RCM) outputs over the area of the Czech Republic. The RCM simulations come from the EU 5th Framework program project PRUDENCE (Prediction of Regional scenarios and Uncertainties for Defining EuropeaN Climate change risks and Effects). Using the analysis of variance we have found that the main source of uncertainty in projected changes of mean seasonal air temperature is the driving global climate model. In case of precipitation changes, the RCM is the largest source of uncertainty in all seasons except for the spring. With the second method, the Reliability Averaging method, we have focused on the uncertainty coming from the RCM itself. The results of both methods showed that the relative contribution of the regional climate model to the uncertainty of simulated mean seasonal air temperature and precipitation changes is largest in summer and smallest in winter.     

Modeling Groundwater Depletion at Regional and Global Scales: Present State and Future Prospects

Except for frozen water in ice and glaciers, groundwater is the world’s largest distributed store of freshwater and has strategic importance to global food and water security. In this paper, the most recent advances quantifying groundwater depletion (GWD) are comprehensively reviewed. This paper critically evaluates the recently advanced modeling approaches estimating GWD at regional and global scales, and the evidence of feedbacks to the Earth system including sea-level rise associated with GWD. Finally, critical challenges and opportunities in the use of groundwater are identified for the adaption to growing food demand and uncertain climate.     

Downscaling IPCC control run and future scenario with focus on the Barents Sea

Global atmosphere-ocean general circulation models are the tool by which projections for climate changes due to radiative forcing scenarios have been produced. Further, regional atmospheric downscaling of the global models may be applied in order to evaluate the details in, e.g., temperature and precipitation patterns. Similarly, detailed regional information is needed in order to assess the implications of future climate change for the marine ecosystems. However, regional results for climate change in the ocean are sparse. We present the results for the circulation and hydrography of the Barents Sea from the ocean component of two global models and from a corresponding pair of regional model configurations. The global models used are the GISS AOM and the NCAR CCSM3. The ROMS ocean model is used for the regional downscaling of these results (ROMS-G and ROMS-N configurations, respectively). This investigation was undertaken in order to shed light on two questions that are essential in the context of regional ocean projections: (1) How should a regional model be set up in order to take advantage of the results from global projections; (2) What limits to quality in the results of regional models are imposed by the quality of global models? We approached the first question by initializing the ocean model in the control simulation by a realistic ocean analysis and specifying air-sea fluxes according to the results from the global models. For the projection simulation, the global models’ oceanic anomalies from their control simulation results were added upon initialization. Regarding the second question, the present set of simulations includes regional downscalings of the present-day climate as well as projected climate change. Thus, we study separately how downscaling changes the results in the control climate case, and how scenario results are changed. For the present-day climate, we find that downscaling reduces the differences in the Barents Sea between the original global models. Furthermore, the downscaled results are closer to observations. On the other hand, the downscaled results from the scenario simulations are significantly different: while the heat transport into the Barents Sea and the salinity distribution change modestly from control to scenario with ROMS-G, in ROMS-N the heat transport is much larger in the scenario simulation, and the water masses become much less saline. The lack of robustness in the results from the scenario simulations leads us to conclude that the results for the regional oceanic response to changes in the radiative forcing depend on the choice of AOGCM and is not settled. Consequently, the effect of climate change on the marine ecosystem of the Barents Sea is anything but certain.     

Interannual variability in glacier contribution to runoff from a high‐elevation Andean catchment: understanding the role of debris cover in glacier hydrology

We present a field‐data rich modelling analysis to reconstruct the climatic forcing, glacier response, and runoff generation from a high‐elevation catchment in central Chile over the period 2000–2015 to provide insights into the differing contributions of debris‐covered and debris‐free glaciers under current and future changing climatic conditions. Model simulations with the physically based glacio‐hydrological model TOPKAPI‐ETH reveal a period of neutral or slightly positive mass balance between 2000 and 2010, followed by a transition to increasingly large annual mass losses, associated with a recent mega drought. Mass losses commence earlier, and are more severe, for a heavily debris‐covered glacier, most likely due to its strong dependence on snow avalanche accumulation, which has declined in recent years. Catchment runoff shows a marked decreasing trend over the study period, but with high interannual variability directly linked to winter snow accumulation, and high contribution from ice melt in dry periods and drought conditions. The study demonstrates the importance of incorporating local‐scale processes such as snow avalanche accumulation and spatially variable debris thickness, in understanding the responses of different glacier types to climate change. We highlight the increased dependency of runoff from high Andean catchments on the diminishing resource of glacier ice during dry years.     

Sensitivity of physical lake processes to climate change within a large Precambrian Shield catchment

Potential future changes in lake physical processes (e.g. stratification and freezing) can be assessed through exploring their sensitivity to climate change, and assessing the current vulnerability of different lake types to plausible changes in meteorological drivers. This study quantifies the impacts of climate change and sensitivity of lake physical processes within a large (5100 km²) Precambrian Shield catchment in south‐central Ontario. Historic regional relationships are established between climate drivers, lake morphology, and lake physical changes through generalized linear modelling (GLM), and are used to quantify likely changes in timing of ice phenology and lake stratification across 72 lakes under a range of future climate models and scenarios. In response to projections of increased temperature (ensemble mean of +3.3 °C), both earlier ice‐off and onset of summer stratification were projected, with later ice‐on and fall turnover compared to the baseline. Process sensitivity to climate change varied by lake type; shallower lakes with a smaller volume (less than 15 m deep and less than 0.05 km³) were more sensitive to processes associated with lake heating (stratification onset and ice‐off), and deeper lakes with a larger surface area (greater than 30 m deep and greater than 1000 ha) were more sensitive to processes associated with lake cooling (fall turnover and ice‐on). These results indicate that whereas small lakes are vulnerable to climate warming because of changes that occur in spring and summer, larger lakes are particularly sensitive during the fall. The findings suggest that lake morphology and associated sensitivity should be considered in the development of sustainable lake management strategies. Copyright © 2016 John Wiley & Sons, Ltd.     

What dominates sea level at the coast: a case study for the Gulf of Guinea

Sea level variations and extreme events are a major threat for coastal zones. This threat is expected to worsen with time because low-lying coastal areas are expected to become more vulnerable to flooding and land loss as sea level rises in response to climate change. Sea level variations in the coastal ocean result from a combination of different processes that act at different spatial and temporal scales. In this study, the relative importance of processes causing coastal sea level variability at different time-scales is evaluated. Contributions from the altimetry-derived sea-level (including the sea level rise due to the ocean warming and land ice loss in response to climate change), dynamical atmospheric forcing induced sea level (surges), wave-induced run-up and set-up, and astronomical tides are estimated from observational datasets and reanalyses. As these processes impact the coast differently, evaluating their importance is essential for assessment of the local coastline vulnerability. A case study is developed in the Gulf of Guinea over the 1993–2012 period. The leading contributors to sea level variability off Cotonou differ depending on the time-scales considered. The trend is largely dominated by processes included in altimetric data and to a lesser extent by swell-waves run-up. The latter dominates interannual variations. Swell-waves run-up and tides dominate subannual variability. Extreme events are due to the conjunction of high tides and large swell run-up, exhibiting a clear seasonal cycle with more events in boreal summer and a trend mostly related to the trend in altimetric-derived sea-level.     

Arctic Climate and Water Change: Model and Observation Relevance for Assessment and Adaptation

The Arctic is subject to growing economic and political interest. Meanwhile, its climate and water systems are in rapid transformation. In this paper, we review and extend a set of studies on climate model results, hydro-climatic change, and hydrological monitoring systems. Results indicate that general circulation model (GCM) projections of drainage basin temperature and precipitation have improved between two model generations. However, some inaccuracies remain for precipitation projections. When considering geographical priorities for monitoring or adaptation efforts, our results indicate that future projections by GCMs and recent observations diverge regarding the basins where temperature and precipitation changes currently are the most pronounced and where they will be so in the future. Regarding late twentieth-century discharge changes in major Arctic rivers, data generally show excess of water relative to precipitation changes. This indicates a possible contribution to sea-level rise of river water that was previously stored in permafrost or groundwater. The river contribution to the increasing Arctic Ocean freshwater inflow is similar in magnitude to the separate contribution from glaciers, which underlines the importance of considering all possible sources of freshwater when assessing sea-level change. We further investigate monitoring systems and find a lack of harmonized water chemistry data, which limits the ability to understand the origin and transport of nutrients, carbon and sediment to the sea. To provide adequate information for research and policy, Arctic hydrological and hydrochemical monitoring needs to be extended, better integrated and made more accessible. Further water-focused data and modeling efforts are required to resolve the source of excess discharge in Arctic rivers. Finally, improvements in climate model parameterizations are needed, in particular for precipitation projections.     

Climate impacts on river flow: projections for the Medway catchment,UK, with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.