Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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 289 records · Page 16

The distribution of galaxies’ gravitational field stemming from their tidal interaction

We calculate the distribution function of astronomical objects’ (like galaxies and/or smooth halos of different kinds) gravitational fields due to their tidal interaction. For that we apply the statistical method of Chandrasekhar, used originally to calculate the famous Holtzmark distribution. We show that in our approach the distribution function is never Gaussian, its form being dictated by the potential of interaction between objects. This calculation permits us to perform a theoretical analysis of the relation between angular momentum and mass (richness) of the galaxy clusters. To do so, we follow the ideas of Catelan and Theuns and Heavens and Peacock. The main difference is that here we reduce the problem to a discrete many-body case, where all physical properties of the system are determined by the interaction potential V(r{sub ij}). The essence of reduction is that we use the multipole (up to quadrupole here) expansion of Newtonian potential so that all hydrodynamic, “extended” characteristics of an object, such as its density mass, are “integrated out,” leaving its “point-like” characteristics, such as mass and quadrupole moment. In that sense we do not distinguish between galaxies and smooth components such as halos. We compare our theoretical results with observational data.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar velocity dispersion and anisotropy of the Milky Way inner halo

We measure the three components of velocity dispersion, σ{sub R}, σ{sub θ}, σ{sub ϕ}, for stars within 6 < R < 30 kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine our measurements with previously published data so that we can more finely sample the stellar halo. We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and covariances assuming only that velocities are normally distributed. The alignment of the velocity ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans equation, the mass of the Milky Way is M(R⩽12 kpc)=1.3×10{sup 11} M{sub ⊙} with an uncertainty of 40%. We also find a region of discontinuity, 15 ≲ R ≲ 25 kpc, where the estimated velocity dispersions and anisotropies diverge from their anticipated values, confirming the break observed by others. We argue that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, such as the Sgr stream. To significantly improve our understanding of halo kinematics will require combining radial velocities with future Gaia proper motions.

79 ASTRONOMY AND ASTROPHYSICS↗

Blob-hole correlation model for edge turbulence and comparisons with NSTX gas puff imaging data [Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data]

We report that gas puff imaging (GPI) observations made in NSTX have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracking of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Finally, other properties of the experimentally observed extrema are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Angular correlation of cosmic neutrinos with ultrahigh-energy cosmic rays and implications for their sources

Cosmic neutrino events detected by the IceCube Neutrino Observatory with energy 0∼> 3 TeV have poor angular resolutions to reveal their origin. Ultrahigh-energy cosmic rays (UHECRs), with better angular resolutions at 0>6 EeV energies, can be used to check if the same astrophysical sources are responsible for producing both neutrinos and UHECRs. We test this hypothesis, with statistical methods which emphasize invariant quantities, by using data from the Pierre Auger Observatory, Telescope Array and past cosmic-ray experiments. We find that the arrival directions of the cosmic neutrinos are correlated with 0≥ 10 EeV UHECR arrival directions at confidence level ≈ 90%. The strength of the correlation decreases with decreasing UHECR energy and no correlation exists at energy 0∼ 6 EeV . A search in astrophysical databases within 3{sup o} of the arrival directions of UHECRs with energy 0≥ 10 EeV, that are correlated with the IceCube cosmic neutrinos, resulted in 18 sources from the Swift-BAT X-ray catalog with redshift z≤ 0.06. We also found 3 objects in the Kühr catalog of radio sources using the same criteria. The sources are dominantly Seyfert galaxies with Cygnus A being the most prominent member. We calculate the required neutrino and UHECR fluxes to produce the observed correlated events, and estimate the corresponding neutrino luminosity (25 TeV–2.2 PeV) and cosmic-ray luminosity (500 TeV–180 EeV), assuming the sources are the ones we found in the Swift-BAT and Kühr catalogs. We compare these luminosities with the X-ray luminosity of the corresponding sources and discuss possibilities of accelerating protons to 0∼> 10 EeV and produce neutrinos in these sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Energetic particle optimization of quasi-axisymmetric stellarator equilibria

Abstract An important goal of stellarator optimization is to achieve good confinement of energetic particles such as, in the case of a reactor, alphas created by deuterium–tritium fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining this improved QA while minimizing the analytical quantity Γ C , which represents the angle between magnetic flux surfaces and contours of J | | , the second adiabatic invariant. This was performed multiple times, resulting in a group of equilibria with significantly reduced energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator have successfully been reduced by optimizing Γ C . The relationship between energetic particle losses and metrics such as QA error ( E q a ) and Γ C in this set of equilibria were examined via statistical methods and a nearly linear relationship between volume-averaged Γ C and prompt particle losses was found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Gaussian process guide for signal regression in magnetic fusion

