Antarctic ice-sheet hysteresis in a three-dimensional hybrid ice-sheet model

Author(s):  
Marisa Montoya ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Javier Blasco ◽  
Ilaria Tabone ◽  
...  

<p>Ice sheets, in particular the Antarctic Ice Sheet (AIS), are considered as potential tipping elements (TEs) of the Earth system. The mechanism underlying tipping is the existence of positive feedbacks leading to self-amplification processes that, once triggered, dominate the dynamics of the system. Positive feedbacks can also lead to hysteresis, with implications for reversibility in the context of long-term future climate change. The main mechanism underlying ice-sheet hysteresis is the positive surface mass balance-elevation feedback.  Marine-based ice sheets, such as the western sector of the AIS, are furthermore subject to specific instability mechanisms that can potentially also lead to hysteresis. Simulations with ice-sheet models have robustly confirmed the presence of different degrees of hysteresis in the evolution of the AIS volume with respect to model parameters and/or climate forcing, suggesting that ice-sheet changes are potentially irreversible on long timescales. Nevertheless, AIS hysteresis is only now becoming a focus of more intensive modeling efforts, including active oceanic forcing in particular. Here, we investigate the hysteresis of the AIS in a three-dimensional hybrid ice-sheet--ice-shelf model with respect to individual atmospheric forcing, ocean forcing and both. The aim is to obtain a probabilistic assessment of the AIS hysteresis and of its critical temperature thresholds by investigating the effect of structural uncertainty, including the representation of ice-sheet dynamics, basal melting and internal feedbacks.</p>

1997 ◽  
Vol 25 ◽  
pp. 259-268 ◽  
Author(s):  
Mikhail Verbitsky ◽  
Barry Saltzman

A three-dimensional (3-D), high-resolution, non-linearly viscous, non-isothermal ice-sheet model is employed to calculate the “present-day” equilibrium regime of the Antarctic ice sheet and its evolution during the last glacial cycle. The model is augmented by an approximate formula for ice-sheet basal temperature, based on a scaling of the thermodynamic equation for the ice flow. Steady-state solutions for both the shape and extent of the areas of basal melting (or freezing) are shown to be in good agreement with those obtained from the solution of the full 3-D thermodynamic equation. The solution for the basal temperature field of the West Antaretie Siple Coast produces areas at the pressure-melting point separated by strips of frozen-to-bed ice, the structure of which is reminiscent of Ice Streams A–E. This configuration appears to be robust, preserving its features in spite of climatic changes during the last glacial cycle. Ice Stream C seems to be more vulnerable to stagnation, switching to a passive mode at least once during the penultimate interglacial. We conjecture that the peculiarities of local topography determine the unique behavior of Ice Stream C: reduced basal stress and, consequently, relatively weak warming due to internal friction and basal sliding is not able to counteract the advective cooling during the periods of increased snowfall rate.


2015 ◽  
Vol 56 (70) ◽  
pp. 63-69 ◽  
Author(s):  
Denis Callens ◽  
Nicolas Thonnard ◽  
Jan T.M. Lenaerts ◽  
Jan M. Van Wessem ◽  
Willem Jan Van de Berg ◽  
...  

AbstractMass changes of polar ice sheets have an important societal impact, because they affect global sea level. Estimating the current mass budget of ice sheets is equivalent to determining the balance between surface mass gain through precipitation and outflow across the grounding line. For the Antarctic ice sheet, grounding line outflow is governed by oceanic processes and outlet glacier dynamics. In this study, we compute the mass budget of major outlet glaciers in the eastern Dronning Maud Land sector of the Antarctic ice sheet using the input/output method. Input is given by recent surface accumulation estimates (SMB) of the whole drainage basin. The outflow at the grounding line is determined from the radar data of a recent airborne survey and satellite-based velocities using a flow model of combined plug flow and simple shear. This approach is an improvement on previous studies, as the ice thickness is measured, rather than being estimated from hydrostatic equilibrium. In line with the general thickening of the ice sheet over this sector, we estimate the regional mass balance in this area at 3.15 ± 8.23 Gt a−1 according to the most recent SMB model results.


2005 ◽  
Vol 18 (13) ◽  
pp. 2194-2198 ◽  
Author(s):  
Mikhail Verbitsky

