Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Material Response”

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 937 records · Page 52

Comparing the Effects of Side Chain Dipole–Dipole Interactions and Hydrogen Bonding on the Mechanical and Electrical Properties of Poly(3-hexylthiophene)

Poly(3-alkylthiophenes) are simple conjugated polymers with good electrical properties, but tend to be brittle, limiting their application. Here, we investigate how side chain modifications can improve their toughness without compromising their charge transport. Specifically, we compare a weak dipole–dipole interaction of ester groups with a stronger hydrogen-bonding interaction of the hydroxyl groups. Two copolymers with targeted 5 and 10 mol % ester-functionalized side chains were synthesized and hydrolyzed in the solid state to introduce hydroxyl groups. Both hydrolyzed polymers became insoluble in most organic solvents, providing a path to layer-by-layer solution processing. All polymers’ optical, thermal, structural, electrical, and mechanical properties were investigated. Optical spectroscopy revealed that functionalization at these low levels did not significantly affect aggregation or π–π stacking. Differential scanning calorimetry indicated reduced crystallinity in the functionalized copolymers compared to that of P3HT. Grazing-incidence wide-angle X-ray diffraction data showed that functionalization led to an increase in lamellar spacing without altering molecular orientation or π–π stacking. Mechanical testing highlighted that at the lower functionalization level, both ester and hydroxyl groups similarly enhanced toughness, suggesting that disorder, rather than the specific nature of the functional group, was responsible for these improvements. At the higher functionalization level, hydroxyl groups significantly increased the elastic modulus while also negatively impacting charge carrier mobility. This study suggests that with stronger hydrogen-bonding groups, only very small amounts can improve P3HT’s toughness. For a wider tunable range, it may be a better strategy to use weak interactions such as dipole–dipole interactions of ester groups.

36 MATERIALS SCIENCE

Physical properties of intergrowth compound Eu 2 CuZn 2 P 3

Here, the physical properties of single-crystalline Eu 2 ⁢CuZn 2 ⁢P 3 were inspected via thermodynamic, transport, and neutron diffraction measurements. Eu 2 ⁢ CuZn 2⁢ P 3 is composed of EuCuP and EuZn 2 ⁢ P 2 motifs stacked in a one-to-one ratio. The magnetic response of Eu 2 ⁢ CuZn 2 ⁢ P 3 is similar to that of EuZn 2 ⁢ P 2 in its A-type antiferromagnetic (semiconducting) state. A Néel temperature of T N = 40.3K is obtained from the specific heat capacity. Single-crystal neutron diffraction reveals a (001) magnetic propagation vector, with primarily in-plane moments coupled along [001] in an up-up-down-down pattern. This spin structure retains the ferromagnetic EuCuP motifs and antiferromagnetic coupling across the EuZn 2 ⁢P 2 motifs. The electrical resistivity is characterized by metallic behavior above 100 K, an increase in ρ upon cooling below ≈ 75K, and a large negative magnetoresistance below T N . First principles calculations evidence a narrow band gap for antiferromagnetic ordering and semimetallic behavior for the ferromagnetic state. This is consistent with the observed effect of magnetic field on the electrical resistivity. Hall effect measurements reveal a field-induced, negative anomalous Hall effect that increases on cooling down to 2 K. Also, an unconventional contribution to the Hall effect dominates the low-field Hall signal, particularly near T N , and is similar in appearance to a topological Hall effect. However, further work is necessary to understand if any of the transport properties have topological origins or if magnetism-induced modifications to the Fermi surface induce these unusual Hall effect signals in this and related materials.

36 MATERIALS SCIENCE

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE

Anomalous Role of Carbon in Pd‐Catalyzed Selective Hydrogenation

Carbonaceous species, including subsurface carbidic carbon and surface carbon, play crucial roles in heterogeneous catalysis. Many reports suggested the importance of subsurface carbon in the selective hydrogenation of alkynes over Pd‐based catalysts. However, the role of surface carbon has been largely overlooked. Here, we demonstrate that subsurface carbon in Pd is not responsible for the selectivity in acetylene hydrogenation. In contrast, the structure of surface carbonaceous species plays a decisive role in hydrogenation selectivity. Electron microscopy and spectroscopy evidence, along with theoretical modelling, reveal that partial graphitization of surface carbonaceous species results in unique spatial confinement of surface reaction intermediates, thus altering the reaction energy landscape in favour of ethylene desorption as opposed to over‐hydrogenation. This mechanism for selectivity control is analogous to enzyme catalysis, where the active centers selectively attract reactants and release products. Similar mechanism may be present in CO/CO 2 hydrogenation and alkane dehydrogenation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Dynamical correlations and nonequilibrium sum rules in photodoped Hubbard ladders

