Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Behavior”

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 541 records · Page 30

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗

In Situ Local Measurement of Austenite Mechanical Stability and Transformation Behavior in Third-Generation Advanced High-Strength Steels

Austenite mechanical stability, i.e., retained austenite volume fraction (RAVF) variation with strain, and transformation behavior were investigated for two third generation advanced high strength steels (3GAHSS) under quasi-static uniaxial tension: a 1200 grade, two-phase medium Mn (10 wt%) TRIP steel, and a 980 grade, three-phase TRIP steel produced with a quenching and partitioning heat treatment. The Mn (10 wt%) TRIP steel deforms inhomogeneously via propagative instabilities (Lüders and Portevin Le Châtelier-like bands) while the 980 grade TRIP steel deforms homogenously up to necking. The dramatically different deformation behaviors of these steels required the development of a new in-situ experimental technique that couples volumetric synchrotron X-ray diffraction measurement of RAVF with surface strain measurement using stereo digital image correlation (DIC) over the beam impingement area. Measurement results with the new technique are compared with those from a more conventional approach wherein strains are measured over the entire gauge region while RAVF measurement is the same as that in the new technique. A determination is made as to the appropriateness of the different measurement techniques in measuring the transformation behaviors for steels with homogeneous and inhomogeneous deformation behaviors. Here, extension of the new in-situ technique to the measurement of austenite transformation under different deformation modes and to higher strain rates is discussed.

36 MATERIALS SCIENCE↗

Assessing Constitutive Models for Prediction of High-Temperature Flow Behavior with a Perspective of Alloy Development

The utility of different constitutive models describing high-temperature flow behavior has been evaluated from the perspective of alloy development. Strain compensated Arrhenius model, modified Johnson–Cook (MJC) model, model D8A and artificial neural network (ANN) have been used to describe flow behavior of different model alloys. These alloys are four grades of SS 316LN with different nitrogen contents ranging from 0.07 to 0.22%. Grades with 0.07%N and 0.22%N have been used to determine suitable material constants of the constitutive equations and also to train the ANN model. While the ANN model has been developed with chemical composition as a direct input, the MJC and D8A models have been amended to incorporate the effect of nitrogen content on flow behavior. The prediction capabilities of all models have been validated using the experimental data obtained from grades containing 0.11%N and 0.14%N. The comparative analysis demonstrates that ‘N-amended D8A’ and ‘N-amended MJC’ are preferable to the ANN model for predicting flow behavior of different grades of 316LN. The work provides detailed insights into the usual statistical error analysis technique and frames five additional criteria which must be considered when a model is analyzed from the perspective of alloy development.

36 MATERIALS SCIENCE↗

Microstructural Characterization and Electrochemical Behavior of AA2014 Al-Cu-Mg-Si Alloy of Various Tempers

AA2014 Al-Cu-Mg-Si alloy is subjected to aging heat treatment, and microstructural features and electrochemical behavior of the AA2014 Al-Cu-Mg-Si alloy of various tempers have been characterized. Hardness measurement of the alloy of different tempers showed typical age hardening behavior. X-ray diffraction and transmission electron microscopy studies revealed the presence of characteristic phases. Differential scanning calorimetry study displayed the sequence of solid-state precipitation and dissolution reactions in the AA2014 alloy of various tempers. Electrochemical behavior was assessed by measurement of variation of open-circuit voltage with time, potentiodynamic polarization and potentiodynamic cyclic polarization of the AA2014 alloy of various tempers in 3.5 wt.% NaCl solution at near neural (pH ~ 7) and at alkaline (pH 12) conditions. Potentiodynamic polarization studies of different tempers in 3.5 wt.% NaCl solution at pH ~ 7 showed that the corrosion potential (E{sub corr}) value has shifted toward more negative (active) direction with increasing aging time. The cyclic polarization in 3.5 wt.% NaCl solution at pH ~ 7 displayed active–passive behavior and also showed a large hysteresis loop indicating high pit growth damage.

