Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpretable models”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Delving into the depths of NGC 3783 with XRISM

We present our study of the X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the Seyfert-1 galaxy NGC 3783. XRISM’s Resolve microcalorimeter has enabled, for the first time, a detailed characterization of the highly ionized outflows in this active galactic nucleus. Our analysis constrains their outflow and turbulent velocities, along with their ionization parameter (ξ) and column density (N H ). The high-resolution Resolve spectrum reveals a distinct series of Fe absorption lines between 6.4 and 7.8 keV, ranging from Fe XVIII to Fe XXVI. At lower energies (1.8−3.3 keV), absorption features from Si, S, and Ar are also detected. Our spectroscopy and photoionization modeling of the time-averaged Resolve spectrum uncovers six outflow components, five of which exhibit relatively narrow absorption lines with outflow velocities ranging from 560 to 1170 km s −1 . In addition, a broad absorption feature is detected, which is consistent with Fe XXVI outflowing at 14 300 km s −1 (0.05 c). The kinetic luminosity of this component is 0.8−3% of the bolometric luminosity. Our analysis of the Resolve spectrum shows that more highly ionized absorption lines are intrinsically broader than those of lower-ionization species, indicating that the turbulent velocity of the six outflow components (ranging from 0 to 3500 km s −1 ) increases with ξ. Furthermore, we find that the column density (N H ) of the outflows generally declines with the ionization parameter up to log ξ = 3.2 but rises beyond this point, suggesting a complex ionization structure. The absorption profile of the Fe XXV resonance line is intriguingly similar to UV absorption lines (Lyα and C IV) observed by the Hubble Space Telescope, from which we infer that the outflows are clumpy in nature. Our XRISM/Resolve results from lower- and higher-ionization regimes support a “hybrid wind” scenario in which the observed outflows have multiple origins and driving mechanisms. We explore various interpretations of our findings within active galactic nucleus wind models.

X-rays: galaxies↗
DOE OSTI2022

A Multi-Phase Modeling Framework Suitable for Dynamic Applications

Under dynamic loading conditions and the associated extreme conditions many metals will undergo phase transformations. The change in crystal structure associated with solid–solid phase transformations can significantly alter the subsequent mechanical response of the material. For the interpretation of experiments involving dynamic loading it is beneficial to have a modeling framework that captures key features of the material response while remaining relatively simple. We introduce a candidate framework and apply it to the metal tin to highlight a range of behaviors that are captured by the model. We also discuss potential extensions to capture additional behaviors that could be important for certain materials and loading scenarios. The model is useful for analysis of results from dynamic experiments and offers a point of departure for more complex model formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Challenges of standard halo models in constraining galaxy properties from cosmic infrared background anisotropies

The halo model, combined with halo occupation distribution (HOD) prescriptions, is widely used to interpret cosmic infrared background (CIB) anisotropies and extract physical information about star-forming galaxies and their connection to large-scale structures. Recent CIB-specific implementations of the halo model have adopted more physical parameterizations. However, the extent to which these models can reliably recover meaningful physical parameters remains uncertain. We assessed whether the current parameterization of CIB halo models is sufficient to recover astrophysical quantities, such as star formation efficiency, η(M h , z), and halo mass at which the peak of star formation efficiency occurs, M max , when fit to mock data. We also assessed whether discrepancies arise from assumptions about galaxy emission (the HOD ingredients) or from more fundamental components in the halo model, such as bias and matter clustering. We fit the M21 CIB HOD model, implemented within the halo model framework, to mock CIB power spectra and star formation rate density (SFRD) data generated from the SIDES-Uchuu simulation, and compared the best-fit parameters to the known simulation inputs. We then repeated the analysis using a simplified version of the simulation (SSU), explicitly designed to match the HOD assumptions. A detailed comparison of model and simulation outputs was carried out to trace the origin of observed discrepancies. While the M21 HOD model provides a good fit to the mock data, it failed to recover the intrinsic parameters accurately, particularly the halo mass at which star formation efficiency peaks. This mismatch persists even when fitting data generated with the same model assumptions. We find strong agreement (within 5%) in the emission-related components (SFRD, emissivity), but observe a scale- and redshift-dependent offset exceeding 20% in the two-halo term of the CIB power spectrum. This likely arises from limitations in the treatment of halo bias and matter clustering within the linear approximation. Additionally, incorporating scatter in the SFR–halo mass relation and the spectral energy distribution (SED) templates significantly affects the shot noise (∼50%), but has only a modest impact (less than 10%) on the clustered component. These results suggest that recovering physical parameters from CIB clustering requires improvements to the cosmological ingredients of the halo model framework, such as adopting scale-dependent halo bias and nonlinear matter power spectra in addition to careful modeling of emission physics.

