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: Modeling Archive

Tang, Guoping; Mills, Richard; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest; Painter, Scott; Thornton, Peter; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn; Lichtner, Peter; Kumar, Jitendra

doi:10.5440/1239799

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: Modeling Archive

Dataset
Other Related Research

Abstract

This Modeling Archive is in support of an NGEE Arctic discussion paper under review and available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model.

Authors:: Tang, Guoping; Mills, Richard; Xu, Xiaofeng; Andre, Ben; Hoffman, Forrest; Painter, Scott; Thornton, Peter; Yuan, Fengming; Bisht, Gautam; Hammond, Glenn; Lichtner, Peter; Kumar, Jitendra

Publication Date:: 2016-03-04

Other Number(s):: NGA076
ngee_3D38B6A9C3DACF8A340B5D4493C65B642018_06_21_111631201

DOE Contract Number:: DE-AC05-00OR22725

Research Org.:: Next Generation Ecosystems Experiment - Arctic, Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (US); NGEE Arctic, Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)

Sponsoring Org.:: USDOE Office of Science (SC), Biological and Environmental Research (BER)

Collaborations:: ORNL

Subject:: 54 Environmental Sciences

Keywords:: Nonnegativity; implicit time stepping; reactive transport model; soil biogeochemistry

OSTI Identifier:: 1239799

DOI:: https://doi.org/10.5440/1239799

Citation Formats


                    Tang, Guoping, Mills, Richard, Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest, Painter, Scott, Thornton, Peter, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn, Lichtner, Peter, and Kumar, Jitendra. 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: Modeling Archive.  United States: N. p., 2016. 
        Web.  doi:10.5440/1239799.

Copy to clipboard


                    Tang, Guoping, Mills, Richard, Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest, Painter, Scott, Thornton, Peter, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn, Lichtner, Peter, & Kumar, Jitendra. 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: Modeling Archive.  United States.  doi:https://doi.org/10.5440/1239799

Copy to clipboard


                    Tang, Guoping, Mills, Richard, Xu, Xiaofeng, Andre, Ben, Hoffman, Forrest, Painter, Scott, Thornton, Peter, Yuan, Fengming, Bisht, Gautam, Hammond, Glenn, Lichtner, Peter, and Kumar, Jitendra. 2016.  
"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: Modeling Archive".  United States.  doi:https://doi.org/10.5440/1239799.  https://www.osti.gov/servlets/purl/1239799. Pub date:Fri Mar 04 00:00:00 EST 2016

Copy to clipboard


                    
@article{osti_1239799,

  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: Modeling Archive},

  author       = {Tang, Guoping and Mills, Richard and Xu, Xiaofeng and Andre, Ben and Hoffman, Forrest and Painter, Scott and Thornton, Peter and Yuan, Fengming and Bisht, Gautam and Hammond, Glenn and Lichtner, Peter and Kumar, Jitendra},

  abstractNote = {This Modeling Archive is in support of an NGEE Arctic discussion paper under review and available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model.},

  doi          = {10.5440/1239799},

  journal      = {},

  number       = ,

  volume       = ,

  place        = {United States},

  year         = {2016},

  month        = {3}

}

Copy to clipboard

Dataset:

View Dataset

DOI: https://doi.org/10.5440/1239799

Save / Share:

Export Metadata

Save to My Library

Similar records in DOE Data Explorer and OSTI.GOV collections:

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

Journal Article Tang, Guoping ; Yuan, Fengming ; Bisht, Gautam ; ... - Geoscientific Model Development (Online)

Here, we explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon–nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration is used to represent the limitationmore »« less
Cited by 9
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

Journal Article Tang, Guoping ; Yuan, Fengming ; Bisht, Gautam ; ... - Geoscientific Model Development (Online)

We explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon–nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration is used to represent the limitation ofmore »« less
Cited by 9
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

Journal Article Tang, Guoping ; Yuan, Fengming ; Bisht, Gautam ; ... - Geoscientific Model Development (Online)

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 atmore »« less
Cited by 9
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

Journal Article Tang, Guoping ; Yuan, Fengming ; Bisht, Gautam ; ... - Geoscientific Model Development (Online)

Here, we explore coupling to a configurable subsurface reactive transport code as a flexible and extensible approach to biogeochemistry in land surface models. A reaction network with the Community Land Model carbon nitrogen (CLM-CN) decomposition, nitrification, denitrification, and plant uptake is used as an example. We implement the reactions in the open-source PFLOTRAN (massively parallel subsurface flow and reactive transport) code and couple it with the CLM. To make the rate formulae designed for use in explicit time stepping in CLMs compatible with the implicit time stepping used in PFLOTRAN, the Monod substrate rate-limiting function with a residual concentration ismore »« less
Cited by 9
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity

Journal Article Tang, Guoping ; Yuan, Fengming ; Bisht, Gautam ; ... - Geoscientific Model Development (Online)

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 atmore »« less
Cited by 9

Similar Records