Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth profiling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Phloem Exudate Protein Profiles during Drought and Recovery Reveal Abiotic Stress Responses in Tomato Vasculature

Drought is the leading cause of agricultural yield loss among all abiotic stresses, and the link between water deficit and phloem protein contents is relatively unexplored. Here we collected phloem exudates from Solanum lycopersicum leaves during periods of drought stress and recovery. Our analysis identified 2558 proteins, the most abundant of which were previously localized to the phloem. Independent of drought, enrichment analysis of the total phloem exudate protein profiles from all samples suggests that the protein content of phloem sap is complex, and includes proteins that function in chaperone systems, branched-chain amino acid synthesis, trehalose metabolism, and RNA silencing. We observed 169 proteins whose abundance changed significantly within the phloem sap, either during drought or recovery. Proteins that became significantly more abundant during drought include members of lipid metabolism, chaperone-mediated protein folding, carboxylic acid metabolism, abscisic acid signaling, cytokinin biosynthesis, and amino acid metabolism. Conversely, proteins involved in lipid signaling, sphingolipid metabolism, cell wall organization, carbohydrate metabolism, and a mitogen-activated protein kinase are decreased during drought. Our experiment has achieved an in-depth profiling of phloem sap protein contents during drought stress and recovery that supports previous findings and provides new evidence that multiple biological processes are involved in drought adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Observation of Quantum Anomalous Hall Effect and Exchange Interaction in Topological Insulator/Antiferromagnet Heterostructure

Integration of a quantum anomalous Hall insulator with a magnetically ordered material provides an additional degree of freedom through which the resulting exotic quantum states can be controlled. Here, an experimental observation is reported of the quantum anomalous Hall e?ect in a magnet-ically-doped topological insulator grown on the antiferromagnetic insulator Cr?O?. The exchange coupling between the two materials is investigated using ?eld-cooling-dependent magnetometry and polarized neutron re?ec-tometry. Both techniques reveal strong interfacial interaction between the antiferromagnetic order of the Cr?O? and the magnetic topological insulator, manifested as an exchange bias when the sample is ?eld-cooled under an out-of-plane magnetic ?eld, and an exchange spring-like magnetic depth pro?le when the system is magnetized within the ?lm plane. These results identify antiferromagnetic insulators as suitable candidates for the manipula-tion of magnetic and topological order in topological insulator ?lms.

Pan, Lei↗
DOE OSTI2020

The Effects of Contaminant-Aging on Decontamination Efficacy for Rapid Remediation of Concrete Surfaces - 20412

Understanding the relationship between decontamination efficacy and contaminant aging is imperative for developing effective remediation strategies following a large-scale contamination event. Non-destructive decontamination of concrete following a contamination incident becomes more difficult with increased time between contamination and decontamination because contaminants may chemically bind to concrete and/or penetrate into the concrete subsurface. In this work, we evaluated the decontamination efficacy of two decontamination methods on concrete samples contaminated with soluble Cesium-137 (Cs-137) and silica particles with 0.5 μm and 2 μm diameters. Concrete samples were aged between 1 and 59 days, with half of the coupons receiving 1 mL of artificial rainwater about once every three days. After aging, coupons were either decontaminated or analyzed to determine contaminant penetration depths. Coupons were decontaminated with 0.1 M Potassium Chloride (KCI) solution applied either by flowing solution across the contaminated coupon face or non-destructive power washing. Contaminant depth profiles were created by removing the top surface of the coupon at least 15 times, measuring the activity in each removed layer, normalizing the total activity removed to the measured decontamination efficacy following the last layer removal, and determining the penetration thickness using the coupon dimensions, bulk density, and mass of material removed for each layer. Pressurized washing effectively removed particles, but not soluble cesium, despite the majority of Cs-137 being located within the first millimeter of the subsurface. Soluble cesium removal by flowing 0.1 M KCI solution across the contaminated surface dropped to less than 12% within the first ten days. This drop in Cs-137 removal efficacy occurred within the first five days for coupons experiencing artificial rainfall events. Flowing 0.1 M KCI solution across the contaminated face was ineffective at removing particles, likely because particles had localized in surface depressions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Topological Antiferromagnetic Van der Waals Phase in Topological Insulator/Ferromagnet Heterostructures Synthesized by a CMOS-Compatible Sputtering Technique

