scholarly journals Simulating the Early Holocene demise of the Laurentide Ice Sheet with BISICLES (public trunk revision 3298)

2020 ◽  
Vol 13 (9) ◽  
pp. 4555-4577
Author(s):  
Ilkka S. O. Matero ◽  
Lauren J. Gregoire ◽  
Ruza F. Ivanovic

Abstract. Simulating the demise of the Laurentide Ice Sheet covering Hudson Bay in the Early Holocene (10–7 ka) is important for understanding the role of accelerated changes in ice sheet topography and melt in the 8.2 ka event, a century long cooling of the Northern Hemisphere by several degrees. Freshwater released from the ice sheet through a surface mass balance instability (known as the saddle collapse) has been suggested as a major forcing for the 8.2 ka event, but the temporal evolution of this pulse has not been constrained. Dynamical ice loss and marine interactions could have significantly accelerated the ice sheet demise, but simulating such processes requires computationally expensive models that are difficult to configure and are often impractical for simulating past ice sheets. Here, we developed an ice sheet model setup for studying the Laurentide Ice Sheet's Hudson Bay saddle collapse and the associated meltwater pulse in unprecedented detail using the BISICLES ice sheet model, an efficient marine ice sheet model of the latest generation which is capable of refinement to kilometre-scale resolutions and higher-order ice flow physics. The setup draws on previous efforts to model the deglaciation of the North American Ice Sheet for initialising the ice sheet temperature, recent ice sheet reconstructions for developing the topography of the region and ice sheet, and output from a general circulation model for a representation of the climatic forcing. The modelled deglaciation is in agreement with the reconstructed extent of the ice sheet, and the associated meltwater pulse has realistic timing. Furthermore, the peak magnitude of the modelled meltwater equivalent (0.07–0.13 Sv) is compatible with geological estimates of freshwater discharge through the Hudson Strait. The results demonstrate that while improved representations of the glacial dynamics and marine interactions are key for correctly simulating the pattern of Early Holocene ice sheet retreat, surface mass balance introduces by far the most uncertainty. The new model configuration presented here provides future opportunities to quantify the range of plausible amplitudes and durations of a Hudson Bay ice saddle collapse meltwater pulse and its role in forcing the 8.2 ka event.

2019 ◽  
Author(s):  
Ilkka S. O. Matero ◽  
Lauren J. Gregoire ◽  
Ruza F. Ivanovic

Abstract. Simulating the demise of the Laurentide Ice Sheet covering the Hudson Bay in the early Holocene (10-7 ka) is important for understanding the role of accelerated changes in ice sheet topography and melt in the 8.2 ka event, a century long cooling of the Northern Hemisphere by several degrees. Freshwater released from the ice sheet through a surface mass balance instability (known as the saddle collapse) has been suggested as a major forcing for the 8.2 ka event, but the temporal evolution of this pulse has not been constrained. Dynamical ice loss and marine interactions could have significantly accelerated the ice sheet demise, but simulating such processes requires computationally expensive models that are difficult to configure and are often impractical for simulating past ice sheets. Here, we developed an ice sheet model setup for studying the Laurentide Ice Sheet’s Hudson Bay saddle collapse and the associated meltwater pulse in unprecedented detail using the BISICLES ice sheet model, an efficient marine ice sheet model of the latest generation, capable of refinement to kilometre-scale resolution and higher-order ice flow physics. The setup draws on previous efforts to model the deglaciation of the North American Ice Sheet for initialising the ice sheet temperature, recent ice sheet reconstructions for developing the topography of the region and ice sheet, and output from a general circulation model for a representation of the climatic forcing. The modelled deglaciation is in agreement with the reconstructed extent of the ice sheet and the associated meltwater pulse has realistic timing. Furthermore,the peak magnitude of the modelled meltwater equivalent (0.07–0.13 Sv) is compatible with geological estimates of freshwater discharge through the Hudson Strait. The results demonstrate that while improved representation of the glacial dynamics and marine interactions are key for correctly simulating the pattern of early Holocene ice sheet retreat, surface mass balance introduces by far the most uncertainty. The new model configuration presented here provides future opportunities to quantify the range of plausible amplitudes and durations of a Hudson Bay ice saddle collapse meltwater pulse and its role in forcing the 8.2 ka event.


