Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process coupling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Tellurium Goes for a Ride on the “Ferrous” Wheel: Interactions of Te(VI) and Te(IV) with Fe(II)-Bearing Minerals

Compared to more abundant metalloids (e.g., As, Se, and Sb), little is known regarding the biogeochemistry of tellurium (Te), a critical metal whose use in CdTe photovoltaic solar cells is driving increasing global demand for this element. In this study, to understand the redox behavior of Te under ferrugenic/sulfidogenic conditions, we examined the potential for reduction of Te(VI) and Te(IV) in aqueous suspensions containing Fe(II) as siderite, vivianite, green rust, magnetite, or mackinawite; Fe(II)-bearing minerals which are often formed during microbial reduction of Fe(III) oxides. In the mackinawite systems, complete reduction of Te(VI) or Te(IV) to Te(0) was observed within 12 h, and near complete reduction (>90% Te(0)) was observed within 2.5 days in the green rust systems. In the presence of siderite, complete reduction of Te(VI) or Te(IV) to Te(0) occurred within 60 days. We observed >60% reduction of Te(VI) to Te(IV) within 2.5 days in the magnetite system, but Te(0) did not form until 120 days; however, complete reduction to Te(0) was observed within 120 days when starting with Te(IV). With vivianite there was >80% reduction of Te(VI) to Te(IV) within 12 h, without further reduction to Te(0) over the duration of the experiment (120 days); in the Te(IV)-vivianite system, there was no evidence of Te(IV) reduction within 120 days. The reduction of Te(VI) and Te(IV) in soils and sediments has been largely attributed to direct reduction by microbes; however, the reduction of Te(VI) and Te(IV) by Fe(II)-bearing minerals suggests that abiotic or coupled biotic-abiotic processes may also play a critical role in Te redox chemistry in natural and engineered environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Novel Cell-Type-Specific Drought-Responsive Proteins in Root Tips of Field-Grown Perennial Switchgrass

The root-tip region of plants, including the root cap, forms the most basal terminal of the root and exhibits a high degree of cellular complexity in terms of morphology, cytological function, and interaction with environmental cues in the soil. Cells in this region follow a developmental trajectory, transitioning from stem cells to meristematic cells, and ultimately to fully differentiated cell types. However, our understanding of root-tip cell-type specific proteomic responses to abiotic stresses, such as drought, particularly under field conditions, remains limited. This study aimed to identify spatially resolved, cell type-specific proteomes in switchgrass (Panicum virgatum) root tips under drought stress. Root tips were collected from seven-year-old, field-grown switchgrass ‘Alamo’ plants excavated under both well-watered and long-term drought conditions. Cell type-specific proteins were identified using laser capture microdissection (LCM) coupled with nanoPOTS (Nanodroplet Processing in One Pot for Trace Samples) and nano-LC-MS proteomics analysis. Five distinct cell types were targeted: (1) cells in the quiescent center and stem cell niche (QuC), (2) protodermal epidermal cells (PEC) in the meristematic zone, (3) epidermal cells in the transition and elongation zones above the root cap (Epi), (4) peripheral root cap cells (PRC), forming 2–3 layers below the PEC and 1–2 layers above the root border cells, and (5) columella root cap cells (Col) comprising of the columella initials and a single underlying layer of cells undergoing active growth. Principal component analysis (PCA) revealed clear separation among the five targeted cell types, confirming distinct proteomic profiles. Proteins predominantly enriched in each cell type were linked to distinct cellular functions, with QuC cells showing involvement in chromosomal behavior, DNA replication, and mitosis—key processes for stem cell niche regulation. Drought stress resulted in alterations of proteostasis, as evidenced by significant decreases in ribosomal proteins and increases in protein synthesis inhibitors. Moreover, drought stress induced unique cell-type–specific proteins involved in phytohormone biosynthesis and signaling pathways, including auxin, cytokinin, and jasmonic acid. In particular, QuC cells were more highly enriched in proteins associated with DNA repair and mitotic processes. Metabolic pathways related to amino acids, carbohydrates, and lipids were differentially affected in a cell-type–dependent manner, whereas general stress-responsive proteins exhibited consistent changes across all five cell types. Overall, this study provides unique spatially resolved, cell-type-specific proteomic profiles in root tips, representing a significant advancement in our understanding of the cellular mechanisms underlying plant responses to drought stress in natural field conditions.

perennial grass↗
DOE OSTI2020

In Situ Identification of Reaction Intermediates and Mechanistic Understandings of Methane Oxidation over Hematite: A Combined Experimental and Theoretical Study