cosmic background radiation↗
DOE OSTI2026

Benchtop Autonomous Electrochemical Characterization System for Combinatorial Thin-Film Solid Oxide Electrodes

The design of materials for electrochemical energy conversion is complicated by a vast search space of candidate materials and multifaceted property requirements: multicarrier conductivity, stability, and catalytic activity are all necessary but rarely intersect. Although self-driving laboratories are rapidly rising to address such material optimization problems, the required infrastructure for integrated, large-scale robotic facilities can be cost-prohibitive. Here we develop and evaluate a closed-loop measurement system for efficient screening of proton-conducting oxide electrodes for ceramic fuel cells and electrolyzers, building on top of an existing benchtop instrument and integrating techniques for rapid impedance measurement and automated analysis. This system exemplifies a “minimum viable” self-driving implementation that can deliver substantial benefits with relatively simple infrastructure. Combinatorial thin-film microelectrode libraries are characterized with a recently developed joint time-domain and frequency-domain impedance measurement technique, which provides an order-of-magnitude acceleration relative to conventional impedance spectroscopy. The distribution of relaxation times is extracted from impedance data and analyzed without human intervention. These results feed an active learning and Bayesian optimization process that learns to predict electrochemical impedance as a function of material composition, measurement temperature, oxygen partial pressure, and electrical bias, which further reduces the screening time by tenfold with optimized experimental sequences. We apply this system to Ba⁡(Co,Fe,Zr,Y)⁢O 3−𝛿 combinatorial libraries and evaluate its effectiveness for learning material property trends and optimizing expensive-to-evaluate properties such as activation energy. This offers insights into key methodological aspects of practical autonomous experimentation, including surrogate model validation, cost-aware acquisition functions, and high-throughput data interpretation. Our results demonstrate the efficacy of the system for rapidly gathering information, but also highlight real-world experimental challenges of thin-film degradation and numerical instability in surrogate models.

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

A Visual Designer of Layer-wise Relevance Propagation Models

Layer-wise Relevance Propagation (LRP) is an emerging and widely-used method for interpreting the prediction results of convolutional neural networks (CNN). LRP developers often select and employ different relevance backpropagation rules and parameters, to compute relevance scores on input images. However, there exists no obvious solution to define a "best" LRP model. A satisfied model is highly reliant on pertinent images and designers' goals. We develop a visual model designer, named as VisLRPDesigner, to overcome the challenges in the design and use of LRP models. Various LRP rules are unified into an integrated framework with an intuitive workflow of parameter setup. VisLRPDesigner thus allows users to interactively configure and compare LRP models. It also facilitates relevance-based visual analysis with two important functions: relevance-based pixel flipping and neuron ablation. Several use cases illustrate the benefits of VisLRPDesigner. The usability and limitation of the visual designer is evaluated by LRP users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A consistent and conservative Phase-Field method for multiphase incompressible flows

In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Inference of experimental radial impurity transport on Alcator C-Mod: Bayesian parameter estimation and model selection

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. Here, we find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Inference of Experimental Radial Impurity Transport on Alcator C-Mod: Bayesian Parameter Estimation and Model Selection.

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. We find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Towards the next generation of Geospatial Artificial Intelligence

Geospatial Artificial Intelligence (GeoAI), as the integration of geospatial studies and AI, has become one of the fastest-developing research directions in spatial data science and geography. This rapid change in the field calls for a deeper understanding of the recent developments and envision where the field is going in the near future. In this work, we provide a quantitative analysis of the GeoAI literature from the spatial, temporal, and semantic aspects. We briefly discuss the history of AI and GeoAI by highlighting some pioneering work. Then we discuss the current landscape of GeoAI by selecting five representative subdomains including remote sensing, urban computing, Earth system science, cartography, and geospatial semantics. Finally, we highlight several unique future research directions of GeoAI which are classified into two groups: GeoAI method development challenges and GeoAI Ethics challenges. Topics include heterogeneity-aware GeoAI, knowledge-guided GeoAI, spatial representation learning, geo-foundation models, fairness-aware GeoAI, privacy-aware GeoAI, as well as interpretable and explainable GeoAI. We hope our review of GeoAI’s past, present, and future is comprehensive and can enlighten the next generation of GeoAI research.

58 GEOSCIENCES↗
DOE OSTI2025

Sequential Infiltration Synthesis of Multinary Metal Oxides: The Role of Precursor Diffusion and Reactivity