Abstract In an earlier paper by Verbitsky and Saltzman, a vertically integrated, high-resolution, nonlinearly viscous, nonisothermal ice sheet model was presented to calculate the “present-day” equilibrium regime of the Antarctic ice sheet. Steady-state solutions for the ice topography and thermodynamics, represented by the extent of the areas of basal melting, were shown to be in good agreement with both observations and results obtained from other three-dimensional thermodynamical equations. The solution for the basal temperature field of the West Antarctic Siple Coast produced areas at the pressure melting point separated by strips of frozen-to-bed ice, the structure of which looks very similar to ice streams A–E. Since the possible response of the Siple Coast basal temperature pattern to global warming and to associated changes in the snowfall rate is not obvious, a special sensitivity study was conducted. Results of such a study suggest that increased precipitation rate and associated intensification of ice advection can effectively “shut down” West Antarctic ice streams.


2018 ◽  
Author(s):  
Javier Blasco ◽  
Ilaria Tabone ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

Abstract. The Antarctic Ice Sheet (AIS) is the largest ice sheet on Earth and hence a major potential contributor to future global sea-level rise. A wealth of studies suggest that increasing oceanic temperatures could cause a collapse of its marine-based western sector, the West Antarctic Ice Sheet, through the mechanism of marine ice-sheet instability, leading to a sea-level increase of 3–5 m. Thus, it is crucial to constrain the sensitivity of the AIS to rapid climate changes. The Last Glacial Period is an ideal benchmark period for this purpose as it was punctuated by abrupt Dansgaard-Oeschger events at millennial timescales. Because their centre of action was in the North Atlantic, where their climate impacts were largest, modelling studies have mainly focused on the millennial-scale evolution of Northern Hemisphere (NH) paleo ice sheets. Sea-level reconstructions attribute the origin of millennial-scale sea-level variations mainly to NH paleo ice sheets, with a minor but not negligible role to the AIS. Here we investigate the AIS response to millennial-scale climate variability for the first time. To this end we use a three-dimensional, thermomechanical hybrid, ice-sheet-shelf model. Different oceanic sensitivities are tested and the sea-level equivalent (SLE) contributions computed. We find that whereas atmospheric variability has no appreciable effect on the AIS, changes in submarine melting rates can have a strong impact on it. We show that in contrast to the widespread assumption that the AIS is a slow reactive and static ice sheet that responds at orbital timescales only, it can lead to ice discharges of almost 15 m of SLE involving substantial grounding line migrations at millennial timescales.


2021 ◽  
Author(s):  
Hannah Chorley

