SEARCH · Engineering Papers

Results for “variance estimation”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1995

On-line estimation of error covariance parameters for atmospheric data assimilation

A simple scheme is presented for on-line estimation of covariance parameters in statistical data assimilation systems. The scheme is based on a maximum-likelihood approach in which estimates are produced on the basis of a single batch of simultaneous observations. Simple-sample covariance estimation is reasonable as long as the number of available observations exceeds the number of tunable parameters by two or three orders of magnitude. Not much is known at present about model error associated with actual forecast systems. Our scheme can be used to estimate some important statistical model error parameters such as regionally averaged variances or characteristic correlation length scales. The advantage of the single-sample approach is that it does not rely on any assumptions about the temporal behavior of the covariance parameters: time-dependent parameter estimates can be continuously adjusted on the basis of current observations. This is of practical importance since it is likely to be the case that both model error and observation error strongly depend on the actual state of the atmosphere. The single-sample estimation scheme can be incorporated into any four-dimensional statistical data assimilation system that involves explicit calculation of forecast error covariances, including optimal interpolation (OI) and the simplified Kalman filter (SKF). The computational cost of the scheme is high but not prohibitive; on-line estimation of one or two covariance parameters in each analysis box of an operational bozed-OI system is currently feasible. A number of numerical experiments performed with an adaptive SKF and an adaptive version of OI, using a linear two-dimensional shallow-water model and artificially generated model error are described. The performance of the nonadaptive versions of these methods turns out to depend rather strongly on correct specification of model error parameters. These parameters are estimated under a variety of conditions, including uniformly distributed model error and time-dependent model error statistics.

Dee, Dick P.↗

NASA NTRSDate not supplied

From Formal Errors Towards Realistic Uncertainties

Evaluation of uncertainties of geodetic parameter estimates is the problem that is not yet solved in a satisfactory way. A direct evaluation of the uncertainties derived from a least square solution is labeled "formal" and is usually biased, sometimes up to an order of magnitude. Customary, the use of formal errors for scientific analysis is discouraged. We claim that the root of the problem is neglecting off-diagonal elements in the variance-covariance matrix of the noise in the data. A careful reconstruction of the full variance-covariance matrix, including the off-diagonal terms greatly improves realism of uncertainty estimates derived from least squares. We processed the dataset of VLBI group delays and built a prior variance-covariance of the atmosphere-driven noise based on analysis of the output of NASA high-resolution numerical weather models. We found that the uncertainties of parameter estimates derived from this least square solution that uses such variance-covariance matrices become much closer to realistic errors. We consider approaches for implementation of this method in routine data analysis of space geodesy data.

Leonid Petrov↗

DOE OSTI2023

Merging flux-variance with surface renewal methods in the roughness sublayer and the atmospheric surface layer

In this article, Two micrometeorological methods utilizing high-frequency sampled air temperature were tested against eddy covariance (EC) sensible heat flux (𝐻) measurements at three sites representing agricultural, agro-forestry, and forestry systems. The two methods cover conventional and newly proposed forms of the flux-variance (FV) and surface renewal (SR) schemes of differing complexities. The sites represent measurements in surface, roughness, and roughness to surface transitional layers. Regression analyzes against EC show that the most reliable FV and SR forms estimate 𝐻 with slopes within ±10% from unity and coefficient of determination 𝑅2 > 0.9 across all the three sites. The best performance of both FV and SR was found at the agricultural site with measurements well within the surface layer, while the worst was found for the tall forest with measurements within the roughness sublayer where its thickness needed to be additionally estimated. The main variable driving 𝐻 in FV is the temperature variance, whereas in SR, it is the geometry of ramp-like structures. Since these structures are also responsible for most of the temperature variance, a novel FV-SR approach emerging from combining the methods is proposed and evaluated against EC measurements and conventional FV and SR schemes. The proposed FV-SR approach requiring only a single fast response thermocouple is potentially independent of calibration and ameliorates some of the theoretical objections that arise when combining ramp statistics with similarity arguments. The combination of methods also provides new insights into the contribution of coherent structures to the temperature variance and its dependence on atmospheric stratification. Other potential utility of the new method is to include it in multi-tool assessments of surface energy fluxes, since a convergence or divergence of the results has a high diagnostic value.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

Rainfall Product Evaluation for the TRMM Ground Validation Program

