SEARCH · Search NASA

Results for “Probability estimates”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Central-place foraging poses variable constraints year-round in a neotropical migrant

Background: “Central-place foragers” are constrained in their habitat selection and foraging range by the frequency with which they need to return to a central place. For example, chick-rearing songbirds that must feed their offspring hourly might be expected to have smaller foraging ranges compared to non-breeding songbirds that return nightly to a roost. Methods: We used GPS units to compare the foraging behaviour of an aerial insectivorous bird, the purple martin (Progne subis), during the breeding season in three regions across North America, as well as the non-breeding season in South America. Specifically, we tested foraging range size and habitat selection. Results: Foraging range did not vary among regions during breeding (14.0 ± 39.2 km 2 ) and was larger during the nonbreeding period (8840 ± 8150 km 2 ). Purple martins strongly preferred aquatic habitats to other available habitats year-round and in the Amazon commuted from night roosts in low productivity sediment-poor water, where risk of predation was probably low, to daytime foraging sites in productive sediment-rich water sites. Conclusions: We provide the first estimates for foraging range size in purple martins and demonstrate foraging preference for aquatic habitats throughout two stages of the annual cycle. Understanding foraging constraints and habitat of aerial insectivores may help plan conservation actions throughout their annual cycle. Future research should quantify foraging behaviour during the post-breeding period and during migration.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Safety Design of Water Infrastructures in A Modern Era

In this chapter, we review the current practice of engineering safety design of large water management infrastructures. Those large infrastructures are designed using probable maximum precipitation (PMP) to withstand extreme conditions, but climate change has been posing challenges on the conventional PMP estimation approach. We discuss multiple limitations of the current approaches, and propose a few possible solutions based on the advancement of atmospheric modeling and statistical analysis of extreme precipitation mechanisms.

Chen, Xiaodong↗

DOE OSTI2022

Abundance and Migration Success of Overshoot Steelhead in the Upper Columbia River

Abstract Summer steelhead Oncorhynchus mykiss may enter freshwater almost a year before spawning and potentially make long migrations (>1,000 km) to interior headwater habitats. However, in response to suboptimal freshwater habitat conditions (e.g., warmer water temperatures), adult summer steelhead may exhibit complex behaviors during upstream migration in the Columbia River basin. Steelhead may migrate upstream of their natal tributary (hereafter, referred to as “overshoot”) and spend days to several months before subsequently migrating downstream (hereafter, referred to as “fallback”) to their natal tributary to spawn. An expansion of an existing Bayesian patch occupancy model, derived from observations of adult steelhead that were PIT-tagged to estimate population-specific abundance upstream of the tagging location, incorporated downstream detection locations to estimate the abundance of overshoot fallbacks. Overshoot steelhead abundance at the tagging location was estimated based on the relationship between the number of known overshoot fallbacks (i.e., the number of steelhead that overshot and successfully migrated downstream to their natal tributary) and their model-estimated abundance. During the study period (2010–2017), the annual mean proportion of overshoot steelhead that successfully migrated downstream of the tagging location (Priest Rapids Dam) was 0.59 (SD = 0.14). The number of dams encountered by overshoot steelhead during their downstream migration was negatively correlated with their downstream migration success probability. Improved downstream passage survival for adult steelhead will increase the abundance of affected populations while reducing potential genetic introgression of upstream populations (i.e., strays). This is the first study to estimate the abundance of overshoot and fallback steelhead, providing the data necessary for scientists to estimate potential conservation benefits of improved downstream survival. For example, surface flow passage routes (e.g., sluiceways and temporary spillway weirs) are very effective in guiding and passing adult steelhead downstream of Columbia River hydroelectric projects and data from this assessment show that changes in dam operations throughout the downstream migration period may maximize conservation benefits.

Murdoch, Andrew R. (ORCID:0000000324827689)↗

DOE OSTI2021

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Quantum Search Approaches to Sampling-Based Motion Planning