<p>The mid-Miocene provides an important example relevant to the response of the East Antarctic Ice Sheet (EAIS) to future anthropogenic climate change. Geological observations and earth system modelling show a broad link between declining carbon dioxide (CO2) concentrations and increasing size and sensitivity of ice sheets in the past. Future projections show CO2 concentrations could reach up to 1000 ppm before the end of the century, with global temperatures 4-5°C warmer - a climate not seen since the mid-Miocene. This time period is therefore becoming increasingly important to understanding future Antarctic Ice Sheet (AIS) response, as CO2 concentrations are already at Pliocene levels (∼400 ppm). An improved, more detailed understanding of the response of the AIS to past climatic variability provides important context for interpreting how the AIS will respond to future climate change under high CO2 scenarios. </p> <p>A dynamic EAIS characterised the mid-Miocene, with major variations in both volume and extent of terrestrial and marine ice sheets. While global climate remained warmer than present-day throughout, this interval was punctuated by an episode of unusual warmth within the Miocene Climatic Optimum (MCO, ∼17-15 Ma). The MCO is one of the warmest intervals since the onset of Antarctic glaciation, with CO2 concentrations of up to 840 ppm during peak warmth and coastal regions characterised by temperate vegetation and mean summer temperatures (MST) of up to ∼10°C. This warmth terminated with major cooling and ice expansion across the mid-Miocene Climate Transition (MMCT, ∼14.8-13.8 Ma). </p> <p>A ∼50 m thick ice-cemented terrestrial glacial sequence was recovered in drill cores from the Friis Hills, McMurdo Dry Valleys in 2016. A chronostratigraphic framework for the cores based on magnetostratigraphy, 40Ar/39Ar isotopic ages, and limited biostratigraphic constraints, revealed 15 sedimentary cycles of the advance and retreat of a temperate alpine glacier system between ∼15.1-13.8 Ma. Each cycle consists of traction tills and moraines deposited during ice advance and intervening glacio-fluvial to glacio-lacustrine lithofacies deposited during ice retreat. This record highlights the influence of increasing glacial-interglacial variability across the MMCT, with till facies becoming progressively thicker, drier and of wider provenance post 14.4 Ma, while interglacial sediments remained similar to those that characterised the late-MCO, sustaining tundra vegetation and MSTs of 6-7°C. </p> <p>An ensemble of model simulations were produced for a recently published mid-Miocene topography and a range of CO2 concentrations, Transantarctic Mountain (TAM) uplift scenarios, and glacial-interglacial orbits in order to better understand the mechanisms driving EAIS variability during the early to mid-Miocene. Sedimentological and palynological data for glacial-interglacial periods of the early to mid-Miocene provide the primary constraint on ice extent and temperature variability. Results of this model-data comparison were used to assess the likely boundary conditions for the MCO and MMCT, and inferred TAM elevations of 300-500 m lower than present-day, modelled CO2 concentrations up to 780 ppm during periods of peak warmth, and a transition to lower CO2 across the MMCT. The onset of marine-based ice advance across the continental shelf was inferred between 280-460 ppm modelled CO2, however, the persistence of a significantly retreated, thick EAIS under even the highest modelled CO2 concentrations is not consistent with proxy data constraints and implies a strong hysteresis effect in the model. The presence of localised tundra vegetation under low CO2 concentrations in the model supports the persistence of higher plants in coastal lowlands post-MMCT, following their extinction at higher elevations after ∼13.8 Ma. </p> <p>Terrestrial, marine, and far-field records were reconciled to assess glacial-interglacial variability and evolution of the EAIS across the mid-Miocene. While 15 cycles were identified within the Friis Hills record, only 7.5 of these are well enough constrained by the age model to be correlated to climate cycles in the δ 18O record, spanning ∼160 ka of the late-MCO and inferring a terrestrial-terminating AIS responding to local insolation controlled by precession. This is consistent with eccentricity modulated precession control implied in other coastal Antarctic and far-field records during the MCO, but results presented here also support a two stepped climatic shift at ∼14.6 and ∼13.8 Ma during the MMCT. This stepwise shift in climatic cooling is attributed to declining CO2, with two boundaries in long-term atmospheric CO2 concentrations crossed during this time: (1) A shift to CO2 concentrations below 460 ppm in the model supported the growth of annual sea-ice and advance of small-scale marine-based ice into the Ross Sea. (2) At 13.8 Ma, a further decline in CO2 concentrations to below 280 ppm supported perennial sea-ice development, limiting the influence of warm, deep-water upwelling, resulting in large-scale marine-based ice advance, ultimately stabilising the AIS. This stepwise mid-Miocene cooling implies threshold behaviour of the AIS during a long-term 200-300 ppm general decline in CO2 proxy records. </p>


2020 ◽  
Author(s):  
Andreas Born ◽  
Alexander Robinson

Abstract. Radar reflections from the interior of the Greenland ice sheet contain a comprehensive archive of past accumulation rates and ice dynamics. Combining this data with dynamic ice sheet models may greatly aid model calibration, improve past and future sea level estimates, and enable insights into past ice sheet dynamics that neither models nor data could achieve alone. Unfortunately, simulating the continental-scale ice sheet stratigraphy represents a major challenge for current ice sheet models. In this study, we present the first three-dimensional ice sheet model that explicitly simulates the Greenland englacial stratigraphy. Individual layers of accumulation are represented on a grid whose vertical axis is time so that they do not exchange mass with each other as the flow of ice deforms them. This isochronal advection scheme is independent from the ice dynamics and only requires modest input data from a host thermomechanical ice-sheet model, making it easy to transfer to a range of models. Using an ensemble of simulations, we show that direct comparison with the dated radiostratigraphy data yields notably more accurate results than selecting simulations based on total ice thickness. We show that the isochronal scheme produces a more reliable simulation of the englacial age profile than traditional age tracers. The interpretation of ice dynamics at different times is possible but limited by uncertainties in the upper and lower boundaries conditions, namely temporal variations in surface mass balance and basal friction.


2021 ◽  
Vol 15 (9) ◽  
pp. 4539-4556
Author(s):  
Andreas Born ◽  
Alexander Robinson

