Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability analysis”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Kinematic Variables and Feature Engineering for Particle Phenomenology

Kinematic variables have been playing an important role in collider phenomenology, as they expedite discoveries of new particles by separating signal events from unwanted background events and allow for measurements of particle properties such as masses, couplings, spins, etc. For the past 10 years, an enormous number of kinematic variables have been designed and proposed, primarily for the experiments at the Large Hadron Collider, allowing for a drastic reduction of high-dimensional experimental data to lower-dimensional observables, from which one can readily extract underlying features of phase space and develop better-optimized data-analysis strategies. We review these recent developments in the area of phase space kinematics, summarizing the new kinematic variables with important phenomenological implications and physics applications. We also review recently proposed analysis methods and techniques specifically designed to leverage the new kinematic variables. As machine learning is nowadays percolating through many fields of particle physics including collider phenomenology, we discuss the interconnection and mutual complementarity of kinematic variables and machine learning techniques. We finally discuss how the utilization of kinematic variables originally developed for colliders can be extended to other high-energy physics experiments including neutrino experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Physics and chemistry from parsimonious representations: image analysis via invariant variational autoencoders

Electron, optical, and scanning probe microscopy methods are generating ever increasing volume of image data containing information on atomic and mesoscale structures and functionalities. This necessitates the development of the machine learning methods for discovery of physical and chemical phenomena from the data, such as manifestations of symmetry breaking phenomena in electron and scanning tunneling microscopy images, or variability of the nanoparticles. Variational autoencoders (VAEs) are emerging as a powerful paradigm for the unsupervised data analysis, allowing to disentangle the factors of variability and discover optimal parsimonious representation. Here, we summarize recent developments in VAEs, covering the basic principles and intuition behind the VAEs. The invariant VAEs are introduced as an approach to accommodate scale and translation invariances present in imaging data and separate known factors of variations from the ones to be discovered. We further describe the opportunities enabled by the control over VAE architecture, including conditional, semi-supervised, and joint VAEs. Several case studies of VAE applications for toy models and experimental datasets in Scanning Transmission Electron Microscopy are discussed, emphasizing the deep connection between VAE and basic physical principles. Python codes and datasets discussed in this article are available at https://github.com/saimani5/VAE-tutorials and can be used by researchers as an application guide when applying these to their own datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Decoding the shift-invariant data: applications for band-excitation scanning probe microscopy *

A shift-invariant variational autoencoder (shift-VAE) is developed as an unsupervised method for the analysis of spectral data in the presence of shifts along the parameter axis, disentangling the physically-relevant shifts from other latent variables. Using synthetic data sets, we show that the shift-VAE latent variables closely match the ground truth parameters. The shift VAE is extended towards the analysis of band-excitation piezoresponse force microscopy data, disentangling the resonance frequency shifts from the peak shape parameters in a model-free unsupervised manner. The extensions of this approach towards denoising of data and model-free dimensionality reduction in imaging and spectroscopic data are further demonstrated. This approach is universal and can also be extended to analysis of x-ray diffraction, photoluminescence, Raman spectra, and other data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Novel Process for Carbon Dioxide Conversion to Fuel

