Calibrating a regional glacier model using post-LIA glacier length changes in the Alps

Author(s):  
Matthias Dusch ◽  
Kurt Nicolussi ◽  
Fabien Maussion

<p>We present an approach to calibrate a regional glacier model based on the well observed period since the mid-19<sup>th</sup> century LIA maximum.<br>We chose 30 glaciers distributed across the entire European Alps with frequent length change observations in that period. These glaciers account for 25% of today's total glacier area in the Alps. We run simulations with the Open Global Glacier Model (OGGM, https://oggm.org) driven by HISTALP (http://www.zamg.ac.at/histalp) gridded climate data. To calibrate the glaciers individually, we vary three model parameters within a reasonable range: (i) a precipitation scaling factor governing average mass-turnover and mass-balance profiles, (ii) the ice creep parameter governing basal sheer stress and the dynamics of ice flow, and (iii) a constant mass balance perturbation applied to the yearly mass-balance. This results in  1365 unique parameter combinations which were tested for all glaciers. We chose individual parameter subsets for every glacier based on objective criteria minimizing the difference between modeled and observed length changes.</p><p>We find that there is no unique parameter combination satisfying our criteria for all glaciers. It is also challenging to identify an ideal parameter combination for each individual glacier, since there is a trade-off between reproducing variability (useful for paleo-climate interpretations) and reproducing observed length change (useful for projections and planing).<br>Furthermore, model and input data uncertainties are variable in time, leading to non-unique optimal parameter sets. Therefore, we rely on an ensemble of simulations consisting of the best runs with respect to multiple statistical measures. Together with a cross-validation procedure, the ensemble produces a probabilistic uncertainty range which can be applied to Holocene glacier reconstructions and future evolution scenarios.</p>

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

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


1995 ◽  
Vol 21 ◽  
pp. 206-212 ◽  
Author(s):  
Wilfried Haeberli ◽  
Martin Hoelzle

A parameterization scheme using simple algorithms for unmeasured glaciers is being applied to glacier inventory data to estimate the basic glaciological characteristics of the inventoried ice bodies and simulate potential climate-change effects on mountain glaciers. For past and potential climate scenarios, glacier changes for assumed mass-balance changes are calculated as step functions between steady-state conditions for time intervals that approximately correspond to the characteristic dynamic response time (a few decades) of the glaciers. In order to test the procedure, a pilot study was carried out in the European Alps where detailed glacier inventories had been compiled around the mid-1970s. Total glacier volume in the Alps is estimated at about 130 km3 for the mid-1970s; strongly negative mass balances are likely to have caused a loss of about 10–20% of this total volume during the decade 1980–90. Backward calculation of glacier-length changes using a mean annual mass balance of 0.25m w.e.a−1 since the end of the “Little Ice Age” around 1850 AD gives considerable scatter but satisfactory overall results as compared with long-term observations. The total loss of Alpine surface ice mass since 1850 can be estimated at about half the original value. An acceleration of this development, with annual mass losses of around 1 m a−1 or more as anticipated from IPCC scenario A for the coming century, could eliminate major parts of the presently existing Alpine ice volume within decades.


1995 ◽  
Vol 21 ◽  
pp. 206-212 ◽  
Author(s):  
Wilfried Haeberli ◽  
Martin Hoelzle

A parameterization scheme using simple algorithms for unmeasured glaciers is being applied to glacier inventory data to estimate the basic glaciological characteristics of the inventoried ice bodies and simulate potential climate-change effects on mountain glaciers. For past and potential climate scenarios, glacier changes for assumed mass-balance changes are calculated as step functions between steady-state conditions for time intervals that approximately correspond to the characteristic dynamic response time (a few decades) of the glaciers. In order to test the procedure, a pilot study was carried out in the European Alps where detailed glacier inventories had been compiled around the mid-1970s. Total glacier volume in the Alps is estimated at about 130 km3 for the mid-1970s; strongly negative mass balances are likely to have caused a loss of about 10–20% of this total volume during the decade 1980–90. Backward calculation of glacier-length changes using a mean annual mass balance of 0.25m w.e.a−1 since the end of the “Little Ice Age” around 1850 AD gives considerable scatter but satisfactory overall results as compared with long-term observations. The total loss of Alpine surface ice mass since 1850 can be estimated at about half the original value. An acceleration of this development, with annual mass losses of around 1 m a−1 or more as anticipated from IPCC scenario A for the coming century, could eliminate major parts of the presently existing Alpine ice volume within decades.