Extracting reliable information from diagnostic data in tokamaks is critical for understanding, analyzing, and controlling the behavior of fusion plasmas and validating models describing that behavior. Recent interest within the fusion community has focused on the use of principled statistical methods, such as Gaussian process regression (GPR), to attempt to develop sharper, more reliable, and more rigorous tools for examining the complex observed behavior in these systems. While GPR is an enormously powerful tool, there is also the danger of drawing fragile, or inconsistent conclusions from naive GPR fits that are not driven by principled treatments. Here we review the fundamental concepts underlying GPR in a way that may be useful for broad-ranging applications in fusion science. We also revisit how GPR is developed for profile fitting in tokamaks. We examine various extensions and targeted modifications applicable to experimental observations in the edge of the DIII-D tokamak. Finally, we discuss best practices for applying GPR to fusion data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A consistent positive association between landscape simplification and insecticide use across the Midwestern US from 1997 through 2012

During 2007, counties across the Midwestern US with relatively high levels of landscape simplification (i.e., widespread replacement of seminatural habitats with cultivated crops) had relatively high crop-pest abundances which, in turn, were associated with relatively high insecticide application. These results suggested a positive relationship between landscape simplification and insecticide use, mediated by landscape effects on crop pests or their natural enemies. A follow-up study, in the same region but using different statistical methods, explored the relationship between landscape simplification and insecticide use between 1987 and 2007, and concluded that the relationship varied substantially in sign and strength across years. Here, we explore this relationship from 1997 through 2012, using a single dataset and two different analytical approaches. We demonstrate that, when using ordinary least squares (OLS) regression, the relationship between landscape simplification and insecticide use is, indeed, quite variable over time. However, the residuals from OLS models show strong spatial autocorrelation, indicating spatial structure in the data not accounted for by explanatory variables, and violating a standard assumption of OLS. When modeled using spatial regression techniques, relationships between landscape simplification and insecticide use were consistently positive between 1997 and 2012, and model fits were dramatically improved. We argue that spatial regression methods are more appropriate for these data, and conclude that there remains compelling correlative support for a link between landscape simplification and insecticide use in the Midwestern US. We discuss the limitations of inference from this and related studies, and suggest improved data collection campaigns for better understanding links between landscape structure, crop-pest pressure, and pest-management practices.

54 ENVIRONMENTAL SCIENCES↗

Granger causal inference for climate change attribution

Abstract Climate change detection and attribution (D&A) is concerned with determining the extent to which anthropogenic activities have influenced specific aspects of the global climate system. D&A fits within the broader field of causal inference, the collection of statistical methods that identify cause and effect relationships. There are a wide variety of methods for making attribution statements, each of which require different types of input data and focus on different types of weather and climate events and each of which are conditional to varying extents. Some methods are based on Pearl causality (direct experimental interference) while others leverage Granger (predictive) causality, and the causal framing provides important context for how the resulting attribution conclusion should be interpreted. However, while Granger-causal attribution analyses have become more common, there is no clear statement of their strengths and weaknesses relative to Pearl-causal attribution and no clear consensus on where and when Granger-causal perspectives are appropriate. In this prospective paper, we provide a formal definition for Granger-based approaches to trend and event attribution and a clear comparison with more traditional methods for assessing the human influence on extreme weather and climate events. Broadly speaking, Granger-causal attribution statements can be constructed quickly from observations and do not require computationally-intesive dynamical experiments. These analyses also enable rapid attribution, which is useful in the aftermath of a severe weather event, and provide multiple lines of evidence for anthropogenic climate change when paired with Pearl-causal attribution. Confidence in attribution statements is increased when different methodologies arrive at similar conclusions. Moving forward, we encourage the D&A community to embrace hybrid approaches to climate change attribution that leverage the strengths of both Granger and Pearl causality.

Risser, Mark D. (ORCID:0000000319561783)↗

Micro on a macroscale: relating microbial-scale soil processes to global ecosystem function

