Self-consistency of the regolith hypothesis for the mid-Pleistocene Transition

Author(s):  
Matthew Drew ◽  
Lev Tarasov

<p>Is the regolith hypothesis consistent with the physics of glacial removal of mechanically weak surface material? </p><p> </p><p>The  mid-Pleistocene transition (MPT) from small 40 kyr glacial cycles to large, abruptly terminating 100 kyr ones represents a major climate system reorganization for which a clear understanding is lacking. A leading mechanism for this transition is a stabilization of ice sheets due to a shift to higher friction substrate. The Pleistocene saw the removal of deformable regolith -- laying bare hard higher-friction bedrock that would help preserve regional ice during warm interstadials. This is the regolith hypothesis. </p><p> </p><p>The removal of regolith by Pleistocene ice sheets remains poorly constrained. To date, only models with a forced change in area of regolith cover or 1D flow line models with simplistic sediment transport have been used to probe the role of regolith in the MPT. It is therefore unclear if the appropriate amount of regolith removal can occur within the time-frame of the MPT.</p><p> </p><p>To properly test the hypothesis, at least three components are required: capable model, observational constraint, and a probe of uncertainties. A capable model must explicitly represent relevant processes in a fully coupled self-consistent manner. We have therefore configured a state of the art 3D glacial systems model (GSM). The GSM incorporates a state-of-the-art fully coupled sediment production/transport model, subglacial hydrology, visco-elastic glacial isostatic adjustment, 3D thermomechanically coupled hybrid shallow ice/shallow shelf ice dynamics, and internal climate solution from an energy balance model. The model generates sediment by quarrying and abrasion, and both subglacial and englacial sediment transport. The subglacial hydrology model employs a linked-cavity system with a flux based switch to tunnel drainage, giving dynamic effective pressure needed for realistic sediment and sliding processes. The coupled model is driven only by prescribed atmospheric CO2 and orbitally derived insolation.</p><p> </p><p>The required observational constraints include present-day regolith distribution and inferred Pleistocene ice volume from proxy records.</p><p> </p><p>The final component is  a large ensemble of full Pleistocene simulations that probe both initial regolith distribution uncertainties and model parametric uncertainties. We present the results of such an ensemble, examining both rates of computed regolith removal and changes in ice volume cycling across the MPT interval.</p>

2014 ◽  
Vol 7 (5) ◽  
pp. 2141-2156 ◽  
Author(s):  
B. de Boer ◽  
P. Stocchi ◽  
R. S. W. van de Wal

Abstract. Relative sea-level variations during the late Pleistocene can only be reconstructed with the knowledge of ice-sheet history. On the other hand, the knowledge of regional and global relative sea-level variations is necessary to learn about the changes in ice volume. Overcoming this problem of circularity demands a fully coupled system where ice sheets and sea level vary consistently in space and time and dynamically affect each other. Here we present results for the past 410 000 years (410 kyr) from the coupling of a set of 3-D ice-sheet-shelf models to a global sea-level model, which is based on the solution of the gravitationally self-consistent sea-level equation. The sea-level model incorporates the glacial isostatic adjustment feedbacks for a Maxwell viscoelastic and rotating Earth model with coastal migration. Ice volume is computed with four 3-D ice-sheet-shelf models for North America, Eurasia, Greenland and Antarctica. Using an inverse approach, ice volume and temperature are derived from a benthic δ18O stacked record. The derived surface-air temperature anomaly is added to the present-day climatology to simulate glacial–interglacial changes in temperature and hence ice volume. The ice-sheet thickness variations are then forwarded to the sea-level model to compute the bedrock deformation, the change in sea-surface height and thus the relative sea-level change. The latter is then forwarded to the ice-sheet models. To quantify the impact of relative sea-level variations on ice-volume evolution, we have performed coupled and uncoupled simulations. The largest differences of ice-sheet thickness change occur at the edges of the ice sheets, where relative sea-level change significantly departs from the ocean-averaged sea-level variations.


2016 ◽  
Vol 10 (4) ◽  
pp. 1753-1769 ◽  
Author(s):  
Johannes Feldmann ◽  
Anders Levermann