In this project, TDA developed a new mixed metal oxide-based sorbent that converts CO2 (captured from coal fired power plants) to CO, which can then be combined with renewable H2 generated by water electrolysis or H2 from steam methane reforming to produce different liquid fuels. TDA’s absorbent-based CO2 conversion process uses a redox process, which splits the catalytic reforming of methane with CO2 reaction into two stages: CO2 reduction to CO and CH4 reforming into H2 and CO which eliminates the equilibrium limitations. The CO produced in the two-stage reactor system can then be further reacted with renewable H2 to produce methanol, naphtha, diesel, or gasoline. We worked with the Advanced Power & Energy Program (APEP) of University of California, Irvine (UCI) to design and develop the liquid fuel synthesis process that is built around this new material. We demonstrated the techno-economic viability of the new sorbent based redox process to convert CO2 into synthesis gas by demonstrating continuous carbon dioxide reduction in a prototype test system for over 585 hours while converting up to 10 kg CO2/day. With the successful completion of the R&D effort, the technology is now ready for a larger pilot-scale demonstration and the technology readiness has been raised from TRL 3 to TRL 5. In collaboration with UCI, we completed a high-fidelity process design and economic analysis. The required sale price (RSP) for gasoline (Case 1 NG-MTG) is $4.91/gal and naphtha and diesel (Case 2 NG-FT) are $4.23/gal and $6.07/gal, respectively, on a 2011 dollar basis. To put these costs in perspective, the California prices in current dollars (with its strict specifications) for regular grade gasoline from last year to current year have varied from a low of $3.10/gal in January 2021 to a high of $5.76/gal in March 2022, while prices for diesel from last year to current year have varied from a low of $3.40/gal in January 2021 to a high of $6.41/gal in May 2022 according to the U.S. Energy Information Agency data. It should be noted that the gasoline and diesel produced by these designs of Case 1 and Case 2 would be of very high quality and both nitrogen and sulfur free. These RSPs are based on a cost of imported electricity of $64/MWh based on the low-end current wind generated electricity cost (Genevieve 2011). This cost is by far the largest component of the variable costs used in computing the RSPs. A sensitivity analysis of these RSPs to the cost of imported electricity shows that the cost of the imported electricity has a significant effect on the RSPs. The life cycle analysis (LCA) shows that the total cradle-to-gate CO2 emissions for both liquid fuels (diesel and gasoline) were negative, indicating that overall more CO2 is consumed than released during production of the fuel from CO2 feed stack for both cases (-296 kgCO2 per MT gasoline for Case 1 and -705 kgCO2 per MT diesel for Case 2). On a cradle-to-grave comparison, the use of diesel produced using TDA’s process (2,457 kgCO2 per MT diesel) would release 37.6% less CO2 compared to petroleum based diesel (3,937 kgCO2 per MT diesel) while the use of gasoline produced using TDA’s process (2,792 kgCO2 per MT gasoline) would release 29.3% less CO2 compared to petroleum based gasoline (3,946 kgCO2 per MT gasoline). With the successful completion of the R&D effort, the technology is now ready for a bench-scale demonstration and the technology readiness has been raised from TRL 3 (Analytical and experimental critical function and/or characteristic proof of concept) to TRL 5 (Laboratory scale similar system validation in relevant environment).

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

A compendium and meta-analysis of flatband potentials for TiO2, ZnO, and SnO2 semiconductors in aqueous media

Semiconductor/electrolyte interfaces are of great interest to numerous scientific fields including renewable energy, (photo)electrochemistry, and energy storage. The semiconductor flatband potential is a key parameter in locating the conduction band minimum or valence band maximum of the semiconductor material in electrolyte. Despite its importance for quantifying the energetic location of the semiconductor bands, literature reports for the same material demonstrate significant variability in the flatband potential. In this compendium and meta-analysis, reported flatband potentials of the common semiconductor materials TiO2, SnO2, and ZnO in aqueous electrolyte were compiled and assessed to quantify the spread in literature flatband potentials as well as determine the factors that lead to the significant spread. For TiO2, SnO2, and ZnO, literature flatband potentials referenced to the reversible hydrogen electrode span a range of nearly 2 V each. Flatband potential tabulations were separated by variables such as the solution pH, the crystalline polymorph, the crystal facet, the morphology, and the dimensions or combinations of these variables to assess the factors that contribute to the observed spread. Important and surprising findings from these categorizations are summarized: (1) Even for the narrowest categorizations, the spread in flatband potential is still large. (2) Flatband potentials of TiO2 and SnO2 follow the expected Nernstian dependence with solution pH. ZnO materials deviate from this Nernstian dependence. (3) In the aggregate, there is no statistically significant difference in the reported flatband potentials of anatase and rutile TiO2. Single crystal tabulations were the only distributions to have statistically significant differences in the flatband potential between anatase and rutile TiO2. (4) Anatase TiO2 materials with a nanotube morphology appear to have a +400 mV difference in mean flatband potential compared to all other morphologies, but we argue that this is likely due to widespread misuse of the Mott–Schottky analysis. Other interesting findings are revealed within the spread of literature flatband potentials, and possible explanations are provided to generate discussion. We also briefly review and discuss common techniques that were used to determine the flatband potential and the pitfalls/criticisms of these techniques. Last, we discuss some ways in which future research on the determination of the flatband potential can be performed to improve the reliability of reported values and the quality of the work. In total, the results from this meta-analysis suggest multiple factors can affect the measured flatband potential and that an abundance of caution should be applied when attempting to quantify the flatband potential of complex or nanostructured systems.

Patel, Milan Y. (ORCID:0000000291948141)↗
DOE OSTI2021

