Bayesian calibration of firn densification models

Firn densification modelling is key to understanding ice sheet mass balance, ice sheet surface elevation change, and the age difference between ice and the air in enclosed air bubbles. This has resulted in the development of many firn models, all relying to a certain degree on parameter calibration...

Full description

Bibliographic Details
Published in:The Cryosphere
Main Authors: Verjans, Vincent, Leeson, Amber A., Nemeth, Christopher, Stevens, C. Max, Kuipers Munneke, Peter, Noël, Brice, van Wessem, Jan Melchior
Format: Article in Journal/Newspaper
Language:English
Published: Copernicus Publications 2020
Subjects:
Online Access:https://doi.org/10.5194/tc-14-3017-2020
https://noa.gwlb.de/receive/cop_mods_00053082
https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00052735/tc-14-3017-2020.pdf
https://tc.copernicus.org/articles/14/3017/2020/tc-14-3017-2020.pdf
Description
Summary:Firn densification modelling is key to understanding ice sheet mass balance, ice sheet surface elevation change, and the age difference between ice and the air in enclosed air bubbles. This has resulted in the development of many firn models, all relying to a certain degree on parameter calibration against observed data. We present a novel Bayesian calibration method for these parameters and apply it to three existing firn models. Using an extensive dataset of firn cores from Greenland and Antarctica, we reach optimal parameter estimates applicable to both ice sheets. We then use these to simulate firn density and evaluate against independent observations. Our simulations show a significant decrease (24 % and 56 %) in observation–model discrepancy for two models and a smaller increase (15 %) for the third. As opposed to current methods, the Bayesian framework allows for robust uncertainty analysis related to parameter values. Based on our results, we review some inherent model assumptions and demonstrate how firn model choice and uncertainties in parameter values cause spread in key model outputs.