Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 883 records · Page 49

Ocean Wave Studies with Applications to Ocean Modeling and Improvement of Satellite Altimeter Measurements

Combining analysis of satellite data (altimeter, scatterometer, high-resolution visible and infrared images, etc.) with mathematical modeling of non-linear wave processes, we investigate various ocean wave fields (on scales from capillary to planetary), their role in ocean dynamics and turbulent transport (of heat and biogeochemical quantities), and their effects on satellite altimeter measuring accuracy. In 1998 my attention was focused on long internal gravity waves (10 to 1000 km), known also as baroclinic inertia-gravity (BIG) waves. We found these waves to be a major factor of altimeter measurements "noise," resulting in a greater uncertainty [up to 10 cm in terms of sea surface height (SSH) amplitude] in the measured SSH signal than that caused by the sea state bias variations (up to 5 cm or so). This effect still remains largely overlooked by the satellite altimeter community. Our studies of BIG waves address not only their influence on altimeter measurements but also their role in global ocean dynamics and in transport and turbulent diffusion of biogeochemical quantities. In particular, in collaboration with Prof Peter Weichman, Caltech, we developed a theory of turbulent diffusion caused by wave motions of most general nature. Applied to the problem of horizontal turbulent diffusion in the ocean, the theory yielded the effective diffusion coefficient as a function of BIG wave parameters obtainable from satellite altimeter data. This effort, begun in 1997, has been successfully completed in 1998. We also developed a theory that relates spatial fluctuations of scalar fields (such as sea surface temperature, chlorophyll concentration, drifting ice concentration, etc.) to statistical characteristics of BIG waves obtainable from altimeter measurements. A manuscript is in the final stages of preparation. In order to verify the theoretical predictions and apply them to observations, we are now analyzing Sea-viewing Wide Field of view Sensor (SeaWiFS) and Field of view Sensor (SeaWiFS) and Advanced Very High-Resolution Radiometer (AVHRR) data on sea surface temperature (SST) and chlorophyll concentration jointly with TOPEX/POSEIDON data on SSH variations.

Glazman, Roman E.↗

Digital zero noise extrapolation for quantum error mitigation

Zero-noise extrapolation (ZNE) is an increasingly popular technique for mitigating errors in noisy quantum computations without using additional quantum resources. We review the fundamentals of ZNE and propose several improvements to noise scaling and extrapolation, the two key components in the technique. We introduce unitary folding and parameterized noise scaling. Furthermore, these are digital noise scaling frameworks, i.e. one can apply them using only gate-level access common to most quantum instruction sets. We also study different extrapolation methods, including a new adaptive protocol that uses a statistical inference framework. Benchmarks of our techniques show error reductions of 18X to 24X over non-mitigated circuits and demonstrate ZNE effectiveness at larger qubit numbers than have been tested previously. In addition to presenting new results, this work is a self-contained introduction to the practical use of ZNE by quantum programmers.

quantum computing↗

Ellicott City Disasters III: Building a Real-Time Statistical Flood Model for Improving Early Warning Systems in Ellicott City, Maryland

As flood events in the United States grow in frequency and intensity, the uses of applied remote sensing analyses are increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP Ellicott City Disasters III project investigated the use of machine learning for applications in flood risk detection to support the improvement of early warning systems in Ellicott City, Maryland. To strengthen the efforts of the Howard County Office of Emergency Management (OEM) in building a more robust flood monitoring system, this term built on the predictive capability of the long-short term memory (LSTM) model created by the second term of this DEVELOP project to create a Sequentially Trained Real-time EstimAted Model (STREAM). Enhancements to the model included the integration of both real-time and predicted weather products from the National Weather Service to increase predictive capacity. These weather products were supplemented by stream gauge data from the OEM as well as real-time radar products. The resultant flood risk model was trained to evaluate input variables and predict stage height in Ellicott City in real time. The model, upgraded to predict stage height up to 8 hours in advance, was incorporated into an online dashboard in a user-friendly interface. The project demonstrated the potential for integration of open data and NASA Earth observations into a flood risk forecasting tool capable of informing real-time decision-making.

Erika Munshi↗

Nonlinear Navigation Analysis for a Europa Lander Mission Concept

Navigation strategies typically use linearized techniques to predict maneuver statistics. This method has been successfully employed on numerous high-energy trajectories, but may not be valid for low-energy ones due to the highly non-linear environment. We explore nonlinear navigational techniques applied to a Europa Lander concept by assessing trajectory sensitivity. The linearity assumption is tested at different phases of the endgame trajectory. Maneuver strategies and placement are considered to improve the delivery to Europa. Several cases of non-linear statistical maneuver Monte Carlo simulations are conducted to optimize maneuvers and compare ∆V statistics against conventional linear simulations.