2010 ◽  
Vol 11 (1) ◽  
pp. 3-25 ◽  
Author(s):  
Sebastian H. Mernild ◽  
Glen E. Liston ◽  
Christopher A. Hiemstra ◽  
Jens H. Christensen

Abstract Fluctuations in the Greenland ice sheet (GrIS) surface mass balance (SMB) and freshwater influx to the surrounding oceans closely follow climate fluctuations and are of considerable importance to the global eustatic sea level rise. A state-of-the-art snow-evolution modeling system (SnowModel) was used to simulate variations in the GrIS melt extent, surface water balance components, changes in SMB, and freshwater influx to the ocean. The simulations are based on the Intergovernmental Panel on Climate Change scenario A1B modeled by the HIRHAM4 regional climate model (RCM) using boundary conditions from the ECHAM5 atmosphere–ocean general circulation model (AOGCM) from 1950 through 2080. In situ meteorological station [Greenland Climate Network (GC-Net) and World Meteorological Organization (WMO) Danish Meteorological Institute (DMI)] observations from inside and outside the GrIS were used to validate and correct RCM output data before they were used as input for SnowModel. Satellite observations and independent SMB studies were used to validate the SnowModel output and confirm the model’s robustness. The authors simulated an ∼90% increase in end-of-summer surface melt extent (0.483 × 106 km2) from 1950 to 2080 and a melt index (above 2000-m elevation) increase of 138% (1.96 × 106 km2 × days). The greatest difference in melt extent occurred in the southern part of the GrIS, and the greatest changes in the number of melt days were seen in the eastern part of the GrIS (∼50%–70%) and were lowest in the west (∼20%–30%). The rate of SMB loss, largely tied to changes in ablation processes, leads to an enhanced average loss of 331 km3 from 1950 to 2080 and an average SMB level of −99 km3 for the period 2070–80. GrIS surface freshwater runoff yielded a eustatic rise in sea level from 0.8 ± 0.1 (1950–59) to 1.9 ± 0.1 mm (2070–80) sea level equivalent (SLE) yr−1. The accumulated GrIS freshwater runoff contribution from surface melting equaled 160-mm SLE from 1950 through 2080.


2008 ◽  
Vol 2 (2) ◽  
pp. 117-129 ◽  
Author(s):  
X. Fettweis ◽  
E. Hanna ◽  
H. Gallée ◽  
P. Huybrechts ◽  
M. Erpicum

Abstract. Results from a regional climate simulation (1970–2006) over the Greenland ice sheet (GrIS) reveals that more than 97% of the interannual variability of the modelled Surface Mass Balance (SMB) can be explained by the GrIS summer temperature anomaly and the GrIS annual precipitation anomaly. This multiple regression is then used to empirically estimate the GrIS SMB since 1900 from climatological time series. The projected SMB changes in the 21st century are investigated with the set of simulations performed with atmosphere-ocean general circulation models (AOGCMs) of the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4). These estimates show that the high surface mass loss rates of recent years are not unprecedented in the GrIS history of the last hundred years. The minimum SMB rate seems to have occurred earlier in the 1930s and corresponds to a zero SMB rate. The AOGCMs project that the SMB rate of the 1930s would be common at the end of 2100. The temperature would be higher than in the 1930s but the increase of accumulation in the 21st century would partly offset the acceleration of surface melt due to the temperature increase. However, these assumptions are based on an empirical multiple regression only validated for recent/current climatic conditions, and the accuracy and time homogeneity of the data sets and AOGCM results used in these estimations constitute a large uncertainty.


2014 ◽  
Vol 8 (2) ◽  
pp. 1453-1477 ◽  
Author(s):  
B. Noël ◽  
X. Fettweis ◽  
W. J. van de Berg ◽  
M. R. van den Broeke ◽  
M. Erpicum

