Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Biogeochemical modeling”

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 307 records · Page 17

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

Simulating physical and biogeochemical dynamics of the Arctic Ocean [Slides]

Ocean physics (currents, mixing, air-sea transfer of heat and gases) and biogeochemistry (carbon and oxygen cycles, ecosystem dynamics, nutrient distributions) are tightly coupled. We used the DOE-developed Energy Exascale Earth System Model (E3SM) to simulate physical and biogeochemical dynamics of the global ocean, with an emphasis on the western North Pacific. In particular, we used the grid refinement capability of E3SM’s ocean component (LANL-developed MPASOcean) to quantify the effect of model resolution on the carbon and oxygen cycles using a pair of century-long simulations that differ only in the western North Pacific grid. Regional grid refinement has a noticeable impact on biogeochemistry. Distribution of important nutrients (such as nitrate) changes the amount of photosynthesis by phytoplankton, which in turn affects air-sea transfer of CO2 and the local carbon cycle.

58 GEOSCIENCES↗

Toward Hybrid Physics-Machine Learning to improve Land Surface Model predictions

A critical challenge for Land Surface Models (LSMs) is to simulate processes at the surface and the subsurface and their feedbacks to the atmosphere. Even using the same climate forcings, different LSMs predict different surface fluxes and soil moisture conditions due to differences in the formulations of individual processes, parameterizations, and representation of spatial heterogeneity. Ultimately, these differences contribute to the LSM prediction errors and uncertainty. This research seeks to address this challenge by coupling physics-based modeling with state-of-the-art machine learning (ML) techniques to describe complex physical and biogeochemical processes and narrow the gap between model predictions and observations.

58 GEOSCIENCES↗

Revealing Hidden Quinones Through Diagnostic MS² Fragmentation of Peptide–Quinone Adducts

Quinones are redox-active components of natural organic matter that mediate electron transfer and influence biogeochemical processes, but many quinones in pyrogenic organic matter (PyOM) remain unresolved because they ionize poorly by mass spectrometry. Here, we present a peptide-tagging approach to improve detection of cysteine-reactive electrophiles in PyOM, with quinones expected to be a dominant subset based on reaction chemistry and selectivity experiments. A cysteine-containing peptide was used to form Michael-addition adducts, enhancing electrospray ionization and enabling untargeted screening by high-performance liquid chromatography-high-resolution tandem mass spectrometry. The method was benchmarked with five quinone standards and applied to extracts from charred plant material as a discovery-level screen for cysteine-reactive targets. We identified 98 quinone-candidate adducts (mean neutral mass ~603 Da), of which more than 70% were not detectable in native MS1 data. Among formula-assigned features, hidden quinone candidates had median (O+N)/C of 0.391 and normalized oxidation state of carbon of -0.281, consistent with relatively low polarity and low oxidation state. These results reveal a previously inaccessible pool of hidden redox-active compounds in PyOM and provide a framework for prioritizing quinone-like electrophiles for confirmation and incorporation into models of fire-driven biogeochemical cycling.

LC-MS/MS↗

Parameterizing Perennial Bioenergy Crops in Version 5 of the Community Land Model Based on Site‐Level Observations in the Central Midwestern United States

Abstract With projected expansion of biofuel production at a global scale, there is a pressing need to develop adequate representation of bioenergy crops in land surface models to help effectively quantify the biogeophysical and biogeochemical effects of its associated land use changes. This study implements two new perennial bioenergy crops, Miscanthus and switchgrass, into the Community Land Model Version 5 based on site‐level observations from the midwestern United States by modifying parameters associated with photosynthesis, phenology, allocation, decomposition, and carbon cost of nitrogen uptake and integrating concomitantly land management practices. Sensitivity analyses indicate that carbon and energy fluxes of the perennial crops are most sensitive to photosynthesis and phenology parameters. Validation of simulated fluxes against site‐level measurements demonstrates that the model is capable of capturing the overall patterns of energy and carbon fluxes, as well as physiological transitions from leaf emergence to senescence. Compared to annual crops, perennial crops feature longer growing season, greater leaf areas, and higher productivity, leading to increased transpiration, lower annual runoff, and larger carbon uptake. The model simulations suggest that with higher CO 2 assimilation rates and lower demands for nutrients and water, high‐yielding perennial crops are promising alternatives of bioenergy feedstocks compared to traditional annual crops not only for mitigating climate change but also for environmental conservation purposes by reducing fertilizer application and therefore alleviating surface‐ and ground‐water contaminations. Although the local‐scale simulations shed light on potential benefits of using perennial grasses as bioenergy feedstocks, quantifying consequences of their plantations at larger scales warrants additional investigation.

