Search NASA⌕ Search

SEARCH · Search NASA

Results for “soil biogeochemistry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Climate warming enhances biodiversity and stability of grassland soil phosphorus-cycling microbial communities

Abstract Climate warming poses significant challenges to global phosphorus sustainability, an essential component of Earth biogeochemistry cycling and water-food-energy nexus. Despite the crucial role of polyphosphate-accumulating organism as key functional microbial agents in phosphorus cycling, the impacts of global climate warming on polyphosphate accumulating organism communities remain largely enigmatic. This study investigates the effects of climate warming on the taxonomic, network, and functional profiles of soil bacterial polyphosphate-accumulating organisms, leveraging fluorescence-activated cell sorting and single-cell Raman spectroscopy. Climate warming enhances both taxonomic and functional biodiversity of polyphosphate-accumulating organisms via biotic interactions and environmental filtering, with observed functionality-biodiversity relationships supporting the functional redundancy theory. Furthermore, polyphosphate-accumulating organism network complexity and stability rise under warming with strengthened positive relationships, supporting stress gradient hypothesis and the belief that complexity begets stability. Finally, polyphosphate-accumulating organisms are significantly correlated to key ecosystem functioning in carbon and phosphorus cycling under warming. Our study suggests that preserving polyphosphate-accumulating organism communities is crucial for maintaining soil ecosystem functioning and sustainable phosphorus management in a warming world and opens avenues for predicting the responses of other functional microbial groups to climate change, beneficially or maliciously.

Environmental Sciences & Ecology↗

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗

Feasibility of Formulating Ecosystem Biogeochemical Models From Established Physical Rules

Abstract To improve the predictive capability of ecosystem biogeochemical models (EBMs), we discuss the feasibility of formulating biogeochemical processes using physical rules that have underpinned the many successes in computational physics and chemistry. We argue that the currently popular empirically based approaches, such as multiplicative empirical response functions and the law of the minimum, will not lead to EBM formulations that can be continuously refined to incorporate improved mechanistic understanding and empirical observations of biogeochemical processes. Instead, we propose that EBM parameterizations, as a lossy data compression problem, can be better formulated using established physical rules widely used in computational physics and chemistry, and different biogeochemical processes can be more robustly integrated within a reactive‐transport framework. Through several examples, we demonstrate how mathematical representations derived from physical rules can improve understanding of relevant biogeochemical processes and enable more effective communication between modelers, observationalists, and experimentalists regarding essential questions, such as what measurements are needed to meaningfully inform models and how can models generate new process‐level hypotheses to test in empirical studies. Finally, while empirical models with more parameters are often less robust, physical rules‐based models can be more robust and show lower predictive equifinality, stemming from their enhanced consistency in representations of processes, interactions and spatial scaling.

54 ENVIRONMENTAL SCIENCES↗

Arbuscular Mycorrhizal Tree Communities Have Greater Soil Fungal Diversity and Relative Abundances of Saprotrophs and Pathogens than Ectomycorrhizal Tree Communities

Trees associating with different mycorrhizas often differ in their effects on litter decomposition, nutrient cycling, soil organic matter (SOM) dynamics, and plant-soil interactions. For example, due to differences between arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) tree leaf and root traits, ECM-associated soil has lower rates of C and N cycling and lower N availability than AM-associated soil. These observations suggest that many groups of nonmycorrhizal fungi should be affected by the mycorrhizal associations of dominant trees through controls on nutrient availability. To test this overarching hypothesis, we explored the influence of predominant forest mycorrhizal type and mineral N availability on soil fungal communities using next-generation amplicon sequencing. Soils from four temperate hardwood forests in southern Indiana, United States, were studied; three forests formed a natural gradient of mycorrhizal dominance (100% AM tree basal area to 100% ECM basal area), while the fourth forest contained a factorial experiment testing long-term N addition in both dominant mycorrhizal types. We found that overall fungal diversity, as well as the diversity and relative abundance of plant pathogenic and saprotrophic fungi, increased with greater AM tree dominance. Additionally, tree community mycorrhizal associations explained more variation in fungal community composition than abiotic variables, including soil depth, SOM content, nitrification rate, and mineral N availability. Our findings suggest that tree mycorrhizal associations may be good predictors of the diversity, composition, and functional potential of soil fungal communities in temperate hardwood forests. These observations help explain differing biogeochemistry and community dynamics found in forest stands dominated by differing mycorrhizal types.

