Submitted as: model description paper 05 Aug 2021

Submitted as: model description paper | 05 Aug 2021

Review status: this preprint is currently under review for the journal GMD.

RADIv1: a non-steady-state early diagenetic model for ocean sediments in Julia and MATLAB/GNU Octave

Olivier Sulpis1,2, Matthew Humphreys3, Monica Wilhelmus4,5, Dustin Carroll5,6, William Berelson7, Dimitris Menemenlis5, Jack Middelburg1, and Jess Adkins8 Olivier Sulpis et al.
  • 1Department of Earth Sciences, Utrecht University, Utrecht, The Netherlands
  • 2Department of Earth and Planetary Sciences, McGill University, Montreal, Canada
  • 3Department of Ocean Systems (OCS), NIOZ Royal Netherlands Institute for Sea Research, Texel, The Netherlands
  • 4Center for Fluid Mechanics, Brown University, Providence, USA
  • 5Jet Propulsion Laboratory, California Institute of Technology, Pasadena, USA
  • 6Moss Landing Marine Laboratories, San José State University, Moss Landing, USA
  • 7Department of Earth Sciences, University of Southern California, Los Angeles, USA
  • 8Geological and Planetary Sciences, California Institute of Technology, Pasadena, USA

Abstract. We introduce a time-dependent, one-dimensional model of early diagenesis that we term RADI, an acronym accounting for the main processes included in the model: chemical Reactions, Advection, molecular and bio-Diffusion, and bio-Irrigation. RADI is targeted for study of deep-sea sediments, in particular those containing calcium carbonates (CaCO3). RADI combines CaCO3 dissolution driven by organic matter degradation with a diffusive boundary layer and integrates state-of-the-art parameterizations of CaCO3 dissolution kinetics in seawater, thus serving as a link between mechanistic surface-reaction modelling and global-scale biogeochemical models. RADI also includes CaCO3 precipitation, providing a continuum between CaCO3 dissolution and precipitation. RADI integrates components rather than individual chemical species for accessibility and is straightforward to compare against measurements. RADI is the first diagenetic model implemented in Julia, a high-performance programming language that is free and open source, and it is also available in MATLAB/GNU Octave. Here, we first describe the scientific background behind RADI and its implementations. Then, we evaluate its performance in three selected locations and explore other potential applications, such as the influence of tides and seasonality on early diagenesis in the deep ocean. RADI is a powerful tool to study the time-transient and steady-state response of the sedimentary system to environmental perturbation, such as deep-sea mining, deoxygenation or acidification events.

Olivier Sulpis et al.

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on gmd-2021-211', Anonymous Referee #1, 06 Sep 2021
  • RC2: 'Comment on gmd-2021-211', David Burdige, 29 Sep 2021

Olivier Sulpis et al.

Olivier Sulpis et al.


Total article views: 539 (including HTML, PDF, and XML)
HTML PDF XML Total Supplement BibTeX EndNote
446 84 9 539 33 0 2
  • HTML: 446
  • PDF: 84
  • XML: 9
  • Total: 539
  • Supplement: 33
  • BibTeX: 0
  • EndNote: 2
Views and downloads (calculated since 05 Aug 2021)
Cumulative views and downloads (calculated since 05 Aug 2021)

Viewed (geographical distribution)

Total article views: 473 (including HTML, PDF, and XML) Thereof 473 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 24 Oct 2021
Short summary
A quarter of the surface of the Earth is covered by marine sediments rich in calcium carbonates, and their dissolution acts as a giant antacid tablet protecting the ocean against human-made acidification caused by massive CO2 emissions. Here, we present a new model of sediment chemistry that incorporates the latest experimental findings on calcium carbonate dissolution kinetics. This model can be used to predict how marine sediments evolve through time in response to environmental perturbations.