In this paper, we present a novel formulation of traditional sampling-based motion planners as database-oracle structures that can be solved via quantum search algorithms. We consider two complementary scenarios: for simpler sparse environments, we formulate the Quantum Full Path Search Algorithm (q-FPS), which creates a superposition of full random path solutions, manipulates probability amplitudes with Quantum Amplitude Amplification (QAA), and quantum measures a single obstacle free full path solution. For dense unstructured environments, we formulate the Quantum Rapidly Exploring Random Tree algorithm, q-RRT, that creates quantum superpositions of possible parent-child connections, manipulates probability amplitudes with QAA, and quantum measures a single reachable state, which is added to a tree. As performance depends on the number of oracle calls and the probability of measuring good quantum states, we quantify how these errors factor into the probabilistic completeness properties of the algorithm. We then numerically estimate the expected number of database solutions to provide an approximation of the optimal number of oracle calls in the algorithm. We compare the q-RRT algorithm with a classical implementation and verify quadratic run-time speedup in the largest connected component of a 2D dense random lattice. We conclude by evaluating a proposed approach to limit the expected number of database solutions and thus limit the optimal number of oracle calls to a given number.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗

DOE OSTI2022

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗

DOE OSTI2023

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Pathways to excitation of atoms with bicircular laser pulses

We study the excitation of the hydrogen atom by bichromatic circularly polarized laser pulses using numerical solutions of the time-dependent Schrödinger equation. The results are in agreement with the selection rules for multiphoton processes in such fields, namely, excited states are populated in which orbital angular momentum and magnetic quantum numbers are either both odd or both even, independent of the relative helicity, peak intensity, and pulse duration of the pulses. For co-rotating pulses the results show that excitation predominantly proceeds to states with magnetic quantum number of the same helicity as the laser pulses. Besides pathways via direct photon absorption from the ground state our results indicate that a transfer of population among the Rydberg states occurs via Λ-type transitions. In the case of counter-rotating pulses the largest excitation probability is found for Rydberg states that differ in magnetic quantum number by Δm=±3. This pattern allows us to estimate how many photons from each of the two bichromatic fields have been absorbed. Finally, we confirm that a population in Rydberg states beyond a maximum orbital angular quantum number is unlikely.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Data-driven Minimum Entropy Control for Stochastic Nonlinear Systems using the Cumulant-Generating Function

Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.

42 ENGINEERING↗

DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗

DOE OSTI2025

A DECADE of dwarfs: first detection of weak lensing around spectroscopically confirmed low-mass galaxies

We present the first detection of weak gravitational lensing around spectroscopically confirmed dwarf galaxies, using the large overlap between DESI DR1 spectroscopic data and DECADE/DES weak lensing catalogs. A clean dwarf galaxy sample with well-defined redshift and stellar mass cuts enables excess surface mass density measurements in two stellar mass bins ($\log \rm{M}_*=[8.2, 9.2]~M_\odot$ and $\log \rm{M}_*=[9.2, 10.2]~M_\odot$), with signal-to-noise ratios of $5.6$ and $12.4$ respectively. This signal-to-noise drops to $4.5$ and $9.2$ respectively for measurements without applying individual inverse probability (IIP) weights, which mitigates fiber incompleteness from DESI's targeting. The measurements are robust against variations in stellar mass estimates, photometric shredding, and lensing calibration systematics. Using a simulation-based modeling framework with stellar mass function priors, we constrain the stellar mass-halo mass relation and find a satellite fraction of $\simeq 0.3$, which is higher than previous photometric studies but $1.5σ$ lower than $Λ$CDM predictions. We find that IIP weights have a significant impact on lensing measurements and can change the inferred $f_{\rm{sat}}$ by a factor of two, highlighting the need for accurate fiber incompleteness corrections for dwarf galaxy samples. Our results open a new observational window into the galaxy-halo connection at low masses, showing that future massively multiplexed spectroscopic observations and weak lensing data will enable stringent tests of galaxy formation models and $Λ$CDM predictions.