Assessment of small mechanical wastewater treatment plants: Relative life cycle environmental impacts of construction and operations

Many slow growing and shrinking rural communities struggle with aging or inadequate wastewater treatment plants (WWTPs), and face challenges in constructing and operating such facilities. Although existing literature has provided insight into the environmental sustainability of large facilities, including both the construction and operational phases, these studies have not examined small, rural facilities treating less than 7,000 m 3 /d (1.8 MGD) of wastewater in adequate depth and breadth. In this study, a detailed inventory of the construction and operational data for 16 case studies of small WWTPs was developed to elucidate their environmental life cycle impacts. Conventional LCA framework was followed. The results show that the environmental impacts of both the construction and operational phases are considerable. Operational impacts are highly related to energy usage. Improving energy efficiency of a plant may reduce the environmental impacts related to operations. Construction impacts can vary considerably between facilities. Process-related factors (e.g., concrete and reinforcing steel used in basins) are typically sized using the design flow; thus much of the variability in construction impacts among plants stems from the non-process related infrastructure. Multiple regression analysis was used as an exploratory tool to identify which non-process related plant aspects contribute to the variable environmental impact of small WWTPs. These factors include aluminum, cast iron, and the capacity utilization ratio (defined as the ratio of average flow to design flow). Furthermore, industry practitioners should consider these factors when aiming to reduce the environmental impacts of a small WWTP related to construction. Scenario sensitivity analyses found that the environmental impact of construction became smaller with longer design life, and the end-of-life consideration does not heavily influence the environmental sustainability of a WWTP.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗
DOE OSTI2022

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗
DOE OSTI2025

Understanding the impacts of inorganic species in woody biomass for preprocessing and pyrolysis–A review

Woody biomass represents an abundant resource for sustainable biofuels, biochemicals, and bioproducts. Technologies for converting woody biomass have been established for decades, and research consistently highlights the critical role of inorganic species and ash plays in feedstock handling and conversion processes, including equipment plugging, corrosion, and catalyst deactivation. A thorough understanding of the variability, transport behavior, and downstream impact of inorganic species in woody biomass is essential for defining feedstock quality specifications and developing effective management strategies for conversion processes. This review compiles critical information in five main sections: 1) inorganic species concentration in woody biomass, based on anatomical fractions and their sources of variability; 2) technique features for quantifying inorganic elemental chemical analysis; 3) impacts of inorganic species on biomass preprocessing; 4) impacts of inorganic species on pyrolysis, and 5) mitigation strategies. Additionally, this review explores future challenges and opportunities in addressing the impacts of inorganic species on biomass quality. These insights aim to support the sustainable development of the biomass-to-bioenergy pipeline and ensure high-quality lignocellulosic feedstocks for efficient downstream conversions. The findings offer valuable guidance to policy makers, industry stakeholders, and researchers in developing effective strategies for managing inorganic species in woody biomass and fostering the sustainable processes for lignocellulosic biorefineries.

09 BIOMASS FUELS↗
DOE OSTI2021

Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors

Cloud drop number concentrations (N d ) over the western North Atlantic Ocean (WNAO) are generally highest during the winter (DJF) and lowest in summer (JJA), in contrast to aerosol proxy variables (aerosoloptical depth, aerosol index, surface aerosol mass concentrations, surface cloud condensation nuclei (CCN) concentrations) that generally peak inspring (MAM) and JJA with minima in DJF. Using aircraft, satellite remote sensing, ground-based in situ measurement data, and reanalysis data, we characterize factors explaining the divergent seasonal cycles and furthermore probe into factors influencing N d on seasonal timescales. The results can be summarized well by features most pronounced in DJF, including features associated with cold-air outbreak (CAO) conditions such as enhanced values of CAO index, planetary boundary layer height (PBLH),low-level liquid cloud fraction, and cloud-top height, in addition to winds aligned with continental outflow. Data sorted into high- and low-N d days in each season, especially in DJF, revealed that all of these conditions were enhanced on the high-N d days, including reduced sea level pressure and stronger wind speeds. Although aerosols may be more abundant in MAM and JJA, the conditions needed to activate those particles into cloud droplets are weaker than in colder months, which is demonstrated by calculations of the strongest (weakest) aerosol indirect effects in DJF (JJA) based on comparing N d to perturbations in four different aerosol proxy variables (total and sulfate aerosol optical depth, aerosol index, surface mass concentration of sulfate). We used three machine learning models and up to 14 input variables to infer about most influential factors related to N d for DJF and JJA, with the best performance obtained with gradient-boosted regression tree (GBRT) analysis. The model results indicated that cloud fraction was the most important input variable, followed by some combination (depending on season) of CAO index and surface mass concentrations of sulfate and organic carbon. Future work is recommended to further understand aspects uncovered here such as impacts of free tropospheric aerosol entrainment on clouds, degree of boundary layer coupling, wet scavenging, and giant CCN effects on aerosol–N d relationships, updraft velocity, and vertical structure of cloud properties such as adiabaticity that impact the satellite estimation of N d .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of historical precipitation interannual variability in CMIP6 over the United States

Interannual precipitation variability profoundly influences society via its effects on agriculture, water resources, infrastructure, and disaster risks. In this study, we use daily in situ precipitation observations from the global historical climatology network-daily (GHCN-D) to assess the ability of 21 Coupled Model Intercomparison Project Phase 6 (CMIP6) models, including the 50-member fifth-generation Canadian Earth System Model single model initial-condition large ensemble (CanESM5_SMILE), to realistically simulate historical interannual precipitation variability trends within 17 regions of the contiguous United States (CONUS). We assess how accurately the CMIP6 simulations align with observational data across annual, summer, and winter periods, focusing on four key hydrometeorological metrics, including interannual precipitation variability, relative interannual precipitation variability (coefficient of variation), annual mean precipitation, and annual wet day frequency. Our findings reveal that CMIP6 ensemble members generally reproduce the spatial patterns of observed trends in annual mean precipitation. In most regions, models agree well with the signs of observed changes in annual mean precipitation, though discrepancies in trend magnitude are evident. Further, observed trends in winter mean precipitation broadly exhibit a spatial pattern similar to that of the observed annual mean. However, analysis of the CanESM5_SMILE shows that trends in precipitation variability may primarily be the result of model-simulated internal variability, suggesting caution in interpreting multi-model single-realization ensemble results. Challenges in accurately simulating interannual precipitation variability underscore the need for ongoing model refinement and validation to enhance climate projections, especially in regions vulnerable to extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Performance of reverse osmosis membrane with large feed pressure fluctuations from a wave-driven desalination system

Wave-driven desalination systems are proposed water treatment systems that involve reverse osmosis of seawater powered directly by wave motion. Such a configuration would result in drastic feed pressure fluctuations. For a technology conventionally operated with a constant feed condition, the effect of these variable pressures on membrane integrity and performance is unknown. Here, experiments were conducted with spiral wound membranes coupled to a system capable of producing feed pressure fluctuations of more than 400 psi. Feed composition included 5, 20, and 35 g/L NaCl, and a synthetic seawater at normal and 1.5x concentration. The variable feed conditions included sine-like pressure waves swings of 200-500 and 500-900 psi with frequencies of 1.25, 7.5, and 12 waves/min, and a model-generated random waveform. Between each wave experiment we performed membrane integrity tests at 650 psi and 25 g/L NaCl feed, which showed a 7.4% drop in the membrane's water permeability coefficient, an 18.4% flux decline, and more than 99% salt rejection over 1770 h of cumulative experimental time. Analysis of permeate samples showed high salt rejection. In general, variable feed pressure had no significant deleterious effect on membrane integrity or performance.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗
DOE OSTI2021

Biogeochemistry of the Antrim Shale Natural Gas Reservoir

The Antrim Shale, located in the Michigan Basin, United States (U.S.), is a major U.S. shale play having produced over 2.5 Trillion Cubic Feet (Tcf) of unconventional shale natural gas as of 2010. The shallow nature of this formation sets it apart from other, more characterized unconventional shale gas plays. The depth of gas production of the Antrim ranges from approximately 150 to 600 m and it is typically vertically drilled, contrary to deeper, horizontally drilled shales. A thorough understanding of the biogeochemistry and microbiology of this complex system will be advantageous for improving well performance, produced water management, and potential biocidal treatment as microbial community composition can vary substantially even among closely spaced wells. In this study, we analyzed produced water collected from nine different wells in the Antrim Shale by investigating the geochemical and microbial community composition of the produced water to gain greater insight into the overall biogeochemistry of this unique shale system. The majority of the wells from this study had high total dissolved solids (TDS) primarily composed of chloride and sodium, averaging 86,804 mg/L with a maximum 116,223 mg/L; however, three of the wells sampled along the northern margin of the basin exhibited significantly lower TDS ranging from 4932 to 6496 mg/L. Overall, our microbial community analysis revealed relatively low abundance within our samples and high variability of the microbial community among the sampled wells. The majority of bacterial sequences were identified within Proteobacteria, Firmicutes, and Actinobacteria phyla and metagenomic sequencing revealed the low presence of Methanobacteriaceae within each sample. We also investigated potential microbial community drivers and found that TDS, sodium, chloride, iodide, bromide, ammonium, potassium, and strontium were significantly correlated with the observed microbial community. The varying geochemical conditions between wells demonstrate different subsurface environmental niches, potentially driving the heterogeneous microbial communities we observed from well to well. This analysis suggests an important relationship between both well location and geochemistry and the observed microbial community that can persist in the reservoir. Continued studies of the Antrim Shale will improve our understanding of the complex interdependencies of this ecosystem.

03 NATURAL GAS↗
DOE OSTI2024

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗
DOE OSTI2024

Atmospheric Variability Drives Anomalies in the Bering Sea Air–Sea Heat Exchange

Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance Statement A long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change.

Hayden, Emily E. [a College of Earth, Ocean, and A↗
DOE OSTI2024

Agent-Based Model of Combined Community- and Jail-Based Take-Home Naloxone Distribution

Importance Opioid-related overdose accounts for almost 80 000 deaths annually across the US. People who use drugs leaving jails are at particularly high risk for opioid-related overdose and may benefit from take-home naloxone (THN) distribution. Objective To estimate the population impact of THN distribution at jail release to reverse opioid-related overdose among people with opioid use disorders. Design, Setting, and Participants This study developed the agent-based Justice-Community Circulation Model (JCCM) to model a synthetic population of individuals with and without a history of opioid use. Epidemiological data from 2014 to 2020 for Cook County, Illinois, were used to identify parameters pertinent to the synthetic population. Twenty-seven experimental scenarios were examined to capture diverse strategies of THN distribution and use. Sensitivity analysis was performed to identify critical mediating and moderating variables associated with population impact and a proxy metric for cost-effectiveness (ie, the direct costs of THN kits distributed per death averted). Data were analyzed between February 2022 and March 2024. Intervention Modeled interventions included 3 THN distribution channels: community facilities and practitioners; jail, at release; and social network or peers of persons released from jail. Main Outcomes and Measures The primary outcome was the percentage of opioid-related overdose deaths averted with THN in the modeled population relative to a baseline scenario with no intervention. Results Take-home naloxone distribution at jail release had the highest median (IQR) percentage of averted deaths at 11.70% (6.57%-15.75%). The probability of bystander presence at an opioid overdose showed the greatest proportional contribution (27.15%) to the variance in deaths averted in persons released from jail. The estimated costs of distributed THN kits were less than $\$$15 000 per averted death in all 27 scenarios. Conclusions and Relevance This study found that THN distribution at jail release is an economical and feasible approach to substantially reducing opioid-related overdose mortality. Training and preparation of proficient and willing bystanders are central factors in reaching the full potential of this intervention.

Tatara, Eric [Argonne National Laboratory (ANL), A↗
DOE OSTI2023

FIND: A Synthetic weather generator to control drought Frequency, Intensity, and Duration

Water systems worldwide are experiencing climate change-induced shifts in drought properties like frequency, intensity, and duration, affecting water security and reliability. To develop and test effective drought preparedness plans, researchers often use synthetic weather generators to create hydrological scenarios that explore drought variability beyond historical records. Existing weather generators typically allow users to adjust streamflow statistics like percentiles or temporal correlation but do not directly control drought properties of frequency, intensity, and duration. To fill this gap, we propose FIND (Frequency, INtensity, and Duration) synthetic weather generator. FIND incorporates a standardized drought index to directly and in dependently control drought frequency, intensity, and duration in generated streamflow time series while preserving observed hydrological variability. Use cases for FIND include i) water systems analysis applications that seek to train and test drought strategies under historical and plausible future drought conditions, and ii) bottom-up vulnerability studies relating system vulnerability outcomes to specific changes in drought properties of frequency, intensity, and duration. Here, we demonstrate FIND’s versatility through three experiments: replicating historically observed drought properties, generating streamflow scenarios for multiple sites preserving correlation between their drought conditions, and generating a set of scenarios with direct and independent changes in drought properties. FIND source code is openly available for applications beyond the scope of this paper.

42 ENGINEERING↗