Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sufficient statistics”

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 361 records · Page 20

Interpreting internal consistency of DES measurements

ABSTRACT Bayesian evidence ratios are widely used to quantify the statistical consistency between different experiments. However, since the evidence ratio is prior dependent, the precise translation between its value and the degree of concordance/discordance requires additional information. The most commonly adopted metric, the Jeffreys scale, can falsely suggest agreement between data sets when priors are chosen to be sufficiently wide. This work examines evidence ratios in a DES-Y1 simulated analysis, focusing on the internal consistency between weak lensing and galaxy clustering. We study two scenarios using simulated data in controlled experiments. First, we calibrate the expected evidence ratio distribution given noise realizations around the best-fitting DES-Y1 ΛCDM cosmology. Secondly, we show the behaviour of evidence ratios for noiseless fiducial data vectors simulated using a modified gravity model, which generates internal tension in the LCDM analysis. We find that the evidence ratio of noise realizations generated at all confidence levels was biased towards agreement and show, with a modified gravity model, that the choice of prior could conceal the discrepancies between weak lensing and galaxy clustering induced by prior effects in unlike cosmological models, concluding that the evidence ratio in a DES-Y1 study is, indeed, biased towards agreement. Boundary effects can also influence conclusions about the inconsistency induced by modified gravity, even in a noiseless data vector simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Monitoring the sea surface with a short pulse radar

A solution is presented for the scattering of short pulses from a stochastic, corrugated surface relative to the sea for the case of a narrow-beam transmitting antenna pointing near nadir. The spectrum of the received power and its time history are calculated and this solution is used to show that a measure of the variance of the surface ordinant can be obtained from the backscattered power. Included explicitly in the analysis is the finite nature of the source and the role of the small-scale wave structure (capillary wave range). It is shown that when sufficiently short pulses are transmitted, one can obtain a measure of the variance of the large scale surface ordinant from either the temporal spacing of the peaks in the returned power or from the envelope of the spectrum of the received power. Assuming an appropriate model for the statistics and spectrum of the surface ordinate, the variance can be used to compute the wind speed and the significant wave height of the surface.

Le Vine, D. M.↗

Energy dependence of light hypernuclei production in heavy-ion collisions from a coalescence and statistical-thermal model perspective

A comparison of light hypernuclei production, from Ultra-relativistic Quantum Molecular Dynamics (UrQMD)+coalescence and the thermal model, in heavy ion collisions over a wide range of beam energies and system sizes is presented. We find that both approaches provide generally similar results, with differences in specific details. Especially the ratios of hypertriton to Λ are affected by both the source radius Δr of the coalescence procedure as well as canonical effects. On the other hand, the double ratio S 3 is almost independent of canonical effects, which is in contrast to coalescence. Thus, both the beam energy dependence and centrality dependence of S 3 can be used to constrain the hypertriton source radius. To do so the currently available data are not yet sufficient. Elliptic flow is shown to be unaffected by the source size of the nuclei and an almost perfect mass scaling of the elliptic flow is observed. Furthermore, our predictions further suggest that the existence of the H-dibaryon(ΛΛ) seems ruled out by ALICE data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The LSST DESC data challenge 1: generation and analysis of synthetic images for next-generation surveys

Data Challenge 1 (DC1) is the first synthetic data set produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST data set. DC1 is comprised of r -band observations of 40 deg 2 to 10 yr LSST depth. In this paper, we present each stage of the simulation and analysis process: (a) generation, by synthesizing sources from cosmological N -body simulations in individual sensor-visit images with different observing conditions; (b) reduction using a development version of the LSST Science Pipelines; and (c) matching to the input cosmological catalogue for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star–galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: (i) survey properties have an impact of 50 per cent of the statistical uncertainty for the scales and models used in DC1; (ii) a selection to eliminate artefacts in the catalogues is necessary to avoid biases in the measured clustering; and (iii) the presence of bright objects has a significant impact (2σ–6σ) in the estimated power spectra at small scales (ℓ > 1200), highlighting the impact of blending in studies at small angular scales in LSST.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Untargeted GC-MS Metabolic Profiling of Anaerobic Gut Fungi Reveals Putative Terpenoids and Strain-Specific Metabolites