36 MATERIALS SCIENCE↗

Dislocation Density-Based Model for Flow Behavior of a Near-α Titanium Alloy Considering Effects of Initial Lamellar Thickness

A dislocation density-based constitutive model that considered the effects of initial lamellar thickness was proposed to predict the flow behavior of a near-α high-temperature Ti-5.4Al-3.7Sn-3.3Zr-0.5Mo-0.4Si alloy. The evolution equations of dislocation density were improved by introducing the influence of the dynamic globularization of lamellar microstructure with different initial thicknesses. The evolution equations and Hall-Petch behavior were coupled into the constitutive model, which was composed of athermal and thermal stresses. According to the established constitutive equations, the variation and contributions of corresponding stress components were analyzed. The calculated dislocation densities were compared with the experimental values obtained by electron backscattered diffraction data to verify the model. Moreover, the flow stresses of the studied alloy with different initial lamellar microstructures under various deformation conditions were predicted. Good agreements between the predicted and experimental results were obtained. Additionally, steady-state flow behavior was further analyzed and the relationship between the calculated steady-state flow stress and the equiaxed α-phase size was quantitatively obtained. The results indicated the validity of the proposed constitutive model for predicting the flow behavior of the alloy with different initial lamellar microstructures over a range of deformation temperatures and strain rates.

36 MATERIALS SCIENCE↗

Microstructure and Work Hardening Behavior of Micro-plasma Arc Welded AISI 316L Sheet Joint

The aim of this research is to investigate the microstructural changes and work hardening behavior of AISI 316L 0.5-mm-thick sheet due to micro-plasma arc welding. AISI 316L similar sheets were welded in single-pass square butt joint configuration. The microstructure in the fusion zone typically contains a variety of complex δ-ferrite-austenitic structure, which significantly increased hardness value compared to the base material because of rapid solidification of the weld pool. Systematic tensile test was performed to investigate strain rate sensitivity of the welded joint under quasi-static loading conditions. The multistage work hardening behavior was determined by Kocks–Mecking (K–M) model and differential Crussard–Jaoul analysis. In both the cases, stage III work hardening behavior and stage IV work hardening behavior were observed. The stress–strain curves of welded specimen at different strain rates were analyzed in terms of Hollomon, Ludwik and Swift equations to determine work hardening exponent values. The welded joints were ruptured in the transition zone of heat affected and fusion zones. The fractographic observation exhibited variation in dimple and void density over strain rate from 0.0001 to 0.0015 s{sup −1}.

36 MATERIALS SCIENCE↗

Mechanical and Electrochemical Behavior of Dual-Phase Steels Having Varying Ferrite–Martensite Volume Fractions

Dual-phase (DP) steels behave differently under corrosive environment depending upon the volume fractions of phases, mainly martensite, ferrite and bainite. The mechanical behavior and corrosion resistance of DP steels with varying martensite (from 32 to 100%) are investigated. The increase in martensite volume fraction of DP steel results in significant improvement in tensile strength and hardness with a subsequent reduction in elongation. Unlike the mechanical behavior, the increase in martensite volume fraction in steel does not exhibit a direct relationship with their corrosion resistance behavior. The corrosion mechanism in DP steels having different martensite and ferrite volume fractions depends on the self-corrosion resistance behavior of both the phases as well as the occurrence of galvanic corrosion between them.

36 MATERIALS SCIENCE↗

Mesoscale Numerical Characterization of Kevlar and Carbon–Kevlar Hybrid Plain-Woven Fabric Compression Behavior

