TY - JOUR
T1 - A model of Greenland ice sheet deglaciation constrained by observations of relative sea level and ice extent
AU - Lecavalier, Benoit S.
AU - Milne, Glenn A.
AU - Simpson, Matthew J.R.
AU - Wake, Leanne
AU - Huybrechts, Philippe
AU - Tarasov, Lev
AU - Kjeldsen, Kristian K.
AU - Funder, Svend
AU - Long, Antony J.
AU - Woodroffe, Sarah
AU - Dyke, Arthur S.
AU - Larsen, Nicolaj K.
N1 - Funding Information:
Support provided by the Natural Science and Engineering Research Council of Canada , the Canada Research Chairs program , the Canadian Foundation for Innovation and the University of Ottawa . We thank Shawn Marshall for discussions that contributed to the results presented in this paper.
PY - 2014/10/15
Y1 - 2014/10/15
N2 - An ice sheet model was constrained to reconstruct the evolution of the Greenland Ice Sheet (GrIS) from the Last Glacial Maximum (LGM) to present to improve our understanding of its response to climate change. The study involved applying a glaciological model in series with a glacial isostatic adjustment and relative sea-level (RSL) model. The model reconstruction builds upon the work of Simpson etal. (2009) through four main extensions: (1) a larger constraint database consisting of RSL and ice extent data; model improvements to the (2) climate and (3) sea-level forcing components; (4) accounting for uncertainties in non-Greenland ice. The research was conducted primarily to address data-model misfits and to quantify inherent model uncertainties with the Earth structure and non-Greenland ice. Our new model (termed Huy3) fits the majority of observations and is characterised by a number of defining features. During the LGM, the ice sheet had an excess of 4.7m ice-equivalent sea-level (IESL), which reached a maximum volume of 5.1m IESL at 16.5calka BP. Modelled retreat of ice from the continental shelf progressed at different rates and timings in different sectors. Southwest and Southeast Greenland began to retreat from the continental shelf by ~16 to 14calka BP, thus responding in part to the Bølling-Allerød warm event (c. 14.5calka BP); subsequently ice at the southern tip of Greenland readvanced during the Younger Dryas cold event. In northern Greenland the ice retreated rapidly from the continental shelf upon the climatic recovery out of the Younger Dryas to present-day conditions. Upon entering the Holocene (11.7calka BP), the ice sheet soon became land-based. During the Holocene Thermal Maximum (HTM; 9-5calka BP), air temperatures across Greenland were marginally higher than those at present and the GrIS margin retreated inland of its present-day southwest position by 40-60km at 4calka BP which produced a deficit volume of 0.16m IESL relative to present. In response to the HTM warmth, our optimal model reconstruction lost mass at a maximum centennial rate of c. 103.4 Gt/yr. Our results suggest that remaining data-model discrepancies are affiliated with missing physics and sub-grid processes of the glaciological model, uncertainties in the climate forcing, lateral Earth structure, and non-Greenland ice (particularly the North American component). Finally, applying the Huy3 Greenland reconstruction with our optimal Earth model we generate present-day uplift rates across Greenland due to past changes in the ocean and ice loads with explicit error bars due to uncertainties in the Earth structure. Present-day uplift rates due to past changes are spatially variable and range from 3.5 to -7 mm/a (including Earth model uncertainty).
AB - An ice sheet model was constrained to reconstruct the evolution of the Greenland Ice Sheet (GrIS) from the Last Glacial Maximum (LGM) to present to improve our understanding of its response to climate change. The study involved applying a glaciological model in series with a glacial isostatic adjustment and relative sea-level (RSL) model. The model reconstruction builds upon the work of Simpson etal. (2009) through four main extensions: (1) a larger constraint database consisting of RSL and ice extent data; model improvements to the (2) climate and (3) sea-level forcing components; (4) accounting for uncertainties in non-Greenland ice. The research was conducted primarily to address data-model misfits and to quantify inherent model uncertainties with the Earth structure and non-Greenland ice. Our new model (termed Huy3) fits the majority of observations and is characterised by a number of defining features. During the LGM, the ice sheet had an excess of 4.7m ice-equivalent sea-level (IESL), which reached a maximum volume of 5.1m IESL at 16.5calka BP. Modelled retreat of ice from the continental shelf progressed at different rates and timings in different sectors. Southwest and Southeast Greenland began to retreat from the continental shelf by ~16 to 14calka BP, thus responding in part to the Bølling-Allerød warm event (c. 14.5calka BP); subsequently ice at the southern tip of Greenland readvanced during the Younger Dryas cold event. In northern Greenland the ice retreated rapidly from the continental shelf upon the climatic recovery out of the Younger Dryas to present-day conditions. Upon entering the Holocene (11.7calka BP), the ice sheet soon became land-based. During the Holocene Thermal Maximum (HTM; 9-5calka BP), air temperatures across Greenland were marginally higher than those at present and the GrIS margin retreated inland of its present-day southwest position by 40-60km at 4calka BP which produced a deficit volume of 0.16m IESL relative to present. In response to the HTM warmth, our optimal model reconstruction lost mass at a maximum centennial rate of c. 103.4 Gt/yr. Our results suggest that remaining data-model discrepancies are affiliated with missing physics and sub-grid processes of the glaciological model, uncertainties in the climate forcing, lateral Earth structure, and non-Greenland ice (particularly the North American component). Finally, applying the Huy3 Greenland reconstruction with our optimal Earth model we generate present-day uplift rates across Greenland due to past changes in the ocean and ice loads with explicit error bars due to uncertainties in the Earth structure. Present-day uplift rates due to past changes are spatially variable and range from 3.5 to -7 mm/a (including Earth model uncertainty).
KW - Deglaciation
KW - Geological data
KW - Glacial isostatic adjustment
KW - Glacialogical model
KW - Greenland ice sheet
KW - Ice sheet reconstruction
KW - Relative sea level
UR - http://www.scopus.com/inward/record.url?scp=84906856541&partnerID=8YFLogxK
U2 - 10.1016/j.quascirev.2014.07.018
DO - 10.1016/j.quascirev.2014.07.018
M3 - Article
AN - SCOPUS:84906856541
SN - 0277-3791
VL - 102
SP - 54
EP - 84
JO - Quaternary Science Reviews
JF - Quaternary Science Reviews
ER -