04 Dec 2023
 | 04 Dec 2023
Status: this preprint is open for discussion.

Description and validation of the ice sheet model Nix v1.0

Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas

Abstract. We present a physical description of the ice-sheet model Nix, an open-source project intended for collaborative development. Nix is a 2D thermomechanical model written in C/C++ that simultaneously solves for the momentum balance equations, mass conservation and temperature evolution. Nix's velocity solver includes a hierarchy of Stokes approximations: Blatter-Pattyn, depth-integrated higher order, shallow-shelf and shallow-ice. The grounding-line position is explicitly solved by a moving coordinate system that avoids further interpolations. The model can be easily forced with any external boundary conditions, including those of stochastic nature. Nix has been verified for standard test problems. Here we show results for a number of benchmark tests from standard intercomparison projects and assess grounding-line migration with an overdeepened bed geometry. Lastly, we further exploit the thermomechanical coupling by designing a suite of experiments where the forcing is a physical variable, unlike previously idealised forcing scenarios where ice temperatures are implicitly fixed via an ice rate factor. Namely, we use atmospheric temperatures and oceanic temperature anomalies to assess model hysteresis behaviour with active thermodynamics. Our results show that hysteresis in an overdeepened bed geometry is similar for atmospheric and oceanic forcings. We find that not only the particular sub-shelf melting parametrisation determines the temperature anomaly at which the ice sheet retreats, but also the particular value of calibrated heat exchange velocities. Notably, the classical hysteresis loop is narrowed for both forcing scenarios (i.e., atmospheric and oceanic) if the ice sheet is thermomechanically active as a results of the internal feedback among ice temperature, stress balance and viscosity. In summary, Nix combines rapid computational capabilities with a Blatter-Pattyn stress balance fully coupled to a thermomechanical solver, not only validating against established benchmarks but also offering a powerful tool for advancing our insight on ice dynamics and grounding line stability.

Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas

Status: open (until 21 Mar 2024)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on egusphere-2023-2690', Anonymous Referee #1, 07 Feb 2024 reply
  • RC2: 'Comment on egusphere-2023-2690', Tijn Berends, 14 Feb 2024 reply
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas


Total article views: 319 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
240 67 12 319 10 5
  • HTML: 240
  • PDF: 67
  • XML: 12
  • Total: 319
  • BibTeX: 10
  • EndNote: 5
Views and downloads (calculated since 04 Dec 2023)
Cumulative views and downloads (calculated since 04 Dec 2023)

Viewed (geographical distribution)

Total article views: 324 (including HTML, PDF, and XML) Thereof 324 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 26 Feb 2024
Short summary
We introduce Nix, an ice-sheet model designed for understanding how large masses of ice behave. Nix as a computer program that simulates the movement and temperature changes in ice sheets. Nix helps us study how ice sheets respond to changes in the atmosphere and ocean. We found that how fast ice melts under the shelves and how heat is exchanged, play a role in determining the future of ice sheets. Nix is a useful tool for learning more about how climate change affects polar ice sheets.