Abstract. During recent summers (2007–2012), several surface melt records were broken over the Greenland Ice Sheet (GrIS). The extreme summer melt resulted in part from a persistent negative phase of the North-Atlantic Oscillation (NAO), favouring warmer than normal conditions over the GrIS. In addition, it has been suggested that significant anomalies in sea ice cover (SIC) and sea surface temperature (SST) may partially explain recent anomalous GrIS surface melt. To assess the impact of 2007–2012 SIC and SST anomalies on GrIS surface mass balance (SMB), a set of sensitivity experiments was carried out with the regional climate model MAR. These simulations suggest that changes in SST and SIC in the seas surrounding Greenland do not significantly impact GrIS SMB, due to the katabatic winds blocking effect. These winds are strong enough to prevent oceanic near-surface air, influenced by SIC and SST variability, from penetrating far inland. Therefore, the ice sheet SMB response is restricted to coastal regions, where katabatic winds are weaker. However, anomalies in SIC and SST could have indirectly affected the surface melt by changing the general circulation in the North Atlantic region, favouring more frequent warm air advection to the GrIS.


1997 ◽  
Vol 25 ◽  
pp. 73-78 ◽  
Author(s):  
Gerhard Krinner ◽  
Christophe Genthon

The Laboratoire de Météorologie Dynamique (LMD) variable-grid atmospheric general circulation model (AGCM) was used in this study for a five-year high-resolution simulation of the Antarctic climate. The horizontal resolution is about 100 km over a large part of the ice sheet. This study focuses on the simulated surface mass balance (precipitation-evaporation sublimation-melt) and on the spatial and temporal variability of snowfall in Antarctica. The simulated annual mean surface mass balance for the whole continent is close to the observed value, and the model simulates well the spatial distribution of the surface mass balance. The annual cycle of snowfall exhibits a clear minimum in summer over the high interior plateau as well as for Antarctica as a whole, in agreement with the observations. In the interior of the continent, the model produces a permanent light background snowfall that accounts for about 5% of the total annual precipitation. The bulk of the snowfall is produced irregularly during periods that generally last only two or three days that are caused by cyclones off the coast.


2008 ◽  
Vol 2 (2) ◽  
pp. 225-254 ◽  
Author(s):  
X. Fettweis ◽  
E. Hanna ◽  
H. Gallée ◽  
P. Huybrechts ◽  
M. Erpicum

Abstract. Results from a regional climate simulation (1970–2006) over the Greenland ice sheet (GrIS) reveals that more than 97% of the interannual variability of the modelled Surface Mass Balance (SMB) can be explained by the GrIS summer temperature anomaly and the GrIS annual precipitation anomaly. This multiple regression is then used to empirically estimate the GrIS SMB since 1900 from climatological time series. The projected SMB changes in the 21st century are investigated with the set of simulations performed with atmosphere-ocean general circulation models (AOGCMs) of the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4). These estimates show that the high surface mass loss rates of recent years are not unprecedented in the GrIS history of the last hundred years. The minimum SMB rate seems to have occurred earlier in the 1930s. The AOGCMs project that the SMB rate of the 1930s would be common at the end of 2100. The temperature would be higher than in the 1930s but the increase of accumulation in the 21st century would partly offset the acceleration of surface melt due to the temperature increase. However, these assumptions are based on an empirical multiple regression only validated for recent/current climatic conditions, and the accuracy and time homogeneity of the data sets and AOGCM results used in these estimations constitute a large uncertainty.


2014 ◽  
Vol 8 (1) ◽  
pp. 181-194 ◽  
Author(s):  
T. L. Edwards ◽  
X. Fettweis ◽  
O. Gagliardini ◽  
F. Gillet-Chaulet ◽  
H. Goelzer ◽  
...  

Abstract. We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR (Modèle Atmosphérique Régional: Fettweis, 2007) regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling ice sheet and climate models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77° N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four "SMB lapse rates", gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CI) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.56 (95% CI: −0.22 to 1.33) kg m−3 a−1 for the north, and 1.91 (1.03 to 2.61) kg m−3 a−1 for the south. Above the ELA, the gradients are much smaller in magnitude: 0.09 (−0.03 to 0.23) kg m−3 a−1 in the north, and 0.07 (−0.07 to 0.59) kg m−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically founded approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections (Edwards et al., 2014).


2015 ◽  
Vol 8 (10) ◽  
pp. 3199-3213 ◽  
Author(s):  
K. Le Morzadec ◽  
L. Tarasov ◽  
M. Morlighem ◽  
H. Seroussi

