Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composition B”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
DOE OSTI2022

Liquid fragility maximum in lithium borate glass‐forming melts related to the local structure

Abstract The structure of liquid lithium pyroborate, Li 4 B 2 O 5 ( J = Li/B = 2), has been measured over a wide temperature range by high‐energy X‐ray diffraction, and compared to that of its glass and borate liquids of other compositions. The results indicate a gradual increase in tetrahedral boron fraction from 3(1)% to 6(1)% during cooling from T = 1271(15) to 721(8) K, consistent with the larger N 4 = 10(1)% found for the glass, and literature 11 B nuclear magnetic resonance measurements. van't Hoff analysis based on a simple boron isomerization reaction BØ 3 O 2 – ⇌ BØO 2 2– yields Δ H = 13(1) kJ mol –1 and Δ S = 40(1) J mol –1 K –1 for the boron coordination change from 4 to 3, which are, respectively, smaller and larger than found for singly charged isomers for J ≤ 1. With these, we extend our model for N 4 ( J , T ), nonbridging oxygen fraction f nbr ( J , T ), configurational heat capacity , and entropy S conf ( J , T ) contributions up to J = 3. A maximum is revealed in at J = 1, and shown semi‐quantitatively to lead to a corresponding maximum in fragility contribution, akin to that observed in the total fragilities by temperature‐modulated differential scanning calorimetry. Lithium is bound to 4.6(2) oxygen in the pyroborate liquid, with 2.7(1) bonds centered around 1.946(8) Å and 1.9(1) around 2.42(1) Å. In the glass, n LiO = 5.4(4), the increase being due to an increase in the number of short Li–O bonds.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Saturn's Rings

The rings are changing before our eyes; structure varies on all timescales and unexpected things have been discovered. Many questions have been answered, but some answers remain elusive. Here we highlight the major ring science progress over the mission to date, and describe new observations planned for Cassini’s final three years. Ring Composition and particle sizes: The rings are nearly all water ice with no other ices – so why are they reddish? The C Ring and Cassini Division are “dirtier” than the more massive B and A Rings, as shown by near- IR and, recently, microwave observations. Particle sizes, from stellar and radio occultation's, vary from place to place. Ring structure, micro and macro: numerous spiral density waves and ubiquitous “self-gravity wakes” reveal processes which fostered planet formation in the solar system and elsewhere. However, big puzzles remain regarding the main ring divisions, the C Ring plateau structures, and the B Ring irregular structure. Moonlets, inside and out, seen and unseen: Two gaps contain sizeable moonlets, but more gaps seem to contain none; even smaller embedded “propeller” objects wander, systematically or randomly, through the A ring. Rubble pile ring moons just outside the rings may escaped from the rings, and the recently discovered “Peggy” may be trying this as we watch. Impact bombardment of the rings: Comet fragments set the rings to rippling on century-timescales, and boulders crash through hourly; meanwhile, the constant hail of in falling Kuiper belt material has a lower mass flux than previously thought. Origin and Age of the Rings: The ring mass and bombardment play key roles. The ring mass is well known everywhere but in the B Ring (where most of it is). New models suggest how tidal breakup of evolving moons may have formed massive ancient rings, of which the current ring is just a shadow. During its last three years, the Cassini tour profile will allow entirely new observations: direct measurement of the still-unknown ring mass; direct in-situ sampling of ring particle composition (targeting the iron- or carbon based red non-icy component); and radar backscattering observations.

Cassini↗
NASA NTRS2020

Capturing an Evolving Nebular Environment: A Petrographic and Geochemical Study of a Type A, B & C CAI

Calcium, Aluminum-rich Inclusions (CAIs) were the first formed solids in our Solar System, with mineral assemblages reflecting the first phases predicted to condense out of a hot nebular gas of Solar composition. Geochemical, textural and crystallographic information in CAIs can be used to constrain the temperature, pressure, and composition (e.g., oxygen fugacity) of the gaseous reservoir(s) from which they formed, as well as any secondary (nebular and parent body) processes they underwent. Coordinated geochemical and textural analyses provide information on nebular conditions (i.e., astrophysical environments and dynamics of nebular gas reservoirs) in which these CAIs formed. In order to better understand the evolution of nebular reservoirs at the time of CAI formation, we analyzed a Type A, B and C CAI using Electron Probe Micro-Analyzer (EPMA) and Electron BackScatter Diffraction (EBSD) at NASA Johnson Space Center (JSC).

Comins, M. B.↗
NASA NTRS1990

Inter- and intra-group compositional variations in Apollo 15 pyroclastic green glass - An electron- and ion-microprobe study

