Search Results
Description and validation of the ice-sheet model Yelmo (version 1.0)
2020, Robinson, Alexander, Alvarez-Solas, Jorge, Montoya, Marisa, Goelzer, Heiko, Greve, Ralf, Ritz, Catherine
We describe the physics and features of the ice-sheet model Yelmo, an open-source project intended for collaborative development. Yelmo is a thermomechanical model, solving for the coupled velocity and temperature solutions of an ice sheet simultaneously. The ice dynamics are currently treated via a “hybrid” approach combining the shallow-ice and shallow-shelf/shelfy-stream approximations, which makes Yelmo an apt choice for studying a wide variety of problems. Yelmo's main innovations lie in its flexible and user-friendly infrastructure, which promotes portability and facilitates long-term development. In particular, all physics subroutines have been designed to be self-contained, so that they can be easily ported from Yelmo to other models, or easily replaced by improved or alternate methods in the future. Furthermore, hard-coded model choices are eschewed, replaced instead with convenient parameter options that allow the model to be adapted easily to different contexts. We show results for different ice-sheet benchmark tests, and we illustrate Yelmo's performance for the Antarctic ice sheet.
Time-scale synchronisation of oscillatory responses can lead to non-monotonous R-tipping
2023, Swierczek-Jereczek, Jan, Robinson, Alexander, Blasco, Javier, Alvarez-Solas, Jorge, Montoya, Marisa
Rate-induced tipping (R-tipping) describes the fact that, for multistable dynamic systems, an abrupt transition can take place not only because of the forcing magnitude, but also because of the forcing rate. In the present work, we demonstrate through the case study of a piecewise-linear oscillator (PLO), that increasing the rate of forcing can make the system tip in some cases but might also prevent it from tipping in others. This counterintuitive effect is further called non-monotonous R-tipping (NMRT) and has already been observed in recent studies. We show that, in the present case, the reason for NMRT is the peak synchronisation of oscillatory responses operating on different time scales. We further illustrate that NMRT can be observed even in the presence of additive white noise of intermediate amplitude. Finally, NMRT is also observed on a van-der-Pol oscillator with an unstable limit cycle, suggesting that this effect is not limited to systems with a discontinuous right-hand side such as the PLO. This insight might be highly valuable, as the current research on tipping elements is shifting from an equilibrium to a dynamic perspective while using models of increasing complexity, in which NMRT might be observed but hard to understand.
Exploring the impact of atmospheric forcing and basal drag on the Antarctic Ice Sheet under Last Glacial Maximum conditions
2021-1-18, Blasco, Javier, Alvarez-Solas, Jorge, Robinson, Alexander, Montoya, Marisa
Little is known about the distribution of ice in the Antarctic Ice Sheet (AIS) during the Last Glacial Maximum (LGM). Whereas marine and terrestrial geological data indicate that the grounded ice advanced to a position close to the continental-shelf break, the total ice volume is unclear. Glacial boundary conditions are potentially important sources of uncertainty, in particular basal friction and climatic boundary conditions. Basal friction exerts a strong control on the large-scale dynamics of the ice sheet and thus affects its size and is not well constrained. Glacial climatic boundary conditions determine the net accumulation and ice temperature and are also poorly known. Here we explore the effect of the uncertainty in both features on the total simulated ice storage of the AIS at the LGM. For this purpose we use a hybrid ice sheet shelf model that is forced with different basal drag choices and glacial background climatic conditions obtained from the LGM ensemble climate simulations of the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). Overall, we find that the spread in the simulated ice volume for the tested basal drag parameterizations is about the same range as for the different general circulation model (GCM) forcings (4 to 6 m sea level equivalent). For a wide range of plausible basal friction configurations, the simulated ice dynamics vary widely but all simulations produce fully extended ice sheets towards the continental-shelf break. More dynamically active ice sheets correspond to lower ice volumes, while they remain consistent with the available constraints on ice extent. Thus, this work points to the possibility of an AIS with very active ice streams during the LGM. In addition, we find that the surface boundary temperature field plays a crucial role in determining the ice extent through its effect on viscosity. For ice sheets of a similar extent and comparable dynamics, we find that the precipitation field determines the total AIS volume. However, precipitation is highly uncertain. Climatic fields simulated by climate models show more precipitation in coastal regions than a spatially uniform anomaly, which can lead to larger ice volumes. Our results strongly support using these paleoclimatic fields to simulate and study the LGM and potentially other time periods like the last interglacial. However, their accuracy must be assessed as well, as differences between climate model forcing lead to a large spread in the simulated ice volume and extension.