A mesoscale finite element model (FEM) of the plain-woven fabric representative elementary volume (REV/unit cell) is developed. The developed FEM is used to predict the nonlinear compaction behavior of monolithic and intraply hybrid fabric made of high-strength carbon and Kevlar fibers. For geometrical modeling of the fabric REV, TexGen, an open-source geometrical modeling software (developed at Nottingham University), is utilized, and for mechanical characterization, FE simulation tool Abaqus{sup ®} is used. Geometrical nonlinearity of the fabric, including periodic boundary conditions and yarn interactions, is incorporated in the REV FE model. Material nonlinearity caused by transversely isotropic behavior of fibrous yarns is implemented through the user material subroutine form in FE simulation along with local orientation assigned to each yarn in the REV. Material properties of fibers are extracted from the literature and supplier data. The compression behavior of monolithic Kevlar and hybrid carbon–Kevlar (carbon warp and Kevlar weft yarns, C–K) fabric is predicted using the proposed FE model. FEM-predicted results were found to have good agreement with the experimentally obtained ones. Results depict that C–K hybrid fabrics have improved compression strength. Hybridization effect is governed by yarns mutual interaction, yarn area at crossovers, and transverse in-plane properties of the interyarn hybridized fabric. The proposed FE model can be used to predict the mechanical behavior of woven fabrics with different weaving patterns, unit cell geometry, material properties, and loading conditions.

36 MATERIALS SCIENCE↗

Tailoring the Properties of Calcium Aluminate Macro-Defect-Free Cements: From Brittle to Ductile Behavior

Cementitious macro-defect-free (MDF) materials can be obtained processing a mixture of a water-soluble polymer, namely polyvinyl alcohol (PVA) and the calcium aluminate cement. Based on process techniques typical of plastics, the MDF products can be obtained in complex forms and, thanks to the reduced water to cement ratio, w/c, their strength is very high compared to standard cements. However, these materials show both the sensitivity typical of water-soluble polymers, presenting a low glass transition temperature (T{sub g}), and the brittle behavior typical of ceramic materials and therefore low toughness. The aim of this work is to modify the mechanical properties of MDF cements by tuning their unacceptable brittle behavior. Water solution of carboxylate styrene butadiene rubber is used as reactive additives, and different mixture formulations are tested changing the rubber-like content both with respect to the fixed water content and PVA amounts. It is shown that the modifications of the organic phase imply small or negligible effects on paste formation of different recipes as detected by the torque–time responses. From flexural and fracture toughness tests, a wide range of responses was observed, and depending on the particular recipe, the mechanical behavior can be tuned from a brittle to a ductile behavior.

36 MATERIALS SCIENCE↗

Enhancing occupant behavior representation for interoperability between building information modeling and building energy modeling

Building Performance Simulation (BPS) has been adopted as an essential tool for designing, operating, and retrofitting buildings to optimize energy efficiency throughout the building life cycle. The Green Building XML (gbXML) schema facilitates seamless data exchange between Building Information Modeling (BIM) and Building Energy Modeling (BEM) software tools. However, limited occupant behavior (OB) representation in BIM often leads to inconsistent and inaccurate energy simulation in BEM software. This paper presents 154 systematic enhancements to the existing occupant behavior XML (obXML) schema v1.3.4, initially developed for standardizing OB representation for BEM, to address existing limitations and improve interoperability with BIM models. The enhancements encompass improved integration with BIM models through extended building representations and system operations, expanded support for advanced OB models with additional environmental parameters and mathematical capabilities, and implementation of a standardized model documentation framework. To facilitate seamless data transformation between gbXML and obXML schemas, we developed a publicly available gb-obXML Schema Converter. Three case studies demonstrate the enhanced schema’s capabilities: representation of building information using a two-story office building model, documentation of a window operation behavior model, and validation of the schema converter’s functionality. The enhanced obXML schema v1.4 enables sophisticated modeling of occupant-building interactions while maintaining consistency with industry-standard BIM schemas. The standardized documentation framework facilitates reproducibility and knowledge sharing in the OB research community, while the schema converter automates the integration of building information into OB simulation workflows. These enhancements establish a foundation for more accurate building performance simulation by supporting sophisticated representation of occupant behavior within the BIM-to-BEM simulation workflows.

Chung, Jihoon↗

Mapping local deformation behavior in single cell metal lattice structures