Abstract. Radar reflections from the interior of the Greenland ice sheet contain a comprehensive archive of past accumulation rates, ice dynamics, and basal melting. Combining these data with dynamic ice sheet models may greatly aid model calibration, improve past and future sea level estimates, and enable insights into past ice sheet dynamics that neither models nor data could achieve alone. Unfortunately, simulating the continental-scale ice sheet stratigraphy represents a major challenge for current ice sheet models. In this study, we present the first three-dimensional ice sheet model that explicitly simulates the Greenland englacial stratigraphy. Individual layers of accumulation are represented on a grid whose vertical axis is time so that they do not exchange mass with each other as the flow of ice deforms them. This isochronal advection scheme does not influence the ice dynamics and only requires modest input data from a host thermomechanical ice sheet model, making it easy to transfer to a range of models. Using an ensemble of simulations, we show that direct comparison with the dated radiostratigraphy data yields notably more accurate results than calibrating simulations based on total ice thickness. We show that the isochronal scheme produces a more reliable simulation of the englacial age profile than traditional age tracers. The interpretation of ice dynamics at different times is possible but limited by uncertainties in the upper and lower boundary conditions, namely temporal variations in surface mass balance and basal friction.


2015 ◽  
Vol 9 (3) ◽  
pp. 881-903 ◽  
Author(s):  
B. de Boer ◽  
A. M. Dolan ◽  
J. Bernales ◽  
E. Gasson ◽  
H. Goelzer ◽  
...  

Abstract. In the context of future climate change, understanding the nature and behaviour of ice sheets during warm intervals in Earth history is of fundamental importance. The late Pliocene warm period (also known as the PRISM interval: 3.264 to 3.025 million years before present) can serve as a potential analogue for projected future climates. Although Pliocene ice locations and extents are still poorly constrained, a significant contribution to sea-level rise should be expected from both the Greenland ice sheet and the West and East Antarctic ice sheets based on palaeo sea-level reconstructions. Here, we present results from simulations of the Antarctic ice sheet by means of an international Pliocene Ice Sheet Modeling Intercomparison Project (PLISMIP-ANT). For the experiments, ice-sheet models including the shallow ice and shelf approximations have been used to simulate the complete Antarctic domain (including grounded and floating ice). We compare the performance of six existing numerical ice-sheet models in simulating modern control and Pliocene ice sheets by a suite of five sensitivity experiments. We include an overview of the different ice-sheet models used and how specific model configurations influence the resulting Pliocene Antarctic ice sheet. The six ice-sheet models simulate a comparable present-day ice sheet, considering the models are set up with their own parameter settings. For the Pliocene, the results demonstrate the difficulty of all six models used here to simulate a significant retreat or re-advance of the East Antarctic ice grounding line, which is thought to have happened during the Pliocene for the Wilkes and Aurora basins. The specific sea-level contribution of the Antarctic ice sheet at this point cannot be conclusively determined, whereas improved grounding line physics could be essential for a correct representation of the migration of the grounding-line of the Antarctic ice sheet during the Pliocene.


2020 ◽  
Author(s):  
Dipayan Choudhury ◽  
Axel Timmermann ◽  
Fabian Schloesser ◽  
Malte Heinemann ◽  
David Pollard

Abstract. It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice-sheets have not been fully understood. Testing theories of ice-ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial build-up from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240–170 ka (thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation/deglaciation events which occurred during MIS 7e-7d-7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere-ocean-vegetation-ice-sheet model (LOVECLIP), we simulate the transient evolution of northern and southern hemisphere ice-sheets during the MIS 7-6 period in response to orbital and greenhouse-gas forcing. For a range of model parameters, the simulations capture the reconstructed evolution of global ice volume reasonably well. It is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary wave – ice-sheet topography feedback can trigger an unabated and unrealistic ice-sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate-ice sheet modeling.


2010 ◽  
Vol 4 (1) ◽  
pp. 1-12 ◽  
Author(s):  
M. Thoma ◽  
K. Grosfeld ◽  
C. Mayer ◽  
F. Pattyn

Abstract. Subglacial lakes in Antarctica influence to a large extent the flow of the ice sheet. In this study we use an idealised lake geometry to study this impact. We employ a) an improved three-dimensional full-Stokes ice flow model with a nonlinear rheology, b) a three-dimensional fluid dynamics model with eddy diffusion to simulate the basal mass balance at the lake-ice interface, and c) a newly developed coupler to exchange boundary conditions between the two individual models. Different boundary conditions are applied over grounded ice and floating ice. This results in significantly increased temperatures within the ice on top of the lake, compared to ice at the same depth outside the lake area. Basal melting of the ice sheet increases this lateral temperature gradient. Upstream the ice flow converges towards the lake and accelerates by about 10% whenever basal melting at the ice-lake boundary is present. Above and downstream of the lake, where the ice flow diverges, a velocity decrease of about 10% is simulated.


Sign in / Sign up

Export Citation Format

Share Document