2020 ◽  
Author(s):  
Luca Carturan ◽  
Aldino Bondesan ◽  
Alberto Carton ◽  
Federico Cazorzi ◽  
Sara Cucchiaro ◽  
...  

<p>The knowledge of past fluctuations of glaciers is the key for understanding their dynamics and their climate-related evolution. Glacier mass balance and length changes are the two metrics normally used for reconstructing past fluctuations series of glaciers. However, length change measurements series are often discontinuous and require validation, whereas mass balance measurements are available for only a few glaciers worldwide and only for the latest decades.</p><p>In the context of glacier reconstructions, other sources of information such as historical-archival, glacio-archaeological and geomorphological data are of critical importance, because they enable the completion and validation of direct measurement series and their extension into the past, providing spatial and temporal constraints.</p><p>A unique source of unexploited historical information dating back to the First World War (WWI, 1915-1918) exists for many glaciers in the Eastern Italian Alps. This information mainly consists of old photos, which however are spread over a multitude of sources, often difficult to access, and in many cases not yet digitized.</p><p>We propose a workflow that enables extracting quantitative information from terrestrial photographs taken during the WWI period, aimed at the reconstruction of glacier area, volume and firn lines by means of the monoplotting technique. This method relies on the availability of high-resolution digital elevation models, which became recently available over wide areas thanks to LiDAR and aerial photogrammetry. This work presents the methods applied, and the results obtained, on several case studies in the Adamello-Presanella, Ortles-Cevedale, Dolomites, and Julian Alps.    </p>


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

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


2020 ◽  
Author(s):  
Frederik Link ◽  
Georg Rümpker

<p>SKS-splitting measurements in the European Alps show an anisotropic fast axes parallel/subparallel relative to the mountain-belt. This indicates a mantle flow with a rotational component according to the orogeny under the assumption that the fast axes directly reflect the flow direction. This might be misleading due to a possible crustal contribution of anisotropy. Therefore, we isolate the crustal anisotropy using an improved receiver function method that accounts for anisotropic and structural properties.</p><p>The analysis for the crustal anisotropy is based on the stacking method proposed by Kaviani & Rümpker (2015). We modify their approach by introducing a time-selective splitting analysis of the crustal Ps- and PpPs-phases. The stacking is performed to the phases after correction of the anisotropic effect according to the model parameters H, the crustal thickness,  κ, the P-to S-wave velocity ratio, a, the percentage of anisotropy and φ, the fast axis orientation.</p><p>The Alps show a considerable Moho-topography due to its mountain root and its complex tectonic history. This can significantly deflect the crustal phases introducing a dominating appearance in the receiver functions. We therefore analyse for a dipping interface (not accounting for anisotropy) and then use an improved model in our analysis to infer the anisotropic properties of the crust.</p><p>Knowing the crustal anisotropic contribution we correct for this effect on the XKS-waveforms to isolate the anisotropy of the mantle. The remaining splitting shows an improved approximation of the flow patterns in the asthenosphere, while complexities might still imply an effect of the lithospheric mantle.</p><p>We apply our approach to stations of the AlpArray network resulting in a detailed distribution of the crustal anisotropy in the European Alps and show first results for the isolated mantle anisotropy from the corrected XKS-waveforms and the crustal anisotropy from the receiver-function analysis.</p>


Author(s):  
Wolfgang Schöner