Using matrix product state techniques we study the nonequilibrium dynamical response of the half-filled Hubbard ladder when subject to an optical pump. Optical pumping offers a way of producing and manipulating new strongly correlated phenomena by suppressing existing magnetic correlations. The ladder allows the effects of pump directionality to be investigated, and compared to a single chain it has strong spin-charge coupling and a fully gapped excitation spectrum, promising different nonequilibrium physics. We compute time-dependent correlations, including the nonequilibrium dynamical structure factors for spin and charge. By deriving a combined spin-charge sum rule that applies both in and out-of-equilibrium, we show that spectral weight is pumped directly from the antiferromagnetic spin response into a low energy 𝜔 ∼ 0 charge response below the Mott gap. The transfer of weight is pump direction dependent: pumping directed along the legs disrupts magnetic correlations more than pumping in the rung-direction, even if the post-pump energy density is similar. The charge correlation length is dramatically enhanced by the pump, whilst the spin correlations are most strongly suppressed at nearest- and next-nearest-neighbor spacings. After the pump the system is in a nonthermal correlated metallic state, with gapless charge excitations and approximately equal spin and charge correlation lengths, emphasizing the importance of treating these degrees of freedom on an equal footing in nonequilibrium systems.

36 MATERIALS SCIENCE

The Kinetic Consequences of Water on Catalytic Methane Pyrolysis

Hydrogen production from biomass and natural gas has emerged as a prominent research area in response to the growing demand for energy from alternative sources that minimize CO 2 emissions. In this study, we investigate the impact of water, which is present in and generated from biomass-derived streams, on carbon nanotube (CNT) growth and hydrogen production during methane decomposition using Ni–Mo/MgO as a catalyst. We reveal here that the role of water on CNT growth is highly complex; its effect depends on the stage of growth at which the water is incorporated. When water is introduced at the beginning of methane decomposition ( t = 0 h), methane conversion rates are negatively impacted. We hypothesize that water inhibits the significant phase changes the Ni–Mo/MgO catalyst undergoes during catalyst carburization. In contrast, the incorporation of a small percentage of water after a stabilization period ( t = 3 h) results in methane conversion rate enhancements that scale with the introduced water partial pressure as water selectively reacts with amorphous carbon deposits that lead to catalyst deactivation, thus prolonging the lifetime of some of the most active sites. Moreover, water incorporation after stabilization significantly reduces the apparent activation energy. Density Functional Theory (DFT) calculations reveal that water preferentially interacts with carbon fragments on the catalyst surface to remove carbon deposits with a barrier lower than that required for methane activation, further supporting its role in cleaning active sites on the catalyst surface. Characterization of the resulting carbon nanotubes reveals the formation of more graphitic materials produced in the presence of water, highlighting the impact of water on nanotube properties. These results provide clarity toward the many ways in which water, or cofeeding of biomass-derived materials, may impact catalytic methane pyrolysis rates.

carbon nanotubes

The effect of heat treatment on the defect evolution in LPBF 316H stainless steel

This work investigates the effect of processing and heat treatment on defect evolution and irradiation response of laser powder bed fusion (LPBF) 316H stainless steel, with comparisons to LPBF 316L and wrought 316L/316H. Using in-situ and ex-situ ion irradiations across a wide parameter space—temperature (300–675 °C), dose (0.2–25 dpa), dose rate (10 -3 –10 -5 dpa/s), and helium co-implantation (20–2500 appm)—we correlated void swelling, dislocation loop evolution, and segregation behavior with pre-irradiation microstructures (as-built, stress-relieved, solution-annealed, and cold-worked). Results show that swelling is strongly controlled by dislocation density: intermediate densities maximize swelling, while solution annealing or cold working reduce susceptibility to levels comparable to wrought alloys. Low dose rates and helium both promote cavity nucleation and lower the incubation barrier, with helium suppressing the role of dislocation density and driving swelling behavior toward wrought-like response. Loop evolution in SA LPBF 316H resembles wrought 316L but shows localized denuded zones near low-angle grain boundaries. STEM-EDS mapping further revealed Ni segregation at voids and sparse Al-rich oxides, without evidence of Ni–Si precipitates. Collectively, these findings identify dislocation density, helium content, and dose rate as the factors governing swelling in LPBF 316H and provide mechanistic datasets for model validation, supporting LAIN and the qualification of AM austenitic steels for nuclear service.

316 stainless steel

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE

High-Sensitivity NO 2 Gas Sensor: Exploiting UV-Enhanced Recovery in a Hexadecafluorinated Iron Phthalocyanine-Reduced Graphene Oxide