The deformation behavior of metal lattice structures is extremely complex and challenging to predict, especially since strain is not uniformly distributed throughout the structure. Understanding and predicting the failure behavior for these types of light-weighting structures is of great interest due to the excellent scaling of stiffness- and strength-to weight ratios they display. Therefore, there is a need to perform simplified experiments that probe unit cell mechanisms. This study reports on high resolution mapping of the heterogeneous structural response of single unit cells to the macro-scale loading condition. Two types of structures, known to show different stress-strain responses, were evaluated using synchrotron radiation micro-tomography while performing in-situ uniaxial compression tests to capture the local micro-strain deformation. These structures included the octet-truss, a stretch-dominated lattice, and the rhombic-dodecahedron, a bend-dominated lattice. The tomographic analysis showed that the stretch- and bend-dominated lattices exhibit different failure mechanisms and that the defects built into the structure cause a heterogeneous localized deformation response. Also shown here is a change in failure mode for stretch-dominated lattices, where there appears to be a transition from buckling to plastic yielding for samples with a relative density between 10 and 20%. In conclusion, the experimental results were also used to inform computational studies designed to predict the mesoscale deformation behavior of lattice structures. Here an equivalent continuum model and a finite element model were used to predict both local strain fields and mechanical behavior of lattices with different topologies.

36 MATERIALS SCIENCE↗

An electronic origin to the oscillatory segregation behavior in Ni-Cr and other BCC defects in FCC metals

Ni-Cr alloys exhibit oscillatory segregation behaviors near low index surfaces, in which the preferred segregation species changes from Ni in the first layer to Cr in the second layer. In many dilute-alloy systems, this oscillatory pattern is attributed to the elastic release of stresses in the local lattice around the segregating solute or impurity atom. These stresses are mostly thought to originate from mismatches in the atomic size of the solute and host atoms. In Ni-Cr alloys, however, an appreciable mismatch in atomic size is not present, leading to questions about the origins of the oscillatory behavior in this alloy. Therefore, using density functional theory, we have modeled the segregation of a single Cr atom in the and surfaces of FCC Ni, an alloy which exhibits this oscillatory behavior. Using Bader charge analysis, we show that the negative energy correlates directly with the amount of charge on the Cr atom. As Ni atoms strip valence charge from the Cr, the Cr contracts slightly in size. The greatest contraction and highest positive charge for the Cr occurs when it is in the second layer of the surface where the system exhibits the oscillating negative segregation energy. We then find that this behavior persists in other alloy systems (Ag-Nb, Cu-Cr, Pt-Nb, and Pt-V), which exhibit similar atomic radii and electronegativity differences between host and solute to Ni-Cr. These represent alloys in which the host metal exhibits an FCC ground-state structure while the solute metal exhibits a BCC ground-state structure.

36 MATERIALS SCIENCE↗

Investigation of rotating mode behavior in BWR out-of-phase limit cycle oscillations – Part 1: Reduced order model

Previous neutronic/thermal-hydraulic (TH) coupled numerical simulations using full-core TRACE/PARCS and SIMULATE-3K BWR models have shown evidence of a specific “rotating mode” behavior (steady rotation of the symmetry line, i.e. constant phase shift of approximately 90° between the first two azimuthal modes) in out-of-phase limit cycle oscillations, regardless of initial conditions and even if the first two azimuthal modes have different natural frequencies. This suggests a nonlinear coupling between these modes; otherwise, the phase shift between these modes would change at a constant rate during the limit cycle. The goal of the present work is to gain further insights on the rotating mode behavior using a simplified mathematical model which contains all of the important physics for this application while providing sufficient flexibility and simplicity to allow for in-depth understanding of the underlying phenomena. This was accomplished using a multi-channel, multi-modal reduced-order model, using a modification of the fixed pressure drop boundary condition to simulate channel coupling via the inlet and outlet plena, in order to destabilize the out-of-phase mode over the in-phase mode. Examination of the time-dependent solution of the nonlinear system showed a clear preference for rotating mode behavior in the four-channel model under stand-alone TH conditions and for conditions with weak neutronic feedback. Furthermore, when neutronic feedback was strengthened (i.e., larger reactivity feedback coefficients), the side-to-side mode (stationary symmetry line) was favored instead. Additional analyses using higher-fidelity numerical modeling, as well as a physical explanation for the rotating behavior seen in both sets of analyses, will be provided in a companion paper (“Part 2”).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Photochemical micromotor of eccentric core in isotropic hollow shell exhibiting multimodal motion behavior

