The Open Global Glacier Model (OGGM) v1.1 Maussion, Fabien; Butenko, Anton; Champollion, Nicolas ...
Geoscientific Model Development,
03/2019, Letnik:
12, Številka:
3
Journal Article
Recenzirano
Odprti dostop
Despite their importance for sea-level rise, seasonal water availability, and
as a source of geohazards, mountain glaciers are one of the few remaining
subsystems of the global climate system for ...which no globally applicable,
open source, community-driven model exists. Here we present the Open Global
Glacier Model (OGGM), developed to provide a modular and open-source
numerical model framework for simulating past and future change of any
glacier in the world. The modeling chain comprises data downloading tools
(glacier outlines, topography, climate, validation data), a preprocessing
module, a mass-balance model, a distributed ice thickness estimation model,
and an ice-flow model. The monthly mass balance is obtained from gridded
climate data and a temperature index melt model. To our knowledge, OGGM is
the first global model to explicitly simulate glacier dynamics: the model
relies on the shallow-ice approximation to compute the depth-integrated flux
of ice along multiple connected flow lines. In this paper, we describe and
illustrate each processing step by applying the model to a selection of
glaciers before running global simulations under idealized climate forcings.
Even without an in-depth calibration, the model shows very realistic
behavior. We are able to reproduce earlier estimates of global glacier volume
by varying the ice dynamical parameters within a range of plausible values.
At the same time, the increased complexity of OGGM compared to other
prevalent global glacier models comes at a reasonable computational cost:
several dozen glaciers can be simulated on a personal computer, whereas
global simulations realized in a supercomputing environment take up to a few
hours per century. Thanks to the modular framework, modules of various
complexity can be added to the code base, which allows for new kinds of model
intercomparison studies in a controlled environment. Future developments will
add new physical processes to the model as well as automated calibration
tools. Extensions or alternative parameterizations can be easily added by the
community thanks to comprehensive documentation. OGGM spans a wide range of
applications, from ice–climate interaction studies at millennial timescales
to estimates of the contribution of glaciers to past and future sea-level
change. It has the potential to become a self-sustained community-driven
model for global and regional glacier evolution.
Damage to speleothems is a common phenomenon in mid-latitude caves, and multiple causes have been proposed. Here we report on one of such type of damage, namely stalagmites that are broken and ...partially sheared near their base but are still in upright position. Such stalagmites occur in the Obir Caves (Austria) associated with cryogenic cave carbonates, demonstrating the former presence of cave ice.
Th dating suggests damage to the speleothems during the Last Glacial Maximum. Numerical modelling combined with laboratory measurements demonstrates that internal deformation within a cave ice body cannot fracture stalagmites, even on a steep slope. Instead, temperature changes lead to thermoelastic stresses within an ice body that reach values equaling to and exceeding the tensile strength of even large stalagmites. Differences in thermal expansion coefficients cause a sharp vertical jump in stress between the stalagmite and the surrounding ice body, and the ice lifts the stalagmite as it expands with increasing temperature. This study refutes the previously accepted model that flow of ice breaks stalagmites, and suggests a link between glacial climate variability and corresponding cooling and warming cycles in the subsurface that weaken and eventually fracture stalagmites due to the opposing thermoelastic properties of calcite and ice.
Non-surface mass balance is non-negligible for glaciers in Iceland. Several Icelandic glaciers are in the neo-volcanic zone where a combination of geothermal activity, volcanic eruptions and ...geothermal heat flux much higher than the global average lead to basal melting close to 150 mm w.e. a−1 for the Mýrdalsjökull ice cap and 75 mm w.e. a−1 for the largest ice cap, Vatnajökull. Energy dissipation in the flow of water and ice is also rather large for the high-precipitation, temperate glaciers of Iceland resulting in internal and basal melting of 20–150 mm w.e. a−1. The total non-surface melting of glaciers in Iceland in 1995–2019 was 45–375 mm w.e. a−1 on average for the main ice caps, and was largest for Mýrdalsjökull, the south side of Vatnajökull and Eyjafjallajökull. Geothermal melting, volcanic eruptions and the energy dissipation in the flow of water and ice, as well as calving, all contribute, and thus these components should be considered in mass-balance studies. For comparison, the average mass balance of glaciers in Iceland since 1995 is −500 to −1500 mm w.e. a−1. The non-surface mass balance corresponds to a total runoff contribution of 2.1 km3 a−1 of water from Iceland.
In the European Alps, the Younger Dryas (YD) was
characterised by the last major glacier advance, with equilibrium line
altitudes being ∼ 220 to 290 m lower than during the Little
Ice Age, and also ...by the development of rock glaciers. Dating of these
geomorphic features, however, is associated with substantial uncertainties,
leading to considerable ambiguities regarding the internal structure of this
stadial, which is the most intensively studied one of the last glacial period. Here,
we provide robust physical evidence based on 230Th-dated cryogenic cave
carbonates (CCCs) from a cave located at 2274 m a.s.l. in the Dolomites of
northern Italy coupled with thermal modelling, indicating that early YD
winters were only moderately cold in this part of the Alps. More precisely,
we find that the mean annual air temperature dropped ≤ 3 ∘C
at the Allerød–YD transition. Our data suggest that autumns and early
winters in the early part of the YD were relatively snow-rich, resulting in
stable winter snow cover. The latter insulated the shallow subsurface in
winter and allowed the cave interior to remain close to the freezing point
(0 ∘C) year-round, promoting CCC formation. The main phase of CCC
precipitation at ∼ 12.2 ka coincided with the mid-YD
transition recorded in other archives across Europe. Based on thermal
modelling we propose that CCC formation at ∼ 12.2 ka was
most likely associated with a slight warming of approximately
+1 ∘C in conjunction with drier autumns and early winters in the
second half of the YD. These changes triggered CCC formation in this Alpine
cave as well as ice glacier retreat and rock glacier expansion across the
Alps.
Glaciers with extensive surface debris cover respond differently to climate forcing than those without supraglacial debris. In order to include debris-covered glaciers in projections of glaciogenic ...runoff and sea level rise and to understand the paleoclimate proxy recorded by such glaciers, it is necessary to understand the manner and timescales over which a supraglacial debris cover develops. Because debris is delivered to the glacier by processes that are heterogeneous in space and time, and these debris inclusions are altered during englacial transport through the glacier system, correctly determining where, when and how much debris is delivered to the glacier surface requires knowledge of englacial transport pathways and deformation. To achieve this, we present a model of englacial debris transport in which we couple an advection scheme to a full-Stokes ice flow model. The model performs well in numerical benchmark tests, and we present both 2-D and 3-D glacier test cases that, for a set of prescribed debris inputs, reproduce the englacial features, deformation thereof and patterns of surface emergence predicted by theory and observations of structural glaciology. In a future step, coupling this model to (i) a debris-aware surface mass balance scheme and (ii) a supraglacial debris transport scheme will enable the co-evolution of debris cover and glacier geometry to be modelled.
A method is described to estimate the thickness of glacier ice using information derived from the measured ice extent, surface topography, surface mass balance, and rate of thinning or thickening of ...the ice column. Shear stress beneath an ice column is assumed to be simply related to ice thickness and surface slope, as for an inclined slab, but this calculation is cast as a linear optimization problem so that a smoothness regularization can be applied. Assignment of bed stress is based on the flow law for ice and a mass balance calculation but must be preceded by delineation of the ice flow drainage basin. Validation of the method is accomplished by comparing thickness estimates to the known thickness generated by a numerical ice dynamics model. Once validated, the method is used to estimate the subglacial topography for all glaciers in western Canada that lie south of 60°N. Adding the present ice volume of each glacier gives the estimated total volume as 2320 km³, equivalent to 5.8 mm of sea level rise. Taking the glaciated area as 26 590 km² gives the average glacier thickness as 87.2 m. A detailed error analysis indicates that systematic errors are likely to increase the estimated sea level rise and when random errors are included the combined result is 6.3 ± 0.6 mmor, expressed as ice volume, 2530 ± 220 km³.
The climate of high midlatitude mountains appears to be warming faster than the global average, but evidence for such elevation-dependent warming (EDW) at higher latitudes is presently scarce. Here, ...we use a comprehensive network of remote meteorological stations, proximal radiosonde measurements, downscaled temperature reanalysis, ice cores, and climate indices to investigate the manifestation and possible drivers of EDW in the St. Elias Mountains in subarctic Yukon, Canada. Linear trend analysis of comprehensively validated annual downscaled North American Regional Reanalysis (NARR) gridded surface air temperatures for the years 1979–2016 indicates a warming rate of 0.028°C a−1 between 5500 and 6000 m above mean sea level (MSL), which is ∼1.6 times larger than the global-average warming rate between 1970 and 2015. The warming rate between 5500 and 6000 m MSL was ∼1.5 times greater than the rate at the 2000–2500 m MSL bin (0.019°C a−1), which is similar to the majority of warming rates estimated worldwide over similar elevation gradients. Accelerated warming since 1979, measured by radiosondes, indicates a maximum rate at 400 hPa (~7010 m MSL). EDW in the St. Elias region therefore appears to be driven by recent warming of the free troposphere. MODIS satellite data show no evidence for an enhanced snow albedo feedback above 2500 m MSL, and declining trends in sulfate aerosols deposited in high-elevation ice cores suggest a modest increase in radiative forcing at these elevations. In contrast, increasing trends in water vapor mixing ratio at the 500-hPa level measured by radiosonde suggest that a longwave radiation vapor feedback is contributing to EDW.
Geothermal heat sources beneath glaciers and ice caps influence local ice-dynamics and mass balance but also control ice surface depression evolution as well as subglacial water reservoir dynamics. ...Resulting jökulhlaups (i.e., glacier lake outburst floods) impose danger to people and infrastructure, especially in Iceland, where they are closely monitored. Due to hundreds of meters of ice, direct measurements of heat source strength and extent are not possible. We present an indirect measurement method which utilizes ice flow simulations and glacier surface data, such as surface mass balance and surface depression evolution. Heat source locations can be inferred accurately to simulation grid scales; heat source strength and spatial distributions are also well quantified. Our methods are applied to the Mýrdalsjökull ice cap in Iceland, where we are able to refine previous heat source estimates.
Sub-meter resolution, stereoscopic satellite images allow for the generation of accurate and high-resolution digital elevation models (DEMs) over glaciers and ice caps. Here, repeated stereo images ...of Drangajökull ice cap (NW Iceland) from Pléiades and WorldView2 (WV2) are combined with in situ estimates of snow density and densification of firn and fresh snow to provide the first estimates of the glacier-wide geodetic winter mass balance obtained from satellite imagery. Statistics in snow- and ice-free areas reveal similar vertical relative accuracy (< 0.5 m) with and without ground control points (GCPs), demonstrating the capability for measuring seasonal snow accumulation. The calculated winter (14 October 2014 to 22 May 2015) mass balance of Drangajökull was 3.33 ± 0.23 m w.e. (meter water equivalent), with ∼ 60 % of the accumulation occurring by February, which is in good agreement with nearby ground observations. On average, the repeated DEMs yield 22 % less elevation change than the length of eight winter snow cores due to (1) the time difference between in situ and satellite observations, (2) firn densification and (3) elevation changes due to ice dynamics. The contributions of these three factors were of similar magnitude. This study demonstrates that seasonal geodetic mass balance can, in many areas, be estimated from sub-meter resolution satellite stereo images.
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.