Monitoring ultralow nitrogen dioxide (NO 2 ) concentrations is crucial for air quality management and public health. However, the existing NO 2 gas sensors have several defects, like high cost and power consumption, and exhibit poor selectivity. This study addresses these challenges by presenting a novel hexadecafluorinated iron phthalocyanine-reduced graphene oxide (FePcF 16 -rGO) covalent hybrid sensor for NO 2 detection. This innovative approach, which overcomes the limitations of fabrication cost, energy efficiency, and gas selectivity, is a significant step forward in gas sensor technology. The sensor demonstrates exceptional sensitivity toward ultralow NO 2 concentrations (15.14% response for 100 ppb) with a rapid 60 s UV light-induced recovery. Additionally, the sensor exhibits high selectivity for NO 2 , achieving a limit of detection (LOD) of 8.59 ppb. This approach paves the way for developing cost-effective, energy-efficient, and miniature NO 2 monitoring devices for improved environmental monitoring and enhanced safety in workplaces where NO 2 exposure is a concern.

36 MATERIALS SCIENCE

Strain-induced superconductivity in RuO 2 (100) thin-films

Recently, superconductivity was discovered with a superconducting transition temperature (T c ) of 2 K in strained (110)-oriented RuO 2 films grown on TiO 2 (110) single-crystal substrates. In this work, we predict and realize superconductivity in strained (100)-oriented RuO 2 thin films grown on TiO 2 (100) single-crystal substrates. We show that while density functional theory predicts the T c of strained RuO 2 (100) films to be even higher than the RuO 2 (110) films, our transport and angle-resolved photoemission spectroscopy measurements find the T c to be about 1 K in strained RuO 2 (100) films grown on TiO 2 (100) substrates. Nonetheless, the thickness dependence of the T c follows a similar trend in both cases. Our scanning SQUID measurements reveal a local superfluid response, indicating a 100 mK inhomogeneity in T c over a 100 μm scale.

Electronic structure

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING

Stratification of fluoride uptake among enamel crystals with age elucidated by atom probe tomography

Dental enamel is subjected to a lifetime of de- and re-mineralization cycles in the oral environment, the cumulative effects of which cause embrittlement with age. However, the understanding of atomic scale mechanisms of dental enamel aging is still at its infancy, particularly regarding where compositional differences occur in the hydroxyapatite nanocrystals and what underlying mechanisms might be responsible. Here, we use atom probe tomography to compare enamel from a young (22 years old) and a senior (56 years old) adult donor tooth. Findings reveal that the concentration of fluorine is elevated in the shells of senior nanocrystals relative to young, with less significant differences between the cores or intergranular phases. It is proposed that the embrittlement of enamel is driven, at least in part, by the infusion of fluorine into the nanocrystals and that the principal mechanism is de- and re-mineralization cycles that preferentially erode and rebuild the nanocrystals shells.

36 MATERIALS SCIENCE

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE

Nondestructive ultrasonic characterization of a triple-weld-bead wire arc additively manufactured ER70S-6 S-curved wall

Understanding build-scale microstructural variation in wire arc additive manufacturing (WAAM) of low-carbon steels is essential for ensuring consistent structure–property relationships throughout large components. Conventional destructive characterization techniques, such as scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), are time-intensive and limited to localized regions, making comprehensive evaluation of large WAAM structures challenging. In this study, a nondestructive ultrasonic approach was employed to characterize a 252 mm tall ER70S-6 S-curved WAAM wall produced using a triple-bead deposition strategy. Optical and SEM analysis revealed a repeating dual-region microstructure consisting of uniform polygonal ferrite at melt pool centers and heterogeneous ferrite with coarse and fine grains near melt pool boundaries, attributed to cyclic thermal conditions. Longitudinal ultrasonic backscatter imaging was used to evaluate the continuity of this periodicity along the full build height. The ultrasonic response exhibited a consistent repeating pattern that correlated with the observed layer-wise microstructural variation. X-ray computed tomography confirmed the absence of detectable porosity, indicating that ultrasonic contrast is primarily governed by grain morphology. Overall, the results indicate that longitudinal backscatter ultrasound is a promising nondestructive characterization technique for validating microstructural variations along the s-curved WAAM wall, with significant potential for microstructure optimization and process control.

36 MATERIALS SCIENCE

Thermal Disorder‐Induced Strain and Carrier Localization Activate Reverse Halide Segregation