Cheng, Yanyan↗

Wellbore cement alteration and roles of CO 2 and shale during underground hydrogen storage

To mitigate climate change and adopt renewable energy, energy storage is crucial and can be done in the form of hydrogen gas (H 2 ). Subsurface geologic reservoirs are positioned to store H 2 on the largest scales for the longest terms of all potential options. However, H 2 injection may boost reactions that consume hydrogen, generate undesired gases, and alter pore structures of geomedia. To explore the extent of H 2 -associated biotic reactions at a near wellbore location, four experiments were conducted under underground storage conditions with wellbore cement cores and, in most instances, shale samples submerged in synthetic formation brine. Post-reaction gas, aqueous, and solid phase samples were analyzed using olfactory screening and, later, gas chromatography (GC-MS), inductively coupled plasma optical emission spectroscopy (ICP-OES), scanning electron microscopy (SEM), and synchrotron micro-scale x-ray fluorescence (μ-XRF). Within a period of 16 weeks, hydrogen sulfide (H 2 S) was generated in systems containing both H 2 and shale. XRF mapping identified a zone enriched in iron(II) and reduced sulfur along the rim of cement cross sections that was largely associated with CO 2 -induced cement carbonation. Shale did not show noticeable alteration, but there is evidence it contributed to the initial inoculation of the system and provided nutrients for microbes via water-rock interactions. Here, this study considers both rock formations and wellbore cement not previously evaluated concurrently. Findings support understanding and modeling of H 2 -associated biogeochemical reactions during underground hydrogen storage.

36 MATERIALS SCIENCE↗

Knowledge-guided graph machine learning for spatially distributed prediction of daily discharge and nitrogen export dynamics

Spatially distributed prediction of streamflow and nitrogen export dynamics is essential for precision management of agricultural watersheds. While temporal deep learning models such as Long Short-Term Memory (LSTM) have shown strong performance at basin scales, their ability to generalize spatially is limited by insufficient representation of spatial dependencies and flow paths, particularly under data-scarce conditions. To address this gap, we propose HydroGraphNet, a knowledge-guided graph machine learning framework that integrates process-based knowledge and explicit spatial learning into temporal modeling. This framework incorporates directed graph topology to encode watershed connectivity and upstream inflows, with mass balance constraints to improve physical consistency. To enhance generalization in sparsely monitored regions, HydroGraphNet is pretrained on synthetic data generated by the SWAT+ (Soil and Water Assessment Tool Plus) model. We evaluated HydroGraphNet in the Upper Sangamon River Basin (44 HUC-12 subwatersheds, 2001–2020) against two LSTM baselines: a lumped basin-level model and a distributed variant. When benchmarked on SWAT+ simulations in pretraining, HydroGraphNet improved test NSEs by 8.9% (discharge) and 13.7% (NO₃–N load) in temporal extrapolation, and by 27.1% and 34.7% in spatial extrapolation, relative to the Lumped LSTM baseline. After fine-tuning with USGS monitoring data, the model achieved mean test NSE (KGE) scores of 0.768 (0.861) for discharge and 0.626 (0.664) for NO₃–N load, substantially outperforming baselines. Attribution analysis further highlighted the importance of upstream inflow representation and graph-based spatial learning in capturing cross-subwatershed dependencies. The model also reproduced seasonal hydrological and biogeochemical patterns consistent with known processes, demonstrating its robustness and process fidelity for spatially distributed prediction. Altogether, HydroGraphNet advances the integration of physical knowledge and spatially explicit learning in hydrological modeling, offering a generalizable framework for distributed modeling to support spatially targeted water quality management in data-scarce watersheds.

54 ENVIRONMENTAL SCIENCES↗

Multiple microbial guilds mediate soil methane cycling along a wetland salinity gradient