59 BASIC BIOLOGICAL SCIENCES↗

Chemistry data from soils and soil incubation experiments from the whole-soil warming experiment at Blodgett Forest, CA, 2018, from: “Metabolic capabilities mute positive response to direct and indirect impacts of warming throughout the soil profile”

This dataset contains soil chemistry data of soils from the whole-soil warming experiment at Blodgett Forest in June 2018, and from laboratory incubation experiments conducted with a subset of those soils. Soil samples were collected at five different depths down to 80 cm from soil plots continuously heated by 4°C above ambient temperature for 4.5 years, and from the respective unheated control plots. Incubation experiments were conducted for 28 days with soils collected between 10-30 cm (“surface soils”) and 60-80 cm (“subsoils”) from heated and unheated field plots, amended in the laboratory with cellobiose (C, carbon amendment), or cellobiose plus ammonium and phosphate (CNP, combined carbon, nitrogen and phosphorus amendment). The dataset includes separate files for: (i) total carbon and nitrogen concentrations of initial field soils; (ii) carbon dioxide efflux rates at multiple times throughout the incubation experiments, and respective cumulative carbon dioxide efflux rates, as well as inorganic nitrogen concentrations at the beginning and end of the incubations, and net nitrogen mineralization rates and carbon use efficiency calculated at the end of the incubation. Analysis of these data has been published in Dove et al., (2021), “Metabolic capabilities mute positive response to direct and indirect impacts of warming throughout the soil profile. Nat. Commun. 12:2089”, https://doi.org/10.1038/s41467-021-22408-5. This research was performed as part of the TES Belowground Biogeochemistry SFA project at Lawrence Berkeley National Laboratory.

54 ENVIRONMENTAL SCIENCES↗

Drought re-routes soil microbial carbon metabolism towards emission of volatile metabolites in an artificial tropical rainforest

Drought impacts on microbial activity can alter soil carbon fate and lead to the loss of stored carbon to the atmosphere as CO 2 and volatile organic compounds (VOCs). Here we examined drought impacts on carbon allocation by soil microbes in the Biosphere 2 artificial tropical rainforest by tracking 13 C from position-specific 13 C-pyruvate into CO 2 and VOCs in parallel with multi-omics. During drought, efflux of 13 C-enriched acetate, acetone and C 4 H 6 O 2 (diacetyl) increased. These changes represent increased production and buildup of intermediate metabolites driven by decreased carbon cycling efficiency. Simultaneously, 13 C-CO 2 efflux decreased, driven by a decrease in microbial activity. However, the microbial carbon allocation to energy gain relative to biosynthesis was unchanged, signifying maintained energy demand for biosynthesis of VOCs and other drought-stress-induced pathways. Overall, while carbon loss to the atmosphere via CO 2 decreased during drought, carbon loss via efflux of VOCs increased, indicating microbially induced shifts in soil carbon fate.

59 BASIC BIOLOGICAL SCIENCES↗

Soil Incubation, Methane Pool Dilution, Edgewater, MD, USA, July 2022 v1.0.1

Isotopic pool dilution is a powerful approach to quantify gross biogeochemical transformation rates, but remains seldom used despite its potential. To facilitate broader implementation of pool dilution methods, we developed a user-friendly R package that optimizes gross production and consumption rates (and optionally fractionation constants as well) based on standard pool dilution time series data. This data package is the soil methane pool dilution data and code for building and using this software, including source code (R) and raw data (.csv). This package features extensive documentation and example analyses, and is easily integrated into analytical pipelines as R software is free and open source. Operation of this package requires R version 2.10 or newer, and we recommend its use in conjunction with packages testthat (3.0.0), knitr, rmarkdown, ggploto2, tibble, and tidyr.

54 ENVIRONMENTAL SCIENCES↗

Characterizing the Atmospheric Mn Cycle and Its Impact on Terrestrial Biogeochemistry