Abstract. The concept of similitude is commonly employed in the fields of fluid dynamics and engineering but rarely used in cryospheric research. Here we apply this method to the problem of ice flow to examine the dynamic similitude of isothermal ice sheets in shallow-shelf approximation against the scaling of their geometry and physical parameters. Carrying out a dimensional analysis of the stress balance we obtain dimensionless numbers that characterize the flow. Requiring that these numbers remain the same under scaling we obtain conditions that relate the geometric scaling factors, the parameters for the ice softness, surface mass balance and basal friction as well as the ice-sheet intrinsic response time to each other. We demonstrate that these scaling laws are the same for both the (two-dimensional) flow-line case and the three-dimensional case. The theoretically predicted ice-sheet scaling behavior agrees with results from numerical simulations that we conduct in flow-line and three-dimensional conceptual setups. We further investigate analytically the implications of geometric scaling of ice sheets for their response time. With this study we provide a framework which, under several assumptions, allows for a fundamental comparison of the ice-dynamic behavior across different scales. It proves to be useful in the design of conceptual numerical model setups and could also be helpful for designing laboratory glacier experiments. The concept might also be applied to real-world systems, e.g., to examine the response times of glaciers, ice streams or ice sheets to climatic perturbations.


2014 ◽  
Vol 7 (3) ◽  
pp. 3505-3544 ◽  
Author(s):  
B. de Boer ◽  
P. Stocchi ◽  
R. S. W. van de Wal

Abstract. Relative sea-level variations during the late Pleistocene cannot be reconstructed regardless of the estimates of ice-volume fluctuations. For the latter, however, the knowledge of regional and global relative sea-level variations is necessary. Overcoming this problem of circularity demands a fully coupled system where ice sheets and sea level vary consistently in space and time and dynamically affect each other. Here we present results for the past 410 000 years (410 kyr) from the coupling of a set of 3-D ice-sheet-shelf models to a global sea-level model based on the solution of gravitationally self-consistent sea-level equation. The sea-level model incorporates all the Glacial Isostatic Adjustment feedbacks for a Maxwell viscoelastic and rotating Earth model with variable coastlines. Ice volume is computed with four 3-D ice-sheet-shelf models for North America, Eurasia, Greenland and Antarctica. With an inverse approach, ice volume and temperature are derived from a benthic δ18O stacked record. The ice-sheet thickness variations are then forwarded to the sea-level model to compute the bedrock deformation, the geoid and the relative sea-level change. The latter are used to generate the new topographies for the next time step, which are forwarded to the ice-sheet models. To quantify the impact of relative sea-level variations on ice-volume evolution, we have performed coupled and uncoupled simulations. The largest differences of ice-sheet thickness change show up in the proximity of the ice-sheets edges, where relative sea-level change significantly departs from the ocean-averaged sea level variation.


2021 ◽  
pp. 1-14
Author(s):  
Guillaume Jouvet ◽  
Guillaume Cordonnier ◽  
Byungsoo Kim ◽  
Martin Lüthi ◽  
Andreas Vieli ◽  
...  

Abstract This paper introduces the Instructed Glacier Model (IGM) – a model that simulates ice dynamics, mass balance and its coupling to predict the evolution of glaciers, icefields or ice sheets. The novelty of IGM is that it models the ice flow by a Convolutional Neural Network, which is trained from data generated with hybrid SIA + SSA or Stokes ice flow models. By doing so, the most computationally demanding model component is substituted by a cheap emulator. Once trained with representative data, we demonstrate that IGM permits to model mountain glaciers up to 1000 × faster than Stokes ones on Central Processing Units (CPU) with fidelity levels above 90% in terms of ice flow solutions leading to nearly identical transient thickness evolution. Switching to the GPU often permits additional significant speed-ups, especially when emulating Stokes dynamics or/and modelling at high spatial resolution. IGM is an open-source Python code which deals with two-dimensional (2-D) gridded input and output data. Together with a companion library of trained ice flow emulators, IGM permits user-friendly, highly efficient and mechanically state-of-the-art glacier and icefields simulations.


2020 ◽  
pp. 1-16
Author(s):  
Tim Hill ◽  
Christine F. Dow ◽  
Eleanor A. Bash ◽  
Luke Copland

Abstract Glacier surficial melt rates are commonly modelled using surface energy balance (SEB) models, with outputs applied to extend point-based mass-balance measurements to regional scales, assess water resource availability, examine supraglacial hydrology and to investigate the relationship between surface melt and ice dynamics. We present an improved SEB model that addresses the primary limitations of existing models by: (1) deriving high-resolution (30 m) surface albedo from Landsat 8 imagery, (2) calculating shadows cast onto the glacier surface by high-relief topography to model incident shortwave radiation, (3) developing an algorithm to map debris sufficiently thick to insulate the glacier surface and (4) presenting a formulation of the SEB model coupled to a subsurface heat conduction model. We drive the model with 6 years of in situ meteorological data from Kaskawulsh Glacier and Nàłùdäy (Lowell) Glacier in the St. Elias Mountains, Yukon, Canada, and validate outputs against in situ measurements. Modelled seasonal melt agrees with observations within 9% across a range of elevations on both glaciers in years with high-quality in situ observations. We recommend applying the model to investigate the impacts of surface melt for individual glaciers when sufficient input data are available.