Abstract. To investigate ice sheet evolution over the timescale of a glacial cycle, 3-D ice sheet models (ISMs) are typically run at "coarse" grid resolutions (10–50 km) that do not resolve individual mountains. This will introduce to-date unquantified errors in sub-grid (SG) transport, accumulation and ablation for regions of rough topography. In the past, synthetic hypsometric curves, a statistical summary of the topography, have been used in ISMs to describe the variability of these processes. However, there has yet to be detailed uncertainty analysis of this approach. We develop a new flow line SG model for embedding in coarse resolution models. A 1 km resolution digital elevation model was used to compute the local hypsometric curve for each coarse grid (CG) cell and to determine local parameters to represent the hypsometric bins' slopes and widths. The 1-D mass transport for the SG model is computed with the shallow ice approximation. We test this model against simulations from the 3-D Ice Sheet System Model (ISSM) run at 1 km grid resolution. Results show that none of the alternative parameterizations explored were able to adequately capture SG surface mass balance and flux processes. Via glacial cycle ensemble results for North America, we quantify the impact of SG model coupling in an ISM. We show that SG process representation and associated parametric uncertainties, related to the exchange of ice between the SG and CG cells, can have significant (up to 35 m eustatic sea level equivalent for the North American ice complex) impact on modelled ice sheet evolution.


Author(s):  
J.M Gregory ◽  
P Huybrechts

Accurate simulation of ice-sheet surface mass balance requires higher spatial resolution than is afforded by typical atmosphere–ocean general circulation models (AOGCMs), owing, in particular, to the need to resolve the narrow and steep margins where the majority of precipitation and ablation occurs. We have developed a method for calculating mass-balance changes by combining ice-sheet average time-series from AOGCM projections for future centuries, both with information from high-resolution climate models run for short periods and with a 20 km ice-sheet mass-balance model. Antarctica contributes negatively to sea level on account of increased accumulation, while Greenland contributes positively because ablation increases more rapidly. The uncertainty in the results is about 20% for Antarctica and 35% for Greenland. Changes in ice-sheet topography and dynamics are not included, but we discuss their possible effects. For an annual- and area-average warming exceeding in Greenland and in the global average, the net surface mass balance of the Greenland ice sheet becomes negative, in which case it is likely that the ice sheet would eventually be eliminated, raising global-average sea level by 7 m.


2013 ◽  
Vol 7 (1) ◽  
pp. 635-674 ◽  
Author(s):  
T. L. Edwards ◽  
X. Fettweis ◽  
O. Gagliardini ◽  
F. Gillet-Chaulet ◽  
H. Goelzer ◽  
...  

Abstract. We present a new parameterisation that relates surface mass balance (SMB: the sum of surface accumulation and surface ablation) to changes in surface elevation of the Greenland ice sheet (GrIS) for the MAR regional climate model. The motivation is to dynamically adjust SMB as the GrIS evolves, allowing us to force ice sheet models with SMB simulated by MAR while incorporating the SMB–elevation feedback, without the substantial technical challenges of coupling the two models. This also allows us to assess the effect of elevation feedback uncertainty on the GrIS contribution to sea level, using multiple global climate and ice sheet models, without the need for additional, expensive MAR simulations. We estimate this relationship separately below and above the equilibrium line altitude (ELA, separating negative and positive SMB) and for regions north and south of 77° N, from a set of MAR simulations in which we alter the ice sheet surface elevation. These give four "SMB lapse rates", gradients that relate SMB changes to elevation changes. We assess uncertainties within a Bayesian framework, estimating probability distributions for each gradient from which we present best estimates and credibility intervals (CIs) that bound 95% of the probability. Below the ELA our gradient estimates are mostly positive, because SMB usually increases with elevation: 0.54 (95% CI: −0.22 to 1.34) kg m−3 a−1 for the north, and 1.89 (1.03 to 2.61) kg m−3 a−1 for the south. Above the ELA the gradients are much smaller: 0.09 (−0.03 to 0.22) kg m−3 a−1 in the north, and 0.06 (−0.07 to 0.56) kg m−3 a−1 in the south, because SMB can either increase or decrease in response to increased elevation. Our statistically based approach allows us to make probabilistic assessments for the effect of elevation feedback uncertainty on sea level projections. In a companion paper we use the best estimates and upper and lower CI bounds in five ice sheet models, and the full probability distributions in another, to adjust simulated SMB from MAR forced by two global climate models for the SRES A1B scenario (Edwards et al., 2013).


Sign in / Sign up

Export Citation Format

Share Document