Articles | Volume 9, issue 10
https://doi.org/10.5194/gmd-9-3729-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/gmd-9-3729-2016
© Author(s) 2016. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Development and technical paper
| 
19 Oct 2016
Development and technical paper |
| 19 Oct 2016
Metos3D: the Marine Ecosystem Toolkit for Optimization and Simulation in 3-D – Part 1: Simulation Package v0.3.2
Jaroslaw Piwonski
CORRESPONDING AUTHOR
Institute for Computer Science and Kiel Marine Science – Centre for Interdisciplinary Marine Science,
Cluster The Future Ocean, Kiel University, 24098 Kiel, Germany
Thomas Slawig
Institute for Computer Science and Kiel Marine Science – Centre for Interdisciplinary Marine Science,
Cluster The Future Ocean, Kiel University, 24098 Kiel, Germany
Related authors
No articles found.
Niko Schmidt, Angelika Humbert, and Thomas Slawig
Geosci. Model Dev., 17, 4943–4959, https://doi.org/10.5194/gmd-17-4943-2024, https://doi.org/10.5194/gmd-17-4943-2024, 2024
Short summary
Short summary
Future sea-level rise is of big significance for coastal regions. The melting and acceleration of glaciers plays a major role in sea-level change. Computer simulation of glaciers costs a lot of computational resources. In this publication, we test a new way of simulating glaciers. This approach produces the same results but has the advantage that it needs much less computation time. As simulations can be obtained with fewer computation resources, higher resolution and physics become affordable.
Maria-Theresia Pelz, Markus Schartau, Christopher J. Somes, Vanessa Lampe, and Thomas Slawig
Geosci. Model Dev., 16, 6609–6634, https://doi.org/10.5194/gmd-16-6609-2023, https://doi.org/10.5194/gmd-16-6609-2023, 2023
Short summary
Short summary
Kernel density estimators (KDE) approximate the probability density of a data set without the assumption of an underlying distribution. We used the solution of the diffusion equation, and a new approximation of the optimal smoothing parameter build on two pilot estimation steps, to construct such a KDE best suited for typical characteristics of geoscientific data. The resulting KDE is insensitive to noise and well resolves multimodal data structures as well as boundary-close data.
Markus Pfeil and Thomas Slawig
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2021-392, https://doi.org/10.5194/gmd-2021-392, 2022
Revised manuscript not accepted
Short summary
Short summary
In investigating the global carbon cycle, shortening the runtime of the simulation of marine ecosystem models is an important issue. We present methods that automatically adjust the time step during the simulation of a steady state using transport matrices. They apply always the time step as large as possible. Two methods reduced the runtime significantly, depending on the complexity of the model. An important property was that small negative concentrations were ignored during the spin-up.
Maria-Theresia Verwega, Christopher J. Somes, Markus Schartau, Robyn Elizabeth Tuerena, Anne Lorrain, Andreas Oschlies, and Thomas Slawig
Earth Syst. Sci. Data, 13, 4861–4880, https://doi.org/10.5194/essd-13-4861-2021, https://doi.org/10.5194/essd-13-4861-2021, 2021
Short summary
Short summary
This work describes a ready-to-use collection of particulate organic carbon stable isotope ratio data sets. It covers the 1960s–2010s and all main oceans, providing meta-information and gridded data. The best coverage exists in Atlantic, Indian and Southern Ocean surface waters during the 1990s. It indicates no major difference between methods and shows decreasing values towards high latitudes, with the lowest in the Southern Ocean, and a long-term decline in all regions but the Southern Ocean.
Alessandro Cotronei and Thomas Slawig
Geosci. Model Dev., 13, 2783–2804, https://doi.org/10.5194/gmd-13-2783-2020, https://doi.org/10.5194/gmd-13-2783-2020, 2020
Short summary
Short summary
We converted the radiation part of the atmospheric model ECHAM to single-precision arithmetic, using a step-by-step change in all modules. A small code portion still requires higher precision. The generated code can be easily changed from double to single precision and vice versa. The quality of the output of the single-precision version is comparable to observational data and the one of the original code. The runtime was reduced by 40 %, and the energy consumption could also be decreased.
Markus Schartau, Philip Wallhead, John Hemmings, Ulrike Löptien, Iris Kriest, Shubham Krishna, Ben A. Ward, Thomas Slawig, and Andreas Oschlies
Biogeosciences, 14, 1647–1701, https://doi.org/10.5194/bg-14-1647-2017, https://doi.org/10.5194/bg-14-1647-2017, 2017
Short summary
Short summary
Plankton models have become an integral part in marine ecosystem and biogeochemical research. These models differ in complexity and in their number of parameters. How values are assigned to parameters is essential. An overview of major methodologies of parameter estimation is provided. Aspects of parameter identification in the literature are diverse. Individual findings could be better synthesized if notation and expertise of the different scientific communities would be reasonably merged.
J. Reimer, M. Schuerch, and T. Slawig
Geosci. Model Dev., 8, 791–804, https://doi.org/10.5194/gmd-8-791-2015, https://doi.org/10.5194/gmd-8-791-2015, 2015
Short summary
Short summary
Model parameters are usually optimized based on measurements. These measurements are often time-consuming or costly. The conditions under which theses measurements are carried out, also called experimental designs, can be optimized so that with minimum effort and cost a maximum accuracy can be achieved. For this, we present different approaches together with their implementation in an MATLAB toolbox. We demonstrate their application to different models for sedimentation in salt marshes.
M. El Jarbi, J. Rückelt, T. Slawig, and A. Oschlies
Biogeosciences, 10, 1169–1182, https://doi.org/10.5194/bg-10-1169-2013, https://doi.org/10.5194/bg-10-1169-2013, 2013
E. Siewertsen, J. Piwonski, and T. Slawig
Geosci. Model Dev., 6, 17–28, https://doi.org/10.5194/gmd-6-17-2013, https://doi.org/10.5194/gmd-6-17-2013, 2013
Related subject area
Biogeosciences
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
Quantifying the role of ozone-caused damage to vegetation in the Earth system: A new parameterization scheme for photosynthetic and stomatal responses
Optimising CH4 simulations from the LPJ-GUESS model v4.1 using an adaptive Markov chain Monte Carlo algorithm
The XSO framework (v0.1) and Phydra library (v0.1) for a flexible, reproducible, and integrated plankton community modeling environment in Python
Exploring the Potential of History Matching for Land Surface Model Calibration
AgriCarbon-EO v1.0.1: large-scale and high-resolution simulation of carbon fluxes by assimilation of Sentinel-2 and Landsat-8 reflectances using a Bayesian approach
SAMM version 1.0: a numerical model for microbial- mediated soil aggregate formation
Using Deep Learning to integrate paleoclimate and global biogeochemistry over Phanerozoic time
A model of the within-population variability of budburst in forest trees
Computationally efficient parameter estimation for high-dimensional ocean biogeochemical models
Radiocarbon analysis reveals underestimation of soil organic carbon persistence in new-generation soil models
The community-centered freshwater biogeochemistry model unified RIVE v1.0: a unified version for water column
EAT v0.9.6: a 1D testbed for physical-biogeochemical data assimilation in natural waters
Observation-based sowing dates and cultivars significantly affect yield and irrigation for some crops in the Community Land Model (CLM5)
The statistical emulators of GGCMI phase 2: responses of year-to-year variation of crop yield to CO2, temperature, water, and nitrogen perturbations
A novel Eulerian model based on central moments to simulate age and reactivity continua interacting with mixing processes
AdaScape 1.0: a coupled modelling tool to investigate the links between tectonics, climate, and biodiversity
An along-track Biogeochemical Argo modelling framework: a case study of model improvements for the Nordic seas
Peatland-VU-NUCOM (PVN 1.0): using dynamic plant functional types to model peatland vegetation, CH4, and CO2 emissions
Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)
SedTrace 1.0: a Julia-based framework for generating and running reactive-transport models of marine sediment diagenesis specializing in trace elements and isotopes
A high-resolution marine mercury model MITgcm-ECCO2-Hg with online biogeochemistry
Improving nitrogen cycling in a land surface model (CLM5) to quantify soil N2O, NO, and NH3 emissions from enhanced rock weathering with croplands
Assimilation of Carbonyl Sulfide (COS) fluxes within the adjoint-based data assimilation system–Nanjing University Carbon Assimilation System (NUCAS v1.0)
Ocean biogeochemistry in the coupled ocean–sea ice–biogeochemistry model FESOM2.1–REcoM3
Forcing the Global Fire Emissions Database burned-area dataset into the Community Land Model version 5.0: impacts on carbon and water fluxes at high latitudes
Modeling boreal forest’s mineral soil and peat C stock dynamics with Yasso07 model coupled with updated moisture modifier
Modeling of non-structural carbohydrate dynamics by the spatially explicit individual-based dynamic global vegetation model SEIB-DGVM (SEIB-DGVM-NSC version 1.0)
MEDFATE 2.9.3: a trait-enabled model to simulate Mediterranean forest function and dynamics at regional scales
Modelling the role of livestock grazing in C and N cycling in grasslands with LPJmL5.0-grazing
Implementation of trait-based ozone plant sensitivity in the Yale Interactive terrestrial Biosphere model v1.0 to assess global vegetation damage
The Permafrost and Organic LayEr module for Forest Models (POLE-FM) 1.0
CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics
Validation of a new spatially explicit process-based model (HETEROFOR) to simulate structurally and compositionally complex forest stands in eastern North America
Global agricultural ammonia emissions simulated with the ORCHIDEE land surface model
ForamEcoGEnIE 2.0: incorporating symbiosis and spine traits into a trait-based global planktic foraminiferal model
FABM-NflexPD 2.0: testing an instantaneous acclimation approach for modeling the implications of phytoplankton eco-physiology for the carbon and nutrient cycles
Evaluating the vegetation–atmosphere coupling strength of ORCHIDEE land surface model (v7266)
Non-Redfieldian carbon model for the Baltic Sea (ERGOM version 1.2) – implementation and budget estimates
Implementation of a new crop phenology and irrigation scheme in the ISBA land surface model using SURFEX_v8.1
Simulating long-term responses of soil organic matter turnover to substrate stoichiometry by abstracting fast and small-scale microbial processes: the Soil Enzyme Steady Allocation Model (SESAM; v3.0)
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.
Fang Li, Zhimin Zhou, Samuel Levis, Stephen Sitch, Felicity Hayes, Zhaozhong Feng, Peter Reich, Zhiyi Zhao, and Yanqing Zhou
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-6, https://doi.org/10.5194/gmd-2024-6, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study introduces a new scheme to simulate 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 response 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.
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
Benjamin Post, Esteban Acevedo-Trejos, Andrew D. Barton, and Agostino Merico
Geosci. Model Dev., 17, 1175–1195, https://doi.org/10.5194/gmd-17-1175-2024, https://doi.org/10.5194/gmd-17-1175-2024, 2024
Short summary
Short summary
Creating computational models of how phytoplankton grows in the ocean is a technical challenge. We developed a new tool set (Xarray-simlab-ODE) for building such models using the programming language Python. We demonstrate the tool set in a library of plankton models (Phydra). Our goal was to allow scientists to develop models quickly, while also allowing the model structures to be changed easily. This allows us to test many different structures of our models to find the most appropriate one.
Nina Raoult, Simon Beylat, James M. Salter, Frédéric Hourdin, Vladislav Bastrikov, Catherine Ottlé, and Philippe Peylin
EGUsphere, https://doi.org/10.5194/egusphere-2023-2996, https://doi.org/10.5194/egusphere-2023-2996, 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. Here we test a technique called “History matching” against more common approaches. This method adapts well to these models, helping better understand how they work and therefore how to make them more realistic.
Taeken Wijmer, Ahmad Al Bitar, Ludovic Arnaud, Remy Fieuzal, and Eric Ceschia
Geosci. Model Dev., 17, 997–1021, https://doi.org/10.5194/gmd-17-997-2024, https://doi.org/10.5194/gmd-17-997-2024, 2024
Short summary
Short summary
Quantification of carbon fluxes of crops is an essential building block for the construction of a monitoring, reporting, and verification approach. We developed an end-to-end platform (AgriCarbon-EO) that assimilates, through a Bayesian approach, high-resolution (10 m) optical remote sensing data into radiative transfer and crop modelling at regional scale (100 x 100 km). Large-scale estimates of carbon flux are validated against in situ flux towers and yield maps and analysed at regional scale.
Moritz Laub, Sergey Blagodatsky, Marijn Van de Broek, Samuel Schlichenmaier, Benjapon Kunlanit, Johan Six, Patma Vityakon, and Georg Cadisch
Geosci. Model Dev., 17, 931–956, https://doi.org/10.5194/gmd-17-931-2024, https://doi.org/10.5194/gmd-17-931-2024, 2024
Short summary
Short summary
To manage soil organic matter (SOM) sustainably, we need a better understanding of the role that soil microbes play in aggregate protection. Here, we propose the SAMM model, which connects soil aggregate formation to microbial growth. We tested it against data from a tropical long-term experiment and show that SAMM effectively represents the microbial growth, SOM, and aggregate dynamics and that it can be used to explore the importance of aggregate formation in SOM stabilization.
Dongyu Zheng, Andrew Merdith, Yves Goddéris, Yannick Donnadieu, Khushboo Gurung, and Benjamin J. W. Mills
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-230, https://doi.org/10.5194/gmd-2023-230, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study uses a deep learning method to enhance the time resolution of paleoclimate simulations to 1 million years. This improved resolution allows a climate-biogeochemical model for more accurate predictions of climate shifts. The method may be critical in developing new fully-continuous methods that are able to translate over a moving continental surface in deep time with high-resolution at a reasonable computational expense.
Jianhong Lin, Daniel Berveiller, Christophe François, Heikki Hänninen, Alexandre Morfin, Gaëlle Vincent, Rui Zhang, Cyrille Rathgeber, and Nicolas Delpierre
Geosci. Model Dev., 17, 865–879, https://doi.org/10.5194/gmd-17-865-2024, https://doi.org/10.5194/gmd-17-865-2024, 2024
Short summary
Short summary
Currently, the high variability of budburst between individual trees is overlooked. The consequences of this neglect when projecting the dynamics and functioning of tree communities are unknown. Here we develop the first process-oriented model to describe the difference in budburst dates between individual trees in plant populations. Beyond budburst, the model framework provides a basis for studying the dynamics of phenological traits under climate change, from the individual to the community.
Skyler Kern, Mary E. McGuinn, Katherine M. Smith, Nadia Pinardi, Kyle E. Niemeyer, Nicole S. Lovenduski, and Peter E. Hamlington
Geosci. Model Dev., 17, 621–649, https://doi.org/10.5194/gmd-17-621-2024, https://doi.org/10.5194/gmd-17-621-2024, 2024
Short summary
Short summary
Computational models are used to simulate the behavior of marine ecosystems. The models often have unknown parameters that need to be calibrated to accurately represent observational data. Here, we propose a novel approach to simultaneously determine a large set of parameters for a one-dimensional model of a marine ecosystem in the surface ocean at two contrasting sites. By utilizing global and local optimization techniques, we estimate many parameters in a computationally efficient manner.
Alexander S. Brunmayr, Frank Hagedorn, Margaux Moreno Duborgel, Luisa I. Minich, and Heather D. Graven
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-242, https://doi.org/10.5194/gmd-2023-242, 2024
Revised manuscript accepted for GMD
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 as 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.
Shuaitao Wang, Vincent Thieu, Gilles Billen, Josette Garnier, Marie Silvestre, Audrey Marescaux, Xingcheng Yan, and Nicolas Flipo
Geosci. Model Dev., 17, 449–476, https://doi.org/10.5194/gmd-17-449-2024, https://doi.org/10.5194/gmd-17-449-2024, 2024
Short summary
Short summary
This paper presents unified RIVE v1.0, a unified version of the freshwater biogeochemistry model RIVE. It harmonizes different RIVE implementations, providing the referenced formalisms for microorganism activities to describe full biogeochemical cycles in the water column (e.g., carbon, nutrients, oxygen). Implemented as open-source projects in Python 3 (pyRIVE 1.0) and ANSI C (C-RIVE 0.32), unified RIVE v1.0 promotes and enhances collaboration among research teams and public services.
Jorn Bruggeman, Karsten Bolding, Lars Nerger, Anna Teruzzi, Simone Spada, Jozef Skákala, and Stefano Ciavatta
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2023-238, https://doi.org/10.5194/gmd-2023-238, 2023
Revised manuscript accepted for GMD
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. Here, we present EAT, a flexible and efficient testbed that allows any scientist to explore and further develop the state-of-the-art in data assimilation.
Sam S. Rabin, William J. Sacks, Danica L. Lombardozzi, Lili Xia, and Alan Robock
Geosci. Model Dev., 16, 7253–7273, https://doi.org/10.5194/gmd-16-7253-2023, https://doi.org/10.5194/gmd-16-7253-2023, 2023
Short summary
Short summary
Climate models can help us simulate how the agricultural system will be affected by and respond to environmental change, but to be trustworthy they must realistically reproduce historical patterns. When farmers plant their crops and what varieties they choose will be important aspects of future adaptation. Here, we improve the crop component of a global model to better simulate observed growing seasons and examine the impacts on simulated crop yields and irrigation demand.
Weihang Liu, Tao Ye, Christoph Müller, Jonas Jägermeyr, James A. Franke, Haynes Stephens, and Shuo Chen
Geosci. Model Dev., 16, 7203–7221, https://doi.org/10.5194/gmd-16-7203-2023, https://doi.org/10.5194/gmd-16-7203-2023, 2023
Short summary
Short summary
We develop a machine-learning-based crop model emulator with the inputs and outputs of multiple global gridded crop model ensemble simulations to capture the year-to-year variation of crop yield under future climate change. The emulator can reproduce the year-to-year variation of simulated yield given by the crop models under CO2, temperature, water, and nitrogen perturbations. Developing this emulator can provide a tool to project future climate change impact in a simple way.
Jurjen Rooze, Heewon Jung, and Hagen Radtke
Geosci. Model Dev., 16, 7107–7121, https://doi.org/10.5194/gmd-16-7107-2023, https://doi.org/10.5194/gmd-16-7107-2023, 2023
Short summary
Short summary
Chemical particles in nature have properties such as age or reactivity. Distributions can describe the properties of chemical concentrations. In nature, they are affected by mixing processes, such as chemical diffusion, burrowing animals, and bottom trawling. We derive equations for simulating the effect of mixing on central moments that describe the distributions. We then demonstrate applications in which these equations are used to model continua in disturbed natural environments.
Esteban Acevedo-Trejos, Jean Braun, Katherine Kravitz, N. Alexia Raharinirina, and Benoît Bovy
Geosci. Model Dev., 16, 6921–6941, https://doi.org/10.5194/gmd-16-6921-2023, https://doi.org/10.5194/gmd-16-6921-2023, 2023
Short summary
Short summary
The interplay of tectonics and climate influences the evolution of life and the patterns of biodiversity we observe on earth's surface. Here we present an adaptive speciation component coupled with a landscape evolution model that captures the essential earth-surface, ecological, and evolutionary processes that lead to the diversification of taxa. We can illustrate with our tool how life and landforms co-evolve to produce distinct biodiversity patterns on geological timescales.
Veli Çağlar Yumruktepe, Erik Askov Mousing, Jerry Tjiputra, and Annette Samuelsen
Geosci. Model Dev., 16, 6875–6897, https://doi.org/10.5194/gmd-16-6875-2023, https://doi.org/10.5194/gmd-16-6875-2023, 2023
Short summary
Short summary
We present an along BGC-Argo track 1D modelling framework. The model physics is constrained by the BGC-Argo temperature and salinity profiles to reduce the uncertainties related to mixed layer dynamics, allowing the evaluation of the biogeochemical formulation and parameterization. We objectively analyse the model with BGC-Argo and satellite data and improve the model biogeochemical dynamics. We present the framework, example cases and routines for model improvement and implementations.
Tanya J. R. Lippmann, Ype van der Velde, Monique M. P. D. Heijmans, Han Dolman, Dimmie M. D. Hendriks, and Ko van Huissteden
Geosci. Model Dev., 16, 6773–6804, https://doi.org/10.5194/gmd-16-6773-2023, https://doi.org/10.5194/gmd-16-6773-2023, 2023
Short summary
Short summary
Vegetation is a critical component of carbon storage in peatlands but an often-overlooked concept in many peatland models. We developed a new model capable of simulating the response of vegetation to changing environments and management regimes. We evaluated the model against observed chamber data collected at two peatland sites. We found that daily air temperature, water level, harvest frequency and height, and vegetation composition drive methane and carbon dioxide emissions.
Chonggang Xu, Bradley Christoffersen, Zachary Robbins, Ryan Knox, Rosie A. Fisher, Rutuja Chitra-Tarak, Martijn Slot, Kurt Solander, Lara Kueppers, Charles Koven, and Nate McDowell
Geosci. Model Dev., 16, 6267–6283, https://doi.org/10.5194/gmd-16-6267-2023, https://doi.org/10.5194/gmd-16-6267-2023, 2023
Short summary
Short summary
We introduce a plant hydrodynamic model for the U.S. Department of Energy (DOE)-sponsored model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES). To better understand this new model system and its functionality in tropical forest ecosystems, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We identified the key parameters that affect the simulated plant hydrodynamics to guide both modeling and field campaign studies.
Jianghui Du
Geosci. Model Dev., 16, 5865–5894, https://doi.org/10.5194/gmd-16-5865-2023, https://doi.org/10.5194/gmd-16-5865-2023, 2023
Short summary
Short summary
Trace elements and isotopes (TEIs) are important tools to study the changes in the ocean environment both today and in the past. However, the behaviors of TEIs in marine sediments are poorly known, limiting our ability to use them in oceanography. Here we present a modeling framework that can be used to generate and run models of the sedimentary cycling of TEIs assisted with advanced numerical tools in the Julia language, lowering the coding barrier for the general user to study marine TEIs.
Siyu Zhu, Peipei Wu, Siyi Zhang, Oliver Jahn, Shu Li, and Yanxu Zhang
Geosci. Model Dev., 16, 5915–5929, https://doi.org/10.5194/gmd-16-5915-2023, https://doi.org/10.5194/gmd-16-5915-2023, 2023
Short summary
Short summary
In this study, we estimate the global biogeochemical cycling of Hg in a state-of-the-art physical-ecosystem ocean model (high-resolution-MITgcm/Hg), providing a more accurate portrayal of surface Hg concentrations in estuarine and coastal areas, strong western boundary flow and upwelling areas, and concentration diffusion as vortex shapes. The high-resolution model can help us better predict the transport and fate of Hg in the ocean and its impact on the global Hg cycle.
Maria Val Martin, Elena Blanc-Betes, Ka Ming Fung, Euripides P. Kantzas, Ilsa B. Kantola, Isabella Chiaravalloti, Lyla L. Taylor, Louisa K. Emmons, William R. Wieder, Noah J. Planavsky, Michael D. Masters, Evan H. DeLucia, Amos P. K. Tai, and David J. Beerling
Geosci. Model Dev., 16, 5783–5801, https://doi.org/10.5194/gmd-16-5783-2023, https://doi.org/10.5194/gmd-16-5783-2023, 2023
Short summary
Short summary
Enhanced rock weathering (ERW) is a CO2 removal strategy that involves applying crushed rocks (e.g., basalt) to agricultural soils. However, unintended processes within the N cycle due to soil pH changes may affect the climate benefits of C sequestration. ERW could drive changes in soil emissions of non-CO2 GHGs (N2O) and trace gases (NO and NH3) that may affect air quality. We present a new improved N cycling scheme for the land model (CLM5) to evaluate ERW effects on soil gas N emissions.
Huajie Zhu, Mousong Wu, Fei Jiang, Michael Vossbeck, Thomas Kaminski, Xiuli Xing, Jun Wang, Weimin Ju, and Jing M. Chen
EGUsphere, https://doi.org/10.5194/egusphere-2023-1955, https://doi.org/10.5194/egusphere-2023-1955, 2023
Short summary
Short summary
In this work, Nanjing University Carbon Assimilation System (NUCAS) was developed. Data assimilation experiments were conducted to demonstrate the robustness and to investigate the feasibility and applicability of NUCAS. The assimilation of COS significantly improved the model performance in gross primary productivity, sensible heat, latent heat and even soil moisture, demonstrating that COS can provide strong constraints on parameters relevant to water, energy and carbon processes.
Özgür Gürses, Laurent Oziel, Onur Karakuş, Dmitry Sidorenko, Christoph Völker, Ying Ye, Moritz Zeising, Martin Butzin, and Judith Hauck
Geosci. Model Dev., 16, 4883–4936, https://doi.org/10.5194/gmd-16-4883-2023, https://doi.org/10.5194/gmd-16-4883-2023, 2023
Short summary
Short summary
This paper assesses the biogeochemical model REcoM3 coupled to the ocean–sea ice model FESOM2.1. The model can be used to simulate the carbon uptake or release of the ocean on timescales of several hundred years. A detailed analysis of the nutrients, ocean productivity, and ecosystem is followed by the carbon cycle. The main conclusion is that the model performs well when simulating the observed mean biogeochemical state and variability and is comparable to other ocean–biogeochemical models.
Hocheol Seo and Yeonjoo Kim
Geosci. Model Dev., 16, 4699–4713, https://doi.org/10.5194/gmd-16-4699-2023, https://doi.org/10.5194/gmd-16-4699-2023, 2023
Short summary
Short summary
Wildfire is a crucial factor in carbon and water fluxes on the Earth system. About 2.1 Pg of carbon is released into the atmosphere by wildfires annually. Because the fire processes are still limitedly represented in land surface models, we forced the daily GFED4 burned area into the land surface model over Alaska and Siberia. The results with the GFED4 burned area significantly improved the simulated carbon emissions and net ecosystem exchange compared to the default simulation.
Boris Ťupek, Aleksi Lehtonen, Alla Yurova, Rose Abramoff, Stefano Manzoni, Bertrand Guenet, Elisa Bruni, Samuli Launiainen, Mikko Peltoniemi, Shoji Hashimoto, Xianglin Tian, Juha Heikkinen, Kari Minkkinen, and Raisa Mäkipää
EGUsphere, https://doi.org/10.5194/egusphere-2023-1523, https://doi.org/10.5194/egusphere-2023-1523, 2023
Short summary
Short summary
Improving representation of peatlands in Earth system models (ESMs) (e.g., in Yasso soil C model coupled with JSBACH model) is rare. We updated Yasso07 moisture dependency for applications on both the mineral soils and peatlands of boreal forest. Our moisture function improved the simulated peatland C stocks and CO2 dynamics. However, due to nature of the measurements used for Bayesian MCMC data assimilation, it contrasts with functions in other ESMs (see graphical abstract).
Hideki Ninomiya, Tomomichi Kato, Lea Végh, and Lan Wu
Geosci. Model Dev., 16, 4155–4170, https://doi.org/10.5194/gmd-16-4155-2023, https://doi.org/10.5194/gmd-16-4155-2023, 2023
Short summary
Short summary
Non-structural carbohydrates (NSCs) play a crucial role in plants to counteract the effects of climate change. We added a new NSC module into the SEIB-DGVM, an individual-based ecosystem model. The simulated NSC levels and their seasonal patterns show a strong agreement with observed NSC data at both point and global scales. The model can be used to simulate the biotic effects resulting from insufficient NSCs, which are otherwise difficult to measure in terrestrial ecosystems globally.
Miquel De Cáceres, Roberto Molowny-Horas, Antoine Cabon, Jordi Martínez-Vilalta, Maurizio Mencuccini, Raúl García-Valdés, Daniel Nadal-Sala, Santiago Sabaté, Nicolas Martin-StPaul, Xavier Morin, Francesco D'Adamo, Enric Batllori, and Aitor Améztegui
Geosci. Model Dev., 16, 3165–3201, https://doi.org/10.5194/gmd-16-3165-2023, https://doi.org/10.5194/gmd-16-3165-2023, 2023
Short summary
Short summary
Regional-level applications of dynamic vegetation models are challenging because they need to accommodate the variation in plant functional diversity. This can be done by estimating parameters from available plant trait databases while adopting alternative solutions for missing data. Here we present the design, parameterization and evaluation of MEDFATE (version 2.9.3), a novel model of forest dynamics for its application over a region in the western Mediterranean Basin.
Jens Heinke, Susanne Rolinski, and Christoph Müller
Geosci. Model Dev., 16, 2455–2475, https://doi.org/10.5194/gmd-16-2455-2023, https://doi.org/10.5194/gmd-16-2455-2023, 2023
Short summary
Short summary
We develop a livestock module for the global vegetation model LPJmL5.0 to simulate the impact of grazing dairy cattle on carbon and nitrogen cycles in grasslands. A novelty of the approach is that it accounts for the effect of feed quality on feed uptake and feed utilization by animals. The portioning of dietary nitrogen into milk, feces, and urine shows very good agreement with estimates obtained from animal trials.
Yimian Ma, Xu Yue, Stephen Sitch, Nadine Unger, Johan Uddling, Lina M. Mercado, Cheng Gong, Zhaozhong Feng, Huiyi Yang, Hao Zhou, Chenguang Tian, Yang Cao, Yadong Lei, Alexander W. Cheesman, Yansen Xu, and Maria Carolina Duran Rojas
Geosci. Model Dev., 16, 2261–2276, https://doi.org/10.5194/gmd-16-2261-2023, https://doi.org/10.5194/gmd-16-2261-2023, 2023
Short summary
Short summary
Plants have been found to respond differently to O3, but the variations in the sensitivities have rarely been explained nor fully implemented in large-scale assessment. This study proposes a new O3 damage scheme with leaf mass per area to unify varied sensitivities for all plant species. Our assessment reveals an O3-induced reduction of 4.8 % in global GPP, with the highest reduction of >10 % for cropland, suggesting an emerging risk of crop yield loss under the threat of O3 pollution.
Winslow D. Hansen, Adrianna Foster, Benjamin Gaglioti, Rupert Seidl, and Werner Rammer
Geosci. Model Dev., 16, 2011–2036, https://doi.org/10.5194/gmd-16-2011-2023, https://doi.org/10.5194/gmd-16-2011-2023, 2023
Short summary
Short summary
Permafrost and the thick soil-surface organic layers that insulate permafrost are important controls of boreal forest dynamics and carbon cycling. However, both are rarely included in process-based vegetation models used to simulate future ecosystem trajectories. To address this challenge, we developed a computationally efficient permafrost and soil organic layer module that operates at fine spatial (1 ha) and temporal (daily) resolutions.
Heewon Jung, Hyun-Seob Song, and Christof Meile
Geosci. Model Dev., 16, 1683–1696, https://doi.org/10.5194/gmd-16-1683-2023, https://doi.org/10.5194/gmd-16-1683-2023, 2023
Short summary
Short summary
Microbial activity responsible for many chemical transformations depends on environmental conditions. These can vary locally, e.g., between poorly connected pores in porous media. We present a modeling framework that resolves such small spatial scales explicitly, accounts for feedback between transport and biogeochemical conditions, and can integrate state-of-the-art representations of microbes in a computationally efficient way, making it broadly applicable in science and engineering use cases.
Arthur Guignabert, Quentin Ponette, Frédéric André, Christian Messier, Philippe Nolet, and Mathieu Jonard
Geosci. Model Dev., 16, 1661–1682, https://doi.org/10.5194/gmd-16-1661-2023, https://doi.org/10.5194/gmd-16-1661-2023, 2023
Short summary
Short summary
Spatially explicit and process-based models are useful to test innovative forestry practices under changing and uncertain conditions. However, their larger use is often limited by the restricted range of species and stand structures they can reliably account for. We therefore calibrated and evaluated such a model, HETEROFOR, for 23 species across southern Québec. Our results showed that the model is robust and can predict accurately both individual tree growth and stand dynamics in this region.
Maureen Beaudor, Nicolas Vuichard, Juliette Lathière, Nikolaos Evangeliou, Martin Van Damme, Lieven Clarisse, and Didier Hauglustaine
Geosci. Model Dev., 16, 1053–1081, https://doi.org/10.5194/gmd-16-1053-2023, https://doi.org/10.5194/gmd-16-1053-2023, 2023
Short summary
Short summary
Ammonia mainly comes from the agricultural sector, and its volatilization relies on environmental variables. Our approach aims at benefiting from an Earth system model framework to estimate it. By doing so, we represent a consistent spatial distribution of the emissions' response to environmental changes.
We greatly improved the seasonal cycle of emissions compared with previous work. In addition, our model includes natural soil emissions (that are rarely represented in modeling approaches).
Rui Ying, Fanny M. Monteiro, Jamie D. Wilson, and Daniela N. Schmidt
Geosci. Model Dev., 16, 813–832, https://doi.org/10.5194/gmd-16-813-2023, https://doi.org/10.5194/gmd-16-813-2023, 2023
Short summary
Short summary
Planktic foraminifera are marine-calcifying zooplankton; their shells are widely used to measure past temperature and productivity. We developed ForamEcoGEnIE 2.0 to simulate the four subgroups of this organism. We found that the relative abundance distribution agrees with marine sediment core-top data and that carbon export and biomass are close to sediment trap and plankton net observations respectively. This model provides the opportunity to study foraminiferal ecology in any geological era.
Onur Kerimoglu, Markus Pahlow, Prima Anugerahanti, and Sherwood Lan Smith
Geosci. Model Dev., 16, 95–108, https://doi.org/10.5194/gmd-16-95-2023, https://doi.org/10.5194/gmd-16-95-2023, 2023
Short summary
Short summary
In classical models that track the changes in the elemental composition of phytoplankton, additional state variables are required for each element resolved. In this study, we show how the behavior of such an explicit model can be approximated using an
instantaneous acclimationapproach, in which the elemental composition of the phytoplankton is assumed to adjust to an optimal value instantaneously. Through rigorous tests, we evaluate the consistency of this scheme.
Yuan Zhang, Devaraju Narayanappa, Philippe Ciais, Wei Li, Daniel Goll, Nicolas Vuichard, Martin G. De Kauwe, Laurent Li, and Fabienne Maignan
Geosci. Model Dev., 15, 9111–9125, https://doi.org/10.5194/gmd-15-9111-2022, https://doi.org/10.5194/gmd-15-9111-2022, 2022
Short summary
Short summary
There are a few studies to examine if current models correctly represented the complex processes of transpiration. Here, we use a coefficient Ω, which indicates if transpiration is mainly controlled by vegetation processes or by turbulence, to evaluate the ORCHIDEE model. We found a good performance of ORCHIDEE, but due to compensation of biases in different processes, we also identified how different factors control Ω and where the model is wrong. Our method is generic to evaluate other models.
Thomas Neumann, Hagen Radtke, Bronwyn Cahill, Martin Schmidt, and Gregor Rehder
Geosci. Model Dev., 15, 8473–8540, https://doi.org/10.5194/gmd-15-8473-2022, https://doi.org/10.5194/gmd-15-8473-2022, 2022
Short summary
Short summary
Marine ecosystem models are usually constrained by the elements nitrogen and phosphorus and consider carbon in organic matter in a fixed ratio. Recent observations show a substantial deviation from the simulated carbon cycle variables. In this study, we present a marine ecosystem model for the Baltic Sea which allows for a flexible uptake ratio for carbon, nitrogen, and phosphorus. With this extension, the model reflects much more reasonable variables of the marine carbon cycle.
Arsène Druel, Simon Munier, Anthony Mucia, Clément Albergel, and Jean-Christophe Calvet
Geosci. Model Dev., 15, 8453–8471, https://doi.org/10.5194/gmd-15-8453-2022, https://doi.org/10.5194/gmd-15-8453-2022, 2022
Short summary
Short summary
Crop phenology and irrigation is implemented into a land surface model able to work at a global scale. A case study is presented over Nebraska (USA). Simulations with and without the new scheme are compared to different satellite-based observations. The model is able to produce a realistic yearly irrigation water amount. The irrigation scheme improves the simulated leaf area index, gross primary productivity, evapotransipiration, and land surface temperature.
Thomas Wutzler, Lin Yu, Marion Schrumpf, and Sönke Zaehle
Geosci. Model Dev., 15, 8377–8393, https://doi.org/10.5194/gmd-15-8377-2022, https://doi.org/10.5194/gmd-15-8377-2022, 2022
Short summary
Short summary
Soil microbes process soil organic matter and affect carbon storage and plant nutrition at the ecosystem scale. We hypothesized that decadal dynamics is constrained by the ratios of elements in litter inputs, microbes, and matter and that microbial community optimizes growth. This allowed the SESAM model to descibe decadal-term carbon sequestration in soils and other biogeochemical processes explicitly accounting for microbial processes but without its problematic fine-scale parameterization.
Cited articles
Balay, S., Gropp, W. D., McInnes, L. C., and Smith, B. F.: Efficient Management of Parallelism in Object Oriented Numerical Software Libraries, in: Modern Software Tools in Scientific Computing, edited by: Arge, E., Bruaset, A. M., and Langtangen, H. P., Birkhäuser Press, Basel, 163–202, 1997.
Balay, S., Brown, J., Buschelman, K., Eijkhout, V., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and Zhang, H.: PETSc Users Manual, Tech. Rep. ANL-95/11 – Revision 3.3, Argonne National Laboratory, Lemont, 2012a.
Balay, S., Buschelman, K., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and Zhang, H.: PETSc Web page, available at: http://www.mcs.anl.gov/petsc/ (last access: 12 July 2013), 2012b.
Bernsen, E., Dijkstra, H. A., and Wubs, F. W.: A method to reduce the spin-up time of ocean models, Ocean Modell., 20, 380–392, https://doi.org/10.1016/j.ocemod.2007.10.008, 2008.
Bryan, K.: Accelerating the Convergence to Equilibrium of Ocean-Climate Models, J. Phys. Oceanogr., 14, 666–673, https://doi.org/10.1175/1520-0485(1984)014<0666:ATCTEO>2.0.CO;2, 1984.
Ciric, L. B.: A Generalization of Banach's Contraction Principle, P. Am. Math. Soc., 45, 267–273, 1974.
Danabasoglu, G., McWilliams, J. C., and Large, W. G.: Approach to Equilibrium in Accelerated Global Oceanic Models, J. Climate, 9, 1092–1110, https://doi.org/10.1175/1520-0442(1996)009<1092:ATEIAG>2.0.CO;2, 1996.
Dennis, J. and Schnabel, R.: Numerical methods for unconstrained optimization and nonlinear equations, Society for Industrial and Applied Mathematics, Philadelphia, 1996.
Dutkiewicz, S., Sokolov, A. P., Scott, J., and Stone, P. H.: A three-dimensional ocean-seaice-carbon cycle model and its coupling to a two-dimensional atmospheric model: Uses in climate change studies, Tech. Rep. 122, MIT Joint Program on the Science and Policy of Global Change, Cambridge, 2005.
Eisenstat, S. C. and Walker, H. F.: Choosing the Forcing Terms in an Inexact Newton Method, SIAM J. Sci. Comput., 17, 16–32, https://doi.org/10.1137/0917003, 1996.
Evans, G. T. and Garçon, V. C.: One-Dimensional Models of Water Column Biogeochemistry, Report of a workshop held in Toulouse, France, November–December 1995, GOFS Report No.23/97, JGOFS Bergen, Norway, 1997.
Fasham, M. J. R. (Ed.): Ocean Biogeochemistry. The Role of the Ocean Carbon Cycle in Global Change, Global Change – The IGBP Series, Springer, Berlin, 2003.
Fennel, K., Losch, M., Schröter, J., and Wenzel, M.: Testing a marine ecosystem model: sensitivity analysis and parameter optimization, J. Marine Syst., 28, 45–63, https://doi.org/10.1016/S0924-7963(00)00083-X, 2001.
Griewank, A. and Walther, A.: Evaluating derivatives: principles and techniques of algorithmic differentiation, Society for Industrial and Applied Mathematics (SIAM), https://doi.org/10.1137/1.9780898717761, 2008.
Kelley, C. T.: Solving nonlinear equations with Newton's method, SIAM, Philadelphia, https://doi.org/10.1137/1.9780898718898, 2003.
Khatiwala, S.: A computational framework for simulation of biogeochemical tracers in the ocean, Global Biogeochem. Cy., 21, https://doi.org/10.1029/2007GB002923, 2007.
Khatiwala, S.: Fast spin up of Ocean biogeochemical models using matrix-free Newton-Krylov, Ocean Modell., 23, 121–129, https://doi.org/10.1016/j.ocemod.2008.05.002, 2008.
Khatiwala, S.: TMM framework web page, available at: http://www.ldeo.columbia.edu/%7Espk/Research/TMM/ (last access: 14 June 2016), 2013.
Khatiwala, S., Visbeck, M., and Cane, M.: Accelerated simulation of passive tracers in ocean circulation models, Ocean Modell., 9, 51–69, 2005.
Knoll, D. and Keyes, D.: Jacobian-free Newton–Krylov methods: a survey of approaches and applications, J. Comput. Phys., 193, 357–397, 2004.
Kriest, I., Khatiwala, S., and Oschlies, A.: Towards an assessment of simple global marine biogeochemical models of different complexity, Prog. Oceanogr., 86, 337–360, https://doi.org/10.1016/j.pocean.2010.05.002, 2010.
Marshall, J., Adcroft, A., Hill, C., Perelman, L., and Heisey, C.: A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers, J. Geophys. Res., 102, 5753–5766, 1997.
McKay, M. D., Beckman, R. J., and Conover, W. J.: Comparison of three methods for selecting values of input variables in the analysis of output from a computer code, Technometrics, 21, 239–245, 1979.
Merlis, T. M. and Khatiwala, S.: Fast dynamical spin-up of ocean general circulation models using Newton–Krylov methods, Ocean Modell., 21, 97–105, 2008.
Paltridge, G. W. and Platt, C. M. R.: Radiative Processes in Meteorology and Climatology, Elsevier, New York, https://doi.org/10.1002/qj.49710343713, 1976.
Saad, Y.: Iterative Methods for Sparse Linear Systems, The Society for Industrial and Applied Mathematics (SIAM), 2003.
Saad, Y. and Schultz, M.: GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems, SIAM J. Sci. Stat. Comput., 7, 856–869, https://doi.org/10.1137/0907058, 1986.
Sarmiento, J. L. and Gruber, N.: Ocean Biogeochemical Dynamics, Princeton University Press, 2006.
Schartau, M. and Oschlies, A.: Simultaneous data-based optimization of a 1d-ecosystem model at three locations in the north Atlantic: Part I – method and parameter estimates, J. Marine Res., 61, 765–793, 2003.
Siewertsen, E., Piwonski, J., and Slawig, T.: Porting marine ecosystem model spin-up using transport matrices to GPUs, Geosci. Model Dev., 6, 17–28, https://doi.org/10.5194/gmd-6-17-2013, 2013.
Stoer, J. and Bulirsch, R.: Introduction to Numerical Analysis, Springer, New York, 3rd Edn., 2002.
Walker, D. W. and Dongarra, J. J.: MPI: A Standard Message Passing Interface, Supercomputer, 12, 56–68, 1996.
Wang, D.: A note on using the accelerated convergence method in climate models, Tellus A, 53, 27–34, https://doi.org/10.1034/j.1600-0870.2001.01134.x, 2001.
Zumbusch, G. W.: Dynamic Load Balancing in a Lightweight Adaptive Parallel Multigrid PDE Solver, in: PPSC, SIAM, Philadelphia, 1999.
Short summary
In order to fundamentally tackle the problem of parameter identification for marine ecosystem models in 3-D, we introduced a general biogeochemical programming interface that fits into the optimization context. Moreover, we implemented a comprehensive parallel solver software for periodic steady states that uses the interface to couple marine ecosystem models to a transport matrix driver. We validated the new implementation using a hierarchy of biogeochemical models.
In order to fundamentally tackle the problem of parameter identification for marine ecosystem...