Background/Objectives: Anaerobic gut fungi (Neocallimastigomycota) are biotechnologically relevant, lignocellulose-degrading microbes with under-explored biosynthetic potential for secondary metabolites. Untargeted metabolomic profiling with gas chromatography–mass spectrometry (GC-MS) was applied to two gut fungal strains, Anaeromyces robustus and Caecomyces churrovis, to establish a foundational metabolomic dataset to identify metabolites and provide insights into gut fungal metabolic capabilities. Methods: Gut fungi were cultured anaerobically in rumen-fluid-based media with a soluble substrate (cellobiose), and metabolites were extracted using the Metabolite, Protein, and Lipid Extraction (MPLEx) method, enabling metabolomic and proteomic analysis from the same cell samples. Samples were derivatized and analyzed via GC-MS, followed by compound identification by spectral matching to reference databases, molecular networking, and statistical analyses. Results: Distinct metabolites were identified between A. robustus and C. churrovis, including 2,3-dihydroxyisovaleric acid produced by A. robustus and maltotriitol, maltotriose, and melibiose produced by C. churrovis. C. churrovis may polymerize maltotriose to form an extracellular polysaccharide, like pullulan. GC-MS profiling potentially captured sufficiently volatile products of proteomically detected, putative non-ribosomal peptide synthetases and polyketide synthases of A. robustus and C. churrovis. The triterpene squalene and triterpenoid tetrahymanol were putatively identified in A. robustus and C. churrovis. Their conserved, predicted biosynthetic genes—squalene synthase and squalene tetrahymanol cyclase—were identified in A. robustus, C. churrovis, and other anaerobic gut fungal genera. Conclusions: This study provides a foundational, untargeted metabolomic dataset to unmask gut fungal metabolic pathways and biosynthetic potential and to prioritize future efforts for compound isolation and identification.

Biochemistry & Molecular Biology↗

Timing the Geminga Pulsar with Gamma-Ray Observations

We present the COS-B/EGRET 1997 ephemeris for the rotation of the Geminga pulsar. This ephemeris is derived from high-energy gamma-ray observations that span 24 yr. The recently obtained accurate position and proper motion are assumed. A cubic ephemeris predicts the rotational phase of Geminga with errors smaller than 50 milliperiods for all existing high-energy gamma-ray observations that span a 24.2 yr timing baseline. The braking index obtained is 17 +/- 1. Further observation is required to ascertain whether this high value truly reflects the rotational energy loss mechanism, or whether it is a manifestation of timing noise. The ephemeris parameters are sufficiently constrained so that timing noise will be the limitation on forward extrapolation. If Geminga continues to rotate without a glitch, as it has for at least 23 yr, we expect this ephemeris to continue to describe the phase, with an error of less than 100 milliperiods, until 2008. Statistically significant timing residuals are detected in the EGRET data that depart from the cubic ephemeris at a level of 30 milliperiods. Although this could simply be an additional manifestation of timing noise, the EGRET timing residuals appear to have a sinusoidal modulation that is consistent with a planet of mass 1.7/sin i solar mass, orbiting Geminga at a radius of 3.3 AU.

Mattox, J. R.↗

The Orbiting Wide-Angle Light Collectors (OWL) Mission for Charged-Particle Astronomy

The space-based OWL mission is designed to perform high-statistics measurements of ultra-high-energy cosmic rays (UHECR) using the Earth's atmosphere as a vast particle calorimeter, furthering the field of charged-particle astronomy. OWL has been developed in formal NASA instrument and mission studies and is comprised of two large telescopes separated by approx.600 km in 1000 km, near-equatorial orbits to stereoscopically image the near-UV air fluorescence emitted by UHECR-induced particle cascades. The High Resolution Fly's Eye (HiRes) Collaboration, and subsequently the Pierre Auger Observatory, recently reported confirmation of the expected Greisen-Zatsepin-Kuzmin (GZK) suppression of the UHECR flux above a few times 10(exp 19) eV. This observation is consistent with the majority of UHECR originating in astrophysical objects and reduces the need to invoke exotic physical processes. Particles observed above the GZK threshold energy must have come from sources within about 100 Mpc from the Earth. The small particle deflection angles expected at UHECR energies, with standard assumptions of extragalactic magnetic fields, are on the order of 1 degree. Thus by observing particles above the GZK threshold with sufficient exposure, there is the potential of identifying and characterizing individual UHECR sources. Auger has reported significant anisotropy in the arrival directions of UHECR at energies above about 6 10(exp 19) eV observed in the South, and a correlation to AGN in the 12th VCV catalog, suggesting that the sources of UHECR are traced by the distribution of luminous matter in the Universe. However, with similar statistics and the same event selection criteria, HiRes observations in the North are consistent with isotropy. Extended observations by Auger-South, by Telescope Array in the North, and possibly by the proposed Auger-North, will further these investigations. However, much greater exposures will be required to fully identify individual sources and measure their cosmic ray spectra. A five-year OWL mission would deliver approximately 10(exp 6)sq km/sr/yr of exposure with full aperture reached at approx. 10(exp 19) eV. The baseline 3 m optical aperture OWL telescopes with 45 degree full field-of-view are easily accomodated on a single conventional launch vehicle. On orbit, the simultaneous viewing of the same volume of atmosphere allows for precise stereo event reconstruction, which is nearly independent of the inclination of the particle track and tolerant of atmospheric conditions. The availability of monocular operation provides increased reliability or can be employed to increase the instantaneous aperture. In stereo mode, OWL is also capable of fully reconstructing horizontal or upward-moving showers and so has unmatched sensitivity to neutrino-induced events. The details of the OWL mission will he presented and its science capabilities will be discussed.