Effective methane utilization for either clean power generation or value-added chemical production has been a subject of growing attention worldwide for decades, yet challenges persist mostly in relation to methane activation under mild conditions. Here, we report hematite, an earth-abundant material, to be highly effective and thermally stable to catalyze methane combustion at low temperatures (<500 °C) with a low light-off temperature of 230 °C and 100% selectivity to CO 2 . The reported performance is impressive and comparable to those of precious-metal-based catalysts, with a low apparent activation energy of 17.60 kcal·mol –1 . Our theoretical analysis shows that the excellent performance stems from a tetra-iron center with an antiferromagnetically coupled iron dimer on the hematite (110) surface, analogous to that of the methanotroph enzyme methane monooxygenase that activates methane at ambient conditions in nature. Isotopic oxygen tracer experiments support a Mars van Krevelen redox mechanism where CH 4 is activated by reaction with a hematite surface oxygen first, followed by a catalytic cycle through a molecular-dioxygen-assisted pathway. Surface studies with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and density functional theory (DFT) calculations reveal the evolution of reaction intermediates from a methoxy CH 3 –O–Fe, to a bridging bidentate formate b-HCOO–Fe, to a monodentate formate m-HCOO–Fe, before CO 2 is eventually formed via a combination of thermal hydrogen-atom transfer (HAT) and proton-coupled electron transfer (PCET) processes. Finally, the elucidation of the reaction mechanism and the intermediate evolutionary profile may allow future development of catalytic syntheses of oxygenated products from CH 4 in gas-phase heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Synthesis, Structure, and Redox Reactivity of Ni Complexes Bearing a Redox and Acid–Base Non-innocent Ligand with Ni II , Ni III , and Ni IV Formal Oxidation States

Here, a series of Ni complexes bearing a redox and acid–base noninnocent tetraamido macrocyclic ligand, H 4 -(TAML-4) {H 4 -(TAML-4) = 15,15-dimethyl-5,8,13,17-tetrahydro-5,8,13,17-tetraaza-dibenzo[a,g]cyclotridecene-6,7,14,16-tetraone}, with formal oxidation states of Ni II , Ni III , and Ni IV were synthesized and characterized structurally and spectroscopically. The X-ray crystallographic analysis of the Ni complexes revealed a square planar geometry, and the [Ni(TAML-4)] complex with the formal oxidation state of Ni IV was characterized to be [Ni III (TAML-4 •+ )] with the oxidation state of the Ni III ion and the one-electron oxidized TAML-4 ligand, TAML-4 •+ . The Ni III oxidation state and the TAML-4 radical cation ligand, TAML-4 •+ , were supported by X-ray absorption spectroscopy and density functional theory calculations. The reversible interconversions between [Ni II (TAML-4)] 2– and [Ni III (TAML-4)] − and between [Ni III (TAML-4)] − and [Ni III (TAML-4 •+ )] were demonstrated in spectroelectrochemical measurements as well as in chemical oxidation and reduction reactions. The reactivities of [Ni III (TAML-4)] − and [Ni III (TAML-4 •+ )] were then investigated in hydride transfer reactions using NADH analogs. Hydride transfer from 9,10-dihydro-10-methylacridine (AcrH 2 ) to [Ni III (TAML-4 •+ )] was found to proceed via electron transfer (ET) from AcrH 2 to [Ni III (TAML-4 •+ )] with no deuterium kinetic isotope effect (k H /k D = 1.0(2)). In contrast, hydride transfer from AcrH 2 to [NiIII(TAML-4)] − proceeded much more slowly via a concerted proton-coupled electron transfer (PCET) process with k H /k D = 7.0(5). In the latter reaction, an electron and a proton were transferred to the Ni III center and the TAML-4 ligand, respectively. The mechanisms of the ET by [Ni III (TAML-4 •+ )] and the concerted PCET by [Ni III (TAML-4)] − were ascribed to the different redox potentials of the Ni complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Seasonal and Morphological Controls on Nitrate Retention in Arctic Deltas

Abstract Estimates of nitrate loading to the Arctic Ocean are limited by the lack of field observations within deltas partly due to logistical constraints. To overcome this limitation, we use a remote sensing framework to estimate retention of nitrate in Arctic deltas. We achieve this by coupling hydrological and biogeochemical process models at the network scale for five major Arctic deltas. Binary masks of delta channels were used to simulate flow direction and magnitude through networks. Models were parameterized using historical and seasonal observations. Simulated nitrate retention ranged from 2.9% to 15% of the incoming load. Retention rates were largest during winter but smallest during spring conditions when increased discharges export large nitrate masses to the coast. Under future climate scenarios, retention rates fall by ∼1%–10%. Arctic deltas have an important effect on the magnitude of nitrate entering Arctic seas and the inclusion of processing in deltas can improve flux estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observations of Offshore Internal Boundary Layers

