Final published version
Licence: CC BY: Creative Commons Attribution 4.0 International License
Research output: Contribution to Journal/Magazine › Journal article › peer-review
Research output: Contribution to Journal/Magazine › Journal article › peer-review
}
TY - JOUR
T1 - Development of physically based liquid water schemes for Greenland firn-densification models
AU - Verjans, Vincent
AU - Leeson, Amber
AU - Stevens, Max
AU - MacFerrin, Michael
AU - Noël, Brice
AU - Van Den Broeke, Michiel
PY - 2019/7/9
Y1 - 2019/7/9
N2 - As surface melt is increasing on the Greenland Ice Sheet (GrIS), quantifying the retention capacity of the firn layer is critical to linking meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation–refreezing process, and more physically based representations of liquid water flow might bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the bucket approach, the Richards equation in a single domain and the Richards equation in a dual domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large (porosity of the upper 15 m firn varies by up to 47 %). The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards equation, and the most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow simulates ice-layer formation more reliably and allows for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties about firn structure (e.g. grain size and shape, presence of ice layers) and about its hydraulic properties, as well as the one-dimensionality of firn models, render the implementation of physically based percolation schemes difficult. Additionally, the performance of firn models is still affected by the various effects affecting the densification process such as microstructural effects, wet snow metamorphism and temperature sensitivity when meltwater is present.
AB - As surface melt is increasing on the Greenland Ice Sheet (GrIS), quantifying the retention capacity of the firn layer is critical to linking meltwater production to meltwater runoff. Firn-densification models have so far relied on empirical approaches to account for the percolation–refreezing process, and more physically based representations of liquid water flow might bring improvements to model performance. Here we implement three types of water percolation schemes into the Community Firn Model: the bucket approach, the Richards equation in a single domain and the Richards equation in a dual domain, which accounts for partitioning between matrix and fast preferential flow. We investigate their impact on firn densification at four locations on the GrIS and compare model results with observations. We find that for all of the flow schemes, significant discrepancies remain with respect to observed firn density, particularly the density variability in depth, and that inter-model differences are large (porosity of the upper 15 m firn varies by up to 47 %). The simple bucket scheme is as efficient in replicating observed density profiles as the single-domain Richards equation, and the most physically detailed dual-domain scheme does not necessarily reach best agreement with observed data. However, we find that the implementation of preferential flow simulates ice-layer formation more reliably and allows for deeper percolation. We also find that the firn model is more sensitive to the choice of densification scheme than to the choice of water percolation scheme. The disagreements with observations and the spread in model results demonstrate that progress towards an accurate description of water flow in firn is necessary. The numerous uncertainties about firn structure (e.g. grain size and shape, presence of ice layers) and about its hydraulic properties, as well as the one-dimensionality of firn models, render the implementation of physically based percolation schemes difficult. Additionally, the performance of firn models is still affected by the various effects affecting the densification process such as microstructural effects, wet snow metamorphism and temperature sensitivity when meltwater is present.
U2 - 10.5194/tc-13-1819-2019
DO - 10.5194/tc-13-1819-2019
M3 - Journal article
VL - 13
SP - 1819
EP - 1842
JO - The Cryosphere
JF - The Cryosphere
IS - 7
ER -