Mitchell, John↗

Fluctuations in heavy ion collisions and global conservation effects

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Infeasibility of Experimental Quantification of Life-Critical Software Reliability

This paper affirms that quantification of life-critical software reliability is infeasible using statistical methods whether applied to standard software or fault-tolerant software. The key assumption of software fault tolerance|separately programmed versions fail independently|is shown to be problematic. This assumption cannot be justified by experimentation in the ultra-reliability region and subjective arguments in its favor are not sufficiently strong to justify it as an axiom. Also, the implications of the recent multi-version software experiments support this affirmation.

Butler, Ricky W.↗

Improving statistical precision in Monte Carlo samples with negative weights via reweighting and uncertainty quantification

High statistical precision is critical for Monte Carlo (MC) samples in high energy physics and is degraded by negatively weighted events. This paper investigates a procedure to learn the relationship between the negative and positive weight distributions of any sample, allowing the reduction of statistical uncertainty by reweighting kinematically equivalent events with the same sign. A robust uncertainty quantification method is required for the practical application of such method. Two methods for the estimation of the reweighting uncertainty are developed: one at the event and another one at the final observable level. The latter method is strongly favored. The gains in statistical precision are then quantified. The method is demonstrated on Sherpa vector boson plus jets samples when using all generated events and when restricted to the signal region of a mock analysis. It is demonstrated to significantly reduce stochastic behavior in sparse MC samples while decreasing the overall uncertainty with a sufficiently well-known reweighting function.

Monte Carlo methods↗

Distributed Acoustic Sensing to Estimate the Permeability

Optical fiber in a borehole can be interrogated with distributed acoustic sensors (DAS) to capture fracture displacements with the potential to map surrounding fracture networks. We designed a laboratory experiment to test the capability of DAS to determine borehole flow characteristics, and we show that for the first time DAS can be used to remotely estimate permeability. Optical fiber was wrapped around a bead filled pipe and the pressure drop and flow velocity were measured to directly calculate permeability. A machine learning model using statistical features from continuous DAS estimated the bulk permeability. Fluid interactions with the permeable material demonstrate insufficient resolution using DAS amplitude-based measurements for estimating pressure drop to infer permeability. Variations in the spectral domain relate DAS measurements to the pressure drop and provide consistent permeability estimates. Resolution with DAS is sufficient to estimate permeability and provides a reliable method to monitor at depth in borehole conditions.

58 GEOSCIENCES↗

Investigation of several aspects of LANDSAT-4 data quality

The Thematic Mapper scene of Sacramento, CA acquired during the TDRSS test was received in TIPS format. Quadrants for both scenes were tested for band-to-band registration using reimplemented block correlation techniques. Summary statistics for band-to-band registrations of TM band combinations for Quadrant 4 of the NE Arkansas scene in TIPS format are tabulated as well as those for Quadrant 1 of the Sacramento scene. The system MTF analysis for the San Francisco scene is completed. The thermal band did not have sufficient contrast for the targets used and was not analyzed.

Wrigley, R. C.↗

Development of a satellite SAR image spectra and altimeter wave height data assimilation system for ERS-1

The applicability of ERS-1 wind and wave data for wave models was studied using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR data. A series of global wave hindcasts is made for the surface stress and surface wind fields by assimilation of scatterometer data for the full 96-day SEASAT and also for two wind field analyses for shorter periods by assimilation with the higher resolution ECMWF T63 model and by subjective analysis methods. It is found that wave models respond very sensitively to inconsistencies in wind field analyses and therefore provide a valuable data validation tool. Comparisons between SEASAT SAR image spectra and theoretical SAR spectra derived from the hindcast wave spectra by Monte Carlo simulations yield good overall agreement for 32 cases representing a wide variety of wave conditions. It is concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra with confidence for wave studies if supported by realistic wave models and theoretical computations of the strongly nonlinear mapping of the wave spectrum into the SAR image spectrum. A closed nonlinear integral expression for this spectral mapping relation is derived which avoids the inherent statistical errors of Monte Carlo computations and may prove to be more efficient numerically.

Hasselmann, Klaus↗

Validation of TRMM Satellite Rainfall Products Over Oklahoma for a Three Year Period (1998-2000)