2018 ◽  
Vol 11 (6) ◽  
pp. 2299-2314 ◽  
Author(s):  
Rubén Banderas ◽  
Jorge Alvarez-Solas ◽  
Alexander Robinson ◽  
Marisa Montoya

Abstract. Offline forcing methods for ice-sheet models often make use of an index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from the Greenland ice-core temperature reconstruction, with present-day climatologies. An important drawback of this approach is that it clearly misrepresents climate variability at millennial timescales. The reason for this is that the spatial glacial–interglacial anomaly field used is associated with orbital climatic variations, while it is scaled following the characteristic time evolution of the index, which includes orbital and millennial-scale climate variability. The spatial patterns of orbital and millennial variability are clearly not the same, as indicated by a wealth of models and data. As a result, this method can be expected to lead to a misrepresentation of climate variability and thus of the past evolution of Northern Hemisphere (NH) ice sheets. Here we illustrate the problems derived from this approach and propose a new offline climate forcing method that attempts to better represent the characteristic pattern of millennial-scale climate variability by including an additional spatial anomaly field associated with this timescale. To this end, three different synthetic transient forcing climatologies are developed for the past 120 kyr following a perturbative approach and are applied to an ice-sheet model. The impact of the climatologies on the paleo-evolution of the NH ice sheets is evaluated. The first method follows the usual index approach in which temperature anomalies relative to the present are calculated by combining a simulated glacial–interglacial climatic anomaly field, interpolated through an index derived from ice-core data, with present-day climatologies. In the second approach the representation of millennial-scale climate variability is improved by incorporating a simulated stadial–interstadial anomaly field. The third is a refinement of the second one in which the amplitudes of both orbital and millennial-scale variations are tuned to provide perfect agreement with a recently published absolute temperature reconstruction over Greenland. The comparison of the three climate forcing methods highlights the tendency of the usual index approach to overestimate the temperature variability over North America and Eurasia at millennial timescales. This leads to a relatively high NH ice-volume variability on these timescales. Through enhanced ablation, this results in too low an ice volume throughout the last glacial period (LGP), below or at the lower end of the uncertainty range of estimations. Improving the representation of millennial-scale variability alone yields an important increase in ice volume in all NH ice sheets but especially in the Fennoscandian Ice Sheet (FIS). Optimizing the amplitude of the temperature anomalies to match the Greenland reconstruction results in a further increase in the simulated ice-sheet volume throughout the LGP. Our new method provides a more realistic representation of orbital and millennial-scale climate variability and improves the transient forcing of ice sheets during the LGP. Interestingly, our new approach underestimates ice-volume variations on millennial timescales as indicated by sea-level records. This suggests that either the origin of the latter is not the NH or that processes not represented in our study, notably variations in oceanic conditions, need to be invoked to explain millennial-scale ice-volume fluctuations. We finally provide here both our derived climate evolution of the LGP using the three methods as well as the resulting ice-sheet configurations. These could be of interest for future studies dealing with the atmospheric or/and oceanic consequences of transient ice-sheet evolution throughout the LGP and as a source of climate input to other ice-sheet models.


1992 ◽  
Vol 338 (1285) ◽  
pp. 235-242 ◽  

The prediction of short-term (100 year) changes in the mass balance of ice sheets and longer-term (1000 years) variations in their ice volumes is important for a range of climatic and environmental models. The Antarctic ice sheet contains between 24 M km 3 and 29 M km 3 of ice, equivalent to a eustatic sea level change of between 60m and 72m. The annual surface accumulation is estimated to be of the order of 2200 Gtonnes, equivalent to a sea level change of 6 mm a -1 . Analysis of the present-day accumulation regime of Antarctica indicates that about 25% ( ca. 500 Gt a -1 ) of snowfall occurs in the Antarctic Peninsula region with an area of only 6.8% of the continent. To date most models have focused upon solving predictive algorithms for the climate-sensitivity of the ice sheet, and assume: (i) surface mass balance is equivalent to accumulation (i.e. no melting, evaporation or deflation); (ii) percentage change in accumulation is proportional to change in saturation mixing ratio above the surface inversion layer; and (iii) there is a linear relation between mean annual surface air tem perature and saturation mixing ratio. For the A ntarctic Peninsula with mountainous terrain containing ice caps, outlet glaciers, valley glaciers and ice shelves, where there can be significant ablation at low levels and distinct climatic regimes, models of the climate response must be more complex. In addition, owing to the high accumulation and flow rates, even short- to medium -term predictions must take account of ice dynamics. Relationships are derived for the mass balance sensitivity and, using a model developed by Hindmarsh, the transient effects of ice dynamics are estimated. It is suggested that for a 2°C rise in mean annual surface tem perature over 40 years, ablation in the A ntarctic Peninsula region would contribute at least 1.0 mm to sea level rise, offsetting the fall of 0.5 mm contributed by increased accum ulation.