Evaluation of the Tropical Rainfall Measuring Mission (TRMM) satellite observations is conducted through a comprehensive Ground Validation (GV) Program. Standardized instantaneous and monthly rainfall products are routinely generated using quality-controlled ground based radar data from four primary GV sites. As part of the TRMM GV program, effort is being made to evaluate these GV products and to determine the uncertainties of the rainfall estimates. The evaluation effort is based on comparison to rain gauge data. The variance between the gauge measurement and the true averaged rain amount within the radar pixel is a limiting factor in the evaluation process. While monthly estimates are relatively simple to evaluate, the evaluation of the instantaneous products are much more of a challenge. Scattegrams of point comparisons between radar and rain gauges are extremely noisy for several reasons (e.g. sample volume discrepancies, timing and navigation mismatches, variability of Z(sub e)-R relationships), and therefore useless for evaluating the estimates. Several alternative methods, such as the analysis of the distribution of rain volume by rain rate as derived from gauge intensities and from reflectivities above the gauge network will be presented. Alternative procedures to increase the accuracy of the estimates and to reduce their uncertainties also will be discussed.

Amitai, E.↗

NASA NTRS2016

How Universal is the Relationship Between Remotely Sensed Vegetation Indices and Crop Leaf Area Index? A Global Assessment

Leaf Area Index (LAI) is a key variable that bridges remote sensing observations to the quantification of agroecosystem processes. In this study, we assessed the universality of the relationships between crop LAI and remotely sensed Vegetation Indices (VIs). We first compiled a global dataset of 1459 in situ quality-controlled crop LAI measurements and collected Landsat satellite images to derive five different VIs including Simple Ratio (SR), Normalized Difference Vegetation Index (NDVI), two versions of the Enhanced Vegetation Index (EVI and EVI2), and Green Chlorophyll Index (CI(sub Green)). Based on this dataset, we developed global LAI-VI relationships for each crop type and VI using symbolic regression and Theil-Sen (TS) robust estimator. Results suggest that the global LAI-VI relationships are statistically significant, crop-specific, and mostly non-linear. These relationships explain more than half of the total variance in ground LAI observations (R2 greater than 0.5), and provide LAI estimates with RMSE below 1.2 m2/m2. Among the five VIs, EVI/EVI2 are the most effective, and the crop-specific LAI-EVI and LAI-EVI2 relationships constructed by TS, are robust when tested by three independent validation datasets of varied spatial scales. While the heterogeneity of agricultural landscapes leads to a diverse set of local LAI-VI relationships, the relationships provided here represent global universality on an average basis, allowing the generation of large-scale spatial-explicit LAI maps. This study contributes to the operationalization of large-area crop modeling and, by extension, has relevance to both fundamental and applied agroecosystem research.

Vegetation Index↗

DOE OSTI2021

Multiscale assessment of land surface phenology from harmonized Landsat 8 and Sentinel-2, PlanetScope, and PhenoCam imagery

As the spatial and temporal resolution of remotely sensed imagery has improved over the last four decades, algorithms for monitoring and mapping seasonal changes in surface properties have evolved rapidly. Most recently, the availability of daily PlanetScope imagery has created new opportunities for monitoring the land surface phenology (LSP) of terrestrial ecosystems at high spatial resolution. However, the quality and value of LSP information from PlanetScope imagery have not been systematically examined. In this paper, we evaluate the character and quality of LSP information derived from PlanetScope by comparing time series of vegetation indices and LSP metrics from PlanetScope to corresponding time series and LSP metrics derived from Harmonized Landsat 8 and Sentinel-2 (HLS) imagery and PhenoCams at six sites that span a diverse range of land cover types and climate. Results show that vegetation index time series from all three data sources show high temporal correlation, and LSP metrics derived from HLS, PlanetScope, and PhenoCam show high agreement with negligible bias. Semi-variograms for phenometrics estimated from PlanetScope imagery indicate that the majority of spatial variance captured in PlanetScope phenometrics occurs well below the spatial resolution HLS imagery. At the same time, LSP metrics from HLS are most strongly correlated with the 50–75% quantiles of 3 m LSP metrics from PlanetScope. This indicates that HLS captures the average phenology at sub-pixel scale captured in PlanetScope imagery. Here, our results represent the first comprehensive comparison of LSP metrics estimated from PlanetScope and publicly available moderate spatial resolution imagery, and provide insights regarding: (1) the quality and character of LSP metrics derived from HLS and PlanetScope; and (2) the relative merits and trade-offs associated with the use of each data source for LSP studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

An error reduction algorithm to improve lidar turbulence estimates for wind energy