The growth of the marine internal boundary layer (MIBL, height h i ) with the shore-normal distance x, is a topic of continuing interest because of its applications in coastal dispersion, offshore wind farm siting, coastal air-sea fluxes and in evaporative ducting. Available data on MIBL are only scarce, given its difficult to measure the variability of coastal winds. During Coupled Air-Sea Processes and Electromagnetic Research (CASPER) campaigns an array of instrumentation was deployed to measure offshore spatial variability and its effect on electromagnetic wave propagation. Meteorological sensors (flux towers and remote sensing) were deployed along the coast of Point Mugu, California, on a research vessel and FLoating Instrument Platform (FLIP) provided surface layer and boundary layer observations. In this article, measurements from multiple remote sensing instruments, such as synchronized triple Doppler lidars, small boat operations with tethered lifting system, and radiosondes provide a holistic view of the MIBL growth and its spatial variability from the coast. A convective and stable MIBL observed during two intensive operating period days showed distinct growth characteristics off the coast of Point-Mugu. During stable stratified atmospheric conditions, an MIBL was observed at least up to 47 km from the coast. The growth of the MIBL within the nearshore adjustment zone was influenced by surrounding atmospheric, oceanographic, and topographic conditions. A new parameterization scheme is developed based on advection-diffusion balance equations accounting for upstream turbulence and compared with h i observations from Doppler lidar and profiles from small boat. An evaluation of existing IBL theories was also conducted.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Monitoring Preferential Flow of Water in Sand Using Thermoacoustics Wave Imaging

Abstract Accurate predictions of fluid flow, mass transport, and reaction rates critically impact the efficiency and reliability of subsurface exploration and sustainable use of subsurface resources. Quantitative dynamical sensing and imaging can play a pivotal role in the ability to make such predictions. Geophysical thermoacoustic technology has the potential to provide the aforementioned capabilities, since it builds upon the principle that electromagnetic and mechanical wave fields can be coupled through a thermodynamic process. In this letter, we present laboratory experiments featuring the efficacy of thermoacoustic imaging in the monitoring of preferential flow of water in porous media. Our laboratory experimental equipment can be readily packaged in a form factor that fits in a borehole, and the use of multiple acoustic transducers—which can be combined with volumetric coding techniques—has the potential to provide quasi‐real‐time imaging (0.5 Hz video rate) of regions in close proximity (a few meters) of an open field well.

Geology↗
DOE OSTI2024

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The environmental controls on efficiency of enhanced rock weathering in soils

Abstract Enhanced rock weathering (ERW) in soils is a promising carbon removal technology, but the realistically achievable efficiency, controlled primarily by in situ weathering rates of the applied rocks, is highly uncertain. Here we explored the impacts of coupled biogeochemical and transport processes and a set of primary environmental and operational controls, using forsterite as a proxy mineral in soils and a multiphase multi-component reactive transport model considering microbe-mediated reactions. For a onetime forsterite application of ~ 16 kg/m 2 , complete weathering within five years can be achieved, giving an equivalent carbon removal rate of ~ 2.3 kgCO 2 /m 2 /yr. However, the rate is highly variable based on site-specific conditions. We showed that the in situ weathering rate can be enhanced by conditions and operations that maintain high CO 2 availability via effective transport of atmospheric CO 2 (e.g. in well-drained soils) and/or sufficient biogenic CO 2 supply (e.g. stimulated plant–microbe processes). Our results further highlight that the effect of increasing surface area on weathering rate can be significant—so that the energy penalty of reducing the grain size may be justified—only when CO 2 supply is nonlimiting. Therefore, for ERW practices to be effective, siting and engineering design (e.g. optimal grain size) need to be co-optimized.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Rigid cationic ligands enable high-efficiency NIR-II photoluminescence in copper( i ) iodide hybrid semiconductors