ABSTRACT Estuarine wetlands harbor considerable carbon stocks, but rising sea levels could affect their ability to sequester soil carbon as well as their potential to emit methane (CH 4 ). While sulfate loading from seawater intrusion may reduce CH 4 production due to the higher energy yield of microbial sulfate reduction, existing studies suggest other factors are likely at play. Our study of 11 wetland complexes spanning a natural salinity and productivity gradient across the San Francisco Bay and Delta found that while CH 4 fluxes generally declined with salinity, they were highest in oligohaline wetlands (ca. 3-ppt salinity). Methanogens and methanogenesis genes were weakly correlated with CH 4 fluxes but alone did not explain the highest rates observed. Taxonomic and functional gene data suggested that other microbial guilds that influence carbon and nitrogen cycling need to be accounted for to better predict CH 4 fluxes at landscape scales. Higher methane production occurring near the freshwater boundary with slight salinization (and sulfate incursion) might result from increased sulfate-reducing fermenter and syntrophic populations, which can produce substrates used by methanogens. Moreover, higher salinities can solubilize ionically bound ammonium abundant in the lower salinity wetland soils examined here, which could inhibit methanotrophs and potentially contribute to greater CH 4 fluxes observed in oligohaline sediments. IMPORTANCE Low-level salinity intrusion could increase CH4 flux in tidal freshwater wetlands, while higher levels of salinization might instead decrease CH4 fluxes. High CH4 emissions in oligohaline sites are concerning because seawater intrusion will cause tidal freshwater wetlands to become oligohaline. Methanogenesis genes alone did not account for landscape patterns of CH4 fluxes, suggesting mechanisms altering methanogenesis, methanotrophy, nitrogen cycling, and ammonium release, and increasing decomposition and syntrophic bacterial populations could contribute to increases in net CH4 flux at oligohaline salinities. Improved understanding of these influences on net CH4 emissions could improve restoration efforts and accounting of carbon sequestration in estuarine wetlands. More pristine reference sites may have older and more abundant organic matter with higher carbon:nitrogen compared to wetlands impacted by agricultural activity and may present different interactions between salinity and CH4. This distinction might be critical for modeling efforts to scale up biogeochemical process interactions in estuarine wetlands.

59 BASIC BIOLOGICAL SCIENCES↗

Calcium is associated with specific soil organic carbon decomposition products at Blodgett Forest Research Center, Georgetown, California as analysed with scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy

This data is from the paper calcium is associated with specific soil organic carbon decomposition products, published in SOIL. DOI: https://doi.org/10.5194/soil-11-381-2025, 2025.This file contains CSVs with spectral data and bulk soil data and there is no specific program required to open this data. The data includes Scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy. data from the measurement of samples from the Whole-soil Warming project, run by the Belowground Biogeochemistry team at Blodgett Forest Research Center, Georgetown, California run by the University of California, Berkeley. It also includes bulk soil chemical properties. The University of California's Blodgett Forest Research Station (Forest) is situated in the Sierra Nevada foothills (1370 m a.s.l.) near Georgetown, California. The samples were collected from here: 38.912013, -120.661469, https://maps.app.goo.gl/291bCJ1zVqUhgktz6. The Forest soils were characterised as Alfisols, which are equivalent to Dystric Cambisols (IUSS Working Group WRB, 2015), and formed in granitic parent materials, in a temperate climate, under thinned, mixed-coniferous forest (Fig. S3; Gaudinski et al., 2009). With these analyses we aimed to answer the question, is calcium associated with a specific type of organic matter enriched in aromatic and phenolic carbon at the microscale in samples from Blodgett Forest Research Center? and how does this specific type of carbon respond to experiments targetted at removing and adding calcium to the soils, specifically cation exchange and incubation after calcium addition? Abstract from the paper can be found below: Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil-types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca-SOC is removed through cation-exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, highlighting a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca, and in turn, this association inhibits dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots on the microscale shifting local decomposition processes and thereby explaining the colocation of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth System Models could improve our understanding, predictions, and management of carbon dynamics in soils, and account for their response to Ca-rich amendments.

54 ENVIRONMENTAL SCIENCES↗

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calcium is associated with specific soil organic carbon decomposition products

Abstract. Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here, we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca–SOC is removed through cation exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, this highlights a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca; in turn, this association can inhibit dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots at the microscale shifting local decomposition processes and thereby explaining the co-location of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth system models could improve our understanding, predictions, and management of carbon dynamics in soils, as well as accounting for their response to Ca-rich amendments.

Rowley, Mike C. (ORCID:0000000224407855)↗

Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions

Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.

54 ENVIRONMENTAL SCIENCES↗

Allometric Scaling of Hyporheic Respiration Across Basins in the Pacific Northwest United States