This study evaluates the level 2 and level 3 rainfall products from the Tropical Rainfall Measuring Mission gridded distribution of 118 high resolution tipping buckets from the Oklahoma Mesonet. The Tropical Rainfall Measuring Mission, a joint satellite mission between NASA and the National Space Development Agency (NASDA) of Japan, was designed to estimate global precipitation between 40 S and 40 N latitude. The TRMM satellite, consisting of two main precipitation sensors, a passive microwave (TMI) and precipitation radar (PR) sensors, was launched in 1997. Although the great advantage of TRMM has been its ability to sample precipitation globally over the tropical oceans in places where ground sensors do not exist, regional points over land still offer the best opportunity for validating these rain estimates. For this study, the gauge data was gridded at I x I degree resolution between 34 and 36 N and -95 and -100 W. The location of Oklahoma is somewhat unique, in that, it is located near the turning point for the satellite. This study investigates effects of temporal sampling on the satellite measurements and the resulting rainfall bias observed between space and ground based sensors. The first part of the study analyzes the monthly rainfall statistics of the 3AI2 (TMI), the 3A25 (PR) and the 31142 (combined sensors) with the rain gauges. It uses regression techniques and a computation of the probability density function (PDF) and cumulative probability density function (CDF) for each sensor to evaluate the effects of temporal sampling on the satellite products and the sufficiency of using two overpasses a day to estimate rainfall on a monthly scale. In further probing this issue, the study also looked at the direct coincidence of ground and satellite measurements by comparing instantaneous comparisons obtained from the level 2 orbital track data, e.g. 2A12, 2A25. Lastly, the effects of temporal sampling were further studied by sub-sampling the gauge data only at overpass times. This methodology provides a quantitative way of separately inferring the fraction of bias due to temporal sampling and the rainfall algorithm. In addition to showing the salient effects of temporal sampling, the study also revisits the historical problem of correlating point gauge estimates with area estimates, though, the problem is ameliorated by averaging and gridding the data to the area resolution scale of the satellite.

Fisher, Brad↗

Redshifts of radio sources in the Million Quasars Catalogue from machine learning

ABSTRACT With the aim of using machine learning techniques to obtain photometric redshifts based upon a source’s radio spectrum alone, we have extracted the radio sources from the Million Quasars Catalogue. Of these, 44 119 have a spectroscopic redshift, required for model validation, and for which photometry could be obtained. Using the radio spectral properties as features, we fail to find a model which can reliably predict the redshifts, although there is the suggestion that the models improve with the size of the training sample. Using the near-infrared–optical–ultraviolet bands magnitudes, we obtain reliable predictions based on the 12 503 radio sources which have all of the required photometry. From the 80:20 training–validation split, this gives only 2501 validation sources, although training the sample upon our previous SDSS model gives comparable results for all 12 503 sources. This makes us confident that SkyMapper, which will survey southern sky in the u, v, g, r, i, z bands, can be used to predict the redshifts of radio sources detected with the Square Kilometre Array. By using machine learning to impute the magnitudes missing from much of the sample, we can predict the redshifts for 32 698 sources, an increase from 28 to 74 per cent of the sample, at the cost of increasing the outlier fraction by a factor of 1.4. While the ‘optical’ band data prove successful, at this stage we cannot rule out the possibility of a radio photometric redshift, given sufficient data which may be necessary to overcome the relatively featureless radio spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗

Integrating Formal Methods and Testing 2002

Traditionally, qualitative program verification methodologies and program testing are studied in separate research communities. None of them alone is powerful and practical enough to provide sufficient confidence in ultra-high reliability assessment when used exclusively. Significant advances can be made by accounting not only tho formal verification and program testing. but also the impact of many other standard V&V techniques, in a unified software reliability assessment framework. The first year of this research resulted in the statistical framework that, given the assumptions on the success of the qualitative V&V and QA procedures, significantly reduces the amount of testing needed to confidently assess reliability at so-called high and ultra-high levels (10-4 or higher). The coming years shall address the methodologies to realistically estimate the impacts of various V&V techniques to system reliability and include the impact of operational risk to reliability assessment. Combine formal correctness verification, process and product metrics, and other standard qualitative software assurance methods with statistical testing with the aim of gaining higher confidence in software reliability assessment for high-assurance applications. B) Quantify the impact of these methods on software reliability. C) Demonstrate that accounting for the effectiveness of these methods reduces the number of tests needed to attain certain confidence level. D) Quantify and justify the reliability estimate for systems developed using various methods.

Cukic, Bojan↗

A hierarchical Bayesian model for background variation in radiation source localization

In this paper, we apply a new model to account for varying background radiation in radiological source localization. We present a hierarchical Bayesian model that simultaneously infers background and source location parameters without requiring separate estimation of the background radiation at each detector location. We employ a simplified photon transport model to reduce the computational expense of Bayesian model calibration. Furthermore, we demonstrate the model accuracy by localizing a cesium-137 source in a simulated city block, and we analyze experimental field measurements with varying background. In both cases, the model provides sufficient fidelity that we can locate the source while simultaneously estimating background radiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