Adaptive motion behavior in response to different environmental stimuli is ubiquitous in biology and enables creatures to achieve a diversity of complex tasks, but not typically observed in synthetic systems. In this report we propose a design of chemically-powered micromotors actively performing multimodal motion behaviors with the external stimulus changing. They have isotropic outer surfaces but inherent inner mass asymmetry, such as eccentric core-in-hollow shell TiO 2 (E-TiO 2 ) microspheres. Their motion behavior can be spontaneously transformed among random Brownian propulsion (stochastic walk), negative phototaxis (moving against incident light), and negative photogravitaxis (moving against gravity) when the light intensity, illumination direction, or fuel concentration change. At a low light intensity and/or low H 2 O 2 fuel concentration, the E-TiO 2 micromotors perform directional movement away from light based on the dominated diffusiophoresis by the photocatalytic reaction over the isotropic shell. With the increase of light intensity or fuel concentration, there are more and more photons or fuels to reach the eccentric core and the contribution of the photocatalytic reaction over it to the diffusiophoresis gradually becomes dominant. In this case, the E-TiO 2 micromotors perform stochastic walks near the substrate due to their Brownian random rotational reorientations. The micromotors are single component, have low density, and can be synthesized in a large scale and at a low cost. This work will facilitate the development of multifunctional micro/nanomotors with varying behaviors and functions depending on environmental cues, e.g., enabling efficient search and delivery operations.

36 MATERIALS SCIENCE↗

Assessment of asymmetrical rheological behavior of cementitious material for 3D printing application

Among the multiplicity of construction methods based on digital concrete technologies, the most common one remains the extrusion-based process. Two opposed strategies have been developed in order to print a concrete structure using extrusion: the first one emphasizes the ease-of-pumping of material which is massively accelerated to ensure the stability of the structure, when leaving the extruder; the second one uses a firm mortar able to sustain the own weight of the structure during the printing, without requiring any accelerator. Both methods have their pros and cons and don't imply the same rheological requirements. Most of the recent literature assumes that both types of material (fluid and firm) follow the same trend of behavior and obey to the Von Mises plasticity criterion, while it has been demonstrated that firm cement-based materials can exhibit a pressure dependent behavior. Moreover, rheological study is mainly based on shear rheometry, while the extrusion-based printing process requires a description of the rheological behavior in compression to predict the global stability of the structure, and in tension to prevent cracks formation. This paper investigates the fresh behavior of cement paste and mortar under different solicitations (shear, compression, tension) for different water to cement ratios. The analysis of these results allows to describe the transition between a ductile fluid material (symmetrical in tension and compression) that obeys to a Von Mises plasticity criterion, and a brittle and firm material (asymmetrical in tension and compression) that requires a pressure dependent plasticity criterion in order to predict its strength under a given solicitation.

36 MATERIALS SCIENCE↗

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗

From observation to replication: machine-learning-driven quantification and replication of fine-scale fish kinematics and behavior