Results are presented on major- and trace-element abundance analyses of Apollo 15 pyroclastic green glasses from groups A, B, C, D, and E, carried out using electron- and ion-microprobe techniques. The diagrams depicting Sr, Zr, Ba, and Nd vs Co variations indicate the presence of a high-Co trend in groups A and D and a low-Co trend in groups B and C. Group-E glasses were found to be significantly enriched in Sr, relative to the other four glass groups. Chemical data of this study were integrated with previous data to evaluate various magmatic processes that have been proposed in the past to explain chemical variations in the lunar green glass. Results of calculations using a source mixing model suggest that the Apollo 15 green glasses represent multiple eruptive events from three chemically distinct but compositionally variable source regions.

Galbreath, K. C.↗
NASA NTRS1994

Ultrasonic techniques for repair of aircraft structures with bonded composite patches

This is a paper on a research and development project to demonstrate a novel ultrasonic process for the field application of boron/epoxy (B/Ep) patches for repair of aircraft structures. The first phase of the project was on process optimization and testing to develop the most practical ultrasonic processing techniques. Accelerated testing and aging behavior of precured B/Ep patches, which were ultrasonically bonded to simulated B-52 wing panel assemblies, were performed by conducting flight-by-flight spectrum loading fatigue tests. The spectrum represented 2340 missions/flights or 30 years of service. The effects of steady-state applied temperature and prior exposure of the B/Ep composite patches were evaluated. Representative experimental results of this phase of the project are presented.

Smith, S. H.↗
DOE OSTI2022

SEI formation mechanisms and Li + dissolution in lithium metal anodes: Impact of the electrolyte composition and the electrolyte-to-anode ratio

The lithium metal battery is one of today's most promising high-energy-density storage devices. Its full-scale implementation depends on solving operational and safety issues intrinsic to the Li metal high reactivity leading to uncontrolled electrolyte decomposition and uneven Li deposition. In this work, we study the spontaneous formation of the solid electrolyte interphase (SEI) upon contact of Li metal with the electrolyte and describe the heterogeneous SEI morphological features. Multiple electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethoxyethane (DME), dimethyl carbonate (DMC), 1,2,2-tetrafluoroethyl-2,2,3,3-tetrafluoropropyl ether (TTE) and bis(2,2,2-trifluoroethyl) ether (BTFE) are used. Findings include the description of the SEI evolution from dispersed LiO, LiS, LiN, and LiF clusters to a continuous and compact inorganic phase in which the LiO and LiF content depend on the presence of fluorine diluents. Further, the role of the DME ether solvent helping the growth of a “wet-SEI” is compared to that of the highly unstable carbonate DMC, which decomposes into complex radical oligomers that might contribute to further electrolyte decomposition. The impact of the electrolyte-to-anode ratio on LiFSI decomposition is highlighted. Finally, we suggest the existence of a critical LiFSI concentration and electrolyte-to-anode ratio that could potentially balance the rate of electrolyte depletion and lithium consumption.

25 ENERGY STORAGE↗
DOE OSTI2024

Understanding Environmental Barrier Coating Lifetimes and Performance for Industrial Gas Turbines

Hydrogen or hydrogen blend fuels are expected to replace natural gas in land-based industrial gas turbines (IGTs) to support a greener power economy. Silicon carbide (SiC) base ceramic matrix composites (CMCs) are considered for replacement of Ni-based superalloys to facilitate future efficiency improvements. SiC CMCs require environmental barrier coatings (EBCs) to mitigate volatilization from high-temperature steam, thus making the EBC lifetime critical information for identifying CMC component lifetimes. Here, the goal of this project is to determine the maximum bond coating temperature underneath the EBC for achieving an IGT component lifetime goal of 25,000 h, which is far greater than current CMC component lifetime requirements for aeroturbine applications. To provide data for the lifetime model, laboratory testing used atmospheric plasma-sprayed rare-earth silicate EBCs on monolithic SiC substrates with an intermediate Si bond coating. Specimens exposed to 1-h thermal cycles in flowing air–steam environments and reaction kinetics were assessed from 700 °C to 1350 °C by measuring the thickness of the thermally grown silica scales. The silica growth and phase transformation appear critical in predicting EBC lifetime and several strategies have been explored to reduce the oxide growth rate and improve EBC durability at elevated temperatures. Advanced characterization using Raman spectroscopy has helped clarify this system.

08 HYDROGEN↗
DOE OSTI2022

Silver Nanoparticles Functionalized Nanosilica Grown over Graphene Oxide for Enhancing Antibacterial Effect

