Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach
<p>To provide estimates of past glacier mass changes over the course of the 20th century, an adequate initial state is required. However, empirical evidence about past glacier states at regional or global scales is largely incomplete, both spatially and temporally, calling for the use of auto...
Main Authors: | , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Copernicus Publications
2019-12-01
|
Series: | The Cryosphere |
Online Access: | https://www.the-cryosphere.net/13/3317/2019/tc-13-3317-2019.pdf |
_version_ | 1819170694867451904 |
---|---|
author | J. Eis F. Maussion B. Marzeion B. Marzeion |
author_facet | J. Eis F. Maussion B. Marzeion B. Marzeion |
author_sort | J. Eis |
collection | DOAJ |
description | <p>To provide estimates of past glacier mass changes over the course of the 20th century, an adequate initial state is required. However, empirical evidence about past glacier states at regional or global scales is largely incomplete, both spatially and temporally, calling for the use of automated numerical methods.
This study presents a new way to initialize the Open Global Glacier Model from past climate information and present-day glacier states. We use synthetic experiments to show that even with these perfectly known but incomplete boundary conditions, the problem of model initialization is an ill-posed inverse problem leading to nonunique solutions, and we propose an ensemble approach as a way forward. The method works as follows: we generate a large set of physically plausible glacier candidates for a given year in the past (e.g., 1850 in the Alps), all of which are then modeled forward to the date of the observed glacier outline and evaluated by comparing the results of the forward runs to the present-day states. We test the approach on 2660 Alpine glaciers and determine error estimates of the method from the synthetic experiments. The results show that the solution is often nonunique, as many of the reconstructed initial states converge towards the observed state in the year of observation. We find that the median state of the best 5 % of all acceptable states is a reasonable best estimate. The accuracy of the method depends on the type of the considered observation for the evaluation (glacier length, area, or geometry). Trying to find past states from only present-day length instead of the full geometry leads to a sharp increase in uncertainty. Our study thus also provides quantitative information on how well the reconstructed initial glacier states are constrained through the limited information available to us. We analyze which glacier characteristics influence the reconstructability of a glacier, and we discuss ways to develop the method further for real-world applications.</p> |
first_indexed | 2024-12-22T19:39:28Z |
format | Article |
id | doaj.art-ee901aaf71b14d4d8736065ab901796b |
institution | Directory Open Access Journal |
issn | 1994-0416 1994-0424 |
language | English |
last_indexed | 2024-12-22T19:39:28Z |
publishDate | 2019-12-01 |
publisher | Copernicus Publications |
record_format | Article |
series | The Cryosphere |
spelling | doaj.art-ee901aaf71b14d4d8736065ab901796b2022-12-21T18:14:54ZengCopernicus PublicationsThe Cryosphere1994-04161994-04242019-12-01133317333710.5194/tc-13-3317-2019Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approachJ. Eis0F. Maussion1B. Marzeion2B. Marzeion3Institute of Geography, University of Bremen, Bremen, GermanyDepartment of Atmospheric and Cryospheric Sciences, University of Innsbruck, Innsbruck, AustriaInstitute of Geography, University of Bremen, Bremen, GermanyMARUM – Center for Marine Environmental Sciences, University of Bremen, Bremen, Germany<p>To provide estimates of past glacier mass changes over the course of the 20th century, an adequate initial state is required. However, empirical evidence about past glacier states at regional or global scales is largely incomplete, both spatially and temporally, calling for the use of automated numerical methods. This study presents a new way to initialize the Open Global Glacier Model from past climate information and present-day glacier states. We use synthetic experiments to show that even with these perfectly known but incomplete boundary conditions, the problem of model initialization is an ill-posed inverse problem leading to nonunique solutions, and we propose an ensemble approach as a way forward. The method works as follows: we generate a large set of physically plausible glacier candidates for a given year in the past (e.g., 1850 in the Alps), all of which are then modeled forward to the date of the observed glacier outline and evaluated by comparing the results of the forward runs to the present-day states. We test the approach on 2660 Alpine glaciers and determine error estimates of the method from the synthetic experiments. The results show that the solution is often nonunique, as many of the reconstructed initial states converge towards the observed state in the year of observation. We find that the median state of the best 5 % of all acceptable states is a reasonable best estimate. The accuracy of the method depends on the type of the considered observation for the evaluation (glacier length, area, or geometry). Trying to find past states from only present-day length instead of the full geometry leads to a sharp increase in uncertainty. Our study thus also provides quantitative information on how well the reconstructed initial glacier states are constrained through the limited information available to us. We analyze which glacier characteristics influence the reconstructability of a glacier, and we discuss ways to develop the method further for real-world applications.</p>https://www.the-cryosphere.net/13/3317/2019/tc-13-3317-2019.pdf |
spellingShingle | J. Eis F. Maussion B. Marzeion B. Marzeion Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach The Cryosphere |
title | Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach |
title_full | Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach |
title_fullStr | Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach |
title_full_unstemmed | Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach |
title_short | Initialization of a global glacier model based on present-day glacier geometry and past climate information: an ensemble approach |
title_sort | initialization of a global glacier model based on present day glacier geometry and past climate information an ensemble approach |
url | https://www.the-cryosphere.net/13/3317/2019/tc-13-3317-2019.pdf |
work_keys_str_mv | AT jeis initializationofaglobalglaciermodelbasedonpresentdayglaciergeometryandpastclimateinformationanensembleapproach AT fmaussion initializationofaglobalglaciermodelbasedonpresentdayglaciergeometryandpastclimateinformationanensembleapproach AT bmarzeion initializationofaglobalglaciermodelbasedonpresentdayglaciergeometryandpastclimateinformationanensembleapproach AT bmarzeion initializationofaglobalglaciermodelbasedonpresentdayglaciergeometryandpastclimateinformationanensembleapproach |