Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility characterization”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

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

Unveiling a pervasive DNA adenine methylation regulatory network in the early-diverging fungus Rhizopus microsporus

Development of the DNA affinity purification and sequencing (DAP-seq) technique has allowed genome-scale studies of transcription factor (TF)-binding sites with high reproducibility. Here, we apply this technique to the human opportunistic pathogen Rhizopus microsporus, a mucoralean fungus belonging to the understudied group of early-diverging fungi. We characterize genome-wide binding sites of 58 TFs encoded by genes regulated through adenine methylation and representing major TF families. This analysis reveals their binding profiles and recognized sequences, expanding and diversifying the catalog of known fungal motifs. By integrating this data with DNA 6-methyladenine profiling, we uncover the extensive direct and indirect impact of this epigenetic modification on the regulation of gene expression. Furthermore, we use the generated data to identify TFs involved in biologically relevant processes such as zinc metabolism and light response. Our work enhances our understanding of regulatory mechanisms in R. microsporus and provides broader insights into gene regulation across the fungal kingdom.

Lax, Carlos [Universidad de Murcia (Spain)] (ORCID↗
DOE OSTI2020

Controlling T c of iridium films using the proximity effect

A superconducting Transition-Edge Sensor (TES) with low-$T_c$ is essential in a high resolution calorimetric detection. With a motivation of developing sensitive calorimeters for applications in cryogenic neutrinoless double beta decay searches, we have been investigating methods to reduce the $T_c$ of an Ir film down to 20 mK. Utilizing the proximity effect between a superconductor and a normal metal, we found two room temperature fabrication recipes of making Ir-based low-$T_c$ films. In the first approach, an Ir film sandwiched between two Au films, a Au/Ir/Au trilayer, has a tunable $T_c$ in the range of 20-100 mK depending on the relative thicknesses. In the second approach, a paramagnetic Pt thin film is used to create Ir/Pt bilayer with a tunable $T_c$ in the same range. We present detailed study of fabrication and characterization of Ir-based low-$T_c$ films, and compare the experimental results to theoretical models. We show that Ir-based films with predictable and reproducible critical temperature can be consistently fabricated for use in large scale detector applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Combined analysis of HPK 3.1 LGADs using a proton beam, beta source, and probe station towards establishing high volume quality control

The upgrades of the CMS and ATLAS experiments for the high luminosity phase of the Large Hadron Collider will employ precision timing detectors based on Low Gain Avalanche Detectors (LGADs). We present a suite of results combining measurements from the Fermilab Test Beam Facility, a beta source telescope, and a probe station, allowing full characterization of the HPK type 3.1 production of LGAD prototypes developed for these detectors. We demonstrate that the LGAD response to high energy test beam particles is accurately reproduced with a beta source. Finally, we further establish that probe station measurements of the gain implant accurately predict the particle response and operating parameters of each sensor, and conclude that the uniformity of the gain implant in this production is sufficient to produce full-sized sensors for the ATLAS and CMS timing detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Limits of Detection for EXAFS Characterization of Heterogeneous Single-Atom Catalysts

Single-atom catalysts (SACs), consisting of individual metal atoms dispersed on a support, attract attention due to their unique reactivity, efficient use of precious metals, and precise chemical tunability. Characterization of the metal species is crucial to substantiate structure–function relationships. Authors often use—and referees often require—X-ray absorption spectroscopy (XAS) data to prove the absence of clustered metal (or metal oxide) structures after pre-treatment and under in situ or operando conditions. However, there has been no critical assessment of the limitations of XAS in substantiating such conclusive statements, which is particularly important given the potential outsized influence of minority catalyst structures in dictating catalytic activity. Here, in this article, we quantitatively assess the detection limits of XAS to identify metal (or metal oxide) clusters in samples containing predominantly single atoms by modeling the extended X-ray absorption fine structure (EXAFS) of mixtures of structures. We identified that a significant fraction of clusters can coexist with SAC active sites (e.g., ~10% metallic Pt or ~40% oxidized Pt clusters in Pt/CeO 2 SACs), while eluding detection via EXAFS with any statistical significance. To generalize these conclusions, a descriptor-based screening of bulk metal oxides using a continuous Cauchy wavelet transform was proposed that suggests certain materials for which differentiating atomically dispersed metal species and metal oxide clusters would be infeasible by EXAFS (e.g., ReO x ). Based on this analysis, we suggest best practices for the study of SACs using EXAFS and provide recommendations to ensure that conclusions do not outpace the evidence used to support them. In this rapidly expanding research area, rigorous characterization will lead to greater understanding of the behavior of SACs and ultimately improved catalytic materials.

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

Successful νp -process in neutrino-driven outflows in core-collapse supernovae

The origin of the solar system abundances of several proton-rich isotopes, especially 92,94 Mo and 96,98 Ru, has been an enduring mystery in nuclear astrophysics. An attractive proposal to solve this problem is the νp-process, which can operate in neutrino-driven outflows in a core-collapse supernova after the shock is launched. Years of detailed studies, however, have cast doubt over the ability of this process to generate sufficiently high absolute and relative amounts of various p-nuclei. The νp-process is also thought to be excluded by arguments based on the long-lived radionuclide 92 Nb. Here, we present explicit calculations, in which both the abundance ratios and the absolute yields of the p-nuclei up to A ≲ 105 are successfully reproduced, even when using the modern (medium enhanced) triple-α reaction rates. The process is also shown to produce the necessary amounts of 92 Nb. The models are characterized by subsonic outflows and by the protoneutron star masses in the ≳ 1.7 M ⊙ range. In conclusion, this suggests that the Mo and Ru p-nuclides observed in the Solar System were made in CCSN explosions characterized by an extended accretion stage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetoelastic transition and magnetocaloric effect in induction melted Fe 100-x Rh x bulk alloys with x = 50, 51

Magnetoelastic transitions (METs) in bulk in nearly equiatomic Fe-Rh alloys produced by arc melting may show poor reproducibility related to insufficient chemical homogeneity and presence of impurity phases in variable concentrations. To better understand the synthesis conditions that reliably yield bulk FeRh materials with reproducible MET characteristics, Fe 100-x Rh x alloys with x = 50, 50.5 and 51 at. % were prepared by induction melting and thermal annealing under identical conditions. The fabricated samples were cut into several slices, followed by characterization of METs in each of the slices using isothermal and isofield magnetization measurements, differential scanning calorimetry, and direct measurements of the magnetocaloric effect. All of the slices exhibit METs between the AFM and FM states, but the transitions are abrupt with nearly the same change of magnetization, ΔM, when x = 50.5 and 51, whereas for the x = 50 alloy the transition spreads over a wide temperature interval and ΔM may fluctuate by as much as 10 % from one specimen to another. A comparison of the magnetocaloric responses of x = 50 and 51 materials is presented. The clearly different effect of the magnetic field on the transition in both directions leads to significant differences in the reversibility and maximum values of the magnetic field-induced entropy and adiabatic temperature changes, as well as average hysteresis losses. In terms of reproducibility, our results suggest that induction melting is a more appropriate technique to prepare these binary alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Activity Characterization for Modeling Behavioral-driven Human Mobility in Platial Networks

The population is increasingly becoming tractable as more and more people carry handheld devices as part of their everyday activities. Recent studies have shown that handheld devices' generated traffic share is now more than 50% of total global online traffic. This has created an unprecedented opportunity for modeling human mobility behavior. For example, aggregate check-ins and dwell time can reveal building level occupancies. However, there are clear limits to accurate modeling (e.g. reproducible, repeatable, and realistic), unless we decipher the underlying reason causing typical mobility patterns. We know that human behavior is a reflection of a set of activities, such as going to the gym or work, and which can be seen as a catalyst for humans to move from one location to another. This work envisions the use of activity characterization for modeling human mobility by introducing a context that maps activities to certain mobility patterns. In the end, we highlight the efficacy of the proposed approach by analyzing the impact of public policies surrounding stay at home order on human mobility.

Thakur, Gautam Malviya↗
DOE OSTI2025

Produced Water DNA Database (PW-DNA): Utilizing KBase to generate an environmental specific curated molecular database

The deep subsurface is estimated to host the majority of Earth’s microbial biomass yet remains one of the most challenging environments to access and study. One common approach to investigate these microbial communities is through the analysis of produced water from subsurface reservoirs, where researchers can assess water and gas chemistry along with molecular (DNA/RNA) sequence data. Advances in high-throughput sequencing have greatly expanded our understanding of these environments and their biotechnological potential. However, further progress requires large-scale, integrative meta-analyses across diverse datasets. To address this need, we developed the Produced Water-DNA (PW-DNA) Database, a curated, publicly available resource that consolidates microbial DNA/RNA sequences, geochemical data, and relevant metadata from in situ hydrocarbon environments such as coal beds, oil reservoirs, and natural gas systems. The PW-DNA database delivers three core benefits to the research community: (1) it improves data sharing by linking environmental microbial datasets with corresponding geochemical parameters, enabling more robust filtering and analysis; (2) it connects with complementary research databases to promote broader dissemination and interoperability; and (3) it supports technological innovation by serving as a resource for identifying microbial trends and exploring genetic potential. While individual studies have highlighted basin-specific microbial communities and functional redundancy in biogeochemical cycling, a comprehensive, system-wide perspective is needed to better understand connectivity and novelty across subsurface ecosystems. By designing the PW-DNA in the KBase platform, we provide a reproducible, visual framework for integrating large-scale genomic and geochemical data, enabling researchers to perform more informed analyses and experimental design. Ultimately, this resource enhances the ability to identify, characterize, and interpret microbial functions across diverse subsurface environments, thereby accelerating discovery in subsurface microbiology and biotechnology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Impact of inoculation practices on microbiota assembly and community stability in a fabricated ecosystem

Studying plant-microbe-soil interactions is challenging due to their high complexity and variability in natural ecosystems. While fabricated ecosystems provide opportunities to recapitulate aspects of these systems in reduced complexity and controlled environments, inoculation can be a significant source of variation. To tackle this, we evaluated how different bacteria inoculation practices and plant harvesting time points affect the reproducibility of a microbial synthetic community (SynCom) in association with the model grass Brachypodium distachyon. We tested three microbial inoculation practices: seed inoculation, transplant inoculation, and seedling inoculation; and two harvesting points: early (14-day-old plants) and late (21 days post-inoculation). We grew our plants and bacterial strains in sterile devices (EcoFABs) and characterized the microbial community from root, rhizosphere, and sand using 16S ribosomal RNA gene sequencing. The results showed that inoculation practices significantly affected the rhizosphere microbial community only when harvesting at an early time point but not at the late stage. As the SynCom showed a persistent association with B. distachyon at 21 days post-inoculation regardless of inoculation practices, we assessed the reproducibility of each inoculation method and found that transplant inoculation showed the highest reproducibility. Moreover, plant biomass was not adversely affected by transplant inoculation treatment. We concluded that bacteria inoculation while transplanting coupled with a later harvesting time point gives the most reproducible microbial community in the EcoFAB-B. distachyon-SynCom fabricated ecosystem and recommend this method as a standardized protocol for use with fabricated ecosystem experimental systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Top-down proteomics

Proteoforms arising from posttranslational modifications, genetic polymorphisms, and RNA splice variants, play a pivotal role as the key drivers in biology. Thus, a comprehensive understanding of proteoforms is essential for unraveling the intricacies of biological systems and bridging the gap between genotype and phenotype. By analyzing whole proteins without digestion, top-down proteomics (TDP) provides a holistic view of the proteome and presents a next-generation approach for deciphering protein function, uncovering disease mechanisms, and advancing precision medicine. This Primer embarks on a journey into the world of TDP by encapsulating its historical context, underlying principles, recent advances, and an outlook on the future of TDP. The experimental section navigates instrumentation, sample preparation, intact protein separation, tandem mass spectrometry techniques, and data collection. Results decipher raw data, visualize intact protein spectra, unravel data analysis, and explain proteoform identification, characterization, and quantitation, as well as statistical analysis. Various applications of TDP spanning the human proteoform project, biomedical, biopharmaceutical, and clinical applications are described. These are complemented by discussions on measurement reproducibility, limitations, and a forward-looking perspective outlining uncharted waters where the field can advance, and potential exciting future applications of TDP.

Roberts, David S.↗
DOE OSTI2024

Benchmarking machine learning interatomic potentials via phonon anharmonicity

Abstract Machine learning approaches have recently emerged as powerful tools to probe structure-property relationships in crystals and molecules. Specifically, machine learning interatomic potentials (MLIPs) can accurately reproduce first-principles data at a cost similar to that of conventional interatomic potential approaches. While MLIPs have been extensively tested across various classes of materials and molecules, a clear characterization of the anharmonic terms encoded in the MLIPs is lacking. Here, we benchmark popular MLIPs using the anharmonic vibrational Hamiltonian of ThO 2 in the fluorite crystal structure, which was constructed from density functional theory (DFT) using our highly accurate and efficient irreducible derivative methods. The anharmonic Hamiltonian was used to generate molecular dynamics (MD) trajectories, which were used to train three classes of MLIPs: Gaussian approximation potentials, artificial neural networks (ANN), and graph neural networks (GNN). The results were assessed by directly comparing phonons and their interactions, as well as phonon linewidths, phonon lineshifts, and thermal conductivity. The models were also trained on a DFT MD dataset, demonstrating good agreement up to fifth-order for the ANN and GNN. Our analysis demonstrates that MLIPs have great potential for accurately characterizing anharmonicity in materials systems at a fraction of the cost of conventional first principles-based approaches.

interatomic potentials↗
DOE OSTI2024

Enabling high-throughput enzyme discovery and engineering with a low-cost, robot-assisted pipeline

Abstract As genomic databases expand and artificial intelligence tools advance, there is a growing demand for efficient characterization of large numbers of proteins. To this end, here we describe a generalizable pipeline for high-throughput protein purification using small-scale expression in E. coli and an affordable liquid-handling robot. This low-cost platform enables the purification of 96 proteins in parallel with minimal waste and is scalable for processing hundreds of proteins weekly per user. We demonstrate the performance of this method with the expression and purification of the leading poly(ethylene terephthalate) hydrolases reported in the literature. Replicate experiments demonstrated reproducibility and enzyme purity and yields (up to 400 µg) sufficient for comprehensive analyses of both thermostability and activity, generating a standardized benchmark dataset for comparing these plastic-degrading enzymes. The cost-effectiveness and ease of implementation of this platform render it broadly applicable to diverse protein characterization challenges in the biological sciences.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improved Performance of Cu(InGa)(SeS) 2 PV Modules Using the Reaction of Metal Precursors. Final Report

This project “Improved Performance of Cu(InGa)(SeS) 2 PV Modules using the Reaction of Metal Precursors” was a partnership led by the Institute of Energy Conversion (IEC) at the University of Delaware with Columbia University and the Molecular Foundry at the Lawrence Berkeley National Laboratory. The aim was to develop pathways to improve Cu(InGa)(SeS) 2 (CIGSS) thin film photovoltaic modules using processes compatible with low manufacturing cost. The CIGSS approach investigated was a two-step process including deposition of metal precursor films following by reaction in hydride gases utilizing IEC’s novel reactor. The process was similar to that under commercial development by the project’s industry partner Stion. When Stion went out of business mid-project the focus changed to a rapid thermal process considered more commercially viable. Approaches to improve the performance of solar cells using the reacted films focused on two material innovations. First, the overall Ga content was increased to increase the operating voltage, which is desirable for scale-up to commercial modules. Second, the processing and performance advantages arising from Ag alloying were investigated. Advanced characterization guided process and material development including control of relative composition gradients. Research on the formation of Cu-Ga-In metal precursors utilized sputtering deposition which is normally used in commercial applications. The work resulted in processes for deposition of precursor stacks with increased relative Ga content and effects of deposition parameters on morphology and phase composition were established. It was shown that the metal precursor films have comparable phase composition and morphology so subsequent reaction follows from the same starting point. The addition of Ag to the metal precursors gave more uniform morphology and improved adhesion of reacted films which enable higher reaction temperature for faster processing. A novel outcome was the discovery of a previously undocumented material phase in sputter-deposited and evaporated Ag-Cu-In-Ga thin films. Hydride gas reaction processes including time-temperature-concentration profiles were developed for different precursor compositions. This enables control of composition profiles to engineer through-film gradients for solar cell optimization with characterization and simulations used to correlate measured film composition profiles to measurements of devices. In particular, the gradient of sulfur at the front of the CIGSS film was found to be critical. The simulations guided process development leading to improved reproducibility of devices improved performance with higher Ga content and higher voltage. With Ag-alloyed precursors, the reaction pathways leading were determined. A significant finding was that Ag-alloying increases the reaction rate to completely convert precursor films to the final chalcopyrite which could enable reduced reaction time to benefit manufacturability. To maintain potential commercial viability, the process under investigation was refocused to a rapid thermal process that could potentially be incorporated into an in-line process for manufacturing. Precursors with different composition were capped with an extra selenium layer and reacted in hydrogen sulfide 5-15 minutes, compared to typically 2 hours in the previous multi-step batch process. Critical RTP parameters were identified to control the reaction. Further optimization would be needed for high efficiency solar cells but pathways to high quality devices with further optimization and improved heating uniformity were developed. The project also developed new optoelectronic characterization approaches with a focus on development and application of spatial- and time-resolved photoluminescence and a custom mapping photoluminescence microscope built. It was shown how critical electronic transport properties strongly depend on the chemical composition of the material and that a wide range of samples show inhomogeneity on a length scale larger than the grains in the films. Additionally, two-photon excitation capability was developed to distinguish bulk vs surface losses. The project advances the state-of-the -art for precursor reaction processes in several ways that could impact manufacturing. This includes validation of approaches to increase voltage and establishment of model-guided control to form optimal composition profiles. The application of process control approaches with knowledge of phase formation and reaction pathways can be critically valuable in designing a large-scale process.

14 SOLAR ENERGY↗
DOE OSTI2023

Origins of Conformational Heterogeneity in Peptoid Helices Formed by Chiral N -1-Phenylethyl Sidechains

N-substituted glycines (polypeptoids) containing chiral hydrophobic sidechains are known to fold into biomimetic alpha helices. These helix formers often produce conformationally heterogeneous structures and are difficult to characterize at a sub-nanometer resolution. Previously, peptoid N-1-phenylethyl (S)-enantiomer sidechains (Nspe) were inferred from various experiments to form right-handed helices and (R)-enantiomers (Nrpe), left-handed helices. Prior computational work for N(s/r)pe oligomers has struggled to reproduce this trend. Herein, quantum mechanics calculations and molecular dynamics simulations are used to understand the origins of this discrepancy. Results from DFT and molecular mechanics calculations on a variety of Nspe and Nrpe oligomers as a function of chain length are in agreement, showing that Nspe and Nrpe prefer left- and right-handed helices, respectively. Additional metadynamics simulations are used to study Nrpe and Nspe oligomers folding in water. These results show that the free-energy driving forces for assembly into a helical backbone configuration are very small (within ~k B T). Lastly, we compare DFT calculations for other experimentally characterized peptoid sidechains, N(r/s)sb, N(r/s)tbe, and N(r/s)npe. In this analysis, we show that peptoid sidechains determined to be more robust experimentally (tbe and npe) have helical preferences opposite the trend seen in less robust assemblies formed by N(r/s)pe and N(r/s)sb chemistries. Furthermore, the more robust tbe and nnpe favor the (S)-enantiomer to right-handed and the (R)-enantiomers to left-handed helices.

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

Optimal adjustment sets for causal query estimation in partially observed biomolecular networks

Abstract Causal query estimation in biomolecular networks commonly selects a ‘valid adjustment set’, i.e. a subset of network variables that eliminates the bias of the estimator. A same query may have multiple valid adjustment sets, each with a different variance. When networks are partially observed, current methods use graph-based criteria to find an adjustment set that minimizes asymptotic variance. Unfortunately, many models that share the same graph topology, and therefore same functional dependencies, may differ in the processes that generate the observational data. In these cases, the topology-based criteria fail to distinguish the variances of the adjustment sets. This deficiency can lead to sub-optimal adjustment sets, and to miss-characterization of the effect of the intervention. We propose an approach for deriving ‘optimal adjustment sets’ that takes into account the nature of the data, bias and finite-sample variance of the estimator, and cost. It empirically learns the data generating processes from historical experimental data, and characterizes the properties of the estimators by simulation. We demonstrate the utility of the proposed approach in four biomolecular Case studies with different topologies and different data generation processes. The implementation and reproducible Case studies are at https://github.com/srtaheri/OptimalAdjustmentSet.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

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

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

54 ENVIRONMENTAL SCIENCES↗