Glaciers are probably the most obvious features of Earth’s changing climate. They enable one to see the effects of a warming or a cooling of the atmosphere by landscape changes on time scales short enough to be perceived or recognized by humans. However, the relationship between a retreating and advancing glacier and the climate is not linear, as glacier flow can filter the direct signal of the climate. Thus, glaciers can advance during periods of warming or, vice versa, retreat during periods of cooling. In fact, it is the mass change of the glacier (i.e., the mass balance) that directly links the glacier reaction to an atmospheric signal. The mechanism-based understanding of the relationship between the changing climate and glacier reaction received important and significant momentum from the science of the Alpine region. This strong momentum from the Alps has to do with the well-established science tradition in Europe in the 19th and beginning of the 20th century, which resulted in a series of important inventions to measure climate and glacier properties. Even at that time, knowledge was gained that is still valid in the early 21st century (e.g., the climate is changing and fluctuating; glacier changes are caused by changing climate; and the ice age was the result of shifting climate). Above all others, Albrecht Penck and Eduard Brückner were the key scientists in this blossoming era of glacier climatology. Interest in a better understanding of the relationship of climate to glaciers was not only driven by curiosity, but also by several impacts of glaciers on human life in the Alps. Investigations of climate–glacier relationships in the Alps began with the expiration of the Little Ice Age (LIA) period when glaciers were particularly large but began to retreat significantly. Observations of post-LIA glacier front positions showed a sharp decline after their maximum extent in about 1850 until the turn of the 19th to 20th centuries, when they began to grow and advance again. They were also forming a prominent moraine around 1920, which was, however, far behind the 1850 extent. Interestingly, climate time series of the post LIA period show a general long-term cooling of summer temperatures and several decades of precipitation deficit in the second half of the 19th century. Thus, the retreat forced by climate changes cannot be simply explained by increasing air temperatures, though calibrated glacier mass balance models are able to simulate this period quite well. Additional effects related to the albedo could be a source for a better understanding. From 1920 onward, the climate moved into a period of warm and high-sunshine summers, which peaked in the 1940s until 1950. Glaciers started again to melt strongly and related discharges of pro-glacial rivers were exceptionally high during this period as glaciers were still quite large and the available energy for melt from radiation was enhanced. With the shift of the Atlantic meridional overturning (AMO), which is an important driver of European climate, into a negative mode in the 1960s, the mass balances of Alpine glaciers experienced more and more positive mass balance years. This finally resulted in a period of advancing glaciers and the development of frontal moraines around 1980 for a large number of glaciers. Thereafter, from 1980 onward, Alpine glaciers moved into an era of continuous negative mass balances and particularly strong retreat. The anthropogenic forcing from greenhouse gases together with global brightening and the increase of anticyclonic weather types in summer moved the climate and thus the mass balances of glaciers into a state far away from equilibrium. Given available scenarios of future climate, this retreat will continue and, even under the optimistic RCP2.6 scenario, glaciers (as derived from model simulations for the future) will not return to an equilibrium mass balance before the end of the 21st century. According to a glacier inventory for the European Alps from Landsat Thematic Mapper scenes of 2003, published by Paul and coworkers in 2011, the total surface of all glaciers and ice patches in the European Alps in 2003 was 2,056 km² (50% in Switzerland, 19% in Italy, 18% in Austria, 13% in France, and <1% in Germany). Generally, the reaction of Alpine glaciers to climate perturbations is rather well understood. For the glaciers of the Alps, important processes of glacier changes are related to the surface energy balance during the ablation season when radiation is the primary source of energy for snow and ice melt. Other ablation processes, such as sublimation and internal and basal ablation, are small compared to surface melt. This specificity enables the use of simple temperature-based models to simulate the mass balance of glaciers sufficiently well. Besides atmospheric forcing of glacier mass balance, glacier flow (which is related to englacial temperature distribution) plays a role, in particular, for observed front position changes of glaciers. Glaciers are continuously adapting their size to the climate, which could work much faster for smaller glaciers compared to large valley glaciers of the Alps having a response time of about 100 years.


Author(s):  
Christoph Kittl ◽  
James Robinson ◽  
Michael J. Raschke ◽  
Arne Olbrich ◽  
Andre Frank ◽  
...  

