Articles | Volume 18, issue 11
https://doi.org/10.5194/gmd-18-3241-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-18-3241-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators
Energy Geosciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA 94720, USA
Benjamin J. Andre
Energy Geosciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA 94720, USA
Jeffrey N. Johnson
Energy Geosciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA 94720, USA
Cohere Consulting, LLC, Seattle, WA 98105, USA
Glenn E. Hammond
Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland, WA 99354, USA
Benjamin N. Sulman
Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
Konstantin Lipnikov
Applied Mathematics and Plasma Physics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
Marcus S. Day
Computational Science Center, National Renewable Energy Laboratory, Golden, CO 80401, USA
James J. Beisman
Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
Daniil Svyatsky
Applied Mathematics and Plasma Physics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
Hang Deng
Department of Energy and Resources Engineering, Peking University, Peking, 100871, China
Peter C. Lichtner
Department of Civil Engineering, University of New Mexico, Albuquerque, NM 87131, USA
Carl I. Steefel
Energy Geosciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA 94720, USA
J. David Moulton
Applied Mathematics and Plasma Physics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
Related authors
Zexuan Xu, Rebecca Serata, Haruko Wainwright, Miles Denham, Sergi Molins, Hansell Gonzalez-Raymat, Konstantin Lipnikov, J. David Moulton, and Carol Eddy-Dilek
Hydrol. Earth Syst. Sci., 26, 755–773, https://doi.org/10.5194/hess-26-755-2022, https://doi.org/10.5194/hess-26-755-2022, 2022
Short summary
Short summary
Climate change could change the groundwater system and threaten water supply. To quantitatively evaluate its impact on water quality, numerical simulations with chemical and reaction processes are required. With the climate projection dataset, we used the newly developed hydrological and chemical model to investigate the movement of contaminants and assist the management of contamination sites.
Michael Nole, Katherine Muller, Glenn Hammond, Xiaoliang He, and Peter Lichtner
EGUsphere, https://doi.org/10.5194/egusphere-2025-1343, https://doi.org/10.5194/egusphere-2025-1343, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Subsurface injection of carbon dioxide (CO2) can be used for a variety of purposes including geologic carbon storage and enhanced oil recovery. Recently, CO2 injection into reactive host rocks has been explored as a way to transform CO2 into dense solid minerals. We present a simulation framework for modeling flow of CO2 due to injection and subsequent reactions that take place to mineralize CO2.
Michael Nole, Jonah Bartrand, Fawz Naim, and Glenn Hammond
Geosci. Model Dev., 18, 1413–1425, https://doi.org/10.5194/gmd-18-1413-2025, https://doi.org/10.5194/gmd-18-1413-2025, 2025
Short summary
Short summary
Safe carbon dioxide (CO2) storage is likely to be critical for mitigating some of the most severe effects of climate change. We present a simulation framework for modeling CO2 storage beneath the seafloor, where CO2 can form a solid. This can aid in permanent CO2 storage for long periods of time. Our models show what a commercial-scale CO2 injection would look like in a marine environment. We discuss what would need to be considered when designing a subsea CO2 injection.
Katherine A. Muller, Peishi Jiang, Glenn Hammond, Tasneem Ahmadullah, Hyun-Seob Song, Ravi Kukkadapu, Nicholas Ward, Madison Bowe, Rosalie K. Chu, Qian Zhao, Vanessa A. Garayburu-Caruso, Alan Roebuck, and Xingyuan Chen
Geosci. Model Dev., 17, 8955–8968, https://doi.org/10.5194/gmd-17-8955-2024, https://doi.org/10.5194/gmd-17-8955-2024, 2024
Short summary
Short summary
The new Lambda-PFLOTRAN workflow incorporates organic matter chemistry into reaction networks to simulate aerobic respiration and biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow in a Jupyter notebook interface that digests raw organic matter chemistry data via Fourier transform ion cyclotron resonance mass spectrometry, develops a representative reaction network, and completes a biogeochemical simulation with the open-source, parallel-reactive-flow, and transport code PFLOTRAN.
Piyoosh Jaysaval, Glenn E. Hammond, and Timothy C. Johnson
Geosci. Model Dev., 16, 961–976, https://doi.org/10.5194/gmd-16-961-2023, https://doi.org/10.5194/gmd-16-961-2023, 2023
Short summary
Short summary
We present a robust and highly scalable implementation of numerical forward modeling and inversion algorithms for geophysical electrical resistivity tomography data. The implementation is publicly available and developed within the framework of PFLOTRAN (http://www.pflotran.org), an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The paper details all the theoretical and implementation aspects of the new capabilities along with test examples.
Utkarsh Mital, Dipankar Dwivedi, James B. Brown, and Carl I. Steefel
Earth Syst. Sci. Data, 14, 4949–4966, https://doi.org/10.5194/essd-14-4949-2022, https://doi.org/10.5194/essd-14-4949-2022, 2022
Short summary
Short summary
We present a new dataset that estimates small-scale variations in precipitation and temperature in mountainous terrain. The dataset is generated using a new machine learning framework that extracts relationships between climate and topography from existing coarse-scale datasets. The generated dataset is shown to capture small-scale variations more reliably than existing datasets and constitutes a valuable resource to model the water cycle in the mountains of Colorado, western United States.
Elchin E. Jafarov, Daniil Svyatsky, Brent Newman, Dylan Harp, David Moulton, and Cathy Wilson
The Cryosphere, 16, 851–862, https://doi.org/10.5194/tc-16-851-2022, https://doi.org/10.5194/tc-16-851-2022, 2022
Short summary
Short summary
Recent research indicates the importance of lateral transport of dissolved carbon in the polygonal tundra, suggesting that the freeze-up period could further promote lateral carbon transport. We conducted subsurface tracer simulations on high-, flat-, and low-centered polygons to test the importance of the freeze–thaw cycle and freeze-up time for tracer mobility. Our findings illustrate the impact of hydraulic and thermal gradients on tracer mobility, as well as of the freeze-up time.
Glenn E. Hammond
Geosci. Model Dev., 15, 1659–1676, https://doi.org/10.5194/gmd-15-1659-2022, https://doi.org/10.5194/gmd-15-1659-2022, 2022
Short summary
Short summary
This paper describes a simplified interface for implementing and testing new chemical reactions within the reactive transport simulator PFLOTRAN. The paper describes the interface, providing example code for the interface. The paper includes several chemical reactions implemented through the interface.
Zexuan Xu, Rebecca Serata, Haruko Wainwright, Miles Denham, Sergi Molins, Hansell Gonzalez-Raymat, Konstantin Lipnikov, J. David Moulton, and Carol Eddy-Dilek
Hydrol. Earth Syst. Sci., 26, 755–773, https://doi.org/10.5194/hess-26-755-2022, https://doi.org/10.5194/hess-26-755-2022, 2022
Short summary
Short summary
Climate change could change the groundwater system and threaten water supply. To quantitatively evaluate its impact on water quality, numerical simulations with chemical and reaction processes are required. With the climate projection dataset, we used the newly developed hydrological and chemical model to investigate the movement of contaminants and assist the management of contamination sites.
Qina Yan, Haruko Wainwright, Baptiste Dafflon, Sebastian Uhlemann, Carl I. Steefel, Nicola Falco, Jeffrey Kwang, and Susan S. Hubbard
Earth Surf. Dynam., 9, 1347–1361, https://doi.org/10.5194/esurf-9-1347-2021, https://doi.org/10.5194/esurf-9-1347-2021, 2021
Short summary
Short summary
We develop a hybrid model to estimate the spatial distribution of the thickness of the soil layer, which also provides estimations of soil transport and soil production rates. We apply this model to two examples of hillslopes in the East River watershed in Colorado and validate the model. The results show that the north-facing (NF) hillslope has a deeper soil layer than the south-facing (SF) hillslope and that the hybrid model provides better accuracy than a machine-learning model.
William R. Wieder, Derek Pierson, Stevan Earl, Kate Lajtha, Sara G. Baer, Ford Ballantyne, Asmeret Asefaw Berhe, Sharon A. Billings, Laurel M. Brigham, Stephany S. Chacon, Jennifer Fraterrigo, Serita D. Frey, Katerina Georgiou, Marie-Anne de Graaff, A. Stuart Grandy, Melannie D. Hartman, Sarah E. Hobbie, Chris Johnson, Jason Kaye, Emily Kyker-Snowman, Marcy E. Litvak, Michelle C. Mack, Avni Malhotra, Jessica A. M. Moore, Knute Nadelhoffer, Craig Rasmussen, Whendee L. Silver, Benjamin N. Sulman, Xanthe Walker, and Samantha Weintraub
Earth Syst. Sci. Data, 13, 1843–1854, https://doi.org/10.5194/essd-13-1843-2021, https://doi.org/10.5194/essd-13-1843-2021, 2021
Short summary
Short summary
Data collected from research networks present opportunities to test theories and develop models about factors responsible for the long-term persistence and vulnerability of soil organic matter (SOM). Here we present the SOils DAta Harmonization database (SoDaH), a flexible database designed to harmonize diverse SOM datasets from multiple research networks.
Yilin Fang, Xingyuan Chen, Jesus Gomez Velez, Xuesong Zhang, Zhuoran Duan, Glenn E. Hammond, Amy E. Goldman, Vanessa A. Garayburu-Caruso, and Emily B. Graham
Geosci. Model Dev., 13, 3553–3569, https://doi.org/10.5194/gmd-13-3553-2020, https://doi.org/10.5194/gmd-13-3553-2020, 2020
Short summary
Short summary
Surface water quality along river corridors can be improved by the area of the stream bed and stream bank in which stream water mixes with shallow groundwater or hyporheic zones (HZs). These zones are ubiquitous and dominated by microorganisms that can process the dissolved nutrients exchanged at this interface of these zones. The modulation of surface water quality can be simulated by connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer representation.
Teamrat A. Ghezzehei, Benjamin Sulman, Chelsea L. Arnold, Nathaniel A. Bogie, and Asmeret Asefaw Berhe
Biogeosciences, 16, 1187–1209, https://doi.org/10.5194/bg-16-1187-2019, https://doi.org/10.5194/bg-16-1187-2019, 2019
Short summary
Short summary
Soil water is a medium from which microbes acquire resources and within which they are able to move. Occupancy and availability of water and oxygen gas in soils are mutually exclusive. In addition, as soil dries the remaining water is held with an increasing degree of adhesive energy, which restricts microbes' ability to extract resources from water. We introduce a mathematical model that describes these interacting effects and organic matter decomposition.
Jennifer W. Harden, Jonathan A. O'Donnell, Katherine A. Heckman, Benjamin N. Sulman, Charles D. Koven, Chien-Lu Ping, and Gary J. Michaelson
SOIL Discuss., https://doi.org/10.5194/soil-2018-41, https://doi.org/10.5194/soil-2018-41, 2019
Revised manuscript not accepted
Short summary
Short summary
We examined changes in soil carbon (C) associated with permafrost thaw, warming, and ecosystem shifts using a space-for-time study. Soil C turnover was estimated for soil C fractions using soil C and radiocarbon data. Observations informed a simple model to track soil C change over time. Both losses and gains of soil C occur in the profile due to shifts in C among density-separated fractions. Thawing initially resulted in C gains to mineral soil and eventually C losses as warming persists.
Gautam Bisht, William J. Riley, Glenn E. Hammond, and David M. Lorenzetti
Geosci. Model Dev., 11, 4085–4102, https://doi.org/10.5194/gmd-11-4085-2018, https://doi.org/10.5194/gmd-11-4085-2018, 2018
Short summary
Short summary
Most existing global land surface models used to study impacts of climate change on water resources routinely use different models for near-surface unsaturated soil and the deeper groundwater table. We developed a model that uses a unified treatment of soil hydrologic processes throughout the entire soil column. Using a calibrated drainage parameter, the new model is able to correctly predict deep water table depth as reported in an observationally constrained global dataset.
Nathaniel W. Chaney, Marjolein H. J. Van Huijgevoort, Elena Shevliakova, Sergey Malyshev, Paul C. D. Milly, Paul P. G. Gauthier, and Benjamin N. Sulman
Hydrol. Earth Syst. Sci., 22, 3311–3330, https://doi.org/10.5194/hess-22-3311-2018, https://doi.org/10.5194/hess-22-3311-2018, 2018
Short summary
Short summary
The petabytes of existing global environmental data provide an invaluable asset to improve the characterization of land heterogeneity in Earth system models. This study introduces a clustering algorithm that summarizes a domain's heterogeneity through spatially interconnected clusters. A series of land model simulations in central California using this approach illustrate the critical role that multi-scale heterogeneity can have on the macroscale water, energy, and carbon cycles.
Guoping Tang, Fengming Yuan, Gautam Bisht, Glenn E. Hammond, Peter C. Lichtner, Jitendra Kumar, Richard T. Mills, Xiaofeng Xu, Ben Andre, Forrest M. Hoffman, Scott L. Painter, and Peter E. Thornton
Geosci. Model Dev., 9, 927–946, https://doi.org/10.5194/gmd-9-927-2016, https://doi.org/10.5194/gmd-9-927-2016, 2016
Short summary
Short summary
We demonstrate that CLM-PFLOTRAN predictions are consistent with CLM4.5 for Arctic, temperate, and tropical sites. A tight relative tolerance may be needed to avoid false convergence when scaling, clipping, or log transformation is used to avoid negative concentration in implicit time stepping and Newton-Raphson methods. The log transformation method is accurate and robust while relaxing relative tolerance or using the clipping or scaling method can result in efficient solutions.
S. Karra, S. L. Painter, and P. C. Lichtner
The Cryosphere, 8, 1935–1950, https://doi.org/10.5194/tc-8-1935-2014, https://doi.org/10.5194/tc-8-1935-2014, 2014
Z. M. Subin, P. C. D. Milly, B. N. Sulman, S. Malyshev, and E. Shevliakova
Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hessd-11-8443-2014, https://doi.org/10.5194/hessd-11-8443-2014, 2014
Preprint withdrawn
Short summary
Short summary
We provide a framework for including the effects of fine-scale hydrology on biogeochemistry in Earth-system models (ESMs). We simulate a representative hillslope in each ESM grid cell. While including the hillslope does not change the average hydrology, it causes greater vegetation and soil carbon to accumulate in lowlands. This is important for understanding how soil carbon might be affected by climate change, particularly in wetlands.
Related subject area
Biogeosciences
Soil nitrous oxide emissions from global land ecosystems and their drivers within the LPJ-GUESS model (v4.1)
Parameterization toolbox for a physical–biogeochemical model compatible with FABM – a case study: the coupled 1D GOTM–ECOSMO E2E for the Sylt–Rømø Bight, North Sea
H2MV (v1.0): global physically constrained deep learning water cycle model with vegetation
NN-TOC v1: global prediction of total organic carbon in marine sediments using deep neural networks
China Wildfire Emission Dataset (ChinaWED v1) for the period 2012–2022
Process-based modeling of solar-induced chlorophyll fluorescence with VISIT-SIF version 1.0
Including the phosphorus cycle into the LPJ-GUESS dynamic global vegetation model (v4.1, r10994) – global patterns and temporal trends of N and P primary production limitation
A comprehensive land-surface vegetation model for multi-stream data assimilation, D&B v1.0
Sources of uncertainty in the SPITFIRE global fire model: development of LPJmL-SPITFIRE1.9 and directions for future improvements
The unicellular NUM v.0.91: a trait-based plankton model evaluated in two contrasting biogeographic provinces
FESOM2.1-REcoM3-MEDUSA2: an ocean–sea ice–biogeochemistry model coupled to a sediment model
Satellite-based modeling of wetland methane emissions on a global scale (SatWetCH4 1.0)
Emulating grid-based forest carbon dynamics using machine learning: an LPJ-GUESS v4.1.1 application
Systematic underestimation of type-specific ecosystem process variability in the Community Land Model v5 over Europe
pyVPRM: A next-generation Vegetation Photosynthesis and Respiration Model for the post-MODIS era
Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling
An improved model for air–sea exchange of elemental mercury in MITgcm-ECCOv4-Hg: the role of surfactants and waves
BOATSv2: new ecological and economic features improve simulations of high seas catch and effort
A dynamical process-based model for quantifying global agricultural ammonia emissions – AMmonia–CLIMate v1.0 (AMCLIM v1.0) – Part 1: Land module for simulating emissions from synthetic fertilizer use
Simulating the drought response of European tree species with the dynamic vegetation model LPJ-GUESS (v4.1, 97c552c5)
Simulating Ips typographus L. outbreak dynamics and their influence on carbon balance estimates with ORCHIDEE r8627
Biological nitrogen fixation of natural and agricultural vegetation simulated with LPJmL 5.7.9
Learning from conceptual models – a study of the emergence of cooperation towards resource protection in a social–ecological system
Development and assessment of the physical-biogeochemical ocean regional model in the Northwest Pacific: NPRT v1.0 (ROMS v3.9–TOPAZ v2.0)
TROLL 4.0: representing water and carbon fluxes, leaf phenology and intraspecific trait variation in a mixed-species individual-based forest dynamics model – Part 1: Model description
The biogeochemical model Biome-BGCMuSo v6.2 provides plausible and accurate simulations of the carbon cycle in central European beech forests
TROLL 4.0: representing water and carbon fluxes, leaf phenology, and intraspecific trait variation in a mixed-species individual-based forest dynamics model – Part 2: Model evaluation for two Amazonian sites
Estimation of above- and below-ground ecosystem parameters for the DVM-DOS-TEM v0.7.0 model using MADS v1.7.3: a synthetic case study
DeepPhenoMem V1.0: deep learning modelling of canopy greenness dynamics accounting for multi-variate meteorological memory effects on vegetation phenology
Impacts of land-use change on biospheric carbon: an oriented benchmark using the ORCHIDEE land surface model
Implementing the iCORAL (version 1.0) coral reef CaCO3 production module in the iLOVECLIM climate model
Assimilation of carbonyl sulfide (COS) fluxes within the adjoint-based data assimilation system – Nanjing University Carbon Assimilation System (NUCAS v1.0)
ML4Fire-XGBv1.0: Improving North American wildfire prediction by integrating a machine-learning fire model in a land surface model
Quantifying the role of ozone-caused damage to vegetation in the Earth system: a new parameterization scheme for photosynthetic and stomatal responses
Radiocarbon analysis reveals underestimation of soil organic carbon persistence in new-generation soil models
Exploring the potential of history matching for land surface model calibration
EAT v1.0.0: a 1D test bed for physical–biogeochemical data assimilation in natural waters
Using deep learning to integrate paleoclimate and global biogeochemistry over the Phanerozoic Eon
Modelling boreal forest's mineral soil and peat C dynamics with the Yasso07 model coupled with the Ricker moisture modifier
Dynamic ecosystem assembly and escaping the “fire trap” in the tropics: insights from FATES_15.0.0
In silico calculation of soil pH by SCEPTER v1.0
Simple process-led algorithms for simulating habitats (SPLASH v.2.0): robust calculations of water and energy fluxes
A global behavioural model of human fire use and management: WHAM! v1.0
Terrestrial Ecosystem Model in R (TEMIR) version 1.0: simulating ecophysiological responses of vegetation to atmospheric chemical and meteorological changes
biospheremetrics v1.0.2: an R package to calculate two complementary terrestrial biosphere integrity indicators – human colonization of the biosphere (BioCol) and risk of ecosystem destabilization (EcoRisk)
Modeling boreal forest soil dynamics with the microbially explicit soil model MIMICS+ (v1.0)
Optimal enzyme allocation leads to the constrained enzyme hypothesis: the Soil Enzyme Steady Allocation Model (SESAM; v3.1)
Implementing a dynamic representation of fire and harvest including subgrid-scale heterogeneity in the tile-based land surface model CLASSIC v1.45
Inferring the tree regeneration niche from inventory data using a dynamic forest model
Optimising CH4 simulations from the LPJ-GUESS model v4.1 using an adaptive Markov chain Monte Carlo algorithm
Jianyong Ma, Almut Arneth, Benjamin Smith, Peter Anthoni, Xu-Ri, Peter Eliasson, David Wårlind, Martin Wittenbrink, and Stefan Olin
Geosci. Model Dev., 18, 3131–3155, https://doi.org/10.5194/gmd-18-3131-2025, https://doi.org/10.5194/gmd-18-3131-2025, 2025
Short summary
Short summary
Nitrous oxide (N2O) is a powerful greenhouse gas mainly released from natural and agricultural soils. This study examines how global soil N2O emissions changed from 1961 to 2020 and identifies key factors driving these changes using an ecological model. The findings highlight croplands as the largest source, with factors like fertilizer use and climate change enhancing emissions. Rising CO2 levels, however, can partially mitigate N2O emissions through increased plant nitrogen uptake.
Hoa Nguyen, Ute Daewel, Neil Banas, and Corinna Schrum
Geosci. Model Dev., 18, 2961–2982, https://doi.org/10.5194/gmd-18-2961-2025, https://doi.org/10.5194/gmd-18-2961-2025, 2025
Short summary
Short summary
Parameterization is key in modeling to reproduce observations well but is often done manually. This study presents a particle-swarm-optimizer-based toolbox for marine ecosystem models, compatible with the Framework for Aquatic Biogeochemical Models, thus enhancing its reusability. Applied to the Sylt ecosystem, the toolbox effectively (1) identified multiple parameter sets that matched observations well, providing different insights into ecosystem dynamics, and (2) optimized model complexity.
Zavud Baghirov, Martin Jung, Markus Reichstein, Marco Körner, and Basil Kraft
Geosci. Model Dev., 18, 2921–2943, https://doi.org/10.5194/gmd-18-2921-2025, https://doi.org/10.5194/gmd-18-2921-2025, 2025
Short summary
Short summary
We use an innovative approach to studying the Earth's water cycle by integrating advanced machine learning techniques with a traditional water cycle model. Our model is designed to learn from observational data, with a particular emphasis on understanding the influence of vegetation on water movement. By closely aligning with real-world observations, our model offers new possibilities for enhancing our understanding of the water cycle and its interactions with vegetation.
Naveenkumar Parameswaran, Everardo González, Ewa Burwicz-Galerne, Malte Braack, and Klaus Wallmann
Geosci. Model Dev., 18, 2521–2544, https://doi.org/10.5194/gmd-18-2521-2025, https://doi.org/10.5194/gmd-18-2521-2025, 2025
Short summary
Short summary
Our research uses deep learning to predict organic carbon stocks in ocean sediments, which is crucial for understanding their role in the global carbon cycle. By analysing over 22 000 samples and various seafloor characteristics, our model gives more accurate results than traditional methods. We estimate that the top 10 cm of ocean sediments hold about 156 Pg of carbon. This work enhances carbon stock estimates and helps plan future sampling strategies to better understand oceanic carbon burial.
Zhengyang Lin, Ling Huang, Hanqin Tian, Anping Chen, and Xuhui Wang
Geosci. Model Dev., 18, 2509–2520, https://doi.org/10.5194/gmd-18-2509-2025, https://doi.org/10.5194/gmd-18-2509-2025, 2025
Short summary
Short summary
The China Wildfire Emission Dataset (ChinaWED v1) estimated wildfire emissions in China during 2012–2022 as 78.13 Tg CO2, 279.47 Gg CH4, and 6.26 Gg N2O annually. Agricultural fires dominated emissions, while forest and grassland emissions decreased. Seasonal peaks occurred in late spring, with hotspots in northeast, southwest, and east China. The model emphasizes the importance of using localized emission factors and high-resolution fire estimates for accurate assessments.
Tatsuya Miyauchi, Makoto Saito, Hibiki M. Noda, Akihiko Ito, Tomomichi Kato, and Tsuneo Matsunaga
Geosci. Model Dev., 18, 2329–2347, https://doi.org/10.5194/gmd-18-2329-2025, https://doi.org/10.5194/gmd-18-2329-2025, 2025
Short summary
Short summary
Solar-induced chlorophyll fluorescence (SIF) is an effective indicator for monitoring photosynthetic activity. This paper introduces VISIT-SIF, a biogeochemical model developed based on the Vegetation Integrative Simulator for Trace gases (VISIT) to represent satellite-observed SIF. Our simulations reproduced the global distribution and seasonal variations in observed SIF. VISIT-SIF helps to improve photosynthetic processes through a combination of biogeochemical modeling and observed SIF.
Mateus Dantas de Paula, Matthew Forrest, David Warlind, João Paulo Darela Filho, Katrin Fleischer, Anja Rammig, and Thomas Hickler
Geosci. Model Dev., 18, 2249–2274, https://doi.org/10.5194/gmd-18-2249-2025, https://doi.org/10.5194/gmd-18-2249-2025, 2025
Short summary
Short summary
Our study maps global nitrogen (N) and phosphorus (P) availability and how they changed from 1901 to 2018. We find that tropical regions are mostly P-limited, while temperate and boreal areas face N limitations. Over time, P limitation increased, especially in the tropics, while N limitation decreased. These shifts are key to understanding global plant growth and carbon storage, highlighting the importance of including P dynamics in ecosystem models.
Wolfgang Knorr, Matthew Williams, Tea Thum, Thomas Kaminski, Michael Voßbeck, Marko Scholze, Tristan Quaife, T. Luke Smallman, Susan C. Steele-Dunne, Mariette Vreugdenhil, Tim Green, Sönke Zaehle, Mika Aurela, Alexandre Bouvet, Emanuel Bueechi, Wouter Dorigo, Tarek S. El-Madany, Mirco Migliavacca, Marika Honkanen, Yann H. Kerr, Anna Kontu, Juha Lemmetyinen, Hannakaisa Lindqvist, Arnaud Mialon, Tuuli Miinalainen, Gaétan Pique, Amanda Ojasalo, Shaun Quegan, Peter J. Rayner, Pablo Reyes-Muñoz, Nemesio Rodríguez-Fernández, Mike Schwank, Jochem Verrelst, Songyan Zhu, Dirk Schüttemeyer, and Matthias Drusch
Geosci. Model Dev., 18, 2137–2159, https://doi.org/10.5194/gmd-18-2137-2025, https://doi.org/10.5194/gmd-18-2137-2025, 2025
Short summary
Short summary
When it comes to climate change, the land surface is where the vast majority of impacts happen. The task of monitoring those impacts across the globe is formidable and must necessarily rely on satellites – at a significant cost: the measurements are only indirect and require comprehensive physical understanding. We have created a comprehensive modelling system that we offer to the research community to explore how satellite data can be better exploited to help us capture the changes that happen on our lands.
Luke Oberhagemann, Maik Billing, Werner von Bloh, Markus Drüke, Matthew Forrest, Simon P. K. Bowring, Jessica Hetzer, Jaime Ribalaygua Batalla, and Kirsten Thonicke
Geosci. Model Dev., 18, 2021–2050, https://doi.org/10.5194/gmd-18-2021-2025, https://doi.org/10.5194/gmd-18-2021-2025, 2025
Short summary
Short summary
Under climate change, the conditions necessary for wildfires to form are occurring more frequently in many parts of the world. To help predict how wildfires will change in future, global fire models are being developed. We analyze and further develop one such model, SPITFIRE. Our work identifies and corrects sources of substantial bias in the model that are important to the global fire modelling field. With this analysis and these developments, we help to provide a basis for future improvements.
Trine Frisbæk Hansen, Donald Eugene Canfield, Ken Haste Andersen, and Christian Jannik Bjerrum
Geosci. Model Dev., 18, 1895–1916, https://doi.org/10.5194/gmd-18-1895-2025, https://doi.org/10.5194/gmd-18-1895-2025, 2025
Short summary
Short summary
We describe and test the size-based Nutrient-Unicellular-Multicellular model, which defines unicellular plankton using a single set of parameters, on a eutrophic and oligotrophic ecosystem. The results demonstrate that both sites can be modeled with similar parameters and robust performance over a wide range of parameters. The study shows that the model is useful for non-experts and applicable for modeling ecosystems with limited data. It holds promise for evolutionary and deep-time climate models.
Ying Ye, Guy Munhoven, Peter Köhler, Martin Butzin, Judith Hauck, Özgür Gürses, and Christoph Völker
Geosci. Model Dev., 18, 977–1000, https://doi.org/10.5194/gmd-18-977-2025, https://doi.org/10.5194/gmd-18-977-2025, 2025
Short summary
Short summary
Many biogeochemistry models assume all material reaching the seafloor is remineralized and returned to solution, which is sufficient for studies on short-term climate change. Under long-term climate change, the carbon storage in sediments slows down carbon cycling and influences feedbacks in the atmosphere–ocean–sediment system. This paper describes the coupling of a sediment model to an ocean biogeochemistry model and presents results under the pre-industrial climate and under CO2 perturbation.
Juliette Bernard, Elodie Salmon, Marielle Saunois, Shushi Peng, Penélope Serrano-Ortiz, Antoine Berchet, Palingamoorthy Gnanamoorthy, Joachim Jansen, and Philippe Ciais
Geosci. Model Dev., 18, 863–883, https://doi.org/10.5194/gmd-18-863-2025, https://doi.org/10.5194/gmd-18-863-2025, 2025
Short summary
Short summary
Despite their importance, uncertainties remain in the evaluation of the drivers of temporal variability of methane emissions from wetlands on a global scale. Here, a simplified global model is developed, taking advantage of advances in remote-sensing data and in situ observations. The model reproduces the large spatial and temporal patterns of emissions, albeit with limitations in the tropics due to data scarcity. This model, while simple, can provide valuable insights into sensitivity analyses.
Carolina Natel, David Martin Belda, Peter Anthoni, Neele Haß, Sam Rabin, and Almut Arneth
EGUsphere, https://doi.org/10.5194/egusphere-2024-4064, https://doi.org/10.5194/egusphere-2024-4064, 2025
Short summary
Short summary
Complex models predict forest carbon responses to future climate change but are slow and computationally intensive, limiting large-scale analyses. We used machine learning to accelerate predictions from the LPJ-GUESS vegetation model. Our emulators, based on random forests and neural networks, achieved 97 % faster simulations. This approach enables rapid exploration of climate mitigation strategies and supports informed policy decisions.
Christian Poppe Terán, Bibi S. Naz, Harry Vereecken, Roland Baatz, Rosie A. Fisher, and Harrie-Jan Hendricks Franssen
Geosci. Model Dev., 18, 287–317, https://doi.org/10.5194/gmd-18-287-2025, https://doi.org/10.5194/gmd-18-287-2025, 2025
Short summary
Short summary
Carbon and water exchanges between the atmosphere and the land surface contribute to water resource availability and climate change mitigation. Land surface models, like the Community Land Model version 5 (CLM5), simulate these. This study finds that CLM5 and other data sets underestimate the magnitudes of and variability in carbon and water exchanges for the most abundant plant functional types compared to observations. It provides essential insights for further research into these processes.
Theo Glauch, Julia Marshall, Christoph Gerbig, Santiago Botía, Michał Gałkowski, Sanam N. Vardag, and André Butz
EGUsphere, https://doi.org/10.5194/egusphere-2024-3692, https://doi.org/10.5194/egusphere-2024-3692, 2025
Short summary
Short summary
The Vegetation Photosynthesis and Respiration Model (VPRM) estimates carbon exchange between the atmosphere and biosphere by modeling gross primary production and respiration using satellite data and weather variables. Our new version, pyVPRM, supports diverse satellite products like Sentinel-2, MODIS, VIIRS and new land cover maps, enabling high spatial and temporal resolution. This improves flux estimates, especially in complex landscapes, and ensures continuity as MODIS nears decommissioning.
Katherine A. Muller, Peishi Jiang, Glenn Hammond, Tasneem Ahmadullah, Hyun-Seob Song, Ravi Kukkadapu, Nicholas Ward, Madison Bowe, Rosalie K. Chu, Qian Zhao, Vanessa A. Garayburu-Caruso, Alan Roebuck, and Xingyuan Chen
Geosci. Model Dev., 17, 8955–8968, https://doi.org/10.5194/gmd-17-8955-2024, https://doi.org/10.5194/gmd-17-8955-2024, 2024
Short summary
Short summary
The new Lambda-PFLOTRAN workflow incorporates organic matter chemistry into reaction networks to simulate aerobic respiration and biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow in a Jupyter notebook interface that digests raw organic matter chemistry data via Fourier transform ion cyclotron resonance mass spectrometry, develops a representative reaction network, and completes a biogeochemical simulation with the open-source, parallel-reactive-flow, and transport code PFLOTRAN.
Ling Li, Peipei Wu, Peng Zhang, Shaojian Huang, and Yanxu Zhang
Geosci. Model Dev., 17, 8683–8695, https://doi.org/10.5194/gmd-17-8683-2024, https://doi.org/10.5194/gmd-17-8683-2024, 2024
Short summary
Short summary
In this study, we incorporate sea surfactants and wave-breaking processes into MITgcm-ECCOv4-Hg. The updated model shows increased fluxes in high-wind-speed and high-wave regions and vice versa, enhancing spatial heterogeneity. It shows that elemental mercury (Hg0) transfer velocity is more sensitive to wind speed. These findings may elucidate the discrepancies in previous estimations and offer insights into global Hg cycling.
Jerome Guiet, Daniele Bianchi, Kim J. N. Scherrer, Ryan F. Heneghan, and Eric D. Galbraith
Geosci. Model Dev., 17, 8421–8454, https://doi.org/10.5194/gmd-17-8421-2024, https://doi.org/10.5194/gmd-17-8421-2024, 2024
Short summary
Short summary
The BiOeconomic mArine Trophic Size-spectrum (BOATSv2) model dynamically simulates global commercial fish populations and their coupling with fishing activity, as emerging from environmental and economic drivers. New features, including separate pelagic and demersal populations, iron limitation, and spatial variation of fishing costs and management, improve the accuracy of high seas fisheries. The updated model code is available to simulate both historical and future scenarios.
Jize Jiang, David S. Stevenson, and Mark A. Sutton
Geosci. Model Dev., 17, 8181–8222, https://doi.org/10.5194/gmd-17-8181-2024, https://doi.org/10.5194/gmd-17-8181-2024, 2024
Short summary
Short summary
A special model called AMmonia–CLIMate (AMCLIM) has been developed to understand and calculate NH3 emissions from fertilizer use and also taking into account how the environment influences these NH3 emissions. It is estimated that about 17 % of applied N in fertilizers was lost due to NH3 emissions. Hot and dry conditions and regions with high-pH soils can expect higher NH3 emissions.
Benjamin Franklin Meyer, João Paulo Darela-Filho, Konstantin Gregor, Allan Buras, Qiao-Lin Gu, Andreas Krause, Daijun Liu, Phillip Papastefanou, Sijeh Asuk, Thorsten E. E. Grams, Christian S. Zang, and Anja Rammig
EGUsphere, https://doi.org/10.5194/egusphere-2024-3352, https://doi.org/10.5194/egusphere-2024-3352, 2024
Short summary
Short summary
Climate change has increased the likelihood of drought events across Europe, potentially threatening European forest carbon sink. Dynamic vegetation models with mechanistic plant hydraulic architecture are needed to model these developments. We evaluate the plant hydraulic architecture version of LPJ-GUESS and show it's capability at capturing species-specific evapotranspiration responses to drought and reproducing flux observations of both gross primary production and evapotranspiration.
Guillaume Marie, Jina Jeong, Hervé Jactel, Gunnar Petter, Maxime Cailleret, Matthew J. McGrath, Vladislav Bastrikov, Josefine Ghattas, Bertrand Guenet, Anne Sofie Lansø, Kim Naudts, Aude Valade, Chao Yue, and Sebastiaan Luyssaert
Geosci. Model Dev., 17, 8023–8047, https://doi.org/10.5194/gmd-17-8023-2024, https://doi.org/10.5194/gmd-17-8023-2024, 2024
Short summary
Short summary
This research looks at how climate change influences forests, and particularly how altered wind and insect activities could make forests emit instead of absorb carbon. We have updated a land surface model called ORCHIDEE to better examine the effect of bark beetles on forest health. Our findings suggest that sudden events, such as insect outbreaks, can dramatically affect carbon storage, offering crucial insights into tackling climate change.
Stephen Björn Wirth, Johanna Braun, Jens Heinke, Sebastian Ostberg, Susanne Rolinski, Sibyll Schaphoff, Fabian Stenzel, Werner von Bloh, Friedhelm Taube, and Christoph Müller
Geosci. Model Dev., 17, 7889–7914, https://doi.org/10.5194/gmd-17-7889-2024, https://doi.org/10.5194/gmd-17-7889-2024, 2024
Short summary
Short summary
We present a new approach to modelling biological nitrogen fixation (BNF) in the Lund–Potsdam–Jena managed Land dynamic global vegetation model. While in the original approach BNF depended on actual evapotranspiration, the new approach considers soil water content and temperature, vertical root distribution, the nitrogen (N) deficit and carbon (C) costs. The new approach improved simulated BNF compared to the scientific literature and the model ability to project future C and N cycle dynamics.
Saeed Harati-Asl, Liliana Perez, and Roberto Molowny-Horas
Geosci. Model Dev., 17, 7423–7443, https://doi.org/10.5194/gmd-17-7423-2024, https://doi.org/10.5194/gmd-17-7423-2024, 2024
Short summary
Short summary
Social–ecological systems are the subject of many sustainability problems. Because of the complexity of these systems, we must be careful when intervening in them; otherwise we may cause irreversible damage. Using computer models, we can gain insight about these complex systems without harming them. In this paper we describe how we connected an ecological model of forest insect infestation with a social model of cooperation and simulated an intervention measure to save a forest from infestation.
Daehyuk Kim, Hyun-Chae Jung, Jae-Hong Moon, and Na-Hyeon Lee
EGUsphere, https://doi.org/10.5194/egusphere-2024-1509, https://doi.org/10.5194/egusphere-2024-1509, 2024
Short summary
Short summary
Physical–biogeochemical ocean global models is difficult to analyze oceanic environmental systems. To accurately understand the physical–biogeochemical processes at the regional scale, physical and biogeochemical models were coupled at a high resolution. The results successfully simulated the seasonal variations of chlorophyll and nutrients, particularly in the marginal seas, which were not captured by global models. The model is an important tool for studying physical–biogeochemical processes.
Isabelle Maréchaux, Fabian Jörg Fischer, Sylvain Schmitt, and Jérôme Chave
EGUsphere, https://doi.org/10.5194/egusphere-2024-3104, https://doi.org/10.5194/egusphere-2024-3104, 2024
Short summary
Short summary
We describe TROLL 4.0, a simulator of forest dynamics that represents trees in a virtual space at one-meter resolution. Tree birth, growth, death and the underlying physiological processes such as carbon assimilation, water transpiration and leaf phenology depend on plant traits that are measured in the field for many individuals and species. The model is thus capable of jointly simulating forest structure, diversity and ecosystem functioning, a major challenge in modelling vegetation dynamics.
Katarína Merganičová, Ján Merganič, Laura Dobor, Roland Hollós, Zoltán Barcza, Dóra Hidy, Zuzana Sitková, Pavel Pavlenda, Hrvoje Marjanovic, Daniel Kurjak, Michal Bošel'a, Doroteja Bitunjac, Maša Zorana Ostrogović Sever, Jiří Novák, Peter Fleischer, and Tomáš Hlásny
Geosci. Model Dev., 17, 7317–7346, https://doi.org/10.5194/gmd-17-7317-2024, https://doi.org/10.5194/gmd-17-7317-2024, 2024
Short summary
Short summary
We developed a multi-objective calibration approach leading to robust parameter values aiming to strike a balance between their local precision and broad applicability. Using the Biome-BGCMuSo model, we tested the calibrated parameter sets for simulating European beech forest dynamics across large environmental gradients. Leveraging data from 87 plots and five European countries, the results demonstrated reasonable local accuracy and plausible large-scale productivity responses.
Sylvain Schmitt, Fabian Fischer, James Ball, Nicolas Barbier, Marion Boisseaux, Damien Bonal, Benoit Burban, Xiuzhi Chen, Géraldine Derroire, Jeremy Lichstein, Daniela Nemetschek, Natalia Restrepo-Coupe, Scott Saleska, Giacomo Sellan, Philippe Verley, Grégoire Vincent, Camille Ziegler, Jérôme Chave, and Isabelle Maréchaux
EGUsphere, https://doi.org/10.5194/egusphere-2024-3106, https://doi.org/10.5194/egusphere-2024-3106, 2024
Short summary
Short summary
We evaluate the capability of TROLL 4.0, a simulator of forest dynamics, to represent tropical forest structure, diversity and functioning in two Amazonian forests. Evaluation data include forest inventories, carbon and water fluxes between the forest and the atmosphere, and leaf area and canopy height from remote-sensing products. The model realistically predicts the structure and composition, and the seasonality of carbon and water fluxes at both sites.
Elchin E. Jafarov, Helene Genet, Velimir V. Vesselinov, Valeria Briones, Aiza Kabeer, Andrew L. Mullen, Benjamin Maglio, Tobey Carman, Ruth Rutter, Joy Clein, Chu-Chun Chang, Dogukan Teber, Trevor Smith, Joshua M. Rady, Christina Schädel, Jennifer D. Watts, Brendan M. Rogers, and Susan M. Natali
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-158, https://doi.org/10.5194/gmd-2024-158, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Thawing permafrost could greatly impact global climate. Our study improves modeling of carbon cycling in Arctic ecosystems. We developed an automated method to fine-tune a model that simulates carbon and nitrogen flows, using computer-generated data. Using computer-generated data, we tested our method and found it enhances accuracy and reduces the time needed for calibration. This work helps make climate predictions more reliable in sensitive permafrost regions.
Guohua Liu, Mirco Migliavacca, Christian Reimers, Basil Kraft, Markus Reichstein, Andrew D. Richardson, Lisa Wingate, Nicolas Delpierre, Hui Yang, and Alexander J. Winkler
Geosci. Model Dev., 17, 6683–6701, https://doi.org/10.5194/gmd-17-6683-2024, https://doi.org/10.5194/gmd-17-6683-2024, 2024
Short summary
Short summary
Our study employs long short-term memory (LSTM) networks to model canopy greenness and phenology, integrating meteorological memory effects. The LSTM model outperforms traditional methods, enhancing accuracy in predicting greenness dynamics and phenological transitions across plant functional types. Highlighting the importance of multi-variate meteorological memory effects, our research pioneers unlock the secrets of vegetation phenology responses to climate change with deep learning techniques.
Thi Lan Anh Dinh, Daniel Goll, Philippe Ciais, and Ronny Lauerwald
Geosci. Model Dev., 17, 6725–6744, https://doi.org/10.5194/gmd-17-6725-2024, https://doi.org/10.5194/gmd-17-6725-2024, 2024
Short summary
Short summary
The study assesses the performance of the dynamic global vegetation model (DGVM) ORCHIDEE in capturing the impact of land-use change on carbon stocks across Europe. Comparisons with observations reveal that the model accurately represents carbon fluxes and stocks. Despite the underestimations in certain land-use conversions, the model describes general trends in soil carbon response to land-use change, aligning with the site observations.
Nathaelle Bouttes, Lester Kwiatkowski, Manon Berger, Victor Brovkin, and Guy Munhoven
Geosci. Model Dev., 17, 6513–6528, https://doi.org/10.5194/gmd-17-6513-2024, https://doi.org/10.5194/gmd-17-6513-2024, 2024
Short summary
Short summary
Coral reefs are crucial for biodiversity, but they also play a role in the carbon cycle on long time scales of a few thousand years. To better simulate the future and past evolution of coral reefs and their effect on the global carbon cycle, hence on atmospheric CO2 concentration, it is necessary to include coral reefs within a climate model. Here we describe the inclusion of coral reef carbonate production in a carbon–climate model and its validation in comparison to existing modern data.
Huajie Zhu, Mousong Wu, Fei Jiang, Michael Vossbeck, Thomas Kaminski, Xiuli Xing, Jun Wang, Weimin Ju, and Jing M. Chen
Geosci. Model Dev., 17, 6337–6363, https://doi.org/10.5194/gmd-17-6337-2024, https://doi.org/10.5194/gmd-17-6337-2024, 2024
Short summary
Short summary
In this work, we developed the Nanjing University Carbon Assimilation System (NUCAS v1.0). Data assimilation experiments were conducted to demonstrate the robustness and investigate the feasibility and applicability of NUCAS. The assimilation of ecosystem carbonyl sulfide (COS) fluxes improved the model performance in gross primary productivity, evapotranspiration, and sensible heat, showing that COS provides constraints on parameters relevant to carbon-, water-, and energy-related processes.
Ye Liu, Huilin Huang, Sing-Chun Wang, Tao Zhang, Donghui Xu, and Yang Chen
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-151, https://doi.org/10.5194/gmd-2024-151, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study integrates machine learning with a land surface model to improve wildfire predictions in North America. Traditional models struggle with accurately simulating burned areas due to simplified processes. By combining the predictive power of machine learning with a land model, our hybrid framework better captures fire dynamics. This approach enhances our understanding of wildfire behavior and aids in developing more effective climate and fire management strategies.
Fang Li, Zhimin Zhou, Samuel Levis, Stephen Sitch, Felicity Hayes, Zhaozhong Feng, Peter B. Reich, Zhiyi Zhao, and Yanqing Zhou
Geosci. Model Dev., 17, 6173–6193, https://doi.org/10.5194/gmd-17-6173-2024, https://doi.org/10.5194/gmd-17-6173-2024, 2024
Short summary
Short summary
A new scheme is developed to model the surface ozone damage to vegetation in regional and global process-based models. Based on 4210 data points from ozone experiments, it accurately reproduces statistically significant linear or nonlinear photosynthetic and stomatal responses to ozone in observations for all vegetation types. It also enables models to implicitly capture the variability in plant ozone tolerance and the shift among species within a vegetation type.
Alexander S. Brunmayr, Frank Hagedorn, Margaux Moreno Duborgel, Luisa I. Minich, and Heather D. Graven
Geosci. Model Dev., 17, 5961–5985, https://doi.org/10.5194/gmd-17-5961-2024, https://doi.org/10.5194/gmd-17-5961-2024, 2024
Short summary
Short summary
A new generation of soil models promises to more accurately predict the carbon cycle in soils under climate change. However, measurements of 14C (the radioactive carbon isotope) in soils reveal that the new soil models face similar problems to the traditional models: they underestimate the residence time of carbon in soils and may therefore overestimate the net uptake of CO2 by the land ecosystem. Proposed solutions include restructuring the models and calibrating model parameters with 14C data.
Nina Raoult, Simon Beylat, James M. Salter, Frédéric Hourdin, Vladislav Bastrikov, Catherine Ottlé, and Philippe Peylin
Geosci. Model Dev., 17, 5779–5801, https://doi.org/10.5194/gmd-17-5779-2024, https://doi.org/10.5194/gmd-17-5779-2024, 2024
Short summary
Short summary
We use computer models to predict how the land surface will respond to climate change. However, these complex models do not always simulate what we observe in real life, limiting their effectiveness. To improve their accuracy, we use sophisticated statistical and computational techniques. We test a technique called history matching against more common approaches. This method adapts well to these models, helping us better understand how they work and therefore how to make them more realistic.
Jorn Bruggeman, Karsten Bolding, Lars Nerger, Anna Teruzzi, Simone Spada, Jozef Skákala, and Stefano Ciavatta
Geosci. Model Dev., 17, 5619–5639, https://doi.org/10.5194/gmd-17-5619-2024, https://doi.org/10.5194/gmd-17-5619-2024, 2024
Short summary
Short summary
To understand and predict the ocean’s capacity for carbon sequestration, its ability to supply food, and its response to climate change, we need the best possible estimate of its physical and biogeochemical properties. This is obtained through data assimilation which blends numerical models and observations. We present the Ensemble and Assimilation Tool (EAT), a flexible and efficient test bed that allows any scientist to explore and further develop the state of the art in data assimilation.
Dongyu Zheng, Andrew S. Merdith, Yves Goddéris, Yannick Donnadieu, Khushboo Gurung, and Benjamin J. W. Mills
Geosci. Model Dev., 17, 5413–5429, https://doi.org/10.5194/gmd-17-5413-2024, https://doi.org/10.5194/gmd-17-5413-2024, 2024
Short summary
Short summary
This study uses a deep learning method to upscale the time resolution of paleoclimate simulations to 1 million years. This improved resolution allows a climate-biogeochemical model to more accurately predict climate shifts. The method may be critical in developing new fully continuous methods that are able to be applied over a moving continental surface in deep time with high resolution at reasonable computational expense.
Boris Ťupek, Aleksi Lehtonen, Alla Yurova, Rose Abramoff, Bertrand Guenet, Elisa Bruni, Samuli Launiainen, Mikko Peltoniemi, Shoji Hashimoto, Xianglin Tian, Juha Heikkinen, Kari Minkkinen, and Raisa Mäkipää
Geosci. Model Dev., 17, 5349–5367, https://doi.org/10.5194/gmd-17-5349-2024, https://doi.org/10.5194/gmd-17-5349-2024, 2024
Short summary
Short summary
Updating the Yasso07 soil C model's dependency on decomposition with a hump-shaped Ricker moisture function improved modelled soil organic C (SOC) stocks in a catena of mineral and organic soils in boreal forest. The Ricker function, set to peak at a rate of 1 and calibrated against SOC and CO2 data using a Bayesian approach, showed a maximum in well-drained soils. Using SOC and CO2 data together with the moisture only from the topsoil humus was crucial for accurate model estimates.
Jacquelyn K. Shuman, Rosie A. Fisher, Charles Koven, Ryan Knox, Lara Kueppers, and Chonggang Xu
Geosci. Model Dev., 17, 4643–4671, https://doi.org/10.5194/gmd-17-4643-2024, https://doi.org/10.5194/gmd-17-4643-2024, 2024
Short summary
Short summary
We adapt a fire behavior and effects module for use in a size-structured vegetation demographic model to test how climate, fire regime, and fire-tolerance plant traits interact to determine the distribution of tropical forests and grasslands. Our model captures the connection between fire disturbance and plant fire-tolerance strategies in determining plant distribution and provides a useful tool for understanding the vulnerability of these areas under changing conditions across the tropics.
Yoshiki Kanzaki, Isabella Chiaravalloti, Shuang Zhang, Noah J. Planavsky, and Christopher T. Reinhard
Geosci. Model Dev., 17, 4515–4532, https://doi.org/10.5194/gmd-17-4515-2024, https://doi.org/10.5194/gmd-17-4515-2024, 2024
Short summary
Short summary
Soil pH is one of the most commonly measured agronomical and biogeochemical indices, mostly reflecting exchangeable acidity. Explicit simulation of both porewater and bulk soil pH is thus crucial to the accurate evaluation of alkalinity required to counteract soil acidification and the resulting capture of anthropogenic carbon dioxide through the enhanced weathering technique. This has been enabled by the updated reactive–transport SCEPTER code and newly developed framework to simulate soil pH.
David Sandoval, Iain Colin Prentice, and Rodolfo L. B. Nóbrega
Geosci. Model Dev., 17, 4229–4309, https://doi.org/10.5194/gmd-17-4229-2024, https://doi.org/10.5194/gmd-17-4229-2024, 2024
Short summary
Short summary
Numerous estimates of water and energy balances depend on empirical equations requiring site-specific calibration, posing risks of "the right answers for the wrong reasons". We introduce novel first-principles formulations to calculate key quantities without requiring local calibration, matching predictions from complex land surface models.
Oliver Perkins, Matthew Kasoar, Apostolos Voulgarakis, Cathy Smith, Jay Mistry, and James D. A. Millington
Geosci. Model Dev., 17, 3993–4016, https://doi.org/10.5194/gmd-17-3993-2024, https://doi.org/10.5194/gmd-17-3993-2024, 2024
Short summary
Short summary
Wildfire is often presented in the media as a danger to human life. Yet globally, millions of people’s livelihoods depend on using fire as a tool. So, patterns of fire emerge from interactions between humans, land use, and climate. This complexity means scientists cannot yet reliably say how fire will be impacted by climate change. So, we developed a new model that represents globally how people use and manage fire. The model reveals the extent and diversity of how humans live with and use fire.
Amos P. K. Tai, David H. Y. Yung, and Timothy Lam
Geosci. Model Dev., 17, 3733–3764, https://doi.org/10.5194/gmd-17-3733-2024, https://doi.org/10.5194/gmd-17-3733-2024, 2024
Short summary
Short summary
We have developed the Terrestrial Ecosystem Model in R (TEMIR), which simulates plant carbon and pollutant uptake and predicts their response to varying atmospheric conditions. This model is designed to couple with an atmospheric chemistry model so that questions related to plant–atmosphere interactions, such as the effects of climate change, rising CO2, and ozone pollution on forest carbon uptake, can be addressed. The model has been well validated with both ground and satellite observations.
Fabian Stenzel, Johanna Braun, Jannes Breier, Karlheinz Erb, Dieter Gerten, Jens Heinke, Sarah Matej, Sebastian Ostberg, Sibyll Schaphoff, and Wolfgang Lucht
Geosci. Model Dev., 17, 3235–3258, https://doi.org/10.5194/gmd-17-3235-2024, https://doi.org/10.5194/gmd-17-3235-2024, 2024
Short summary
Short summary
We provide an R package to compute two biosphere integrity metrics that can be applied to simulations of vegetation growth from the dynamic global vegetation model LPJmL. The pressure metric BioCol indicates that we humans modify and extract > 20 % of the potential preindustrial natural biomass production. The ecosystems state metric EcoRisk shows a high risk of ecosystem destabilization in many regions as a result of climate change and land, water, and fertilizer use.
Elin Ristorp Aas, Heleen A. de Wit, and Terje K. Berntsen
Geosci. Model Dev., 17, 2929–2959, https://doi.org/10.5194/gmd-17-2929-2024, https://doi.org/10.5194/gmd-17-2929-2024, 2024
Short summary
Short summary
By including microbial processes in soil models, we learn how the soil system interacts with its environment and responds to climate change. We present a soil process model, MIMICS+, which is able to reproduce carbon stocks found in boreal forest soils better than a conventional land model. With the model we also find that when adding nitrogen, the relationship between soil microbes changes notably. Coupling the model to a vegetation model will allow for further study of these mechanisms.
Thomas Wutzler, Christian Reimers, Bernhard Ahrens, and Marion Schrumpf
Geosci. Model Dev., 17, 2705–2725, https://doi.org/10.5194/gmd-17-2705-2024, https://doi.org/10.5194/gmd-17-2705-2024, 2024
Short summary
Short summary
Soil microbes provide a strong link for elemental fluxes in the earth system. The SESAM model applies an optimality assumption to model those linkages and their adaptation. We found that a previous heuristic description was a special case of a newly developed more rigorous description. The finding of new behaviour at low microbial biomass led us to formulate the constrained enzyme hypothesis. We now can better describe how microbially mediated linkages of elemental fluxes adapt across decades.
Salvatore R. Curasi, Joe R. Melton, Elyn R. Humphreys, Txomin Hermosilla, and Michael A. Wulder
Geosci. Model Dev., 17, 2683–2704, https://doi.org/10.5194/gmd-17-2683-2024, https://doi.org/10.5194/gmd-17-2683-2024, 2024
Short summary
Short summary
Canadian forests are responding to fire, harvest, and climate change. Models need to quantify these processes and their carbon and energy cycling impacts. We develop a scheme that, based on satellite records, represents fire, harvest, and the sparsely vegetated areas that these processes generate. We evaluate model performance and demonstrate the impacts of disturbance on carbon and energy cycling. This work has implications for land surface modeling and assessing Canada’s terrestrial C cycle.
Yannek Käber, Florian Hartig, and Harald Bugmann
Geosci. Model Dev., 17, 2727–2753, https://doi.org/10.5194/gmd-17-2727-2024, https://doi.org/10.5194/gmd-17-2727-2024, 2024
Short summary
Short summary
Many forest models include detailed mechanisms of forest growth and mortality, but regeneration is often simplified. Testing and improving forest regeneration models is challenging. We address this issue by exploring how forest inventories from unmanaged European forests can be used to improve such models. We find that competition for light among trees is captured by the model, unknown model components can be informed by forest inventory data, and climatic effects are challenging to capture.
Jalisha T. Kallingal, Johan Lindström, Paul A. Miller, Janne Rinne, Maarit Raivonen, and Marko Scholze
Geosci. Model Dev., 17, 2299–2324, https://doi.org/10.5194/gmd-17-2299-2024, https://doi.org/10.5194/gmd-17-2299-2024, 2024
Short summary
Short summary
By unlocking the mysteries of CH4 emissions from wetlands, our work improved the accuracy of the LPJ-GUESS vegetation model using Bayesian statistics. Via assimilation of long-term real data from a wetland, we significantly enhanced CH4 emission predictions. This advancement helps us better understand wetland contributions to atmospheric CH4, which are crucial for addressing climate change. Our method offers a promising tool for refining global climate models and guiding conservation efforts
Cited articles
Andre, B., Molins, S., Johnson, J., and Steefel, C.: Alquimia, U.S. Department of Energy Office of Scientific and Technical Information, https://doi.org/10.11578/DC.20210416.49, 2013. a, b
Andre, B., Molins, S., Johnson, J., and Steefel, C. I.: Alquimia, GitHub, https://github.com/LBL-EESA/alquimia-dev (last access: 25 May 2025), 2015. a
Balos, C. J., Day, M., Esclapez, L., Felden, A. M., Gardner, D. J., Hassanaly, M., Reynolds, D. R., Rood, J. S., Sexton, J. M., Wimer, N. T., and Woodward, C. S.: SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations, Int. J. High Perform. C., 39, 123–146, https://doi.org/10.1177/10943420241280060, 2025. a
Bao, C., Li, L., Shi, Y., and Duffy, C.: Understanding watershed hydrogeochemistry: 1. Development of RT?Flux?PIHM, Water Resour. Res., 53, 2328–2345, https://doi.org/10.1002/2016WR018934, 2017. a
Bartlett, R., Demeshko, I., Gamblin, T., Hammond, G., Heroux, M. A., Johnson, J., Klinvex, A., Li, X., McInnes, L. C., Moulton, J. D., Osei-Kuffuor, D., Sarich, J., Smith, B., Willenbring, J., and Yang, U. M.: xSDK Foundations: Toward an Extreme-scale Scientific Software Development Kit, Supercomputing Frontiers and Innovations, 4, 69–82, https://doi.org/10.14529/jsfi170104, 2017. a, b
Beisman, J., Maxwell, R., Navarre-Sitchler, A., Steefel, C., and Molins, S.: ParCrunchFlow: an efficient, parallel reactive transport simulation tool for physically and chemically heterogeneous saturated subsurface environments, Comput. Geosci., 19, 403–422, https://doi.org/10.1007/s10596-015-9475-x, 2015. a, b, c
Bethke, C. M.: The Geochemist's Workbench® Release 17 ChemPlugin™ User's Guide, Aqueous Solutions, LLC, Champaign, Illinois, https://www.gwb.com/pdf/GWB/ChemPluginUsersGuide.pdf (last access: 24 May 2025), 2024. a
CFFI: CFFI documentation – CFFI 1.15.1 documentation, https://cffi.readthedocs.io/en/latest/ (last access: 24 May 2025), 2023. a
Chang, E., Zavarin, M., Beverly, L., and Wainwright, H.: A chemistry-informed hybrid machine learning approach to predict metal adsorption onto mineral surfaces, Appl. Geochem., 155, 105731, https://doi.org/10.1016/j.apgeochem.2023.105731, 2023. a
Charlton, S. R. and Parkhurst, D. L.: Modules based on the geochemical model PHREEQC for use in scripting and programming languages, Comput. Geosci., 37, 1653–1663, https://doi.org/10.1016/j.cageo.2011.02.005, 2011. a
Coon, E., Svyatsky, D., Jan, A., Kikinzon, E., Berndt, M., Atchley, A., Harp, D., Manzini, G., Shelef, E., Lipnikov, K., Garimella, R., Xu, C., Moulton, D., Karra, S., Painter, S., Jafarov, E., and Molins, S.: Advanced Terrestrial Simulator, U.S. Department of Energy Office of Scientific and Technical Information, https://doi.org/10.11578/DC.20190911.1, 2019. a, b
De Lucia, M., Kühn, M., Lindemann, A., Lübke, M., and Schnor, B.: POET (v0.1): speedup of many-core parallel reactive transport simulations with fast DHT lookups, Geosci. Model Dev., 14, 7391–7409, https://doi.org/10.5194/gmd-14-7391-2021, 2021. a
Dubey, A., Almgren, A., Bell, J., Berzins, M., Brandt, S., Bryan, G., Colella, P., Graves, D., Lijewski, M., L”offler, F., O'Shea, B., Schnetter, E., Van Straalen, B., and Weide, K.: A survey of high level frameworks in block-structured adaptive mesh refinement packages, J. Parallel and Distr. Com., 74, 3217–3227, https://doi.org/10.1016/j.jpdc.2014.07.001, 2014. a, b, c
GitHub: GitHub Actions documentation, https://docs.github.com/en/actions (last access: 24 May 2025), 2024. a
Greskowiak, J., Gwo, J., Jacques, D., Yin, J., and Mayer, K. U.: A benchmark for multi-rate surface complexation and 1D dual-domain multi-component reactive transport of U(VI), Comput. Geosci., 19, 585–597, https://doi.org/10.1007/s10596-014-9457-4, 2015. a
Gupta, A. D., Lake, L. W., Pope, G. A., Sepehrnoori, K. T. U., and King, M. J. B. R.: High-Resolution Monotonic Schemes for Reservoir Fluid Flow Simulation, In Situ, (United States), 15, https://www.osti.gov/scitech/biblio/5799660 (last access: 24 May 2025), 1991. a
Hammond, G. E.: The PFLOTRAN Reaction Sandbox, Geosci. Model Dev., 15, 1659–1676, https://doi.org/10.5194/gmd-15-1659-2022, 2022. a, b, c, d
Hammond, G. E., Lichtner, P. C., and Mills, R. T.: Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN, Water Resour. Res., 50, 208–228, https://doi.org/10.1002/2012WR013483, 2014. a
Jara, D., de Dreuzy, J.-R., and Cochepin, B.: TReacLab: An object-oriented implementation of non-intrusive splitting methods to couple independent transport and geochemical software, Comput. Geosci., 109, 281–294, https://doi.org/10.1016/j.cageo.2017.09.005, 2017. a
Jaysaval, P., Hammond, G. E., and Johnson, T. C.: Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN, Geosci. Model Dev., 16, 961–976, https://doi.org/10.5194/gmd-16-961-2023, 2023. a
LaFond-Hudson, S. and Sulman, B.: Modeling strategies and data needs for representing coastal wetland vegetation in land surface models, New Phytol., 238, 938–951, https://doi.org/10.1111/nph.18760, 2023. a, b, c
Leal, A. M. M., Kyas, S., Kulik, D. A., and Saar, M. O.: Accelerating Reactive Transport Modeling: On-Demand Machine Learning Algorithm for Chemical Equilibrium Calculations, Transport Porous Med., 133, 161–204, https://doi.org/10.1007/s11242-020-01412-1, 2020. a
Li, L., Maher, K., Navarre-Sitchler, A., Druhan, J., Meile, C., Lawrence, C., Moore, J., Perdrial, J., Sullivan, P., Thompson, A., Jin, L., Bolton, E. W., Brantley, S. L., Dietrich, W. E., Mayer, K. U., Steefel, C. I., Valocchi, A., Zachara, J., Kocar, B., Mcintosh, J., Bao, C., Tutolo, B. M., Beisman, J., Kumar, M., and Sonnenthal, E.: Expanding the role of reactive transport models in critical zone processes, Earth-Sci. Rev., 165, 280–301, https://doi.org/10.1016/j.earscirev.2016.09.001, 2017. a, b
Lichtner, P. C.: The Quasi-Stationary State Approximation to Fluid/Rock Reaction: Local Equilibrium Revisited, in: Diffusion, Atomic Ordering, and Mass Transport: Selected Topics in Geochemistry, edited by: Ganguly, J., Advances in Physical Geochemistry, Springer US, New York, NY, 452–560, ISBN 978-1-4613-9019-0, https://doi.org/10.1007/978-1-4613-9019-0_13, 1991. a
Lichtner, P. C., Hammond, G. E., Lu, C., Karra, S., Bisht, G., Andre, B. N. C. F. A. R., Lawrence Berkeley National Lab. (LBNL), B., Mills, R. I. C., Univ. of Tennessee, K., and Kumar, J.: Pflotran User Manual: A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes, Tech. Rep. LA-UR–15-20403, Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States), Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States), Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States), OFM Research, Redmond, WA, USA, https://doi.org/10.2172/1168703, 2015. a
Maxwell, R. M. and Miller, N. L.: Development of a Coupled Land Surface and Groundwater Model, J. Hydrometeorol., 6, 233–247, https://doi.org/10.1175/JHM422.1, 2005. a
Maxwell, R. M., Putti, M., Meyerhoff, S., Delfs, J.-O., Ferguson, I. M., Ivanov, V., Kim, J., Kolditz, O., Kollet, S. J., Kumar, M., Lopez, S., Niu, J., Paniconi, C., Park, Y.-J., Phanikumar, M. S., Shen, C., Sudicky, E. A., and Sulis, M.: Surface-subsurface model intercomparison: A first set of benchmark results to diagnose integrated hydrology and feedbacks, Water Resour. Res., 50, 1531–1549, https://doi.org/10.1002/2013WR013725, 2014. a
Mayer, K. U., Frind, E. O., and Blowes, D. W.: Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions, Water Resour. Res., 38, 1174, https://doi.org/10.1029/2001WR000862, 2002. a
Molins, S.: Reactive Interfaces in Direct Numerical Simulation of Pore-Scale Processes, Rev. Mineral. Geochem., 80, 461–481, https://doi.org/10.2138/rmg.2015.80.14, 2015. a
Molins, S., Trebotich, D., Steefel, C. I., and Shen, C.: An investigation of the effect of pore scale flow on average geochemical reaction rates using direct numerical simulation, Water Resour. Res., 48, https://doi.org/10.1029/2011WR011404, 2012. a, b
Molins, S., Soulaine, C., Prasianakis, N. I., Abbasi, A., Poncet, P., Ladd, A. J. C., Starchenko, V., Roman, S., Trebotich, D., Tchelepi, H. A., and Steefel, C. I.: Simulation of mineral dissolution at the pore scale with evolving fluid-solid interfaces: review of approaches and benchmark problem set, Comput. Geosci., 25, 1285–1318, https://doi.org/10.1007/s10596-019-09903-x, 2020. a, b
Molins, S., Andre, B., Johnson, J., Hammond, G., Sulman, B., Lipnikov, K., Day, M., Beisman, J., Svyatskiy, D., Deng, H., Lichtner, P., Steefel, C., and Moulton, D.: Alquimia: A generic interface to biogeochemical codes – A tool for interoperable development, prototyping and benchmarking for multiphysics simulators, Zenodo [code], https://doi.org/10.5281/zenodo.11414442, 2024a. a
Molins, S., Trebotich, D., and Steefel, C. I.: Approaches for the simulation of coupled processes in evolving fractured porous media enabled by exascale computing, Comput. Sci. Eng., 26, 33–42, https://doi.org/10.1109/MCSE.2024.3403983, 2024b. a
Nardi, A., Idiart, A., Trinchero, P., de Vries, L. M., and Molinero, J.: Interface COMSOL-PHREEQC (iCP), an efficient numerical framework for the solution of coupled multiphysics and geochemistry, Comput. Geosci., 69, 10–21, https://doi.org/10.1016/j.cageo.2014.04.011, 2014. a
OpenFOAM: Open-source Field Operation And Manipulation, https://www.openfoam.com/ (last access: 24 May 2025), 2022. a
Parkhurst, D. L. and Wissmeier, L.: PhreeqcRM: A reaction module for transport simulators based on the geochemical model PHREEQC, Adv. Water Resour., 83, 176–189, https://doi.org/10.1016/j.advwatres.2015.06.001, 2015. a
Parkhurst, D. L., Kipp, K. L., and Charlton, S. R.: PHAST Version 2 – A program for simulating groundwater flow, solute transport, and multicomponent geochemical reactions, in: U.S. Geological Survey Techniques and Methods, 6-A35, 235, https://doi.org/10.3133/tm6A35, 2010. a
Prommer, H. and Post, V.: PHT3D: A ReactiveMulticomponent Transport Model for Saturated Porous Media, User's Manual v2.10, http://www.pht3d.org (last access: 24 May 2025), 2010. a
Simunek, J., Jacques, D., Langergraber, G., Bradford, S. A., Šejna, M., and Genuchten, M. T. v.: Numerical Modeling of Contaminant Transport Using HYDRUS and its Specialized Modules, J. Indian I. Sci., 93, 265–284, https://journal.iisc.ac.in/index.php/iisc/article/view/1224 (last access: 24 May 2025), 2013. a
Sonnenthal, E., Spycher, N., Xu, T., and Zheng, L.: TOUGHREACT V4.12-OMP and TReactMech V1.0 Geochemical and Reactive-Transport User Guide, https://escholarship.org/uc/item/8945d2c1 (last access: 24 May 2025), 2021. a
Sphinx: Welcome – Sphinx documentation, https://www.sphinx-doc.org/en/master/ (last access: 24 May 2025), 2024. a
Steefel, C. I., DePaolo, D. J., and Lichtner, P. C.: Reactive transport modeling: An essential tool and a new research approach for the Earth sciences, Earth Planet. Sc. Lett., 240, 539–558, https://doi.org/10.1016/j.epsl.2005.09.017, 2005. a
Steefel, C. I., Appelo, C. A. J., Arora, B., Jacques, D., Kalbacher, T., Kolditz, O., Lagneau, V., Lichtner, P. C., Mayer, K. U., Meeussen, J. C. L., Molins, S., Moulton, D., Shao, H., Šimůnek, J., Spycher, N., Yabusaki, S. B., and Yeh, G. T.: Reactive transport codes for subsurface environmental simulation, Comput. Geosci., 19, 445–478, https://doi.org/10.1007/s10596-014-9443-x, 2015. a, b, c, d, e, f, g, h, i
Sulman, B., Yuan, F., O'Meara, T., Graham, D., Gu, B., Herndon, E., and Zheng, J.: Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive, ESS-Dive [code], https://doi.org/10.5440/1814844, 2020. a, b
Sulman, B., Wang, J., LaFond-Hudson, S., O'Meara, T., Yuan, F., Molins, S., Forbrich, I., Cardon, Z., and Giblin, A.: Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN, ESS-Dive [code], https://doi.org/10.15485/1991625, 2023. a
Sulman, B. N., Yuan, F., O'Meara, T., Gu, B., Herndon, E. M., Zheng, J., Thornton, P. E., and Graham, D. E.: Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil, J. Geophys. Res.-Biogeo., 127, e2021JG006662, https://doi.org/10.1029/2021JG006662, 2022. a, b, c
Sulman, B. N., Wang, J., LaFond-Hudson, S., O'Meara, T., Yuan, F., Molins, S., Hammond, G. E., Forbrich, I., Cardon, Z., and Giblin, A.: Integrating tide-driven wetland soil redox and biogeochemical interactions into a land surface model, J. Adv. Model. Earth Sy., 16, e2023MS004002, https://doi.org/10.1029/2023MS004002, 2024. a, b, c, d, e, f
Tang, G., Yuan, F., Bisht, G., Hammond, G. E., Lichtner, P. C., Kumar, J., Mills, R. T., Xu, X., Andre, B., Hoffman, F. M., Painter, S. L., and Thornton, P. E.: Addressing numerical challenges in introducing a reactive transport code into a land surface model: a biogeochemical modeling proof-of-concept with CLM–PFLOTRAN 1.0, Geosci. Model Dev., 9, 927–946, https://doi.org/10.5194/gmd-9-927-2016, 2016. a
van der Lee, J., De Windt, L., Lagneau, V., and Goblet, P.: Module-oriented modeling of reactive transport with HYTEC, Comput. Geosci., 29, 265–275, https://doi.org/10.1016/S0098-3004(03)00004-9, 2003. a
Wang, J., O'Meara, T., LaFond-Hudson, S., He, S., Maiti, K., Ward, E. J., and Sulman, B. N.: Subsurface Redox Interactions Regulate Ebullitive Methane Flux in Heterogeneous Mississippi River Deltaic Wetland, J. Adv. Model. Earth Sy., 16, e2023MS003762, https://doi.org/10.1029/2023MS003762, 2024. a, b, c
Zhang, Q., Deng, H., Dong, Y., Molins, S., Li, X., and Steefel, C.: Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro-Continuum Reactive Transport Model, Water Resour. Res., 58, e2021WR030578, https://doi.org/10.1029/2021WR030578, 2022. a
Zhang, Q., Dong, Y., Molins, S., and Deng, H.: The Impacts of Micro-Porosity and Mineralogical Texture on Fractured Rock Alteration, Water Resour. Res., 60, e2023WR036266, https://doi.org/10.1029/2023WR036266, 2024. a
Short summary
Developing scientific software and making sure it functions properly requires a significant effort. As we advance our understanding of natural systems, however, there is the need to develop yet more complex models and codes. In this work, we present a piece of software that facilitates this work, specifically with regard to reactive processes. Existing tried-and-true codes are made available via this new interface, freeing up resources to focus on the new aspects of the problems at hand.
Developing scientific software and making sure it functions properly requires a significant...