BER launches Environmental System Science Program. Visit our new website under construction!

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 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.

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 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.

BER Program Manager
Daniel Stover
Terrestrial Ecosystem Science, SC-23.1 (301-903-0289)

Principal Investigator
William J. Riley
Lawrence Berkeley National Laboratory
Berkeley, CA 94720; 510-495-2223

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

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

Topic Areas:

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

Division: SC-33.1 Earth 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

Mar 23, 2021
Molecular Connections from Plants to Fungi to Ants
Lipids transfer energy and serve as an inter-kingdom communication tool in leaf-cutter ants&rsqu [more...]

Mar 19, 2021
Microbes Use Ancient Metabolism to Cycle Phosphorus
Microbial cycling of phosphorus through reduction-oxidation reactions is older and more widespre [more...]

Feb 22, 2021
Warming Soil Means Stronger Microbe Networks
Soil warming leads to more complex, larger, and more connected networks of microbes in those soi [more...]

Jan 27, 2021
Labeling the Thale Cress Metabolites
New data pipeline identifies metabolites following heavy isotope labeling.

Analysis [more...]

Aug 31, 2020
Novel Bacterial Clade Reveals Origin of Form I Rubisco

  • All plant biomass is sourced from the carbon-fixing enzyme Rub [more...]

List all highlights (possible long download time)