Abstract The role of manganese (Mn) in ecosystem carbon (C) biogeochemical cycling is gaining increasing attention. While soil Mn is mainly derived from bedrock, atmospheric deposition could be a major source of Mn to surface soils, with implications for soil C cycling. However, quantification of the atmospheric Mn cycle, which comprises emissions from natural (desert dust, sea salts, volcanoes, primary biogenic particles, and wildfires) and anthropogenic sources (e.g., industrialization and land‐use change due to agriculture), transport, and deposition, remains uncertain. Here, we use compiled emission data sets for each identified source to model and quantify the atmospheric Mn cycle by combining an atmospheric model and in situ atmospheric concentration measurements. We estimated global emissions of atmospheric Mn in aerosols (<10 μm in aerodynamic diameter) to be 1,400 Gg Mn year −1 . Approximately 31% of the emissions come from anthropogenic sources. Deposition of the anthropogenic Mn shortened Mn “pseudo” turnover times in 1‐m‐thick surface soils (ranging from 1,000 to over 10,000,000 years) by 1–2 orders of magnitude in industrialized regions. Such anthropogenic Mn inputs boosted the Mn‐to‐N ratio of the atmospheric deposition in non‐desert dominated regions (between 5 × 10 −5 and 0.02) across industrialized areas, but that was still lower than soil Mn‐to‐N ratio by 1–3 orders of magnitude. Correlation analysis revealed a negative relationship between Mn deposition and topsoil C density across temperate and (sub)tropical forests, consisting with atmospheric Mn deposition enhancing carbon respiration as seen in in situ biogeochemical studies.

Environmental Sciences & Ecology↗

A restructured and updated global soil respiration database (SRDB-V5)

Field-measured soil respiration (R S , the soil-to-atmosphere CO 2 flux) observations were compiled into a global soil respiration database (SRDB) a decade ago, a resource that has been widely used by the biogeochemistry community to advance our understanding of R S dynamics. Novel carbon cycle science questions require updated and augmented global information with better interoperability among datasets. Here, we restructured and updated the global R S database to version SRDB-V5. The updated version has all previous fields revised for consistency and simplicity, and it has several new fields to include ancillary information (e.g., R S measurement time, collar insertion depth, collar area). The new SRDB-V5 includes published papers through 2017 (800 independent studies), where total observations increased from 6633 in SRDB-V4 to 10 366 in SRDB-V5. The SRDB-V5 features more R S data published in the Russian and Chinese scientific literature and has an improved global spatio-temporal coverage and improved global climate space representation. We also restructured the database so that it has stronger interoperability with other datasets related to carbon cycle science. For instance, linking SRDB-V5 with an hourly timescale global soil respiration database (HGRsD) and a community database for continuous soil respiration (COSORE) enables researchers to explore new questions. The updated SRDB-V5 aims to be a data framework for the scientific community to share seasonal to annual field R S measurements, and it provides opportunities for the biogeochemistry community to better understand the spatial and temporal variability in R S , its components, and the overall carbon cycle.

54 ENVIRONMENTAL SCIENCES↗

1000 Soils Pilot Dataset, version 8, May 2025

This record hosts data generated by the 1000 Soils Pilot. Data will be updated as more become available. Please see the most recent data upload for current data. A beta visualization tool is available for some data types at https://shinyproxy.emsl.pnnl.gov/app/1000soils. Please submit any suggestions or comments through the 'contact' tab. We are actively working to improve visualizations and value all feedback. Data completed include: Geochemistry, texture, respiration, and enzyme activities FTICR-MS organic matter chemistry Microbial biomass C and N TOC/TDN of water-extractable OM X-ray computed tomography (derived metrics available here, raw data available upon request) Metagenomes; a variety of data formats are available upon request Soil hydraulic properties Data in progress: LC-MS/MS in development, timeline TBD, inquire for status 1000S_processed_BGC_summary.csv contains all available biogeochemical data; microbial biomass C and N; and TOC/TDN of water-extractable OM; and 1000S_Tomography.xslx contains a summary of data generated via X-ray computed tomography. icr_v2_corems2.csv contains FTICR-MS data processed by CoreMS version 2. These data are merged by formula across instrument runs to enable cross-sample comparisons. Technical replicates are merged by retaining peaks present in 2 out of 3 replicates. 1000Soils_Metadata_Site_Mastersheet_v1.csv contains site information. Soil Hydraulics_corrected_02042025.xlsx contains soil hydraulics information. Readme File_v4.xlsx is the readme file. Please contact the MONet project (monet.emsl@pnnl.gov) or Emily Graham (emily.graham@pnnl.gov) with questions. The following file and all raw data are available upon request: icr_by_mass_for_single_sample_analysis_only.csv contains FTICR-MS data processed by CoreMS and is intended for usage in the calculation of biochemical transformations within samples only. These data are not acceptable for cross-sample comparison of masses because they are from multiple instrument runs. For more information, please see: https://www.emsl.pnnl.gov/monet and https://sc-data.emsl.pnnl.gov/monet Acknowledgment: Soil data were provided by the Molecular Observation Network (MONet) at the Environmental Molecular Sciences Laboratory (https://ror.org/04rc0xn13), a DOE Office of Science user facility sponsored by the Biological and Environmental Research program under Contract No. DE-AC05-76RL01830. The work (proposal: 10.46936/10.25585/60008970) conducted by the U.S. Department of Energy, Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science user facility, is supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. The Molecular Observation Network (MONet) database is an open, FAIR, and publicly available compilation of the molecular and microstructural properties of soil. Data in the MONet open science database can be found at https://sc-data.emsl.pnnl.gov/.