Abstract Hyporheic zones regulate biogeochemical processes in streams and rivers, but high spatiotemporal heterogeneity makes it difficult to predict how these processes scale from individual reaches to river basins. Recent work applying allometric scaling (i.e., power‐law relationships between size and function) to river networks provides a new paradigm for understanding cumulative hyporheic biogeochemical processes. We used previously published model predictions of reach‐scale hyporheic aerobic respiration to explore patterns in allometric scaling across two climatically divergent basins with differing characteristics in the Pacific Northwest, United States. In the model, hydrologic exchange fluxes (HEFs) regulate hyporheic respiration, so we examined how HEFs might influence allometric scaling of respiration. We found consistent scaling behaviors where HEFs were either very low or very high, but differences between basins when HEFs were moderate. Our findings provide initial model‐generated hypotheses for factors influencing allometric scaling of hyporheic respiration. These hypotheses can be used to optimize new data generation efforts aimed at developing predictive understanding of allometries that can, in turn, be used to scale biogeochemical dynamics across watersheds.

59 BASIC BIOLOGICAL SCIENCES↗

Data from: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N 2 O) fluxes from US croplands

Nitrous oxide (N₂O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N₂O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N2O fluxes from US cropland. Trained and validated on approximately 12,000 N2O chamber measurements at 17 U.S. Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N2O at both training (R² = 0.84, RMSE = 16.4 g N ha⁻¹ d⁻¹) and held-out testing sites (R² = 0.84, RMSE = 6.2 g N ha⁻¹ d⁻¹). Analyses identified six dominant N₂O drivers: soil organic carbon (SOC), NH₄⁺, NO₃⁻, water-filled pore space (WFPS), soil temperature, and biomass production. Wet, warm soils produced large N₂O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N₂O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.

agricultural sciences↗

Coupled machine learning–ecosystem ensemble models substantially improve predictions of nitrous oxide (N 2 O) fluxes from US croplands

Nitrous oxide (N 2 O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N 2 O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high-emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N 2 O fluxes from US cropland. Trained and validated on ~12,000 N 2 O chamber measurements at 17 US Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N 2 O at both training (R 2 = 0.84, RMSE = 16.4 g N ha −1 d −1 ) and held-out testing sites (R 2 = 0.84, RMSE = 6.2 g N ha −1 d −1 ). Analyses identified six dominant N 2 O drivers: soil organic carbon (SOC), NH 4 + , NO 3 - , water-filled pore space, temperature, and aboveground biomass production. Wet, warm soils produced large N 2 O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N 2 O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.

AI↗

Mountainous Floodplain Connectivity in Response to Hydrological Transitions

In mountainous watersheds, floodplain sediments are typically characterized by gravel bed layers capped by an overlying soil unit that serves as a hotspot for biogeochemical reactivity. However, the influence of soil biogeochemistry on gravel bed underflow composition remains unclear, especially during hydrological transitions that alter the vertical connectivity between overlaying soils and the underlying gravel bed. This study investigates these dynamics by measuring hydraulic gradients and water compositions over three hydrological years in a typical mountainous, low-order stream floodplain in the Upper Colorado River Basin. Results indicate that the timing of hydrological conditions strongly influences the vertical exchanges that control water quality. Specifically, during flooding events such as beaver ponding, that induce downward flushing of the soil, anoxic conditions prevalent in the biogeochemically active soil are transferred downstream via gravel bed underflow. Conversely, snowmelt and drought conditions increase oxic conditions in the gravel bed due to diminished hydrological connectivity with the overlying soil. To compare water quality response to hydrological transitions across similar floodplain environments, we propose a conceptual model that quantifies the inundation-induced flushing of soil porewater to measure solute exchange efficiency with the gravel bed solute convergence efficiency (SCE). This model provides a framework for quantifying biogeochemical processes in hydrological underflow systems, which is critical for water and elemental budgets in these globally important mountainous ecosystems.

58 GEOSCIENCES↗

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗

2016 International Land Model Benchmarking (ILAMB) Workshop Report

As earth system models (ESMs) become increasingly complex, there is a growing need for comprehensive and multi-faceted evaluation of model projections. To advance understanding of terrestrial biogeochemical processes and their interactions with hydrology and climate under conditions of increasing atmospheric carbon dioxide, new analysis methods are required that use observations to constrain model predictions, inform model development, and identify needed measurements and field experiments. Better representations of biogeochemistryclimate feedbacks and ecosystem processes in these models are essential for reducing the acknowledged substantial uncertainties in 21st century climate change projections.

earth system models (ESMs)↗