Understanding volatile metal–organic precursor transport in polymer matrices during sequential infiltration synthesis (SIS) is critical to the design of hybrid materials of tailored composition and uniform structure. Here, we investigate the diffusion and reaction behavior of diethylzinc (DEZ) in poly(methyl methacrylate) (PMMA) matrices with variable indium and zinc incorporation. Unlike conventional SIS where precursors adduct or react directly with the polymer backbone, DEZ shows minimal interaction with PMMA. However, DEZ may interact strongly with In–OH or Zn–OH moieties that are previously incorporated into a PMMA thick film. Using spectroscopic ellipsometry, X-ray photoelectron spectroscopy depth profiling, and EDX line scans, we quantify Zn incorporation as a function of exposure time, polymer thickness, and previous infiltration steps. The Zn distribution evolves from surface-localized to deep within the ∼1 μm thick film with increasing DEZ exposures up to 30 min. A reaction-diffusion model incorporating a hindering factor to capture diffusivity decay due to product accumulation reproduces the Zn depth profile and yields physically meaningful parameters to guide future growth. Furthermore, this study provides a framework for interpreting metal–organic precursor infiltration behavior and highlights the utility of model-guided SIS process design.

diethylzinc↗
DOE OSTI2023

The star-formation history in the last 10 billion years from CIB cross-correlations

ABSTRACT The cosmic infrared background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range z ≲ 2 to measure the bias-weighted star-formation rate (SFR) density 〈bρSFR〉 as a function of time in a model independent way. This quantity is complementary to direct measurements of the SFR density ρSFR, giving a higher weight to more massive haloes, and thus provides additional information to constrain the physical properties of star formation. Using cross-correlations of the CIB with galaxies from the DESI Legacy Survey and the extended Baryon Oscillation Spectroscopic Survey, we obtain high signal-to-noise ratio measurements of 〈bρSFR〉, which we then use to place constraints on halo-based models of the star-formation history. We fit halo-based SFR models to our data and compare the recovered ρSFR with direct measurements of this quantity. We find a qualitatively good agreement between both independent data sets, although the details depend on the specific halo model assumed. This constitutes a useful robustness test for the physical interpretation of the CIB, and reinforces the role of CIB maps as valuable astrophysical probes of the large-scale structure. We report our measurements of 〈bρSFR〉 as well as a thorough account of their statistical uncertainties, which can be used to constrain star-formation models in combination with other data.

Astronomy & Astrophysics↗
DOE OSTI2024

Photoproduction of K0K0-bar with the GlueX Experiment

In this dissertation, we study photoproduction of K0 ?K 0 through the KSKLp and KSKSp ?nal states with the GlueX Phase-I (GlueX-I) data set. We measure the spin-density matrix elements (SDMEs) and di?erential cross section of ?(1020) ! KSKL to better understand photoproduction of light vector mesons. A mass independent Partial Wave Analysis (PWA) of the KSKL and KSKS meson spectrum below 2 GeV is also undertaken to study the meson spectrum with JPC = even++ and odd??. The ?(1020) di?erential cross section is measured at E = 8:2 ? 8:8 GeV and ?t = 0:15 ? 1:00 GeV2. The di?erential cross section is well described by an exponential decay with slope 4:44?0:01 GeV?2 and the integrated cross section is determined to be 295:7?0:4 nb, only statistical uncertainties are quoted. Both measurements are consistent and far more precise than the previous measurement by Ballam et al. [1]. The ?(1020) SDMEs were measured in nine bins of ?t in the same range. At low ?t, we ?nd the data were consistent with s-channel helicity conservation, SCHC, i. e. the only non-zero SDMEs were ?1 1?1 = ?Im(?2 1?1) = 1=2. At higher ?t, the SDMEs deviate from SCHC and the measurements indicate this is due to natural parity exchange since we observe that the SDMEs are consistent with zero unnatural exchange. We also ?find that the contribution from helicity double-flip amplitudes is consistent with zero. The measured SDMEs are in poor agreement with theoretical predictions put forward by JPAC [35]. These measurements will serve as input to re?ne models of production processes, which will be essential for the interpretation of possible signals of exotic mesons in GlueX. Analysis of the spectrum above the ?(1020) indicates the presence of at least three particles at ? 1:50, ? 1:75, and 2:20 GeV. We employ simple parametrizations of the resonance line shape based on relativistic Breit-Wigner functions considering cases with and without interference. The resonance at ? 1:50 GeV is not identifi?ed as any specifi?c resonance, but could be due to interference between the ?(1450) and !(1420). Our models favor the resonance at ~1:75 GeV as the X(1750) rather than the ?(1680). We ?find that introducing a third Breit-Winger into the models improves the ?t quality with the ?2/ndf going from 1.75 to 1.59 and 1.41 to 1.18 for models without and with interference respectively. A PWA of the KSKL spectrum indicates that the spectrum is consistent with exclusively spin-1 contribution up to ? 1:6 GeV. Above ? 1:6 GeV we ?find evidence for a small spin-3 contribution. We also ?nd that the spectrum predominantly positive reflectivity up to ? 1:6 GeV, after which the spectrum becomes a nearly equal mix of both. The KSKS system shows a rich spectrum with multiple resonance structures but is more statistically limited than the KSKL system. Modelling the line shape we ?find evidence at greater than 10? level in favor of a resonance at ? 1:75 GeV. The model parameters suggest this resonance is the f0(1710). The PWA of the KSKS system was inconclusive. However, using a small set of amplitudes suggests that the spin-2 contribution to the spectrum is primarily in the range 1:2?1:6 GeV, where f2(1270), a2(1320) and f02 (1525) are expected.