Hernandez, Sonia↗

A comparison of model validation approaches for echo state networks using climate model replicates

As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. Here, in this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.

54 ENVIRONMENTAL SCIENCES↗

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗

Atomic resolution coherent x-ray imaging with physics-based phase retrieval

Coherent x-ray imaging and scattering from accelerator based sources such as synchrotrons continue to impact biology, medicine, technology, and materials science. Many synchrotrons around the world are currently undergoing major upgrades to increase their available coherent x-ray flux by approximately two orders of magnitude. The improvement of synchrotrons may enable imaging of materials in operando at the atomic scale which may revolutionize battery and catalysis technologies. Current algorithms used for phase retrieval in coherent x-ray imaging are based on the projection onto sets method. These traditional iterative phase retrieval methods will become more computationally expensive as they push towards atomic resolution and may struggle to converge. Additionally, these methods do not incorporate physical information that may additionally constrain the solution. In this work, we present an algorithm which incorporates molecular dynamics into Bragg coherent diffraction imaging (BCDI). This algorithm, which we call PRAMMol (Phase Retrieval with Atomic Modeling and Molecular Dynamics) combines statistical techniques with molecular dynamics to solve the phase retrieval problem. We present several examples where our algorithm is applied to simulated coherent diffraction from 3D crystals and show convergence to the correct solution at the atomic scale.

47 OTHER INSTRUMENTATION↗

Multiscale Flow for robust and optimal cosmological analysis

We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as machine learning–based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of the muon decay-in-flight in the τ - → μ - ν ¯ μ ν τ decay to measure the Michel parameter ξ '

We present the first measurement of the Michel parameter ξ' in the $\tau$ - → μ¯ $\overline{v}$ μ v $\tau$ decay using the full data sample of 988 fb -1 collected by the Belle detector operating at the KEKB asymmetric energy e + ⁢e - collider. The method is based on the reconstruction of the μ¯ → e⁢¯$\overline{v}$ e v μ decay-in-flight in the Belle central drift chamber and relies on the correlation between muon spin and its daughter electron momentum. We study the main sources of the background that can imitate the signal decay, such as kaon and pion decays-in-flight and charged particle scattering on the detector material. Highly efficient methods of their suppression are developed and applied to select 165 signal-candidate events. We obtain ξ' = 0.22±0.94±0.42 where the first uncertainty is statistical, and the second one is systematic. The result is in agreement with the Standard Model prediction of ξ' = 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A constrained discrete optimal guidance strategy for low-thrust space flight

A low-thrust guidance algorithm suitable for operational use is formulated. Using linear stochastic control theory, a linear feedback control law is obtained for a minimum terminal miss performance criterion. The control variables are restricted to constant values over specified time intervals and are explicitly limited. Both fixed and variable time-of-arrival guidance are considered. Guidance equations are presented in deterministic (operational) and in statistical form; and computational procedures necessary for implementation are outlined. The performance of the guidance algorithm is evaluated by applying it to the approach phase of the 1980 rendezvous mission with the comet Encke.

Jacobson, R. A.↗

Strong scintillations in astrophysics. I - The Markov approximation, its validity and application to angular broadening

The Markov approximation to the propagation of waves in an extended, irregular medium is discussed in an astrophysical context. A new derivation is presented which is simple and which shows that the assumption of Gaussian statistics used by previous authors is irrelevant. We discuss the relevance of the approximation and show that it may apply in many situations of interest, including interstellar scintillations of pulsar signals. The approximation does not require the assumption of weak scattering or Gaussian correlation functions. The Markov equation for the angular spectrum is particularly simple, and solutions are discussed for typical turbulence spectra. It is found that the equation for the angular spectrum is very nearly that used by previous authors, and the present discussion shows that these results are much more general than previously thought. A possible observational test for distinguishing between Gaussian and power-law interstellar density spectra is discussed.

Lee, L. C.↗

Thematic mapper data quality and performance assessment in renewable resource/agricultural remote sensing

A "quick look" investigation of the initial LANDSAT-4, thematic mapper (TM) scene received from Goddard Space Flight Center was performed to gain early insight into the characteristics of TM data. The initial scene, containing only the first four bands of the seven bands recorded by the TM, was acquired over the Detroit, Michigan, area on July 20, 1982. It yielded abundant information for scientific investigation. A wide variety of studies were conducted to assess all aspects of TM data. They ranged from manual analyses of image products to detect obvious optical, electronic, or mechanical defects to detailed machine analyses of the digital data content for evaluation of spectral separability of vegetative/nonvegetative classes. These studies were applied to several segments extracted from the full scene. No attempt was made to perform end-to-end statistical evaluations. However, the output of these studies do identify a degree of positive performance from the TM and its potential for advancing state-of-the-art crop inventory and condition assessment technology.