Breaking time-reversal symmetry by introducing magnetic order, thereby opening a gap in the topological surface state bands, is essential for realizing useful topological properties such as the quantum anomalous Hall and axion insulator states. In this work, a novel topological antiferromagnetic (AFM) phase is created at the interface of a sputtered, c-axis-oriented, topological insulator/ferromagnet heterostructure—Bi 2 Te 3 /Ni 80 Fe 20 because of diffusion of Ni in Bi 2 Te 3 (Ni-Bi 2 Te 3 ). The AFM property of the Ni-Bi2Te3 interfacial layer is established by observation of spontaneous exchange bias in the magnetic hysteresis loop and compensated moments in the depth profile of the magnetization using polarized neutron reflectometry. Analysis of the structural and chemical properties of the Ni-Bi2Te3 layer is carried out using selected-area electron diffraction, electron energy loss spectroscopy, and X-ray photoelectron spectroscopy. These studies, in parallel with first-principles calculations, indicate a solid-state chemical reaction that leads to the formation of Ni=Te bonds and the presence of topological antiferromagnetic (AFM) compound NiBi 2 Te 4 in the Ni-Bi 2 Te 3 interface layer. The Neél temperature of the Ni-Bi 2 Te 3 layer is ≈ 63 K, which is higher than that of typical magnetic topological insulators (MTIs). The presented results provide a pathway toward industrial complementary metal-oxide-semiconductor (CMOS)-process-compatible sputtered-MTI heterostructures, leading to novel materials for topological quantum devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unraveling the Stable Cathode Electrolyte Interface in all Solid‐State Thin‐Film Battery Operating at 5 V

Abstract Spinel‐type LiNi 0.5 Mn 1.5 O 4 (LNMO) is one of the most promising 5 V‐class cathode materials for Li‐ion batteries that can achieve high energy density and low production costs. However, in liquid electrolyte cells, the high voltage causes continuous cell degradation through the oxidative decomposition of carbonate‐based liquid electrolytes. In contrast, some solid‐state electrolytes have a wide electrochemical stability range and can withstand the required oxidative potential. In this work, a thin‐film battery consisting of an LNMO cathode with a solid lithium phosphorus oxynitride (LiPON) electrolyte is tested and their interface before and after cycling is characterized. With Li metal as the anode, this system can deliver stable performance for 600 cycles with an average Coulombic efficiency >99%. Neutron depth profiling indicates a slight overlithiated layer at the interface prior to cycling, a result that is consistent with the excess charge capacity measured during the first cycle. Cryogenic electron microscopy further reveals intimate contact between LNMO and LiPON without noticeable structure and chemical composition evolution after extended cycling, demonstrating the superior stability of LiPON against a high voltage cathode. Consequently, design guidelines are proposed for interface engineering that can accelerate the commercialization of a high voltage cell with solid or liquid electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2026

Combining ToF‐SIMS and Multivariate Analysis to Resolve Active Sites on Ni‐Based HER Catalysts

Unambiguous identification of active sites in heterogeneous catalysis remains a major challenge, particularly for materials with ultrathin, chemically mixed surface layers. Here, we demonstrate a generalizable approach that combines time-of-flight secondary ion mass spectrometry (ToF-SIMS) with multivariate statistical analysis (principal component analysis [PCA] and multivariate curve resolution [MCR]) to resolve catalytically relevant motifs at the nanoscale. Using Ni electrodes as a model system, PCA distinguished hydroxide-enriched domains from oxide- and metal-rich regions, while MCR decomposed depth profiles and 3D images into hydroxide, oxide, and metallic layers with nanometer resolution. A unique secondary-ion fragment, NiO 3 H 3 − (m/z 108.94), emerged as a marker of hydroxide-rich environments and correlated with hydrogen evolution reaction (HER) activity across a series of Ni electrodes. Complementary density functional theory (DFT) calculations revealed that Ni(OH) 2 clusters adjacent to metallic Ni offer the most favorable water dissociation energetics, establishing the structural origin of the marker. While illustrated here for Ni-based HER, this workflow provides a broadly applicable framework to isolate and rank near-surface patterns that govern catalytic activity, thereby extending ToF-SIMS from a qualitative probe to a predictive tool for active site identification.