Linera, Gabriel Rodriguez↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

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

A Mass‐Conserving‐Perceptron for Machine‐Learning‐Based Modeling of Geoscientific Systems

Although decades of effort have been devoted to building Physical-Conceptual (PC) models for predicting the time-series evolution of geoscientific systems, recent work shows that Machine Learning (ML) based Gated Recurrent Neural Network technology can be used to develop models that are much more accurate. However, the difficulty of extracting physical understanding from ML-based models complicates their utility for enhancing scientific knowledge regarding system structure and function. Here, we propose a physically interpretable Mass-Conserving-Perceptron (MCP) as a way to bridge the gap between PC-based and ML-based modeling approaches. The MCP exploits the inherent isomorphism between the directed graph structures underlying both PC models and GRNNs to explicitly represent the mass-conserving nature of physical processes while enabling the functional nature of such processes to be directly learned (in an interpretable manner) from available data using off-the-shelf ML technology. As a proof of concept, we investigate the functional expressivity (capacity) of the MCP, explore its ability to parsimoniously represent the rainfall-runoff (RR) dynamics of the Leaf River Basin, and demonstrate its utility for scientific hypothesis testing. To conclude, we discuss extensions of the concept to enable ML-based physical-conceptual representation of the coupled nature of mass-energy-information flows through geoscientific systems.

58 GEOSCIENCES↗
DOE OSTI2021

Full Vorticity Budget of the Arabian Sea from a 0.1° Ocean Model: Sverdrup Dynamics, Rossby Waves, and Nonlinear Eddy Effects

Abstract The Arabian Sea, influenced by the Indian monsoon, has many unique features, including its basin-scale seasonally reversing surface circulation and the Great Whirl, a seasonal anticyclonic system appearing during the southwest monsoon close to the western boundary. To establish a comprehensive dynamical picture of the Arabian Sea, we utilize numerical model output and design a full vorticity budget that includes a fully decomposed nonlinear term. The ocean general circulation model has 0.1° resolution and is mesoscale eddy-resolving in the region. In the western boundary current system, we highlight the role of nonlinear eddies in the life cycle of the Great Whirl. The nonlinear eddy term is of leading-order importance in this feature’s vorticity balance. Specifically, it contributes to the Great Whirl’s persistence in boreal fall after the weakening of the southwesterly winds. In the open ocean, Sverdrup dynamics and annual Rossby waves are found to dominate the vorticity balance; the latter is considered as a key factor in the formation of the Great Whirl and the seasonal reversal of the western boundary current. In addition, we discuss different forms of vertically integrated vorticity equations in the model and argue that the bottom pressure torque term can be interpreted analogously as friction in the western boundary and vortex stretching in the open ocean.

58 GEOSCIENCES↗
DOE OSTI2022

Development of Cr, Se, U, Sb, and Te isotopes as indicators of redox reactions, contaminant fate, and contaminant transport in aqueous systems: A review

Cr, Se, U, Sb, and Te are toxic, redox-active elements that are more mobile and environmentally problematic in their oxidized forms, and less mobile and bioavailable in their reduced forms. This chapter reviews the development of Cr, Se, U, Sb, and Te isotope ratio measurements as new indicators of redox reactions and contaminant migration. Reliable analytical methods exist, but are still evolving. Understanding of isotopic fractionation induced by various (bio)geochemical processes has been explored in dozens of publications, yet is far from complete: Reduction reactions, the major driver of isotopic variation, have been relatively well studied. However, the magnitude of fractionation is variable and the systematics of that variation are still being explored. Isotopic fractionation induced by oxidation reactions is not well understood. Non-redox reactions, which involve smaller changes in bonding of these elements, tend to induce less isotopic fractionation, but can nonetheless cause significant isotopic shifts. Field applications of Cr, Se, U isotope ratios have demonstrated that they are useful as indicators of reduction in natural systems. A few studies suggest they are also useful as indicators of oxidation and contaminant sources. The physical and chemical complexity of groundwater systems hinders accurate quantitative interpretation of Cr, Se, U isotope data using simple models. Numerical models have been developed that capture the behavior of complex, coupled systems and enable the most effective extraction of information from field data sets.

58 GEOSCIENCES↗