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...

Full description

Bibliographic Details
Main Authors: J. Eis, F. Maussion, B. Marzeion
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&thinsp;% 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&thinsp;% 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