HER active sites↗
DOE OSTI2021

Hydraulic redistribution by deeply rooted grasses and its ecohydrologic implications in the southern Great Plains of North America

Perennial bioenergy crops with deep (>1 m) rooting systems, such as switchgrass (Panicum virgatum L.), are hypothesized to increase carbon storage in deep soil. Deeply rooted plants may also affect soil hydrology by accessing deep soil water for transpiration, which can affect soil water content and infiltration in deep soil layers, thereby affecting groundwater recharge. Using stable H and O isotope (δ 2 H and δ 18 O) and 3 H values, we studied the soil water conditions at 20–30 cm intervals to depths of 2.4–3.6 m in paired fields of switchgrass and shallow rooted crops at three sites in the southern Great Plains of North America. We found that soil under switchgrass had consistently higher soil water content than nearby soil under shallow-rooted annual crops by a margin of 15%–100%. Soil water content and isotopic depth profiles indicated that hydraulic redistribution of deep soil water by switchgrass roots explained these observed soil water differences. To our knowledge, these are the first observations of hydraulic redistribution in deeply rooted grasses, and complement earlier observations of dynamic soil water fluxes under shallow-rooted grasses. Hydraulic redistribution by switchgrass may be a strategy for drought avoidance, wherein the plant may actively prevent water limitation. Furthermore, this raises the possibility that deeply rooted grasses may be used to passively subsidize soil water to more shallow-rooted species in inter-cropping arrangements.

3H↗
DOE OSTI2020

Analysis of the Stable Interphase Responsible for the Excellent Electrochemical Performance of Graphite Electrodes in Sodium-Ion Batteries

Considerable efforts have been exerted to understand the formation and properties of the solid electrolyte interphase (SEI) in sodium ion batteries. However, the puzzling existence and role of SEI behind the huge volume changes of the graphite electrodes need to be answered. Herein, the reason of how ether-derived SEI maintains excellent reversibility despite the huge volume changes during cycling is unraveled. Theoretical simulations and Fourier-transform infrared spectroscopy demonstrate the formation mechanism of an SEI between the graphite anode and electrolyte. Furthermore, the high mechanical tolerance of the ether-derived SEI is confirmed in atomic force microscopy. A depth profile of X-ray photoelectron spectroscopy points to a multilayer structure of the ether-derived SEI. The outer layer comprises organics (sodium alkoxide), while the inorganics (Na 2 CO 3 , NaF) in interior region are mixed with some organics. Notably, the presence of organics ensures the adaptability of the SEI to the volume expansion of graphite during cycling, and the concentrated distribution of inorganics improves the Young's modulus (resistance to deformation). Therefore, the graphite anode exhibits high cycle stability (96.6% capacity retention ratio at 1 A g –1 over 860 cycles) and efficiency (≈99.5%).

25 ENERGY STORAGE↗
DOE OSTI2021

Migration Kinetics of Surface Ions in Oxygen-Deficient Perovskite During Topotactic Transitions

Oxygen diffusivity and surface exchange kinetics underpin the ionic, electronic, and catalytic functionalities of complex multivalent oxides. Towards understanding and controlling the kinetics of oxygen transport in emerging technologies, it is highly desirable to reveal the underlying lattice dynamics and ionic activities related to oxygen variation. In this study, the evolution of oxygen content is identified in real-time during the progress of a topotactic phase transition in La 0.7 Sr 0.3 MnO 3-δ epitaxial thin films, both at the surface and throughout the bulk. Using polarized neutron reflectometry, a quantitative depth profile of the oxygen content gradient is achieved, which, alongside atomic-resolution scanning transmission electron microscopy, uniquely reveals the formation of a novel structural phase near the surface. Surface-sensitive X-ray spectroscopies further confirm a significant change of the electronic structure accompanying the transition. The anisotropic features of this novel phase enable a distinct oxygen diffusion pathway in contrast to conventional observation of oxygen motion at moderate temperatures. The results provide insights furthering the design of solid oxygen ion conductors within the framework of topotactic phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impact of Salt Chemistry and Interlayer Architecture on Solid Electrolyte Interphase Composition and Performance