biogeochemistry↗

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem

Background: Viruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics’ seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes. Results: We identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages ( Leviviricetes ), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a “seed-bank” viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: “Kill-The-Winner” which proposes that lytic phages are actively infecting abundant bacteria, and “Piggyback-The-Persistent” which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth. Conclusion: Taken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗

A framework for soil microbial ecology in urban ecosystems

Nearly all ecosystems host diverse microbiomes that support vital ecosystem processes. At the same time, these ecosystems and their microbiomes are increasingly altered by human activities, particularly in highly managed urban environments. While microbial ecologists are beginning to understand the drivers of microbial assembly and the link between community structure and function in many ecosystems, few of these advances have been applied to urban ecosystems. In this synthesis, we review research on the urban soil microbiome and develop a framework to integrate soil microbial communities with urban ecosystem function. We identify disturbance, altered resources, and heterogeneity as key drivers through which human activities including urban development affect soils and their resident microorganisms. Steep environmental gradients in many urban systems present a unique opportunity to address fundamental questions in microbial ecology, such as how microbes respond to stress and how biogeochemical rates relate to microbial diversity and composition. Soil microbiomes in cities also provide ecosystem services and harms, making it crucial to understand how human activity drives those functions and the consequences for environmental and human health. We argue that much-needed integration across disturbance ecology, urban ecology, and microbial ecology will help generate practical and equitable strategies for managing ecosystem benefits in cities where most humans now live.

54 ENVIRONMENTAL SCIENCES↗

Microbial dormancy and its impacts on northern temperate and boreal terrestrial ecosystem carbon budget

A large amount of soil carbon in northern temperate and boreal regions could be emitted as greenhouse gases in a warming future. However, lacking detailed microbial processes such as microbial dormancy in current biogeochemistry models might have biased the quantification of the regional carbon dynamics. Here the effect of microbial dormancy was incorporated into a biogeochemistry model to improve the quantification for the last century and this century. Compared with the previous model without considering the microbial dormancy, the new model estimated the regional soils stored 75.9 Pg more C in the terrestrial ecosystems during the last century and will store 50.4 and 125.2 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, in this century. This study highlights the importance of the representation of microbial dormancy in earth system models to adequately quantify the carbon dynamics in the northern temperate and boreal natural terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Synthetic Soil Aggregates: Bioprinted Habitats for High-Throughput Microbial Metaphenomics

The dynamics of microbial processes are difficult to study in natural soil, owing to the small spatial scales on which microorganisms operate and to the opacity and chemical complexity of the soil habitat. To circumvent these challenges, we have created a 3D-bioprinted habitat that mimics aspects of natural soil aggregates while providing a chemically defined and translucent alternative culturing method for soil microorganisms. Our Synthetic Soil Aggregates (SSAs) retain the porosity, permeability, and patchy resource distribution of natural soil aggregates—parameters that are expected to influence emergent microbial community interactions. We demonstrate the printability and viability of several different microorganisms within SSAs and show how the SSAs can be integrated into a multi-omics workflow for single SSA resolution genomics, metabolomics, proteomics, lipidomics, and biogeochemical assays. We study the impact of the structured habitat on the distribution of a model co-culture microbial community and find that it is significantly different from the spatial organization of the same community in liquid culture, indicating a potential for SSAs to reproduce naturally occurring emergent community phenotypes. The SSAs have the potential as a tool to help researchers quantify microbial scale processes in situ and achieve high-resolution data from the interplay between environmental properties and microbial ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