Submitted as: model description paper 09 Oct 2020

Submitted as: model description paper | 09 Oct 2020

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

Model of Early Diagenesis in the Upper Sediment with Adaptable complexity – MEDUSA (v. 2): a time-dependent biogeochemical sediment module for Earth System Models, process analysis and teaching

Guy Munhoven Guy Munhoven
  • Dépt. d’Astrophysique, Géophysique et Océanographie, Université de Liège, 4000 Liège, Belgium

Abstract. MEDUSA is a time-dependent one-dimensional numerical model of coupled early diagenetic processes in the surface seafloor sediment. In the vertical, the sediment is subdivided into two different zones. Solids (biogenic, mineral, etc.) raining down from the surface of the ocean are collected by the reactive mixed layer at the top. This is where chemical reactions take place. Solids are transported by bioturbation and advection, solutes by diffusion and bioirrigation. The classical coupled time-dependent early diagenesis equations (advection-diffusion reaction equations) are used to describe the evolutions of the solid and solute components here. Solids that get transported deeper than the bottom boundary of the reactive mixed layer enter the second zone underneath, where reactions and mixing are neglected. Gradually as solid material gets transferred here from the overlying reactive layer, it is buried and preserved in a stack of layers that make up a synthetic sediment core.

MEDUSA has been extensively modified since its first release from 2007. The composition of the two phases, the processes (chemical reactions) and chemical equilibria between solutes are not fixed any more, but get assembled from a set of XML based description files that are processed by a code generator to produce the required Fortran code. 1D, 2D and 2D × 2D interfaces have been introduced to facilitate the coupling to common grid configurations and material compositions used in biogeochemical models. MEDUSA can also be run in parallel computing environments using the Message Passing Interface (MPI).

Guy Munhoven

Status: final response (author comments only)
Status: final response (author comments only)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Login for authors/topical editors] [Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement

Guy Munhoven


Total article views: 353 (including HTML, PDF, and XML)
HTML PDF XML Total Supplement BibTeX EndNote
261 84 8 353 35 11 10
  • HTML: 261
  • PDF: 84
  • XML: 8
  • Total: 353
  • Supplement: 35
  • BibTeX: 11
  • EndNote: 10
Views and downloads (calculated since 09 Oct 2020)
Cumulative views and downloads (calculated since 09 Oct 2020)

Viewed (geographical distribution)

Total article views: 287 (including HTML, PDF, and XML) Thereof 284 with geography defined and 3 with unknown origin.
Country # Views %
  • 1
Latest update: 15 Apr 2021
Short summary
Sea-floor sediments play an important role for biogeochemical cycling of elements (e.g., carbon, silicon, nutrients,...) in the ocean. Realistic sediment modules are, however, not yet commonly used in global ocean biogeochemical models. Here we present MEDUSA, a model of the processes taking place in the surface sea-floor sediments and which control the interaction between the sediments and the ocean. MEDUSA can be configured to meet the exact needs of any given ocean biogeochemical model.