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
Abstract
Reactive transport codes (e.g., PFLOTRAN) are increasingly used to improve the representation of biogeochemical processes in terrestrial ecosystem models (e.g., the Community Land Model, CLM). As CLM and PFLOTRAN use explicit and implicit time stepping, implementation of CLM biogeochemical reactions in PFLOTRAN can result in negative concentration, which is not physical and can cause numerical instability and errors. The objective of this work is to address the nonnegativity challenge to obtain accurate, efficient, and robust solutions. We illustrate the implementation of a reaction network with the CLM-CN decomposition, nitrification, denitrification, and plant nitrogen uptake reactions and test the implementation at arctic, temperate, and tropical sites. We examine use of scaling back the update during each iteration (SU), log transformation (LT), and downregulating the reaction rate to account for reactant availability limitation to enforce nonnegativity. Both SU and LT guarantee nonnegativity but with implications. When a very small scaling factor occurs due to either consumption or numerical overshoot, and the iterations are deemed converged because of too small an update, SU can introduce excessive numerical error. LT involves multiplication of the Jacobian matrix by the concentration vector, which increases the condition number, decreases the time step size, and increases the computationalmore »
- Authors:
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States); Sandia National Lab. (SNL-CA), Livermore, CA (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC); Work for Others (WFO)
- OSTI Identifier:
- 1240298
- Alternate Identifier(s):
- OSTI ID: 1240537; OSTI ID: 1261540; OSTI ID: 1266840; OSTI ID: 1327768
- Grant/Contract Number:
- AC05-00OR22725; NGEE-Arctic; AC04-94AL85000
- Resource Type:
- Published Article
- Journal Name:
- Geoscientific Model Development (Online)
- Additional Journal Information:
- Journal Name: Geoscientific Model Development (Online) Journal Volume: 9 Journal Issue: 3; Journal ID: ISSN 1991-9603
- Publisher:
- Copernicus Publications, EGU
- Country of Publication:
- Germany
- Language:
- English
- Subject:
- 58 GEOSCIENCES
Citation Formats
Tang, Guoping, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn E., Lichtner, Peter C., Kumar, Jitendra, Mills, Richard T., Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest M., Painter, Scott L., and Thornton, Peter 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. Germany: N. p., 2016.
Web. doi:10.5194/gmd-9-927-2016.
Tang, Guoping, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn E., Lichtner, Peter C., Kumar, Jitendra, Mills, Richard T., Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest M., Painter, Scott L., & Thornton, Peter 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. Germany. https://doi.org/10.5194/gmd-9-927-2016
Tang, Guoping, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn E., Lichtner, Peter C., Kumar, Jitendra, Mills, Richard T., Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest M., Painter, Scott L., and Thornton, Peter E. Fri .
"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". Germany. https://doi.org/10.5194/gmd-9-927-2016.
@article{osti_1240298,
title = {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},
author = {Tang, Guoping and Yuan, Fengming and Bisht, Gautam and Hammond, Glenn E. and Lichtner, Peter C. and Kumar, Jitendra and Mills, Richard T. and Xu, Xiaofeng and Andre, Ben and Hoffman, Forrest M. and Painter, Scott L. and Thornton, Peter E.},
abstractNote = {Reactive transport codes (e.g., PFLOTRAN) are increasingly used to improve the representation of biogeochemical processes in terrestrial ecosystem models (e.g., the Community Land Model, CLM). As CLM and PFLOTRAN use explicit and implicit time stepping, implementation of CLM biogeochemical reactions in PFLOTRAN can result in negative concentration, which is not physical and can cause numerical instability and errors. The objective of this work is to address the nonnegativity challenge to obtain accurate, efficient, and robust solutions. We illustrate the implementation of a reaction network with the CLM-CN decomposition, nitrification, denitrification, and plant nitrogen uptake reactions and test the implementation at arctic, temperate, and tropical sites. We examine use of scaling back the update during each iteration (SU), log transformation (LT), and downregulating the reaction rate to account for reactant availability limitation to enforce nonnegativity. Both SU and LT guarantee nonnegativity but with implications. When a very small scaling factor occurs due to either consumption or numerical overshoot, and the iterations are deemed converged because of too small an update, SU can introduce excessive numerical error. LT involves multiplication of the Jacobian matrix by the concentration vector, which increases the condition number, decreases the time step size, and increases the computational cost. Neither SU nor SE prevents zero concentration. When the concentration is close to machine precision or 0, a small positive update stops all reactions for SU, and LT can fail due to a singular Jacobian matrix. The consumption rate has to be downregulated such that the solution to the mathematical representation is positive. A first-order rate downregulates consumption and is nonnegative, and adding a residual concentration makes it positive. For zero-order rate or when the reaction rate is not a function of a reactant, representing the availability limitation of each reactant with a Monod substrate limiting function provides a smooth transition between a zero-order rate when the reactant is abundant and first-order rate when the reactant becomes limiting. When the half saturation is small, marching through the transition may require small time step sizes to resolve the sharp change within a small range of concentration values. Our results from simple tests and CLM-PFLOTRAN simulations caution against use of SU and indicate that accurate, stable, and relatively efficient solutions can be achieved with LT and downregulation with Monod substrate limiting function and residual concentration.},
doi = {10.5194/gmd-9-927-2016},
journal = {Geoscientific Model Development (Online)},
number = 3,
volume = 9,
place = {Germany},
year = {Fri Mar 04 00:00:00 EST 2016},
month = {Fri Mar 04 00:00:00 EST 2016}
}
https://doi.org/10.5194/gmd-9-927-2016
Web of Science
Works referenced in this record:
Investigating soil moisture–climate interactions in a changing climate: A review
journal, May 2010
- Seneviratne, Sonia I.; Corti, Thierry; Davin, Edouard L.
- Earth-Science Reviews, Vol. 99, Issue 3-4
Long residence times of rapidly decomposable soil organic matter: application of a multi-phase, multi-component, and vertically resolved model (BAMS1) to soil carbon dynamics
journal, January 2014
- Riley, W. J.; Maggi, F.; Kleber, M.
- Geoscientific Model Development, Vol. 7, Issue 4
Generalized model for NO x and N 2 O emissions from soils
journal, August 2001
- Parton, W. J.; Holland, E. A.; Del Grosso, S. J.
- Journal of Geophysical Research: Atmospheres, Vol. 106, Issue D15
Reactive transport codes for subsurface environmental simulation
journal, September 2014
- Steefel, C. I.; Appelo, C. A. J.; Arora, B.
- Computational Geosciences, Vol. 19, Issue 3
Preference for different inorganic nitrogen forms among plant functional types and species of the Patagonian steppe
journal, June 2013
- Gherardi, Laureano A.; Sala, Osvaldo E.; Yahdjian, Laura
- Oecologia, Vol. 173, Issue 3
Uptake of nitrate, ammonium and glycine by plants of Tasmanian wet eucalypt forests
journal, March 2007
- Warren, C. R.; Adams, P. R.
- Tree Physiology, Vol. 27, Issue 3
Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN: Evaluating the Parallel Performance of Pflotran
journal, January 2014
- Hammond, G. E.; Lichtner, P. C.; Mills, R. T.
- Water Resources Research, Vol. 50, Issue 1
Competition between roots and microorganisms for nitrogen: mechanisms and ecological relevance
journal, March 2013
- Kuzyakov, Yakov; Xu, Xingliang
- New Phytologist, Vol. 198, Issue 3
Growth and mineral nutrition of non-mycorrhizal and mycorrhizal Norway spruce (Picea abies) seedlings grown in semi-hydroponic sand culture. I. Growth and mineral nutrient uptake in plants supplied with different forms of nitrogen
journal, July 1996
- Eltrop, Ludger; Marschner, Horst
- New Phytologist, Vol. 133, Issue 3
PFLOTRAN User Manual: A Massively Parallel Reactive Flow and Transport Model for Describing Surface and Subsurface Processes
report, January 2015
- Lichtner, Peter; Hammond, Glenn; Lu, Chuan
CLM4-BeTR, a generic biogeochemical transport and reaction module for CLM4: model development, evaluation, and application
journal, January 2013
- Tang, J. Y.; Riley, W. J.; Koven, C. D.
- Geoscientific Model Development, Vol. 6, Issue 1
Effect of 7 yr of experimental drought on vegetation dynamics and biomass storage of an eastern Amazonian rainforest
journal, June 2010
- da Costa, Antonio Carlos Lola; Galbraith, David; Almeida, Samuel
- New Phytologist, Vol. 187, Issue 3
Application of a coupled ecosystem-chemical equilibrium model, DayCent-Chem, to stream and soil chemistry in a Rocky Mountain watershed
journal, January 2007
- Hartman, Melannie D.; Baron, Jill S.; Ojima, Dennis S.
- Ecological Modelling, Vol. 200, Issue 3-4
Soil carbon and nitrogen mineralization: Theory and models across scales
journal, July 2009
- Manzoni, Stefano; Porporato, Amilcare
- Soil Biology and Biochemistry, Vol. 41, Issue 7
Uptake of ammonium and nitrate by heterotrophic bacteria and phytoplankton in the sub-Arctic Pacific
journal, February 1998
- Kirchman, David L.; Wheeler, Patricia A.
- Deep Sea Research Part I: Oceanographic Research Papers, Vol. 45, Issue 2-3
Extreme Oxygen Sensitivity in Methanogenic Archaebacteria
journal, May 1985
- Jarrell, Ken F.
- BioScience, Vol. 35, Issue 5
Technical Note: A generic law-of-the-minimum flux limiter for simulating substrate limitation in biogeochemical models
journal, January 2016
- Tang, J. Y.; Riley, W. J.
- Biogeosciences, Vol. 13, Issue 3
Nitrogen uptake by heterotrophic bacteria and phytoplankton in the nitrate-rich Thames estuary
journal, January 2000
- Middelburg, Jj; Nieuwenhuize, J.
- Marine Ecology Progress Series, Vol. 203
The Rumen Microbial Ecosystem
journal, November 1975
- Hungate, R. E.
- Annual Review of Ecology and Systematics, Vol. 6, Issue 1
Confronting model predictions of carbon fluxes with measurements of Amazon forests subjected to experimental drought
journal, July 2013
- Powell, Thomas L.; Galbraith, David R.; Christoffersen, Bradley O.
- New Phytologist, Vol. 200, Issue 2
Estimated change in tundra ecosystem function near Barrow, Alaska between 1972 and 2010
journal, January 2012
- Lara, M. J.; Villarreal, S.; Johnson, D. R.
- Environmental Research Letters, Vol. 7, Issue 1
Ecosystem model spin-up: Estimating steady state conditions in a coupled terrestrial carbon and nitrogen cycle model
journal, November 2005
- Thornton, Peter E.; Rosenbloom, Nan A.
- Ecological Modelling, Vol. 189, Issue 1-2
Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4)
journal, October 2012
- Bonan, Gordon B.; Hartman, Melannie D.; Parton, William J.
- Global Change Biology, Vol. 19, Issue 3
Kinetics of nitrate and ammonium absorption and accompanying H+ fluxes in roots of Lolium perenne L. and N2-fixing Trifolium repens L.
journal, September 1997
- Hogh-Jensen, H.; Wollenweber, B.; Schjoerring, J. K.
- Plant, Cell and Environment, Vol. 20, Issue 9
Nitrogen uptake by Eucalyptus regnans and Acacia spp. - preferences, resource overlap and energetic costs
journal, January 2009
- Pfautsch, S.; Rennenberg, H.; Bell, T. L.
- Tree Physiology, Vol. 29, Issue 3
Spatial and temporal patterns of active layer thickness at Circumpolar Active Layer Monitoring (CALM) sites in northern Alaska, 1995–2000
journal, January 2003
- Hinkel, K. M.
- Journal of Geophysical Research, Vol. 108, Issue D2
Oak Forest Carbon and Water Simulations: Model Intercomparisons and Evaluations Against Independent data
journal, August 2004
- Hanson, P. J.; Amthor, J. S.; Wullschleger, S. D.
- Ecological Monographs, Vol. 74, Issue 3
Generalized model for N 2 and N 2 O production from nitrification and denitrification
journal, September 1996
- Parton, W. J.; Mosier, A. R.; Ojima, D. S.
- Global Biogeochemical Cycles, Vol. 10, Issue 3
A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components
journal, January 1989
- Yeh, G. T.; Tripathi, V. S.
- Water Resources Research, Vol. 25, Issue 1
Combined effects of short term rainfall patterns and soil texture on soil nitrogen cycling — A modeling analysis
journal, March 2010
- Gu, Chuanhui; Riley, William J.
- Journal of Contaminant Hydrology, Vol. 112, Issue 1-4
Persistence of soil organic matter as an ecosystem property
journal, October 2011
- Schmidt, Michael W. I.; Torn, Margaret S.; Abiven, Samuel
- Nature, Vol. 478, Issue 7367
Modeling the Exchanges of Energy, Water, and Carbon Between Continents and the Atmosphere
journal, January 1997
- Sellers, P. J.
- Science, Vol. 275, Issue 5299
Climate Change 2013 – The Physical Science Basis
book, March 2014
- Intergovernmental Panel on Climate Change,
Modeling the Production of and Competition for Hydrogen in a Dechlorinating Culture
journal, August 1998
- Fennell, Donna E.; Gossett, James M.
- Environmental Science & Technology, Vol. 32, Issue 16
Nitrogen uptake by heterotrophic bacteria and phytoplankton in Arctic surface waters
journal, January 2007
- Fouilland, E.; Gosselin, M.; Rivkin, R. B.
- Journal of Plankton Research, Vol. 29, Issue 4
A generic biogeochemical module for Earth system models: Next Generation BioGeoChemical Module (NGBGC), version 1.0
journal, January 2013
- Fang, Y.; Huang, M.; Liu, C.
- Geoscientific Model Development, Vol. 6, Issue 6
Soil nitrogen form and plant nitrogen uptake along a boreal forest productivity gradient
journal, September 2001
- Nordin, Annika; Högberg, Peter; Näsholm, Torgny
- Oecologia, Vol. 129, Issue 1
Modelling changes in nitrogen cycling to sustain increases in forest productivity under elevated atmospheric CO 2 and contrasting site conditions
journal, January 2013
- Grant, R. F.
- Biogeosciences, Vol. 10, Issue 11
The uptake of inorganic nutrients by heterotrophic bacteria
journal, September 1994
- Kirchman, D. L.
- Microbial Ecology, Vol. 28, Issue 2
Interspecies differences in the preference of ammonium and nitrate in vascular plants
journal, January 1995
- Falkengren-Grerup, Ursula
- Oecologia, Vol. 102, Issue 3
A mechanistic treatment of the dominant soil nitrogen cycling processes: Model development, testing, and application: NITROGEN CYCLE MODELING
journal, April 2008
- Maggi, F.; Gu, C.; Riley, W. J.
- Journal of Geophysical Research: Biogeosciences, Vol. 113, Issue G2
Non-negative solutions of ODEs
journal, November 2005
- Shampine, L. F.; Thompson, S.; Kierzenka, J. A.
- Applied Mathematics and Computation, Vol. 170, Issue 1
Microbial uptake of dissolved organic and inorganic nitrogen in Randers Fjord
journal, November 2004
- Veuger, Bart; Middelburg, Jack J.; Boschker, Henricus T. S.
- Estuarine, Coastal and Shelf Science, Vol. 61, Issue 3