Reversible lithium metal anodes require interphases that are chemically stable, mechanically robust, and compatible with scalable processing. Here, we use slot-die-coated, salt-rich PEO interlayers to design SEI nanolayers on Cu current collectors and isolate how salt anion chemistry controls interphase composition and lithium plating/stripping behavior in architectures relevant to Li plating and anode-free cells. Mesoscale modeling shows that dendrite initiation at SEI defects is governed by a tradeoff among interphase thickness, ionic conductivity, and stiffness through defect-driven current focusing, motivating controlled formation from thin precursor layers. Using UHV-transferred XPS and ToF-SIMS depth profiling after controlled Li deposition, we find that fluorinated sulfonylimide salts (LiTFSI, LiFSI) form mixed inorganic/organic interphases containing LiF and PEO-derived alkoxides, with LiFSI producing a more LiF-dominant Li-facing surface than LiTFSI under identical conditions. In contrast, LiNO 3 -containing interlayers yield a comparatively thick nitrate-derived interphase that reduces upon Li contact to an N-rich inorganic layer (Li–N/N 3− , LiN x O y ). In Cu| |Li cells, these designed interlayers reduce interfacial resistance and improve Coulombic efficiency (CE) and critical current density (CCD) relative to uncoated controls, with LiFSI- and LiNO 3 -based interlayers providing the highest CE/CCD. Overall, the results demonstrate a manufacturing-relevant approach to engineer SEI nanolayers by salt chemistry, providing pathways to improve reversible lithium metal anodes while reducing reliance on LiTFSI.

anion chemistry↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗
DOE OSTI2023

Natural attenuation of uranium in a fluvial Wetland: Importance of hydrology and speciation

A nuclear fuel fabrication facility released 43,500 kg of uranium into a riparian wetland located on the Savannah River Site between 1955 and 1988. Studies were undertaken to evaluate hydrological and geochemical processes influencing uranium accumulation in the wetland. Gamma-radiation-mapping surveys were conducted by systematically walking over the contaminated wetland with backpacks equipped with global positioning systems and NaI gamma detectors. Based on maps compiled from >700,000 gamma spectra and eight sediment uranium depth profiles, it was determined that 94% of the released uranium remained in the wetland. The uranium in the wetland is concentrated in five multi-hectare areas along the stream, accounting for ~11% of the land area adjacent to the stream. While land type (upland or wetland) and topography provided a reasonable first approximation of where much of the uranium was deposited, hydrological watershed modeling revealed that the stream velocity was especially slow through many of the hot spots. Here, using autoradiography combined with SEM/EDX measurements of contaminated sediments, surprisingly few hot particles were detected. Instead, uranium was evenly distributed throughout the sampled sediment, suggesting that dissolved uranium had bound to sediment particles that became suspended and later deposited in low energy (low flow velocity) portions of the stream. EXAFS suggested that U atoms were present as individual ions in disordered complexes within the sediment. Furthermore, linear combination analyses suggested that the predominant component of the U(VI) was adsorbed to sediment minerals (~70%) and a minor component (~30%) was associated with organic matter phases. Furthermore, these studies show that wetlands can be extraordinarily effective at binding and retaining uranium, thereby providing a natural barrier to the transport of uranium out of a watershed. However, significant anthropogenic or climatic changes to wetlands, such as those associated with flooding, forest fires, or land use, may disrupt the complex hydrological and biogeochemical balance necessary to maintain long-term immobilization of uranium.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Effect of thermal treatments on electrochemical behavior of binder jetted 17-4 PH stainless steel

