U.S. Department of Energy Office of Biological and Environmental Research

BER Research Highlights


An Improved Numerical Method for Solving Depth-Resolved Biogeochemical Models
Published: February 20, 2018
Posted: April 10, 2018

A method of alternating characteristics with application to advection-dominated environmental systems.

The Science
Scientists at Lawrence Berkeley National Lab (LBNL) propose a numerical integration method, termed the method of alternating characteristics (MAC), to efficiently and accurately solve systems of partial differential equations that arise in modeling environmental processes. They highlight the advantages of the MAC with emphasis on advection-dominated environmental systems with biogeochemical reactions.

The Impact
The proposed method is uniquely suited for solving depth-resolved models of advection-dominated environmental systems with biogeochemical reactions and offers advantages in performance over other numerical integration schemes that often require considerable computational resources.

Summary
Here, LBNL scientists present a numerical integration method for solving systems of partial differential equations (PDEs) that arise in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The proposed method, termed the method of alternating characteristics (MAC), then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.

Contacts (BER PM)
Daniel Stover
Terrestrial Ecosystem Science, SC-23.1
Daniel.Stover@science.doe.gov (301-903-0289)

(PI Contact)
William J. Riley
Lawrence Berkeley National Laboratory
wjriley@lbl.gov; 510-495-2223

Funding
This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research, Terrestrial Ecosystem Science Program, under contract number DE-AC02-05CH11231. K.G. acknowledges support from the U.S. Department of Energy, Office of Science, Office of Workforce Development for Teachers and Scientists, Office of Science Graduate Student Research (SCGSR) program. The SCGSR program is administered by the Oak Ridge Institute for Science and Education (ORISE) for the DOE. ORISE is managed by ORAU under contract number DE-SC0014664.

Publications
Georgiou, K., J. Harte, A. Mesbah, and W. J. Riley. "A method of alternating characteristics with application to advection-dominated environmental systems." Computational Geosciences, (2018). [DOI:10.1007/s10596-018-9729-5]

Topic Areas:

  • Research Area: Earth and Environmental Systems Modeling
  • Research Area: Terrestrial Ecosystem Science

Division: SC-23.1 Climate and Environmental Sciences Division, BER

 

BER supports basic research and scientific user facilities to advance DOE missions in energy and environment. More about BER

Recent Highlights

May 10, 2019
Quantifying Decision Uncertainty in Water Management via a Coupled Agent-Based Model
Considering risk perception can improve the representation of human decision-making processes in age [more...]

May 09, 2019
Projecting Global Urban Area Growth Through 2100 Based on Historical Time Series Data and Future Scenarios
Study provides country-specific urban area growth models and the first dataset on country-level urba [more...]

May 05, 2019
Calibrating Building Energy Demand Models to Refine Long-Term Energy Planning
A new, flexible calibration approach improved model accuracy in capturing year-to-year changes in bu [more...]

May 03, 2019
Calibration and Uncertainty Analysis of Demeter for Better Downscaling of Global Land Use and Land Cover Projections
Researchers improved the Demeter model’s performance by calibrating key parameters and establi [more...]

Apr 22, 2019
Representation of U.S. Warm Temperature Extremes in Global Climate Model Ensembles
Representation of warm temperature events varies considerably among global climate models, which has [more...]

List all highlights (possible long download time)