2020 ◽  
Author(s):  
Matthias Scheiter ◽  
Marius Schaefer ◽  
Eduardo Flández ◽  
Deniz Bozkurt ◽  
Ralf Greve

Abstract. Glaciers and ice caps are thinning and retreating along the entire Andes ridge, and drivers of this mass loss vary between the different climate zones. The southern part of the Andes (Wet Andes) has the highest abundance of glaciers in number and size, and a proper understanding of ice dynamics is important to assess their evolution. In this contribution, we apply the ice sheet model SICOPOLIS to the Mocho-Choshuenco ice cap in the Chilean Lake District (40° S, 72° W, Wet Andes) to reproduce its current state and to project its evolution until the end of the 21st century under different global warming scenarios. First, we create a model spin-up using surface mass balance data observed on the south-eastern catchment, extrapolating them to the whole ice cap using an exposition-dependent parameterization. This spin-up is able to reproduce the most important present-day glacier features. Based on the spin-up, we then run the model 80 years into the future, forced by projected surface temperature anomalies from different global circulation models under different radiative pathway scenarios to obtain estimates of the ice cap's state by the end of the 21st century. The mean projected ice volume losses are 25 ± 19 % (RCP2.6), 64 ± 14 % (RCP4.5) and 94 ± 3 % (RCP8.5) with respect to the ice volume estimated by radio-echo sounding data from 2013. We estimate the uncertainty of our projections based on the spread of the results when forcing with different global climate models and on the uncertainty associated with the variation of the equilibrium line altitude with temperature change. Considering our results, we project an considerable deglaciation of the Chilean Lake District by the end of the 21st century.


Author(s):  
Antragama Ewa Abbas ◽  
◽  
Wirawan Agahari ◽  
Montijn van de Ven ◽  
Anneke Zuiderwijk ◽  
...  

Data marketplaces are expected to play a crucial role in tomorrow’s data economy but hardly achieve commercial exploitation. Currently, there is no clear understanding of the knowledge gaps in data marketplace research, especially neglected research topics that may contribute to advancing data marketplaces towards commercialization. This study provides an overview of the state of the art of data marketplace research. We employ a Systematic Literature Review (SLR) approach and structure our analysis using the Service-TechnologyOrganization-Finance (STOF) model. We find that the extant data marketplace literature is primarily dominated by technical research, such as discussions about computational pricing and architecture. To move past the first stage of the platform’s lifecycle (i.e., platform design) to the second stage (i.e., platform adoption), we call for empirical research in non-technological areas, such as customer expected value and market segmentation.


2018 ◽  
Author(s):  
Niall Gandy ◽  
Lauren J. Gregoire ◽  
Jeremy C. Ely ◽  
Christopher D. Clark ◽  
David M. Hodgson ◽  
...  

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. A valuable case to examine these processes is the last deglaciation of the British-Irish Ice Sheet. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British-Irish Ice Sheet during the last deglaciation, is well constrained, with abundant empirical data which could be used to inform, validate and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We simulate retreat from the shelf edge under constant "warm" surface mass balance and subshelf melt, to isolate the role of internal ice dynamics from external forcings. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at a "pinning-point" between mainland Scotland and the Isle of Lewis. At this stage, the presence of ice shelves became a major control on deglaciation, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch Strait produces an acceleration in retreat, leading to a "collapsed" state, even when the climate returns to the initial "cold" conditions. Our simulations demonstrate the importance of the Marine Ice Sheet Instability and ice shelf buttressing during the deglaciation of parts of the British-Irish Ice Sheet. Thus, geological data could be used to constrain these processes in ice sheet models used for projecting the future of our contemporary ice sheets.


Sign in / Sign up

Export Citation Format

Share Document