Binder jetted 17–4 PH stainless steel was post-processed to relate heat treatment, microstructure, and corrosion in 3.5 wt% NaCl. Specimens were sintered at 1380 or 1400 °C, solution-annealed at 1055 °C for 1 h, and aged at 482 °C for 1 h. Here, as-sintered parts showed α′-martensitic matrix with a δ-ferrite network and Cu-rich precipitates in ferrite; inclusions (MnS, NbC) promoted localized attack. solutionizing redistributed elements and reduced ferrite, while aging generated coherent Cu nano-precipitates. Corrosion resistance was highly sensitive to post-processing in which aged specimens exhibited the lowest corrosion current density and formed a thicker, stable Cr 2 O 3 -rich passive film, whereas sintered specimens degraded most. Pitting potential depended on sintering temperature and microstructure, with 1400 °C sintering yielding more positive pitting potentials and the best overall performance after aging. Although the 1400-solutionized condition showed a relatively noble pitting response versus 1380-solutionized, it displayed unstable corrosion kinetics attributed to an imperfect passive film linked to higher NbC density. XPS depth profiles corroborated these trends, showing thicker, more continuous Cr 2 O 3 in aged states and discontinuous/thinning oxides in less resistant conditions. Practically, high-temperature sintering (∼1400 °C) followed by solutionizing and aging is recommended, with further gains expected from reducing NbC/MnS populations and porosity via powder and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multifunctional transitional metal-based phosphide nanoparticles towards improved polysulfide confinement and redox kinetics for highly stable lithium-sulfur batteries

The shuttle effect and the sluggish redox kinetics of lithium polysulfides (LiPSs) are the major issues impeding the practical applications of lithium-sulfur batteries (LSBs). Herein, a highly-efficient Ni 2 P electrocatalyst supported on N, P co-doped graphene (Ni 2 P@NPG) is developed via a simple “recrystallization-self-assembly” method to address the above issues. Here, the ultrafine Ni 2 P nanoparticles ensure abundant adsorption-diffusion-conversion interfaces for accelerating LiPSs transformation and Li 2 S deposition, which extremely decreases the accumulation of LiPSs in the electrolyte and therefore prevents the migration of LiPSs. Their superior catalytic performance is demonstrated by reduced Gibbs free energy changes of rate-limiting step based on the systematic theoretical calculations and the reduced shuttle effect is tested by the three-dimensional reconstructions of Raman depth profiles. Benefiting from these synergistic effects, the LSBs with Ni 2 P@NPG modified separators present a superior cycling performance with an average capacity decay rate of 0.048 % per cycle at 1C around the 400 cycles and a high-rate capacity of 731 mAh/g at 2C. Even with a high-sulfur loading of 3.53 mg cm –2 , the cell can still contain a reversible capacity of 809 mAh/g at 0.2C with a remarkable columbic efficiency of 98.4 %.

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

Parabolic oxidation kinetics of a plutonium alloy at room temperature

Auger electron spectroscopy (AES) was used to characterize the oxide scale formed on an ~0.5 wt. % Ga plutonium metal alloy. Analysis of changes to the characteristic Auger transitions for Pu and O during sputter depth profiles were used in combination with previously measured sputter yields to determine the thickness of the oxide scale and provide insight into its chemistry. From this data, the parabolic kinetics for the Pu-Ga alloy were determined as a function of relative humidity, RH = 14–92%, and oxygen partial pressure, p(O 2 ) = 31–157 Torr, at 25°C. For each relative humidity investigated, the parabolic rate constant was found to be independent of oxygen partial pressure. Furthermore, power-law fits to the relative humidity effect on the rate constant was found to have a power of n = 1.265, indicating the diffusion species which rate-limits the reaction isn’t molecular water, which subsequent dissociates at the oxide-metal interface, or hydroxyl from dissociation at the gas-surface interface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Study the existing form of copper ( p-type oxide/segregation ) and its release mechanism from the passive film of Ti-7Cu alloy

For this work, the existing form of copper and its release mechanism from the passive film of Ti-7Cu alloy was studied by Quantitative analysis of depth profile and series of electrochemical tests. Results were interpreted within the Point Defect Model-II. Copper does not form any p-type (cuprous or cupric) oxide but preferentially oxidises at the metal/film interface and segregated as Cu interstitials ($Cu^+_i$) and become substituted as $Cu^{x'}_{Ti}$ in the oxygen-deficient TiO 2-x barrier layer having n-type semiconducting character. Eventually, the segregated species of Cu preferentially dissolute through the film as Cu ions into the solution at film/solution interface via the Esaki-tunnelling effect.

Siddiqui, Muhammad Ali↗
DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