Skip to main content
Publications | Persons | Institutes | Projects
[ report an error in this record ]basket (0): add | show Print this page

Reactive transport in surface sediments: 2. Media: an object-oriented problem-solving environment for early diagenesis
Meysman, F.J.R.; Middelburg, J.J.; Herman, P.M.J.; Heip, C.H.R. (2003). Reactive transport in surface sediments: 2. Media: an object-oriented problem-solving environment for early diagenesis. Comput. Geosci. 29(3): 301-318
In: Computers and Geosciences. Elsevier Science: Oxford; New York. ISSN 0098-3004; e-ISSN 1873-7803, more
Peer reviewed article  

Available in  Authors 

Keywords
    Cycles > Chemical cycles > Geochemical cycle > Biogeochemical cycle
    Equilibrium
    Modelling
    Physics > Mechanics > Kinetics
    Sedimentation > Diagenesis

Authors  Top 
  • Meysman, F.J.R., correspondent, more
  • Middelburg, J.J., more
  • Herman, P.M.J., more
  • Heip, C.H.R., more

Abstract
    The MEDIA (Modelling Early DIAgenesis) software package comprises a flexible and extensible software system that provides problem-solving assistance for simulating 1D reactive transport in surface sediments. MEDIA allows multiple diagenetic models to be built by extending a model template with new model components from a toolbox of available objects (elements, species, parameters, reactions). A detailed review is given of the transport and reaction components available for model construction. Upon assemblage, the model is channelled to the numerical subunit of the MEDIA package. Via a canonical transformation, the user-defined mixed kinetic-equilibrium model is rearranged into a proper set of differential algebraic equations (DAE), for which both steady state and transient solutions can be calculated. Steady-state profiles are obtained either directly using a Newton-Raphson method, or alternatively, as the asymptotic result of a dynamic simulation. Dynamic simulations involve a global implicit procedure based on the stiff- ODE solver package VODE, employing a direct substitution approach to reduce the number of equations in the DAE system. Verification of the MEDIA code was accomplished (1) by comparison with analytic models and (2) by emulating the model formulation and output of the existing diagenetic model code STEADYSED. As an example application, a diagenetic model was constructed to analyse an extensive dataset collected from a marine sediment in the Santa Barbara Basin (California). The different pathways of organic matter mineralization were modelled, and the coupling among the biogeochemical cycles of C, O, N, S, Mn and Fe was investigated. Depth profiles of both porewater and solid-phase constituents could be reproduced with great accuracy.

All data in the Integrated Marine Information System (IMIS) is subject to the VLIZ privacy policy Top | Authors