Abstract. Remote-sensing devices such as lidars are currently being investigated as alternatives to cup anemometers on meteorological towers for the measurement of wind speed and direction. Although lidars can measure mean wind speeds at heights spanning an entire turbine rotor disk and can be easily moved from one location to another, they measure different values of turbulence than an instrument on a tower. Current methods for improving lidar turbulence estimates include the use of analytical turbulence models and expensive scanning lidars. While these methods provide accurate results in a research setting, they cannot be easily applied to smaller, vertically profiling lidars in locations where high-resolution sonic anemometer data are not available. Thus, there is clearly a need for a turbulence error reduction model that is simpler and more easily applicable to lidars that are used in the wind energy industry. In this work, a new turbulence error reduction algorithm for lidars is described. The Lidar Turbulence Error Reduction Algorithm, L-TERRA, can be applied using only data from a stand-alone vertically profiling lidar and requires minimal training with meteorological tower data. The basis of L-TERRA is a series of physics-based corrections that are applied to the lidar data to mitigate errors from instrument noise, volume averaging, and variance contamination. These corrections are applied in conjunction with a trained machine-learning model to improve turbulence estimates from a vertically profiling WINDCUBE v2 lidar. The lessons learned from creating the L-TERRA model for a WINDCUBE v2 lidar can also be applied to other lidar devices. L-TERRA was tested on data from two sites in the Southern Plains region of the United States. The physics-based corrections in L-TERRA brought regression line slopes much closer to 1 at both sites and significantly reduced the sensitivity of lidar turbulence errors to atmospheric stability. The accuracy of machine-learning methods in L-TERRA was highly dependent on the input variables and training dataset used, suggesting that machine learning may not be the best technique for reducing lidar turbulence intensity (TI) error. Future work will include the use of a lidar simulator to better understand how different factors affect lidar turbulence error and to determine how these errors can be reduced using information from a stand-alone lidar.

17 WIND ENERGY↗

NASA NTRS2023

Testing Cosmic Microwave Background Anomalies in E-mode Polarization with Current and Future Data

In this paper, we explore the power of the cosmic microwave background (CMB) polarization (E-mode) data to corroborate four potential anomalies in CMB temperature data: the lack of large angular-scale correlations, the alignment of the quadrupole and octupole (Q–O), the point-parity asymmetry, and the hemispherical power asymmetry. We use CMB simulations with noise representative of three experiments—the Planck satellite, the Cosmology Large Angular Scale Surveyor (CLASS), and the LiteBIRD satellite—to test how current and future data constrain the anomalies. We find the correlation coefficients ρ between temperature and E-mode estimators to be less than 0.1, except for the point-parity asymmetry (ρ = 0.17 for cosmic-variance-limited simulations), confirming that E-modes provide a check on the anomalies that is largely independent of temperature data. Compared to Planck component-separated CMB data (smica), the putative LiteBIRD survey would reduce errors on E-mode anomaly estimators by factors of ∼3 for hemispherical power asymmetry and point-parity asymmetry, and by ∼26 for lack of large-scale correlation. The improvement in Q–O alignment is not obvious due to large cosmic variance, but we found the ability to pin down the estimator value will be improved by a factor ≳100. Improvements with CLASS are intermediate to these.

Cosmic microwave background radiation↗

DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

Validity conditions for moment closure approximations in stochastic chemical kinetics

Approximations based on moment-closure (MA) are commonly used to obtain estimates of the mean molecule numbers and of the variance of fluctuations in the number of molecules of chemical systems. The advantage of this approach is that it can be far less computationally expensive than exact stochastic simulations of the chemical master equation. Here, we numerically study the conditions under which the MA equations yield results reflecting the true stochastic dynamics of the system. We show that for bistable and oscillatory chemical systems with deterministic initial conditions, the solution of the MA equations can be interpreted as a valid approximation to the true moments of the chemical master equation, only when the steady-state mean molecule numbers obtained from the chemical master equation fall within a certain finite range. The same validity criterion for monostable systems implies that the steady-state mean molecule numbers obtained from the chemical master equation must be above a certain threshold. For mean molecule numbers outside of this range of validity, the MA equations lead to either qualitatively wrong oscillatory dynamics or to unphysical predictions such as negative variances in the molecule numbers or multiple steady-state moments of the stationary distribution as the initial conditions are varied. Our results clarify the range of validity of the MA approach and show that pitfalls in the interpretation of the results can only be overcome through the systematic comparison of the solutions of the MA equations of a certain order with those of higher orders.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗

DOE OSTI2022

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Dimensionality Reduction Library v 0.2