Near-infrared (NIR) luminescent materials are pivotal for advanced optoelectronic and biomedical applications, yet attaining efficient emission in the NIR-II region (950-1400 nm) remains challenging. Here, we introduce a ligand cationization strategy for designing copper(i) iodide-organic hybrid materials that emit in the NIR-II region (920-1120 nm) with PLQYs up to 8.58%. By incorporating rigid cationic ligands with CuI modules, we synergistically achieve bandgap narrowing (to 1.51 eV) and structural rigidification via ionic-dative bonding, effectively suppressing non-radiative decay while extending emission beyond 1100 nm. Coupled with solution processability-enabled by the successful synthesis of nanometer-sized nanoparticles in various shapes-and excellent thermal stability (≥210 °C), this work establishes ligand cationization as a universal approach for designing efficient NIR-II emitters.

Chen, Jingwen↗
DOE OSTI2020

Virtual Hawking Radiation

We consider the effects of off-shell Hawking radiation on scattering processes involving black holes coupled to quantum fields. The focus here is to the case of gravitational scattering of a scalar field mediated by the exchange of virtual Hawking gravitons from a four-dimensional Schwarzschild black hole. Our result is obtained in the context of a worldline effective field theory for the black hole, and is valid in the semi-classical limit where the Schwarzschild radius rs is larger than the Planck length 1/m Pl . In addition, we assume that four-momentum exchange q is smaller than r s –1 and that the incoming particle has energy larger then the black hole's Hawking temperature. The inelastic cross section we obtain is a new, leading order quantum gravity effect, arising at the same order in q 2 /$m$ $^{2}_{Pl}$ as the well understood one-loop graviton vacuum polarization corrections to gravitational scattering between massive particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗
DOE OSTI2022

Theoretical Analysis of Critical Conditions for Crack Formation and Propagation, and Optimal Operation of SOECs

A theoretical analysis on crack formation and propagation was performed based on the coupling between the electrochemical process, classical elasticity, and fracture mechanics. The chemical potential of oxygen, thus oxygen partial pressure, at the oxygen electrode-electrolyte interface ( μ O 2 OE∣El ) was investigated as a function of transport properties, electrolyte thickness and operating conditions (e.g., steam concentration, constant current, and constant voltage). Our analysis shows that: a lower ionic area specific resistance (ASR), r i O E , and a higher electronic ASR ( r e O E ) of the oxygen electrode/electrolyte interface are in favor of suppressing crack formation. The μ O 2 O E ∣ E l , thus local pO 2 , are sensitive towards the operating parameters under galvanostatic or potentiostatic electrolysis. Constant current density electrolysis provides better robustness, especially at a high current density with a high steam content. While constant voltage electrolysis leads to greater variations of μ O 2 O E ∣ E l . Constant current electrolysis, however, is not suitable for an unstable oxygen electrode because μ O 2 O E ∣ E l can reach a very high value with a gradually increased r i O E . A crack may only occur under certain conditions when p O 2 T P B > p c r .

42 ENGINEERING↗
DOE OSTI2026

HIGRAD/FIRETEC

This is the public-facing repository for the LANL-developed and maintained HIGRAD/FIRETEC software, a CFD software designed to simulate the coupled fire-atmospheric processes at a landscape-scale of resolution. Distribution is possible under OSS license O5086.

Josephson, Alexander [Los Alamos National Laborato↗
DOE OSTI2025

Redox poise in R. rubrum phototrophic growth drives large-scale changes in macromolecular pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H 2 ), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H 2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Subsurface microbial community structure shifts along the geological features of the Central American Volcanic Arc

Subduction of the Cocos and Nazca oceanic plates beneath the Caribbean plate drives the upward movement of deep fluids enriched in carbon, nitrogen, sulfur, and iron along the Central American Volcanic Arc (CAVA). These compounds fuel diverse subsurface microbial communities that in turn alter the distribution, redox state, and isotopic composition of these compounds. Microbial community structure and functions vary according to deep fluid delivery across the arc, but less is known about how microbial communities differ along the axis of a convergent margin as geological features (e.g., extent of volcanism and subduction geometry) shift. Here, we investigate changes in bacterial 16S rRNA gene amplicons and geochemical analysis of deeply-sourced seeps along the southern CAVA, where subduction of the Cocos Ridge alters the geological setting. We find shifts in community composition along the convergent margin, with communities in similar geological settings clustering together independently of the proximity of sample sites. Microbial community composition correlates with geological variables such as host rock type, maturity of hydrothermal fluid and slab depth along different segments of the CAVA. This reveals tight coupling between deep Earth processes and subsurface microbial activity, controlling community distribution, structure and composition along a convergent margin.

Science & Technology - Other Topics↗
DOE OSTI2025

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.

54 ENVIRONMENTAL SCIENCES↗