The reversal of halide ions is studied under various conditions. However, the underlying mechanism of heat-induced reversal remains unclear. This work finds that dynamic disorder-induced localization of self-trapped polarons and thermal disorder-induced strain (TDIS) can be co-acting drivers of reverse segregation. Localization of polarons results in an order of magnitude decrease in excess carrier density (polaron population), causing a reduced impact of the light-induced strain (LIS – responsible for segregation) on the perovskite framework. Meanwhile, exposing the lattice to TDIS exceeding the LIS can eliminate the photoexcitation-induced strain gradient, as thermal fluctuations of the lattice can mask the LIS strain. Under continuous 0.1 W cm -2 illumination (upon segregation), the strain disorder is estimated to be 0.14%, while at 80 °C under dark conditions, the strain is 0.23%. However, in situ heating of the segregated film to 80 °C under continuous illumination (upon reversal) increases the total strain disorder to 0.25%, where TDIS is likely to have a dominant contribution. Therefore, the contribution of entropy to the system's free energy is likely to dominate, respectively. Various temperature-dependent in situ measurements and simulations further support the results. These findings highlight the importance of strain homogenization for designing stable perovskites under real-world operating conditions.

36 MATERIALS SCIENCE

From Detector Layout to Signal Analysis: Geometry Optimization and Neutron-Gamma Tagging in Plastic Scintillator Detectors

Neutrinos are elementary particles with many properties still unknown. Their masses so far have only upper and lower limits. Still, due to neutrino oscillations, it is clear that they are not massless, as stated by the Standard Model of Elementary Particles. Neutrinos are also present in the Universe in vast amounts, but they rarely interact with the surrounding matter. Their abundance makes them very interesting for many theories beyond the Standard Model, e.g., dark matter searches and CP symmetry violation in the leptonic sector, which could be (partially) responsible for the observed matter-antimatter asymmetry in today's Universe.\\ The Deep Underground Neutrino Experiment (DUNE) is a next-generation accelerator-based neutrino oscillation experiment that will study neutrinos with unprecedented precision and may answer many open questions. DUNE uses a powerful neutrino beam from Fermilab. It consists of a Near Detector (ND) complex to measure neutrinos before oscillat ion, and a Far Detector (FD) complex $1300\ \mathrm{km}$ away to measure them after oscillation. As part of the ND complex, measurements are also possible with different angles to the neutrino beam. This enables excellent control of systematic uncertainties of e.g., neutrino cross section measurements.\\ One detector in the near detector complex is The Muon Spectrometer (TMS), an extension of a Liquid Argon (LAr) detector that measures the charge and momentum of muons produced in neutrino interactions within the LAr. The design of this detector, which consists of alternating layers of steel and plastic scintillator bars, must be optimized for the expected muon energies.\\ In this thesis, a study of the optimal module orientation plan is presented, which is necessary for the physics performance of the near detector complex and, by extension, DUNE. As part of this study, the event reconstruction was also developed and improved. Simulated muons are then reconstructed, and the performan ce of different module orientation plans is tested.\\ As a second part, a study of neutron and gamma tagging using a Pulse Shape Discriminating (PSD) plastic scintillator is presented. The properties of this material allow particle differentiation based on the temporal distribution of emitted light. A novel approach to using the individual light signals was successfully tested using data from a small, local test setup.

Nehm, Asa [Mainz U.]

In Situ Insights into Enhanced Cooperative Ligand Exchange Kinetics via Solvent-Induced Restacking in a 2D Metal–Organic Framework

Understanding the reaction kinetics at catalytically active sites is crucial for integrating catalytic two-dimensional (2D) materials into industrial processes. This study focuses on in situ observation of ligand exchange kinetics and solvent-assisted structural restacking transition in the 2D paddle wheel-based MOF [Cu 2 (dttc) 2 ] n (DUT-134(Cu), dttc = dithieno[3,2-b:2′,3′- d]thiophene-2,6-dicarboxylate). The ligand exchange process, involving the replacement of dimethylformamide (DMF) with nitriles such as acetonitrile (ACN), pentanenitrile, and heptanenitrile, was investigated using advanced in situ characterization techniques with high temporal resolution, including powder X-ray diffraction and Raman spectroscopy. The larger analytes exhibited reduced exchange rates, consistent with enhanced steric hindrance and greater diffusion constraints. Interestingly, the study revealed that the exchange of DMF with ACN induces a structural transition to higher symmetry within few seconds, a transition from AB to AA stacking mode of the layers, and a widening of the interlayer distance. Crucially, this structural transition dramatically accelerates the solvent exchange process through cooperative effects, offering critical advantages for catalytic applications. Notably, the reverse exchange from ACN to DMF proceeds more slowly and does not reverse the structural changes, but a new phase is formed with preserved AA stacking. By isotope labeling of linker molecules in combination with two complementary theoretical vibrational simulation methods, the precise assignment of Raman bands and the vibrational modes associated with the ligand exchange process could be achieved. These pioneering insights into the dynamic behavior of 2D MOFs, coupled with ligand exchange, establish a highly promising and transformative approach to achieving enhanced tunability and responsiveness in future catalytic applications.

Layers

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator