Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariates”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine learning-assisted identification of potential sources of bias in measurements of prompt-fission neutron spectra

Unrecognized sources of uncertainty (USU) can bias the reported mean and/or covariance of experimental nuclear data. These biases, in turn, can propagate through evaluated nuclear data to application simulations or may poorly inform nuclear theory that is fitted to the experimental data. Such unknown sources of bias must be tied to the inherent physical constituents of the measurements such as the characteristics of a detector response or a background reduction technique. Here, in this article, a sparse Bayesian learning model is used to support experts in their efforts to identify and characterize USU in experimental prompt fission neutron spectra (PFNS) for spontaneous fissioning of 252 Cf by linking observed biases to features of the measurement system. Three different bias components were found. The first acts as a verification case for the algorithm as it identifies a bias coming from a well-known source related to the use of 6 Li in the neutron detection system. The second two cases demonstrate how this method can benefit the evaluation of experimental nuclear data by identifying, quantifying, and relating unknown biases to potential causes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Identification of over ten thousand candidate structured RNAs in viruses and phages

Structured RNAs play crucial roles in viruses, exerting influence over both viral and host gene expression. However, the extensive diversity of structured RNAs and their ability to act in cis or trans positions pose challenges for predicting and assigning their functions. While comparative genomics approaches have successfully predicted candidate structured RNAs in microbes on a large scale, similar efforts for viruses have been lacking. In this study, we screened over 5 million DNA and RNA viral sequences, resulting in the prediction of 10,006 novel candidate structured RNAs. These predictions are widely distributed across taxonomy and ecosystem. We found transcriptional evidence for 206 of these candidate structured RNAs in the human fecal microbiome. These candidate RNAs exhibited evidence of nucleotide covariation, indicative of selective pressure maintaining the predicted secondary structures. Our analysis revealed a diverse repertoire of candidate structured RNAs, encompassing a substantial number of putative tRNAs or tRNA-like structures, Rho-independent transcription terminators, and potentially cis-regulatory structures consistently positioned upstream of genes. In summary, our findings shed light on the extensive diversity of structured RNAs in viruses, offering a valuable resource for further investigations into their functional roles and implications in viral gene expression and pave the way for a deeper understanding of the intricate interplay between viruses and their hosts at the molecular level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Intensive and extensive margins of the peak load: Measuring adaptation with mixed frequency panel data

In this work we investigate the response of daily electricity peak load to daily maximum temperatures across states in Europe and India. We propose a method that decomposes short- from medium/long-run effects, retains the high frequency nature of the load-weather covariation and treats economic growth as a modulating factor. By simultaneously exploiting variation in unexpected daily weather anomalies and decade-long climatic changes in each location we decompose transitory - intensive margin - adjustments from permanent - extensive margin - adjustments. We find that the shocks over the long-run differ substantially from the short-run dynamics. Furthermore, we find evidence that per capita income modulates the adjustments over the short- and long-run. We project that in response to climate change around 2050 the peak load may increase by up to 20%-30% in Southern Europe and in several states in India, depending on the degree of warming and the evolution of socio-economic conditions. Even with a limited scope to two world regions, we identify that the structure of the economy and differences in future income growth matter in shaping the adaptation to climate change. Our decomposition allows to identify how future weather anomalies can further amplify the relative increase associated to the shift in the climate norm. Assuming that the interannual variability of maximum temperatures follows the distribution observed in the past, we find a doubling of the impacts of climate change during the summer in Europe. Uncertainty around the distribution of future weather anomalies may lead to further unexpected peak load amplifications. Our results have important policy implications for power systems’ generation capacity, transmission and storage, as we show that the challenges to accommodate the peak load in days with extreme temperatures may substantially increase already around mid-century.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Environmental exposure to industrial air pollution is associated with decreased male fertility

Objective: To understand how chronic exposure to industrial air pollution is associated with male fertility through semen parameters. Design: Retrospective cohort study. Subjects: Men in the Subfertility, Health and Assisted Reproduction cohort who underwent a semen analysis 2005-2017 with ≥1 measured semen parameter (N=21,563). Intervention(s): Residential histories for each man were constructed using locations from administrative records linked through the Utah Population Database. Industrial facilities with air emissions of nine endocrine disrupting compound chemical classes were identified from the Environmental Protection Agency Risk-Screening Environmental Indicators microdata. Chemical levels were linked with residential histories for the 5 years prior to each semen analysis. Main Outcome Measures: Semen analyses were classified as azoospermic or oligozoospermic (< 15 M/mL) using World Health Organization cutoffs for concentration. Bulk semen parameters such as concentration, total count, ejaculate volume, total motility, total motile count, and total progressive motile count were also measured. Multivariable regression models with robust standard errors were used to associate exposure quartiles for each of the nine chemical classes with each semen parameter, adjusting for age, race, and ethnicity, as well as neighborhood socioeconomic disadvantage. Results: After adjustment for demographic covariates, several chemical classes were associated with azoospermia and decreased total motility and volume. For exposure in the 4th relative to 1st quartile, significant associations were observed for acrylonitrile (β total motility = -0.87 pp), aromatic hydrocarbons (odds ratio [OR]azoospermia = 1.53; β volume = -0.14 mL), dioxins (OR azoospermia = 1.31; β volume = -0.09 mL; β total motility = -2.65 pp), heavy metals (β total motility = -2.78pp), organic solvents (OR azoospermia = 1.75; β volume = -0.10 mL), organochlorines (OR azoospermia = 2.09; β volume = -0.12 mL), phthalates (OR azoospermia = 1.44; β volume = -0.09 mL; β total motility = -1.21 pp), and silver particles (OR azoospermia = 1.64; β volume = -0.11 mL). All semen parameters significantly decreased with increasing socioeconomic disadvantage. Men who lived in the most disadvantaged areas had concentration, volume, and total motility of 6.70 M/mL, 0.13 mL, and 1.79 pp lower, respectively. Count, motile count, and total progressive motile count all decreased by 30–34 M.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2023

Profiling overstory survival trends following varying thinning and burning disturbance regimes in a mixed pine-hardwood forest in the US South

Prescribed burning as a silvicultural tool may be effective for achieving various management objectives, yet some practitioners have concerns about inadvertently increasing overstory tree mortality. Using extensive data from a long-term, ongoing study in northcentral Alabama, USA, that systematically applies a variety of thinning and prescribed burning disturbances on mixed pine-hardwood stands, we assess survival trends for different groups of trees. The primary research question is whether more frequent prescribed burns adversely affects overstory survival versus infrequent or no burns; secondary questions explore whether overstory survival trends differ based on thinning level, species group, or size class. Interest is in broad groups of trees, not individual tree survival or mortality; consequently, survival analysis methods are used, including nonparametric Kaplan-Meier (KM) techniques for examining single grouping factors, and parametric accelerated failure time models for analyzing simultaneous effects of multiple covariates. Leveraging relatively recent methodological advances, the statistical techniques are adjusted for interval censored data. KM survival curves showed that more frequent prescribed fires did not result in differing survival trends. Also, trees in unthinned stands had the lowest 14-year survival compared to thinned stands, oaks had the highest survival compared to pines and the others species group, and the smallest size class of overstory trees had the lowest survival. In conclusion, these results support managing transitional mixed pine-hardwoods using thinning and multiple prescribed fires to restore specific species composition and structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sulfur isotope fractionation between aqueous and carbonate-associated sulfate in abiotic calcite and aragonite

Sulfate (SO 4 2– ) incorporated into calcium carbonate minerals enables measurements of sulfur (S) isotope ratios in carbonate rocks. This Carbonate Associated Sulfate (CAS) in marine carbonate minerals is thought to faithfully represent the S isotope composition of the seawater sulfate incorporated into the mineral, with little or no S isotope fractionation in the process. However, comparison between different calcifying species reveals both positive and negative S isotope fractionation between CAS and seawater sulfate, and a large range of S isotope ratios can be found within a single rock sample, depending on the component measured. To better understand the isotopic effects associated with sulfate incorporation into carbonate minerals, we precipitated inorganic calcite and aragonite over a range covering more than two orders of magnitude of sulfate concentration and precipitation rate. Coupled measurements of CAS concentration, S isotope composition and X-ray absorption near-edge spectra (XANES) permit characterization and explanation of the observed dependence of S isotope fractionation between CAS and aqueous sulfate (CAS-SO 4 2– isotope fractionation) on sulfate concentration and precipitation rate. In aragonite, the CAS-SO 4 2– isotope fractionation is 1.0 ± 0.3‰ and independent of the sulfate (and CAS) concentration. In contrast, we note the CAS-SO 4 2– isotope fractionation in calcite covaries strongly with the sulfate concentration and weakly with the precipitation rate, between values of 1.3 ± 0.1 and 3.1 ± 0.6‰. We suggest that the correlation between aqueous sulfate concentration and CAS-SO 4 2– isotope fractionation in calcite reflects a dependence of the equilibrium S isotope fractionation on the concentration of CAS, through the effect of the sulfate impurity on the carbonate mineral’s energetic state.

58 GEOSCIENCES↗
DOE OSTI2021

Origin of the isotopic composition of natural perchlorate: Experimental results for the impact of reaction pathway and initial ClO x reactant

Natural perchlorate (ClO 4 - ) exists in many places on Earth, in lunar regolith, meteorites, and on the surface of Mars. Terrestrial natural ClO 4 - has widely variable Cl and O stable isotopic compositions (δ 37 Cl, δ 18 O, Δ 17 O). The δ 18 O and Δ 17 O values of ClO 4 - from the most hyper-arid locations co-vary. ClO 4 - from less arid areas has relatively little 17 O excess and poor Δ 17 O-δ 18 O correlation. ClO 4 - from the Atacama Desert has unusually low δ 37 Cl (<-10‰) and exhibits a positive correlation between δ 37 Cl and δ 18 O, while the δ 37 Cl of ClO 4 - from all other locations varies between -5 and +7‰ with no δ 37 Cl-δ 18 O covariation. To evaluate the impact of different precursors (ClO x ) and reaction pathways on the isotopic composition of ClO 4 - , we measured the isotopic composition of ClO 4 - produced in the laboratory by UV or O 3 mediated aqueous oxidation of Cl-, OCl-, ClO2-, and ClO2° as well as O 3 mediated oxidation of dry NaCl. ClO x oxidation in aqueous or dry systems enriched in O 3 produced ClO 4 - with Δ 17 O values that generally increased with the number of O atoms required and included evidence that the site-specific 17 O anomaly in O 3 was preferentially transferred to ClO 4 - . Based on the inferred number of O atoms sourced from O 3 , and known Cl and O reaction pathways, it appears that ClO 2 ° and ClO 3 * were required intermediates in the production of ClO 4 - in the O 3 experiments. ClO x aqueous oxidation by UV irradiation produced ClO 4 - with a large range of δ 18 O values and little or no 17 O anomaly. ClO 3 - was produced to a much greater extent than ClO 4 - in all experiments except dry oxidation of NaCl by O 3 . The isotopic composition of ClO 3 - was distinct from that of ClO 4 - produced from the same initial reactants. Combined results of O 3 and UV mediated reactions largely bracketed the range of natural ClO 4 - δ 18 O and Δ 17 O values as well as δ 37 Cl values of non-Atacama natural samples, but no conditions produced the low δ 37 Cl values of Atacama ClO 4 - . Finally, our results indicate that variation in production mechanisms, possibly combined with isotopically variable precursors, could be responsible for much of the observed isotopic variation in natural ClO 4 - and ClO 3 - .

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

Multi-site evaluation of stratified and balanced sampling of soil organic carbon stocks in agricultural fields

Estimating soil organic carbon (SOC) stocks in agricultural fields is essential for environmental and agronomic research, management, and policy. Stratified sampling is a classic strategy for estimating mean soil properties, and has recently been codified in SOC monitoring protocols. However, for the specific task of estimating the SOC stock of an agricultural field, concrete guidance is needed for which covariates to stratify on and how much stratification can improve estimation efficiency. It is also unknown how stratified sampling of SOC stocks compares to modern alternatives, notably doubly balanced sampling. To address these gaps, we collected high-density (average of 7 samples ha -1 ) and deep (average of 75 cm) measurements of SOC stocks at eight commercial fields under maize-soybean production in two US Midwestern states. We combined these measurements with a Bayesian geostatistical model to evaluate stratified and balanced sampling strategies that use a set of readily-available geographic, topographic, spectroscopic, and soil survey data. We examined the number of samples needed to achieve a given level of SOC stock estimation accuracy. While stratified sampling using these variables enables an average sample size reduction of 17% (95% CI, 11% to 23%) compared to simple random sampling, doubly balanced sampling is consistently more efficient, reducing sample sizes by 32% (95% CI, 25% to 37%). The data most important to these efficiency gains are a remotely-sensed SOC index, SSURGO estimates of SOC stocks, and the topographic wetness index. We conclude that in order to meet the urgent challenge of climate change, SOC stocks in agricultural fields could be more efficiently estimated by taking advantage of this readily-available data, especially with doubly balanced sampling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Regional-scale soil carbon predictions can be enhanced by transferring global-scale soil–environment relationships

Accurate modelling and mapping soil organic carbon are crucial for supporting soil health restoration and climate change mitigation at both regional and global scales. However, regional soil predictions often suffer from data scarcity and high prediction uncertainty. Utilizing a pre-trained global-to-regional soil carbon predictive model can be a potential solution to address this challenge. Despite its promise, how to construct and apply the global-scale model to enhance regional-scale soil carbon mapping remains largely unexplored. Here, we propose the Global Soil Carbon Pre-trained Model (GSoilCPM), a deep-learning-based domain adaptative model, to enhance regional-scale soil carbon predictions. Based on large amount of environmental covariate data and 106,167 soil samples across the globe, we verify our hypothesis of the effectiveness of this 'global-to-regional' modelling strategy. The pre-trained model can be then transferred and fine-tuned to bridge the regional- and global-scale soil–environment relationships. We applied and validated this modelling strategy in four regional-scale study areas, three in the Northern Hemisphere and one in the Southern Hemisphere, each with distinct environmental background. Compared to traditional modelling approaches as a baseline, four case studies all demonstrated significant improvement in prediction accuracy across diverse environments and varying data availabilities. The average percentage improvement across all regions is 10.93% (absolute values decreased by 1.20 g kg−1 averagely) in MAE and 29.04% (absolute values increased by 0.10 averagely) in CCC. The applicability and future horizons of using GSoilCPM were further discussed. We further reveal that regions with fewer soil samples or lower baseline accuracy benefit more from the pre-trained global model. Our findings highlight the advantages of leveraging the generalized knowledge from global models to enhance specifically localized soil modelling, positioning a potential paradigm shift in digital soil mapping, and far-reaching implications for soil monitoring and land management.

Deep learning↗
DOE OSTI2021

A time-varying vulnerability index for COVID-19 in New Mexico, USA using generalized propensity scores

The coronavirus disease (COVID-19) pandemic has highlighted systemic inequities in the United States and resulted in a larger burden of negative social outcomes for marginalized communities. New Mexico, a state in the southwestern US, has a unique population with a large racial minority population and a high rate of poverty that may make communities more vulnerable to negative social outcomes from COVID-19. To identify which communities may be at the highest relative risk, we created a county-level vulnerability index. After the first COVID-19 case was reported in New Mexico on March 11, 2020, we fit a generalized propensity score model that incorporates sociodemographic factors to predict county-level viral exposure and thus, the generic risk to negative social outcomes such as unemployment or mental health impacts. We used four static sociodemographic covariates important for the state of New Mexico—population, poverty, household size, and minority population—and weekly cumulative case counts to iteratively run our model each week and normalize the exposure score to create a time-varying vulnerability index. We found the relative vulnerability between counties varied in the first eight weeks from the initial COVID-19 case before stabilizing. This framework for creating a location-specific vulnerability index in response to an ongoing disaster may be used as a quick, deployable metric to inform health policy decisions such as allocating state resources to the county level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Robust PCA-Deep Belief Network Surrogate Model for Distribution System Topology Identification with DERs

With the expansion of distribution networks and increased penetration of distributed energy resources (DERs), it is becoming increasingly important to obtain accurate distribution network topology in real-time. In this paper, a robust principal component analysis coupled deep belief network (PCA-DBN) surrogate model is proposed for distribution system topology identification. It integrates the benefits of robust feature extraction from PCA to deal with data quality issues and filter out noise, and the strength of DBN in capturing the nonlinear relationship between voltage amplitudes and the binary states of switchable connections. This also significantly reduces the DBN training complexity without loss of accuracy. It is shown that the widely used standard deviation of voltage drop and the voltage covariance matrix features yield less accuracy as compared to that of the voltage amplitudes in presence of high penetration of DERs and ZIP loads. Comparison results with other alternatives, such as the random forest (RF), multi-output regression (MOR) and the traditional DBN methods demonstrate that the proposed method can achieve a much higher topology identification accuracy while maintaining robustness to missing data and measurement noise under various penetration levels of DERs.

deep belief network↗
DOE OSTI2022

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Understanding Competing Risks

Highlights: • Clinicians and researchers need to be cognizant of competing risks (ie, the occurrence of an event that precludes future development of the outcome of interest). • Censoring is appropriate when a patient being censored at a particular time (eg, owing to loss to follow-up) is deemed to have approximately the same future risk of developing the outcome of interest as patients still being followed, and when the reason for censoring is not related to events of interest. • Censoring is not appropriate in many scenarios; for example, a dead patient can never develop a future cancer recurrence or late effect. • Kaplan-Meier (KM) analysis is commonly used for a time-to-event analysis when the primary outcome (eg, overall survival) has no meaningful competing risk and censoring is appropriate. Cox proportional hazards modeling can be used for multivariable analysis to examine the association between covariates and the primary outcome. • In the presence of competing events, incorrect use of the KM analysis will likely lead to an overestimation of the risk of the outcome of interest and biased estimates from both univariable and multivariable analyses. Therefore, KM usually should not be used to determine local recurrence risk because humanity is always at risk for dying, which is a competing risk. • A competing risk analysis is therefore most appropriate when there are competing events to the primary outcome of interest. The most commonly used multivariable model in these situations is the Fine-Gray model. • Many fundamental questions in radiation oncology—local and regional recurrence, treatment-related toxicity—are best addressed using competing risk methodologies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Adjusting solar-induced fluorescence to nadir-viewing provides a better proxy for GPP

Solar-induced fluorescence (SIF) provides key information for inferring terrestrial gross primary productivity (GPP). However, there is little research on analyzing the performance of nadir-adjusted SIF (SIF nadir ) on GPP, compared to total SIF emitted by all leaves (SIF total ) for reducing the viewing angle effects and estimating GPP. Besides, there have been controversial opinions on whether SIF total is better correlated to GPP than raw SIF observations (SIF obs ). Here we systematically analyzed the relationship between raw/nadir/total SIF (i.e., SIF obs , SIF nadir and SIF total ) and GPP and investigated the underlying mechanism, using multi-angular field measurements and eddy covariance data in wheat and corn crops at sub-daily scale. We further compared the performance of SIF obs , SIF nadir and SIF total in estimating GPP using the TROPOspheric Monitoring Instrument (TROPOMI) SIF and concurrent AmeriFlux measurements at daily scale. Furthermore, the results indicate that diurnal SIF nadir has stronger correlations to GPP than SIF obs for field measurements, with the increase of mean coefficient of determination (R 2 ) by 0.05–0.07 for far-red band and 0.11–0.20 for red band. SIF nadir shows comparable performance with SIF total for both far-red and red bands. Although the viewing angle effects dominate the difference in estimating GPP between SIF nadir and SIF obs , the correlation between light use efficiency (LUE) and fesc further determines the different performance of SIF nadir and SIF total in estimating GPP. TROPOMI-based analysis further confirms that SIF nadir overall has higher correlations to AmeriFlux GPP than SIF obs for different plant functional types and shows similar performance with SIF total . Compared to SIF total , the estimation of SIF nadir independent of viewing angle effects does not require any canopy structure parameters, and thus offers promising potential for reliably estimating regional and global GPP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Updated neutrino mass constraints from galaxy clustering and CMB lensing-galaxy cross-correlation measurements

In this work, we revisit cosmological constraints on the sum of the neutrino masses Σm v from a combination of full-shape BOSS galaxy clustering [P(k)] data and measurements of the cross-correlation between Planck Cosmic Microwave Background (CMB) lensing convergence and BOSS galaxy overdensity maps [$C_{l}^{kg}$], using a simple but theoretically motivated model for the scale-dependent galaxy bias in auto- and cross-correlation measurements. We improve upon earlier related work in several respects, particularly through a more accurate treatment of the correlation and covariance between P(k) and $C_{l}^{kg}$ measurements. When combining these measurements with Planck CMB data, we find a 95% confidence level upper limit of Σm v sum of neutrino masses < 0.14 eV, while slightly weaker limits are obtained when including small-scale ACTPol CMB data, in agreement with our expectations. We confirm earlier findings that (once combined with CMB data) the full-shape information content is comparable to the geometrical information content in the reconstructed BAO peaks given the precision of current galaxy clustering data, discuss the physical significance of our inferred bias and shot noise parameters, and perform a number of robustness tests on our underlying model. While the inclusion of $C_{l}^{kg}$ measurements does not currently appear to lead to substantial improvements in the resulting Σm v of neutrino mass constraints, we expect the converse to be true for near-future galaxy clustering measurements, whose shape information content will eventually supersede the geometrical one.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Global LAnd Surface Satellite (GLASS) evapotranspiration product Version 5.0: Algorithm development and preliminary validation

An accurate estimation of spatially and temporally continuous global terrestrial evapotranspiration (ET) is essential in the assessment of surface energy, water and carbon cycles. The Global LAnd Surface Satellite (GLASS) ET product Version 4.0 (v4.0) based on the Bayesian model averaging (BMA) method was generated to estimate global terrestrial ET. However, certain uncertainty for the GLASS ET product v4.0 limits its application. In this study, we introduced the deep neural networks (DNN) merging framework to improve terrestrial ET estimation for GLASS ET product Version 5.0 (v5.0) generation by integrating five satellite-derived ET products [Moderate Resolution Imaging Spectroradiometer (MODIS) ET product (MOD16), Shuttleworth–Wallace dual-source ET product (SW), Priestley–Taylor-based ET product (PT-JPL), modified satellite-based Priestley–Taylor ET product (MS-PT) and simple hybrid ET product (SIM)]. We compared the performance of DNN method against other merging methods, including GLASS ET algorithm v4.0 (BMA), the gradient boosting regression tree (GBRT) method and the random forest (RF) method, based on 195 global eddy covariance (EC) flux towers covering observations from 2000 through 2015. Validations indicated that the DNN had the highest accuracy among four merging methods across different land cover types, yielding the highest average determination coefficients (R 2 , 0.62), root-mean-squared-error (RMSE, 24.1 W/m 2 ) and Kling–Gupta efficiency (KGE, 0.77) with a of 99% confidence interval. Compared with GLASS ET algorithm v4.0, the DNN improved on the R 2 by approximately 7% (p < 0.01) and the KGE by 10%. Based on the DNN, we then generated 8-day GLASS ET product v5.0 globally with a 1 km spatial resolution from 2001 to 2015 driven by GLASS vegetation and surface net radiation (R n ) datasets and Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA2) datasets. Finally, this global terrestrial ET product provides a valuable dataset for monitoring regional and global water resources and environmental changes.

54 ENVIRONMENTAL SCIENCES↗