Long-term quantification of fish behavior is essential for aquatic ecology, wildlife telemetry, and biomechanical device development. However, the observation duration required to obtain reliable behavioral and kinematic metrics remains unclear, and few tools exist to physically reproduce natural swimming motion for controlled experimentation. We address these challenges by developing a generalizable framework that models behavioral reliability (Spearman–Brown reliability index) as a function of observation duration and derives metric-specific monitoring thresholds. Using juvenile white sturgeon (Acipenser transmontanus) as a case study, we demonstrate that the minimum duration needed for reliable estimates varies substantially across kinematic features: to exceed a reliability of 0.8, total distance traveled requires 12 days, average curvature (mm?¹) 15 days, tail-beat frequency (Hz) 8 days, and average speed (body length/s) 17 days. We further bridge digital analysis and physical testing by developing a hardware-in-the-loop simulator that reconstructs machine-learning-derived swimming kinematics with high fidelity (correlation coefficient 0.98–0.99, RMSE 1.22–1.27 mm over a 5-minute segment). This platform enables realistic, repeatable motion stimuli for evaluating aquatic sensing technologies and bio-integrated devices under controlled conditions. Together, these contributions provide a scalable approach for designing long-term behavioral studies and a data-driven connection between ecological observation and robotic experimentation.

Hwang, SungJoo↗

Rare earth element behavior during groundwater–seawater mixing along the Kona Coast of Hawaii

Groundwater and seawater samples were collected from nearshore wells and offshore along the Kona Coast of the Big Island of Hawaii to investigate rare earth element (REE) behavior in local subterranean estuaries. Previously we saw that submarine groundwater discharge (SGD) is the predominant flux of terrestrial waters to the coastal ocean along the arid Kona Coast of Hawaii. Groundwater and seawater samples were filtered through 0.45 μm and 0.02 μm pore-size filters to evaluate the importance of colloidal and soluble (i.e., truly dissolved ionic species and/or low molecular weight [LMW] colloids) fractions of the REEs in the local subterranean estuaries. Mixing experiments using groundwater collected immediately down gradient from a wastewater treatment facility (WWTF) proximal to the Kaloko-Hanokohau National Historic Park, and more “pristine” groundwater from a well constructed in a lava tube at Kiholo Bay, were conducted with local seawater to study the effect of solution composition (i.e., pH, salinity) on the concentrations and fractionation behavior of the REEs as groundwater mixes with seawater in Kona Coast subterranean estuaries. The mixed waters were also filtered through 0.45 or 0.02 μm filters to ascertain the behavior of colloidal and soluble fractions of the REEs across the salinity gradient in each mixing experiment. Concentrations of the REEs were statistically identical (two-tailed Student t-test, 95% confidence) between the sequentially filtered sample aliquots, indicating that the REEs occur as dissolved ionic species and/or LMW colloids in Kona Coast groundwaters. The mixing experiments revealed that the REEs are released to solution from suspended particles or colloids when Kona Coast groundwater waters mix with local seawater. The order of release that accompanies increasing pH and salinity follows light REE (LREE) > middle REE (MREE) > heavy REE (HREE). Release of REEs in the mixing experiments is driven by decreases in the free metal ion activity in solution and the concomitant increase in the amount of each REE that occurs in solution as dicarbonato complexes [i.e., Ln(CO 3 ) 2 - ] as pH increases across the salinity gradient. Input-normalized REE patterns of Kona Coast groundwater and coastal seawater are nearly identical and relatively flat compared to North Pacific seawater, indicating that SGD is the chief source of these trace elements to the ocean along the Kona Coast. Additionally, REE concentrations of the coastal seawater are between 10 and 50 times higher than previously reported open-ocean seawater values from the North Pacific, further demonstrating the importance of SGD fluxes of REEs to these coastal waters. Altogether, these observations indicate that large-scale removal of REEs, which characterizes the behavior of REEs in the low salinity reaches of many surface estuaries, is not a feature of the subterranean estuary along the Kona Coast. A large positive gadolinium (Gd) anomaly characterizes groundwater from the vicinity of the WWTF. The positive Gd anomaly can be traced to the coastal ocean, providing further evidence of the impact of SGD on the coastal waters. Estimates of the SGD fluxes of the REEs to the coastal ocean along the Kona Coast (i.e., 1.3–2.6 mmol Nd day -1 ) are similar to recent estimates of SGD fluxes of REEs along Florida’s east coast and to Rhode Island Sound, all of which points to the importance of SGD as significant flux of REEs to the coastal ocean.

58 GEOSCIENCES↗