ABSTRACT Soil microorganisms play a key role in driving major biogeochemical cycles and in global responses to climate change. However, understanding and predicting the behavior and function of these microorganisms remains a grand challenge for soil ecology due in part to the microscale complexity of soils. It is becoming increasingly clear that understanding the microbial perspective is vital to accurately predicting global processes. Here, we discuss the microbial perspective including the microbial habitat as it relates to measurement and modeling of ecosystem processes. We argue that clearly defining and quantifying the size, distribution and sphere of influence of microhabitats is crucial to managing microbial activity at the ecosystem scale. This can be achieved using controlled and hierarchical sampling designs. Model microbial systems can provide key data needed to integrate microhabitats into ecosystem models, while adapting soil sampling schemes and statistical methods can allow us to collect microbially-focused data. Quantifying soil processes, like biogeochemical cycles, from a microbial perspective will allow us to more accurately predict soil functions and address long-standing unknowns in soil ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗

Getting more from heterogeneous HIV-1 surveillance data in a high immigration country: estimation of incidence and undiagnosed population size using multiple biomarkers

Most HIV infections originate from individuals who are undiagnosed and unaware of their infection. Estimation of this quantity from surveillance data is hard because there is incomplete knowledge about (i) the time between infection and diagnosis (TI) for the general population, and (ii) the time between immigration and diagnosis for foreign-born persons. We developed a new statistical method for estimating the incidence of HIV-1 and the number of undiagnosed people living with HIV (PLHIV), based on dynamic modelling of heterogeneous HIV-1 surveillance data. The methods consist of a Bayesian non-linear mixed effects model using multiple biomarkers to estimate TI of HIV-1-positive individuals, and a novel incidence estimator which distinguishes between endogenous and exogenous infections by modelling explicitly the probability that a foreign-born person was infected either before or after immigration. The incidence estimator allows for direct calculation of the number of undiagnosed persons. The new methodology is illustrated combining heterogeneous surveillance data from Sweden between 2003 and 2015. A leave-one-out cross-validation study showed that the multiple-biomarker model was more accurate than single biomarkers (mean absolute error 1.01 vs ≥1.95). We estimate that 816 [95% credible interval (CI) 775-865] PLHIV were undiagnosed in 2015, representing a proportion of 10.8% (95% CI 10.3-11.4%) of all PLHIV. Here, the proposed methodology will enhance the utility of standard surveillance data streams and will be useful to monitor progress towards and compliance with the 90–90-90 UNAIDS target.

60 APPLIED LIFE SCIENCES↗

Realistic galaxy image simulation via score-based generative models

ABSTRACT We show that a denoising diffusion probabilistic model (DDPM), a class of score-based generative model, can be used to produce realistic mock images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument (DESI) grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real data set. We quantify the similarity by borrowing from the deep generative learning literature, using the ‘Fréchet inception distance’ to test for subjective and morphological similarity. We also introduce the ‘synthetic galaxy distance’ metric to compare the emergent physical properties (such as total magnitude, colour, and half-light radius) of a ground truth parent and synthesized child data set. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as adversarial networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate inpainting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we ‘DESI-fy’ cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

79 ASTRONOMY AND ASTROPHYSICS↗

Forecasting constraints on the mean free path of ionizing photons at z ≥ 5.4 from the Lyman-α forest flux autocorrelation function