Erickson, J. D.↗

Texture classification using autoregressive filtering

A general theory of image texture models is proposed and its applicability to the problem of scene segmentation using texture classification is discussed. An algorithm, based on half-plane autoregressive filtering, which optimally utilizes second order statistics to discriminate between texture classes represented by arbitrary wide sense stationary random fields is described. Empirical results of applying this algorithm to natural and sysnthesized scenes are presented and future research is outlined.

Lawton, W. M.↗

Context dependent anti-aliasing image reconstruction

Image Reconstruction has been mostly confined to context free linear processes; the traditional continuum interpretation of digital array data uses a linear interpolator with or without an enhancement filter. Here, anti-aliasing context dependent interpretation techniques are investigated for image reconstruction. Pattern classification is applied to each neighborhood to assign it a context class; a different interpolation/filter is applied to neighborhoods of differing context. It is shown how the context dependent interpolation is computed through ensemble average statistics using high resolution training imagery from which the lower resolution image array data is obtained (simulation). A quadratic least squares (LS) context-free image quality model is described from which the context dependent interpolation coefficients are derived. It is shown how ensembles of high-resolution images can be used to capture the a priori special character of different context classes. As a consequence, a priori information such as the translational invariance of edges along the edge direction, edge discontinuity, and the character of corners is captured and can be used to interpret image array data with greater spatial resolution than would be expected by the Nyquist limit. A Gibb-like artifact associated with this super-resolution is discussed. More realistic context dependent image quality models are needed and a suggestion is made for using a quality model which now is finding application in data compression.

Beaudet, Paul R.↗

Validation of the AMSU-B Bias Corrections Based on Satellite Measurements from SSM/T-2

The NOAA-15 Advanced Microwave Sounding Unit-B (AMSU-B) was designed in the same spirit as the Special Sensor Microwave Water Vapor Profiler (SSM/T-2) on board the DMSP F11-14 satellites, to perform remote sensing of spatial and temporal variations in mid and upper troposphere humidity. While the SSM/T-2 instruments have a 48 km spatial resolution at nadir and 28 beam positions per scan, AMSU-B provides an improvement with a 16 km spatial resolution at nadir and 90 beam positions per scan. The AMSU-B instrument, though, has been experiencing radio frequency interference (RFI) contamination from the NOAA-15 transmitters whose effect is dependent upon channel, geographic location, and current spacecraft antenna configuration. This has lead to large cross-track biases reaching as high as 100 Kelvin for channel 17 (150 GHz) and 50 Kelvin for channel 19 (183 +/-3 GHz). NOAA-NESDIS has recently provided a series of bias corrections for AMSU-B data starting from March, 1999. These corrections are available for each of the five channels, for every third field of view, and for three cycles within an eight second period. There is also a quality indicator in each data record to indicate whether or not the bias corrections should be applied. As a precursor to performing retrievals of mid and upper troposphere humidity, a validation study is performed by statistically analyzing the differences between the F14 SSM/T-2 and the bias corrected AMSU-B brightness temperatures for three months in the spring of 1999.

Kolodner, Marc A.↗

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗

Transmission measurement based on STM observation to detect the penetration depth of low-energy heavy ions in botanic samples

The penetration depth of low-energy heavy ions in botanic samples was detected with a new transmission measurement. In the measurement, highly oriented pyrolytic graphite (HOPG) pieces were placed behind the botanic samples with certain thickness. During the irradiation of heavy ions with energy of tens of keV, the energetic particles transmitted from those samples were received by the HOPG pieces. After irradiation, scanning tunneling microscope (STM) was applied to observe protrusion-like damage induced by these transmitted ions on the surface of the HOPG. The statistical average number density of protrusions and the minimum transmission rate of the low-energy heavy ions can be obtained. The detection efficiency of the new method for low-energy heavy ions was about 0.1-1 and the background in the measurement can be reduced to as low as 1.0 x 10(8) protrusions/cm2. With this method, the penetration depth of the energetic particles was detected to be no less than 60 micrometers in kidney bean slices when the slices were irradiated by 100 keVAr+ ion at the fluence of 5 x 10(16) ions/cm2. c2002 Elsevier Science Ltd. All rights reserved.

DNA Damage↗