Abstract Purpose The purpose of this study was to examine the length change patterns of the native medial structures of the knee and determine the effect on graft length change patterns for different tibial and femoral attachment points for previously described medial reconstructions. Methods Eight cadaveric knee specimens were prepared by removing the skin and subcutaneous fat. The sartorius fascia was divided to allow clear identification of the medial ligamentous structures. Knees were then mounted in a custom-made rig and the quadriceps muscle and the iliotibial tract were loaded, using cables and hanging weights. Threads were mounted between tibial and femoral pins positioned in the anterior, middle, and posterior parts of the attachment sites of the native superficial medial collateral ligament (sMCL) and posterior oblique ligament (POL). Pins were also placed at the attachment sites relating to two commonly used medial reconstructions (Bosworth/Lind and LaPrade). Length changes between the tibiofemoral pin combinations were measured using a rotary encoder as the knee was flexed through an arc of 0–120°. Results With knee flexion, the anterior fibres of the sMCL tightened (increased in length 7.4% ± 2.9%) whilst the posterior fibres slackened (decreased in length 8.3% ± 3.1%). All fibre regions of the POL displayed a uniform lengthening of approximately 25% between 0 and 120° knee flexion. The most isometric tibiofemoral combination was between pins placed representing the middle fibres of the sMCL (Length change = 5.4% ± 2.1% with knee flexion). The simulated sMCL reconstruction that produced the least length change was the Lind/Bosworth reconstruction with the tibial attachment at the insertion of the semitendinosus and the femoral attachment in the posterior part of the native sMCL attachment side (5.4 ± 2.2%). This appeared more isometric than using the attachment positions described for the LaPrade reconstruction (10.0 ± 4.8%). Conclusion The complex behaviour of the native MCL could not be imitated by a single point-to-point combination and surgeons should be aware that small changes in the femoral MCL graft attachment position will significantly effect graft length change patterns. Reconstructing the sMCL with a semitendinosus autograft, left attached distally to its tibial insertion, would appear to have a minimal effect on length change compared to detaching it and using the native tibial attachment site. A POL graft must always be tensioned near extension to avoid capturing the knee or graft failure.


Water ◽  
2021 ◽  
Vol 13 (11) ◽  
pp. 1484
Author(s):  
Dagmar Dlouhá ◽  
Viktor Dubovský ◽  
Lukáš Pospíšil

We present an approach for the calibration of simplified evaporation model parameters based on the optimization of parameters against the most complex model for evaporation estimation, i.e., the Penman–Monteith equation. This model computes the evaporation from several input quantities, such as air temperature, wind speed, heat storage, net radiation etc. However, sometimes all these values are not available, therefore we must use simplified models. Our interest in free water surface evaporation is given by the need for ongoing hydric reclamation of the former Ležáky–Most quarry, i.e., the ongoing restoration of the land that has been mined to a natural and economically usable state. For emerging pit lakes, the prediction of evaporation and the level of water plays a crucial role. We examine the methodology on several popular models and standard statistical measures. The presented approach can be applied in a general model calibration process subject to any theoretical or measured evaporation.


1993 ◽  
Vol 70 (2) ◽  
pp. 549-558 ◽  
Author(s):  
R. Hallworth ◽  
B. N. Evans ◽  
P. Dallos

1. The microchamber method was used to examine the motile responses of isolated guinea pig outer hair cells to electrical stimulation. In the microchamber method, an isolated cell is drawn partway into a suction pipette and stimulated transcellularly. The relative position of the cell in the microchamber is referred to as the exclusion fraction. 2. The length changes of the included and excluded segments were compared for constant sinusoidal stimulus amplitude as functions of the exclusion fraction. Both included and excluded segments showed maximal responses when the cell was excluded approximately halfway. Both segments showed smaller or absent responses when the cell was almost fully excluded or almost fully included. 3. When the cell was near to, but not at, the maximum exclusion, the included segment response amplitude was zero, whereas the excluded segment response amplitude was nonzero. In contrast, when the cell was nearly fully included, the excluded segment response amplitude was zero, but the included segment response amplitude was still detectable. A simple model of outer hair cell motility based on these results suggests that the cell has finite-resistance terminations and that the motors are restricted to a region above the nucleus and below its ciliated apex (cuticular plate). 4. The function describing length change as a function of command voltage was measured for each segment as the exclusion fraction was varied. The functions were similar at midrange exclusions (i.e., when the segments were about equal length), showing nonlinearity and saturability. The functions were strikingly different when the segment lengths were different. The effects of exclusion on the voltage to length-change functions suggested that the nonlinearity and saturability are local properties of the motility mechanism. 5. The diameter changes of both segments were examined. The segment diameter changes were always antiphasic to the length changes. This finding implies that the motility mechanism has an active antiphasic diameter component. The diameter change amplitude was a monotonically increasing function of exclusion for the included segment, and a decreasing function for the excluded segment. 6. The voltage to length-change and voltage to diameter-change functions were measured for the same cell and exclusion fraction. The voltage to diameter-change function was smaller in amplitude than the voltage to length-change function. The functions were of opposite polarity to each other, but were otherwise similar in character. Thus it is likely that the same motor mechanism is responsible for both axial and diameter deformations.


Sign in / Sign up

Export Citation Format

Share Document