Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity
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...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , , , , , , , , , , |
Language: | unknown |
Published: |
2022
|
Subjects: | |
Online Access: | http://www.osti.gov/servlets/purl/1240537 https://www.osti.gov/biblio/1240537 https://doi.org/10.5194/gmd-9-927-2016 |
id |
ftosti:oai:osti.gov:1240537 |
---|---|
record_format |
openpolar |
spelling |
ftosti:oai:osti.gov:1240537 2023-07-30T04:02:08+02:00 Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity Tang, Guoping Yuan, Fengming Bisht, Gautam Hammond, Glenn E. Lichtner, Peter C. Collier, Nathaniel O. Kumar, Jitendra Mills, Richard T. Xu, Xiaofeng Andre, Ben Hoffman, Forrest M. Painter, Scott L. Thornton, Peter E. 2022-05-23 application/pdf http://www.osti.gov/servlets/purl/1240537 https://www.osti.gov/biblio/1240537 https://doi.org/10.5194/gmd-9-927-2016 unknown http://www.osti.gov/servlets/purl/1240537 https://www.osti.gov/biblio/1240537 https://doi.org/10.5194/gmd-9-927-2016 doi:10.5194/gmd-9-927-2016 58 GEOSCIENCES 2022 ftosti https://doi.org/10.5194/gmd-9-927-2016 2023-07-11T09:05:10Z 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 ... Other/Unknown Material Arctic SciTec Connect (Office of Scientific and Technical Information - OSTI, U.S. Department of Energy) Arctic Geoscientific Model Development 9 3 927 946 |
institution |
Open Polar |
collection |
SciTec Connect (Office of Scientific and Technical Information - OSTI, U.S. Department of Energy) |
op_collection_id |
ftosti |
language |
unknown |
topic |
58 GEOSCIENCES |
spellingShingle |
58 GEOSCIENCES Tang, Guoping Yuan, Fengming Bisht, Gautam Hammond, Glenn E. Lichtner, Peter C. Collier, Nathaniel O. Kumar, Jitendra Mills, Richard T. Xu, Xiaofeng Andre, Ben Hoffman, Forrest M. Painter, Scott L. Thornton, Peter E. Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
topic_facet |
58 GEOSCIENCES |
description |
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 ... |
author |
Tang, Guoping Yuan, Fengming Bisht, Gautam Hammond, Glenn E. Lichtner, Peter C. Collier, Nathaniel O. Kumar, Jitendra Mills, Richard T. Xu, Xiaofeng Andre, Ben Hoffman, Forrest M. Painter, Scott L. Thornton, Peter E. |
author_facet |
Tang, Guoping Yuan, Fengming Bisht, Gautam Hammond, Glenn E. Lichtner, Peter C. Collier, Nathaniel O. Kumar, Jitendra Mills, Richard T. Xu, Xiaofeng Andre, Ben Hoffman, Forrest M. Painter, Scott L. Thornton, Peter E. |
author_sort |
Tang, Guoping |
title |
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
title_short |
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
title_full |
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
title_fullStr |
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
title_full_unstemmed |
Coupling a Reactive Transport Code with a Global Land Surface Model for Mechanistic Biogeochemistry Representation: 1. Addressing the Challenge of Nonnegativity |
title_sort |
coupling a reactive transport code with a global land surface model for mechanistic biogeochemistry representation: 1. addressing the challenge of nonnegativity |
publishDate |
2022 |
url |
http://www.osti.gov/servlets/purl/1240537 https://www.osti.gov/biblio/1240537 https://doi.org/10.5194/gmd-9-927-2016 |
geographic |
Arctic |
geographic_facet |
Arctic |
genre |
Arctic |
genre_facet |
Arctic |
op_relation |
http://www.osti.gov/servlets/purl/1240537 https://www.osti.gov/biblio/1240537 https://doi.org/10.5194/gmd-9-927-2016 doi:10.5194/gmd-9-927-2016 |
op_doi |
https://doi.org/10.5194/gmd-9-927-2016 |
container_title |
Geoscientific Model Development |
container_volume |
9 |
container_issue |
3 |
container_start_page |
927 |
op_container_end_page |
946 |
_version_ |
1772812865147764736 |