The continuous growth of multidrug-resistant bacteria due to the overuse of antibiotics and antibacterial agents poses a threat to human health. Silver nanoparticles, silica-based materials, and graphene-based materials have become potential antibacterial candidates. In this study, we developed an effective method of enhancing the antibacterial property of graphene oxide (GO) by growing nanosilica (NS) of approximately 50 nm on the graphene oxide (GO) surface. The structures and compositions of the materials were characterized through powdered X-ray diffraction (P-XRD), transmission electron microscopy (TEM), scanning electron microscopy coupled with energy dispersive X-ray spectroscopy (SEM-EDS), ultraviolet–visible spectroscopy (UV–VIS), dynamic light scattering (DLS), Raman spectroscopy (RM), Fourier-transform infrared spectroscopy (FTIR), Brunauer–Emmet–Teller (BET) surface area, and pore size determination. The silver nanoparticles (AgNPs) with an average diameter of 26 nm were functionalized on the nanosilica (NS) surface. The composite contained approximately 3% of silver nanoparticles. The silver nanoparticles on nanosilica supported over graphene oxide (GO/NS/AgNPs) exhibited a 7-log reduction of Escherichia coli and a 5.2-log reduction of Bacillus subtilis within one hour of exposure. Both GO/NS and GO/NS/AgNPs exhibited substantial antimicrobial effects against E. coli and B. subtilis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Distribution of Infiltration in the 216-U-10, 216-B-3 Pond, and 216-T-4 Pond Systems 1944-1997

A modeling system to simulate fate and transport of contaminant plumes in groundwater beneath the Central Plateau to support the Cumulative Impact Evaluation (CIE) and Composite Analysis (CA) is being developed. The CIE will evaluate the cumulative effect of multiple Central Plateau radiological and chemical sources on groundwater quality. This will enable cleanup decisions to be evaluated in the context of all Central Plateau waste sites and existing groundwater conditions. The purpose of this Environmental Calculation File (ECF) is to partition the infiltration of wastewater and contaminant inventory between the main pond and each expansion lobe of the 216-B-3 Pond system, and to partition the infiltration of wastewater and contaminant inventory between the ditches and ponds of the 216-U-10 Pond and 216-T-4 Pond systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2014

Three-Dimensional Imaging and Numerical Reconstruction of Graphite/Epoxy Composite Microstructure Based on Ultra-High Resolution X-Ray Computed Tomography

A combined experimental and computational study aimed at high-resolution 3D imaging, visualization, and numerical reconstruction of fiber-reinforced polymer microstructures at the fiber length scale is presented. To this end, a sample of graphite/epoxy composite was imaged at sub-micron resolution using a 3D X-ray computed tomography microscope. Next, a novel segmentation algorithm was developed, based on concepts adopted from computer vision and multi-target tracking, to detect and estimate, with high accuracy, the position of individual fibers in a volume of the imaged composite. In the current implementation, the segmentation algorithm was based on Global Nearest Neighbor data-association architecture, a Kalman filter estimator, and several novel algorithms for virtualfiber stitching, smoothing, and overlap removal. The segmentation algorithm was used on a sub-volume of the imaged composite, detecting 508 individual fibers. The segmentation data were qualitatively compared to the tomographic data, demonstrating high accuracy of the numerical reconstruction. Moreover, the data were used to quantify a) the relative distribution of individual-fiber cross sections within the imaged sub-volume, and b) the local fiber misorientation relative to the global fiber axis. Finally, the segmentation data were converted using commercially available finite element (FE) software to generate a detailed FE mesh of the composite volume. The methodology described herein demonstrates the feasibility of realizing an FE-based, virtual-testing framework for graphite/fiber composites at the constituent level.

Czabaj, M. W.↗
DOE OSTI2021

Composite fermion nonlinear sigma models

We study particle-hole symmetry at the integer quantum Hall plateau transition using composite fermion mean-field theory. Because this theory implicitly includes some electron-electron interactions, it also has applications to certain fractional quantum Hall plateau transitions. Previous work [P. Kumar et al., Phys. Rev. B 100, 235124 (2019)] using this approach showed that the diffusive quantum criticality of this transition is described by a nonlinear sigma model with topological θ = π term. This result, which holds for both the Dirac and Halperin, Lee, and Read composite fermion theories, signifies an emergent particle-hole (reflection) symmetry of the integer (fractional) quantum Hall transition. Here we consider the stability of this result to various particle-hole symmetry-violating perturbations. In the presence of quenched disorder that preserves particle-hole symmetry, we find that finite longitudinal conductivity at this transition requires the vanishing of a symmetry-violating composite fermion effective mass, which if present would generally lead to θ ≠ π and a corresponding violation of particle-hole symmetric electrical transport σ xy ≠ $\frac1{2}$$\frac{e^2}{h}$ . When the disorder does not preserve particle-hole symmetry, we find that θ can vary continuously within the diffusive regime. Our results call for further study of the universality of the quantum Hall plateau transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The Post‐Perovskite Transition in Fe‐ and Al‐Bearing Bridgmanite: Effects on Seismic Observables

