SEARCH · Search NASA
Results for “Reactive Transport 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.
Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments
Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.
High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor
Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.
Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling
Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.
BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale
Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.
Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model
Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.
Effects of Heterogeneity in Reactive Transport Models for Carbon Sequestration in the Tuscaloosa Group, Southwest Alabama
Effects of Heterogeneity in Reactive Transport Models for Carbon Sequestration in the Tuscaloosa Group, Southwest Alabama
Reactive transport modeling for supporting climate resilience at groundwater contamination sites
Abstract. Climate resilience is an emerging issue at contaminated sites and hazardous waste sites, since projected climate shifts (e.g., increased/decreased precipitation) and extreme events (e.g., flooding, drought) could affect ongoing remediation or closure strategies. In this study, we develop a reactive transport model (Amanzi) for radionuclides (uranium, tritium, and others) and evaluate how different scenarios under climate change will influence the contaminant plume conditions and groundwater well concentrations. We demonstrate our approach using a two-dimensional (2D) reactive transport model for the Savannah River Site F-Area, including mineral reaction and sorption processes. Different recharge scenarios are considered by perturbing the infiltration rate from the base case as well as considering cap-failure and climate projection scenarios. We also evaluate the uranium and nitrate concentration ratios between scenarios and the base case to isolate the sorption effects with changing recharge rates. The modeling results indicate that the competing effects of dilution and remobilization significantly influence pH, thus changing the sorption of uranium. At the maximum concentration on the breakthrough curve, higher aqueous uranium concentration implies that sorption is reduced with lower pH due to remobilization. To better evaluate the climate change impacts in the future, we develop the workflow to include the downscaled CMIP5 (Coupled Model Intercomparison Project) climate projection data in the reactive transport model and evaluate how residual contamination evolves through 2100 under four climate Representative Concentration Pathway (RCP) scenarios. The integration of climate modeling data and hydrogeochemistry models enables us to quantify the climate change impacts, assess which impacts need to be planned for, and therefore assist climate resiliency efforts and help guide site management.
Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations
Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.
Understanding Contaminant Migration within a Dynamic River Corridor through Field Experiments and Reactive Transport Modeling
The behavior of a persistent uranium plume within an extended river corridor at the DOE Hanford site is dominantly controlled by river stage fluctuations in the adjacent Columbia River. The plume behavior is further complicated by substantial heterogeneity in physical and geochemical properties of the host aquifer sediments. Multi-scale field and laboratory experiments and reactive transport modeling were integrated to understand the complex plume behavior influenced by highly variable hydrologic and geochemical conditions in time and space. In this paper, we (1) describe multiple data sets from field-scale uranium adsorption and desorption experiments performed at our experimental well-field, (2) develop a reactive transport model that incorporates hydrologic and geochemical heterogeneities characterized from multi-scale and multi-type datasets and a surface complexation reaction network based on laboratory studies, and (3) compare the modeling and observation results to provide insights on how to refine the conceptual model and reduce prediction uncertainties. The experimental results revealed significant spatial variability in uranium adsorption/desorption behavior, while modeling demonstrated that ambient hydrologic and geochemical conditions and heterogeneities in sediment physical and chemical properties both contributed to complex plume behavior and its persistence. This research underscores the great challenges in adequately characterizing this type of site to model the reactive transport processes over scales of 10 m or more. Our analysis provides important insights into the characterization, understanding, modeling, and remediation of groundwater contaminant plumes influenced by dynamic surface water and groundwater interactions.
Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios
Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.
Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles
In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.
BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale
Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.
Reactive-transport modeling of fracture flow to quantify the changes in flow pathways from matrix thermal contraction and mineral precipitation and dissolution; influence of grid resolution
A series of reactive-transport models of Enhanced Geothermal Systems (EGS) were constructed using the reactive transport code PFLOTRAN to examine the effect of matrix thermal contraction and mineral dissolution/precipitation on fracture flow in the context of grid cell size and model complexity. It was found that for thermal drawdown at production well, the impact of fracture zone grid cell size is negligible.
ATS (Advanced Terrestrial Simulator) integrated hydrology and reactive transport model output in Copper Creek, Colorado
This dataset is generated using the ATS (Advanced Terrestrial Simulator) model at Copper Creek, Colorado, the largest catchment in the East River watershed. ATS is an integrated hydrology and reactive transport model to simulate the Concentration-Discharge (C-Q) relationship, and is used to quantify the geochemical export and understand how watersheds respond to climate disturbances. Key hydrogeochemical variables, including discharge, precipitation, evapotranspiration, and geochemical concentrations, were simulated and output at a daily frequency at the Copper Creek outlet. A group of hydrology and geochemical model variables over the spatial domain are output at a bi-monthly frequency. The model used 1-km resolution daily frequency Daymet meteorological forcing, and mineral composition from previous studies. The simulations are performed on NERSC. The datasets include the input files for running the ATS (Advanced Terrestrial Simulator) simulation, including the mesh file, input script and geochemical model files. The output files from the ATS (Advanced Terrestrial Simulator) are the daily time-series discharge and geochemical concentration simulated at the outlet of Copper Creek, and the bi-monthly spatial output over the simulation domain.
Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching
Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.
2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results
This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.
Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications
This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.