To, Chun-Hao [Chicago U., Astron. Astrophys. Ctr.;↗

DOE OSTI2024

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

DOE OSTI2020

Toward Complete Characterization: Prospects for Directly Imaging Transiting Exoplanets

High-contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of “standard planets” must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semimajor axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ∼8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at the V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Machine learning and shallow groundwater chemistry to identify geothermal prospects in the Great Basin, USA

This study discovers various geothermal prospects in the Great Basin, USA based on shallow groundwater chemical (geochemical) data. The geochemical data are expected to include hidden (latent) information that is a proxy for geothermal prospectivity. We processed the sparse geochemical data in the Great Basin at 14,341 locations including 18 attributes. Next, a non-negative matrix factorization with customized k-means clustering is applied to the geochemical data matrix that automatically finds three hidden geothermal signatures representing modestly, moderately, and highly confident geothermal prospects. The algorithm also evaluated the probability of occurrence of these types of resources through the studied region. There is a consistency between regional geothermal prospectivity as estimated by our ML methodology and the traditional play fairway analysis conducted over a portion of the study area. We also identify the dominant data attributes associated with each signature. Finally, our ML analyses allow us to reconstruct attributes from sparse into continuous over the study domain. The predicted continuous attributes can be used for future detailed geothermal explorations in the Great Basin.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Equatorial waves triggering extreme rainfall and floods in southwest Sulawesi, Indonesia

On the basis of detailed analysis of a case study and long-term climatology, it is shown that equatorial waves and their interactions serve as precursors for extreme rain and flood events in the central Maritime Continent region of southwest Sulawesi, Indonesia. Meteorological conditions on January 22, 2019, leading to heavy rainfall and devastating flooding in this area were studied. It is shown that a convectively coupled Kelvin wave (CCKW) and a convectively coupled Rossby wave (CCERW) embedded within the larger-scale envelope of the Madden-Julian Oscillation (MJO) enhanced convective phase, contributed to the onset of a mesoscale convective system which developed over the Java Sea. Low-Level convergence from the CCKW forced mesoscale convective organization and orographic ascent of moist air over the slopes of southwest Sulawesi. Climatological analysis shows that 92 % of December-January-February floods and 76% of extreme rain events in this region were immediately preceded by positive low-level westerly wind anomalies. It is estimated that both CCKWs and CCERWs propagating over Sulawesi double the chance of floods and extreme rain event development, which the probability of such hazardous events occurring during their combined activity is eight times greater than on a random day. While the MJO is a key component shaping tropical atmospheric variability, it is shown that its usefulness as a single factor for extreme weather-driven hazard prediction is limited.

Latos, Beata↗

DOE OSTI2024

Systematic Effects in Galaxy-Galaxy Lensing with DESI

The Dark Energy Spectroscopic Instrument (DESI) survey will measure spectroscopic redshifts for millions of galaxies across roughly $14,000 \, \mathrm{deg}^2$ of the sky. Cross-correlating targets in the DESI survey with complementary imaging surveys allows us to measure and analyze shear distortions caused by gravitational lensing in unprecedented detail. In this work, we analyze a series of mock catalogs with ray-traced gravitational lensing and increasing sophistication to estimate systematic effects on galaxy-galaxy lensing estimators such as the tangential shear $\gamma_{\mathrm{t}}$ and the excess surface density $\Delta\Sigma$ . We employ mock catalogs tailored to the specific imaging surveys overlapping with the DESI survey: the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey, and the Kilo-Degree Survey (KiDS). Among others, we find that fiber incompleteness can have significant effects on galaxy-galaxy lensing estimators but can be corrected effectively by up-weighting DESI targets with fibers by the inverse of the fiber assignment probability. Similarly, we show that intrinsic alignment and lens magnification are expected to be statistically significant given the precision forecasted for the DESI year-1 data set. Our study informs several analysis choices for upcoming cross-correlation studies of DESI with DES, HSC, and KiDS.

79 ASTRONOMY AND ASTROPHYSICS↗