Dimensionality Reduction Library is a C++ library for dimensionality reduction. In the context of this library, dimensionality reduction is considered to consist of 1)estimation of the intrinsic dimensionality using sampled data, 2) Finding maps that reduce the dimensionality of data (forward map) or increase the dimensionality of data (reverse map) and 3) mapping arbitray coordinates to high and low dimensionalities. The library is intended to provide a consistent interface to multiple dimensionality reduction algorithms with an efficient C++ interface that runs efficiently on multicore architectures. A few routines have been optimized with an option for GPU acceleration or distributed computation. Currently the library offers intrinsic dimensionality estimation using point-PCA reconstruction error and/ residual variance. The following dimensionality reduction methods have been implemented: Principal Component Analysis Multidimensional Scaling Locally Linear Embedding IsoMap Autoencoder Neutral Networks An executable is also supplied that can be built to allow for command-line access to the library routines. A description for an application of the library for molecular structure analysis has been published.

BROWN, WILLIAM↗

DOE OSTI2018

Tree Ring Chronology Indexes and Reconstructions of Precipitation in Central Iowa, USA (1984) (NDP-002)

Tree core samples (4 mm in diameter) were extracted from the trunks of white oak (Quercus alba) at three sites in central Iowa (Duvick Back Woods, Ledges State Park, and Pammel). At least 60 trees were sampled at each site, and at least two cores were taken from each tree. The growth rings of each core were dated by calendar year and measured; the measurements were then transformed into dimensionless ring-width indices and correlated with annual precipitation. Data were collected for the years 1680 through 1979. Each tree ring was characterized by the site, year, tree-ring-width index, number of core samples, decade year, and the annual reconstructed precipitation estimate. These data have more than 50% of their variance in common with the known annual statewide average precipitation for Iowa and serve as useful indicators of the precipitation and drought history of the region for the past 300 years. The data are in two files: tree-ring-chronology data (8 kB) and the annual reconstructed precipitation data for central Iowa (2 kB).For access to the data files, click this link to the CDIAC data transition website: http://cdiac.ess-dive.lbl.gov/ndps/ndp002.html

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Scalable Combinatorial Tools for Health Disparities Research

Despite staggering investments made in unraveling the human genome, current estimates suggest that as much as 90% of the variance in cancer and chronic diseases can be attributed to factors outside an individual’s genetic endowment, particularly to environmental exposures experienced across his or her life course. New analytical approaches are clearly required as investigators turn to complicated systems theory and ecological, place-based and life-history perspectives in order to understand more clearly the relationships between social determinants, environmental exposures and health disparities. While traditional data analysis techniques remain foundational to health disparities research, they are easily overwhelmed by the ever-increasing size and heterogeneity of available data needed to illuminate latent gene x environment interactions. This has prompted the adaptation and application of scalable combinatorial methods, many from genome science research, to the study of population health. Most of these powerful tools are algorithmically sophisticated, highly automated and mathematically abstract. Their utility motivates the main theme of this paper, which is to describe real applications of innovative transdisciplinary models and analyses in an effort to help move the research community closer toward identifying the causal mechanisms and associated environmental contexts underlying health disparities. The public health exposome is used as a contemporary focus for addressing the complex nature of this subject.

Environmental Sciences & Ecology↗

NASA NTRS1975

Status of input design for aircraft parameter identification

Results are presented on the design of aircraft inputs (i.e. elevator, rudder and aeleron deflection time histories) to identify aircraft stability and control derivatives from flight test data. The problem is first reduced to an optimization problem with differential and integral constraints. The criteria used are either expressed in terms of the Cramer-Rao lower bound on the covariance matrix of the parameter estimates or in terms of the maximum prediction error variance. Both time-domain longitudinal and lateral dynamics of C-8 and Jet Star aircrafts and comparison with doublet type inputs are made.

Mehra, R. K.↗

NASA NTRS1979

Estimating regional heat flux from scanning radiometer data

Method using matrix approach to determine best linear estimate of total flux consistent with condition of minimum variance proves useful in meteorological studies. It is useful for researchers studying options in scanner design such as scan pattern, scan rate, and radiometer field-of-view.

Bess, T. D.↗

NASA NTRS1978

A method for determining the weak statistical stationarity of a random process

A method for determining the weak statistical stationarity of a random process is presented. The core of this testing procedure consists of generating an equivalent ensemble which approximates a true ensemble. Formation of an equivalent ensemble is accomplished through segmenting a sufficiently long time history of a random process into equal, finite, and statistically independent sample records. The weak statistical stationarity is ascertained based on the time invariance of the equivalent-ensemble averages. Comparison of these averages with their corresponding time averages over a single sample record leads to a heuristic estimate of the ergodicity of a random process. Specific variance tests are introduced for evaluating the statistical independence of the sample records, the time invariance of the equivalent-ensemble autocorrelations, and the ergodicity. Examination and substantiation of these procedures were conducted utilizing turbulent velocity signals.

Sadeh, W. Z.↗