Abstract The primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al 3+ , Fe 3+ )‐, (Al 3+ )‐, (Fe 2+ )‐, and (Fe 3+ )‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe 3+ ‐, Al 3+ ‐, or (Al 3+ , Fe 3+ )‐alloying increases the transition pressure, while Fe 2+ ‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO 3 , and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.

Geochemistry & Geophysics↗
DOE OSTI2022

Extreme fast charge aging: Correlation between electrode scale and heterogeneous degradation in Ni-rich layered cathodes

Extreme fast charging (XFC) is a key requirement for the adoption of battery-based electric vehicles by the transportation sector. However, XFC has been shown to accelerate degradation, causing the capacity, life, and safety of batteries to deteriorate. There are no systematic studies in the open literature regarding aging modes in Ni-rich Li y Ni 0.a Mn 0.b Co 0.c O 2 (NMCabc) cathodes caused by fast charging. Herein, we report the effects of cathode composition and electrode loading in pouch cells containing NMC532, 622 and 811 paired with graphite and cycled under XFC conditions. The relative anisotropic volume change in the unit cell increases with Ni content in low-loading cells, while it levels up for all three NMC cathodes in high-loading cells because of substantial Li plating. The amounts of lithium plating and heterogeneity on the anode, determined by crystallographic phase quantification, were found to be correlated with electrode loading and cathode heterogeneity. The concentrations of the transition metals deposited on the anodes depend on NMCabc composition in a complex way. More particle cracking and surface degradation was found in NMC811. Here, the findings in this work provide a new understanding of the failure mechanisms and their practical implications for compositional tuning of future high-Ni NMCabc cathode materials during fast charging.

25 ENERGY STORAGE↗
DOE OSTI2022

Turbulent Flux Measurements of the Near‐Surface and Residual‐Layer Small Particle Events

Abstract According to recent field studies, almost half of the New Particle Formation (NPF) events occur aloft, in a residual layer, near the top of the boundary layer. Therefore, measurements of the meteorological parameters, precursor gas concentrations, and aerosol loadings conducted at the ground level are often not representative of the conditions where the NPFs take place. This paper presents new measurements obtained during the Turbulent Flux Measurements of the Residual Layer Nucleation Particles, conducted at the Southern Great Plains research site. Vertical turbulent fluxes of 3–10 nm‐sized particles were measured using a sonic anemometer and two condensation particle counters with nominal cutoff diameters of 3 nm and 10 nm mounted at the top of the 10‐m telescoping tower. Aerosol number size distribution (5–300 nm) was determined through the ground‐based Scanning Mobility Particle Sizers. The size selected (15–50 nm) particle hygroscopicity was derived with the Humidified Tandem Differential Mobility Analyzer. The ground level observations were supplemented by vertically‐resolved measurements of horizontal and vertical wind speeds and aerosol backscatter. The data analysis suggests that (a) turbulent flux measurements of 3–10 nm particles can distinguish between near‐surface and residual‐layer small particle events; (b) sub‐50 nm particles had a hygroscopicity value of 0.2, suggesting that organic compounds dominate atmospheric nanoparticle chemical composition at the site; and (c) current methodologies are inadequate for estimating the dry deposition velocity of sub‐10 nm particles because it is not feasible to measure particle concentration very near the surface, in the diffusion sublayer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Train, Inform, Borrow, or Combine? Approaches to Process–Guided Deep Learning for Groundwater–Influenced Stream Temperature Prediction

Although groundwater discharge is a critical stream temperature control process, it is not explicitly represented in many stream temperature models, an omission that may reduce predictive accuracy, hinder management of aquatic habitat, and decrease user confidence. We assessed the performance of a previously-described process-guided deep learning model of stream temperature in the Delaware River Basin (USA). We found lower accuracy (root mean square error [RMSE] of 1.71 versus 1.35°C) and stronger seasonal bias (absolute mean monthly bias of 1.06 vs. 0.68°C) for reaches primarily influenced by deep groundwater as compared to atmospheric conditions. We then tested four approaches for improving groundwater process representation: (a) a custom loss function leveraging the unique patterns of air and water temperature coupling characteristic of different temperature drivers, (b) inclusion of additional groundwater-relevant catchment attributes, (c) incorporation of additional process model outputs, and (d) a composite model. The custom loss function and the additional attributes significantly improved the predictive accuracy in groundwater-dominated reaches (RMSE of 1.37 and 1.26°C) and reduced the seasonal bias (absolute mean monthly bias of 0.44 and 0.48°C), but neither approach could identify holdout groundwater reaches. Variable importance analysis indicates the custom loss function nudges the model to use the existing inputs more efficiently, whereas with the added features the model relies on a broader suite of inputs. This analysis is a substantial step toward more accurately representing groundwater discharge processes in stream temperature models and will improve predictive accuracy and inform habitat management.

54 ENVIRONMENTAL SCIENCES↗