ABSTRACT Fluctuations in Lyman-α (Ly α) forest transmission towards high-z quasars are partially sourced from spatial fluctuations in the ultraviolet background, the level of which are set by the mean free path of ionizing photons (λmfp). The autocorrelation function of Ly α forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to λmfp. Recent measurements at z ∼ 6 suggest a rapid evolution of λmfp at z > 5.0 which would leave a signature in the evolution of the autocorrelation function. For this forecast, we model mock Ly α forest data with properties similar to the XQR-30 extended data set at 5.4 ≤ z ≤ 6.0. At each z, we investigate 100 mock data sets and an ideal case where mock data matches model values of the autocorrelation function. For ideal data with λmfp = 9.0 cMpc at z = 6.0, we recover $\lambda _{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of λmfp from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $\sim 40~{{\ \rm per\ cent}}$ reduction in the error bars over all z. The distribution of mock values of the autocorrelation function in this work is highly non-Gaussian for high-z, which should caution work with other statistics of the high-z Ly α forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements.

Astronomy & Astrophysics↗

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗

A steep slope and small scatter for the high-mass end of the L–σ relation at z ~ 0.55

We measure the intrinsic relation between velocity dispersion (σ) and luminosity (L) for massive, luminous red galaxies at redshift z ~ 0.55. Here, we achieve unprecedented precision by using a sample of 600 000 galaxies with spectra from the Baryon Oscillation Spectroscopic Survey of the third Sloan Digital Sky Survey (SDSS-III), covering a range of stellar masses M* ≳ 10 11 M ⊙ . We deconvolve the effects of photometric errors, limited spectroscopic signal-to-noise ratio, and red–blue galaxy confusion using a novel hierarchical Bayesian formalism that is generally applicable to any combination of photometric and spectroscopic observables. For an L–σ relation of the form L ∝ σ β , we find β = 7.8 ± 1.1 for σ corrected to the effective radius, and a very small intrinsic scatter of s = 0.047 ± 0.004 in log10σ at fixed L. No significant redshift evolution is found for these parameters. The evolution of the zero-point within the redshift range considered is consistent with the passive evolution of a galaxy population that formed at redshift z = 2–3, assuming single stellar populations. An analysis of previously reported results seems to indicate that the passively evolved high-mass L–σ relation at z ~ 0.55 is consistent with the one measured at z = 0.1. Finally, our results, in combination with those presented in the LF work of Montero-Dorta et al., provide a detailed description of the high-mass end of the red sequence (RS) at z ~ 0.55. This characterization, in the light of previous literature, suggest that the high-mass RS distribution corresponds to the ‘core’ elliptical population.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A direct measurement of the high-mass end of the velocity dispersion function at z ~ 0.55 from SDSS-III/BOSS

When two galaxies that are distant from one another (and also distant from Earth) happen to lie along a single line of sight in the sky, the resulting phenomenon is known as a “gravitational lens.” The gravity of the more nearby galaxy warps the image of the more distant galaxy into multiple images or complete rings (know as “Einstein rings” since the quantitative description of the gravitational lensing effect relies on Einstein’s theory of gravity.) Strong gravitational lens systems have multiple scientific applications. If the more distant galaxy happens to contain a time-varying quasar (bright emission powered by a supermassive black hole at the galaxy’s center) or supernova explosion, the time delay between multiple images can be used as a probe of the expansion rate of the universe (and other cosmological parameters.) Forecasting the incidence of gravitational lenses in future large-scale sky surveys relies on quantifying the population of potential lens galaxies in the universe in terms of their abundance and their lensing efficiency. The lensing efficiency is most directly correlated with the galaxy’s “velocity dispersion:” the characteristic speed with which stars in the galaxy are orbiting under the influence of the galaxy’s overall gravitational field. This paper uses previous results quantifying the combined demographics of galaxies in brightness and velocity dispersion to compute the demographics of massive “elliptical” galaxies in velocity dispersion alone, thereby providing the essential ingredient for forecasting the expected incidence of strong gravitational lensing by these types of galaxies in future sky surveys such as DESI and LSST. These results are also applicable to the association of massive galaxies with their associated dark-matter “halos,” which is an essential ingredient for the most accurate and informative extraction of cosmological parameters from the data sets produced by large-scale surveys of the universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Debiasing cosmic gravitational wave sirens

ABSTRACT Accurate estimation of the Hubble constant, and other cosmological parameters, from distances measured by cosmic gravitational wave sirens requires sufficient allowance for the dark energy evolution. We demonstrate how model-independent statistical methods, specifically Gaussian process regression, can remove bias in the reconstruction of H(z), and can be combined to model independently with supernova distances. This allows stringent tests of both H0 and Λ cold dark matter, and can detect unrecognized systematics. We also quantify the redshift systematic control necessary for the use of dark sirens, showing that it must approach spectroscopic precision to avoid significant bias.

Keeley, Ryan E.↗

Intercomparison of Commercially Available Active Radon Measurement Devices in a Discovered Radon Chamber

An unventilated 70 m 3 underground space, with 2.1 m high ceilings and an inner set of concrete walls to damp vibrations, was discovered to have relatively stable radon levels of ~1,080 Bq m -3 over a 1-y period. The addition of radium sources was not required to achieve high and somewhat variable radon levels, as these arise naturally in the space. Several active radon monitoring devices designed for homeowners, radon screeners, radon mitigation professionals, and researchers were tested in the space. These devices were variable in both cost and intended user sophistication. Data were collected at specific time intervals that varied for each device. A cellphone was deployed with an available camera application to periodically capture photographs of the displays of devices not enabled for automatic temporal recording. Attempts were made to alter the radon concentration using fans, radium-laden objects, and radon impermeable sheets. Statistical methods were employed to compare the device's performance under the minor radon transients produced by the addition of 226 Ra sources during the test period. The evaluation of the overall accuracy, noise, and response to transients revealed highly affordable devices that may be useful for future research applications.While unsuitable for precise calibrations at stable radon levels, the space was determined to be adequate as a radon chamber for general research purposes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