SEARCH · Engineering Papers

Results for “covariance and correlation”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates

Herein, the use of many control variates is proposed as a method to accelerate the second- and third-order Monte Carlo (MC) many-body perturbation (MC-MP2 and MC-MP3) calculations. A control variate is an exactly integrable function that is strongly correlated or anti-correlated with the target function to be integrated by the MC method. Evaluating both integrals as well as their covariances in the same MC run, one can effect a mutual cancellation of the statistical uncertainties and biases in the MC integrations, thereby accelerating its convergence considerably. Six and thirty-six control variates, whose integrals are known a priori, are generated for MC-MP2 and MC-MP3, respectively, by systematically replacing one or more two-electron-integral vertexes of certain configurations by zero-valued overlap-integral vertexes in their Goldstone diagrams. The variances and co- variances of these control variates are computed at a marginal cost, enhancing the overall efficiency of the MC-MP2 and MC-MP3 calculations by a factor of up to 14 and 20, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Event-by-event neutron–photon multiplicity correlations in 252Cf(sf)

Excited nuclear fragments are emitted during nuclear fission. The de-excitation of these fission fragments takes place as sequential emission of neutrons followed by photons. A correlation between neutron and photon multiplicities accompanying fission is thus expected. Fission event generators based on established statistical nuclear physics models predict a negative event-by-event correlation in neutron-photon multiplicity. A survey of published experimental results of an event-by-event covariance between the neutron - photon multiplicities emitted following the spontaneous fission of 252 Cf is presented. Analytic unfolding expressions are developed in this work to determine the bias introduced by background sources, particle misclassification, pulse pileup, and inelastic photon production. The published experimental data are re-analyzed using these unfolding techniques and are found to be in qualitative agreement with the predictions of model-based calculations. Finally, in particular, we have concluded that there exists a significant event-by-event neutron-photon emission competition following the spontaneous fission of 252 Cf.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1984

Autocovariance Computer

Laser-velocimeter covariance processor calculates autocovariance and cross-covariance functions for turbulent flow field, based on Poissonsampled measurements in time from laser velocimeter. Hardware implementation of correlation technique to laser velocimeter flow-field diagnostic system fast enough for online production applications.

Meyers, J. F.↗

NASA NTRS1997

The Influence of Correlated Crustal Signals in Modelling the Main Geomagnetic Field

Algorithms used in geomagnetic main-field modelling have for the most part treated the noise in the field measurements as if it were white. A major component of the noise consists of the field due to magnetization in the crust and it has been realized for some time that such signals are highly correlated at satellite altitude. Hence approximation by white noise, while of undoubted utility, is of unknown validity. In this paper we study two plausible statistical models for the crustal magnetization, in which the magnetization is a realization of a stationary, isotropic, random process. At a typical satellite altitude the associated fields exhibit significant correlation over ranges as great as 15 deg. or more, which introduces off-diagonal elements into the covariance matrix, elements that have usually been neglected in modelling procedures. Dealing with a full covariance matrix for a large data set would present a formidable computational challenge, but fortunately most of the entries in the covariance matrix are so small that they can be replaced by zeros. The resultant matrix comprises only about 3 per cent non-zero entries and thus we can take advantage of efficient sparse matrix techniques to solve the numerical system. We construct several main-field models based on vertical-component data from a selected 5 deg. by 5 deg. data set derived from the Magsat mission. Models with and without off-diagonal terms are compared.

Rygaard-Hjalsted, C.↗

NASA NTRS1993

The Effects of Noise on Polarimetric SAR Data

Polarimetric SAR data can provide a great deal of information about the scattering behavior of the surface under observation. Polarimetric SAR systems often measure the scattering matrices of the areas under observation in linear polarizations (H and V). From the scattering matrix commonly used forms such as the covariance matrix and the Stokes matrix can be easily derived. Other measures derived from polarimetric SAR data include correlation coefficients between scattering matrix terms and the mode and variance of phase differences between scattering matrix terms. The effects of additive system noise on these measurements is not often considered in the literature on this subject. In this paper, the effects of additive system noise on measurements derived from polarimetric SAR data will be examined. It will be shown how first-order noise effects can be removed and how second-order noise effects can be reduced for some measurements...

Freeman, A.↗

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 NTRS1998

The Importance of Semi-Major Axis Knowledge in the Determination of Near-Circular Orbits

Modem orbit determination has mostly been accomplished using Cartesian coordinates. This usage has carried over in recent years to the use of GPS for satellite orbit determination. The unprecedented positioning accuracy of GPS has tended to focus attention more on the system's capability to locate the spacecraft's location at a particular epoch than on its accuracy in determination of the orbit, per se. As is well-known, the latter depends on a coordinated knowledge of position, velocity, and the correlation between their errors. Failure to determine a properly coordinated position/velocity state vector at a given epoch can lead to an epoch state that does not propagate well, and/or may not be usable for the execution of orbit adjustment maneuvers. For the quite common case of near-circular orbits, the degree to which position and velocity estimates are properly coordinated is largely captured by the error in semi-major axis (SMA) they jointly produce. Figure 1 depicts the relationships among radius error, speed error, and their correlation which exist for a typical low altitude Earth orbit. Two familiar consequences are the relationship Figure 1 shows are the following: (1) downrange position error grows at the per orbit rate of 3(pi) times the SMA error; (2) a velocity change imparted to the orbit will have an error of (pi) divided by the orbit period times the SMA error. A less familiar consequence occurs in the problem of initializing the covariance matrix for a sequential orbit determination filter. An initial covariance consistent with orbital dynamics should be used if the covariance is to propagate well. Properly accounting for the SMA error of the initial state in the construction of the initial covariance accomplishes half of this objective, by specifying the partition of the covariance corresponding to down-track position and radial velocity errors. The remainder of the in-plane covariance partition may be specified in terms of the flight path angle error of the initial state. Figure 2 illustrates the effect of properly and not properly initializing a covariance. This figure was produced by propagating the covariance shown on the plot, without process noise, in a circular low Earth orbit whose period is 5828.5 seconds. The upper subplot, in which the proper relationships among position, velocity, and their correlation has been used, shows overall error growth, in terms of the standard deviations of the inertial position coordinates, of about half of the lower subplot, whose initial covariance was based on other considerations.

Carpenter, J. Russell↗

DOE OSTI2021

Systematic and statistical uncertainties of the hilbert-transform based high-precision FID frequency extraction method

Pulsed nuclear magnetic resonance (NMR) is widely used in high-precision magnetic field measurements. The absolute value of the magnetic field is determined from the precession frequency of nuclear magnetic moments. The Hilbert transform is widely used to extract the phase function from the observed free induction decay (FID) signal and then its frequency. In this paper, a detailed implementation of a Hilbert-transform based FID frequency extraction method is described. How artifacts and noise level in the FID signal affect the extracted phase function are derived analytically. A method of mitigating the artifacts in the extracted phase function of an FID is discussed. Correlations between noises of the phase function samples are studied for different noise spectra. We discovered that the error covariance matrix for the extracted phase function is nearly singular and improper for constructing the $\chi^2$ used in the fitting routine. A down-sampling method for fixing the singular covariance matrix has been developed, so that the minimum $\chi^2$-fit yields properly the statistical uncertainty of the extracted frequency. Other practical methods of obtaining the statistical uncertainty are also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Primordial non-Gaussianity from the completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey – I: Catalogue preparation and systematic mitigation

ABSTRACT We investigate the large-scale clustering of the final spectroscopic sample of quasars from the recently completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). The sample contains 343 708 objects in the redshift range 0.8 < z < 2.2 and 72 667 objects with redshifts 2.2 < z < 3.5, covering an effective area of $4699\, {\rm deg}^{2}$. We develop a neural network-based approach to mitigate spurious fluctuations in the density field caused by spatial variations in the quality of the imaging data used to select targets for follow-up spectroscopy. Simulations are used with the same angular and radial distributions as the real data to estimate covariance matrices, perform error analyses, and assess residual systematic uncertainties. We measure the mean density contrast and cross-correlations of the eBOSS quasars against maps of potential sources of imaging systematics to address algorithm effectiveness, finding that the neural network-based approach outperforms standard linear regression. Stellar density is one of the most important sources of spurious fluctuations, and a new template constructed using data from the Gaia spacecraft provides the best match to the observed quasar clustering. The end-product from this work is a new value-added quasar catalogue with the improved weights to correct for non-linear imaging systematic effects, which will be made public. Our quasar catalogue is used to measure the local-type primordial non-Gaussianity in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction [DI]) and how resistance to one drug impacts resistance to the other (collateral effects [CE]). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug-sensitive (ancestral) and drug-resistant (mutant) populations. We show that evolved resistance to the component drugs – and in turn, the adaptation of growth rate – is governed by a Price equation whose covariance terms encode geometric features of both the two-drug-response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multidrug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally explicit selective outcomes.

Gjini, Erida↗

NASA NTRS1991

Characterization of frequency standard instability by estimation of their covariance matrix

The popular 3-cornered hat method used for evaluating the noise contributions of individual frequency standards is revisited. This method is used in several cases, but sometimes the results are not consistent because one or more estimated clock variances turn out to be negative. Different causes of this unacceptable result have been conjectured: among them one regards the hypothesis of uncorrelated clocks, essential in this method. Since recently realistic cases of correlation between clocks, mainly due to the environmental conditions, have been observed, this paper proposes an entirely revisited version of the 3-cornered hat method which permits to evaluate the individual variances and also the possible covariances between clocks, by relaxing the hypothesis of uncorrelation. The uncertainty and the lack of contemporaneity of the measurement series are assumed to be negligible. The lack of the uncorrelation hypothesis calls for a more general mathematical model leading to an underdetermined linear system. The estimates of the (co)variances of the measurement series us well us those of the individual clocks are introduced by means of the scalar product of the related time series and arranged in the respective covariance matrices S and R. Since covariance matrix is positive definite by definition, the problem consists in estimating the unknown R, subject to the constraint of positive definiteness, from the known S. Unfortunately, this constraint is not sufficient to estimate R. Therefore a suitable optimization criterion is proposed, which assures the positive definiteness of R and, at the same time, minimizes the global correlation among clocks. Examples of frequency instability measurements processed by the "classical" 3-cornered hat method and the here-revisited method are presented showing that the solutions are identical only when the uncorrelation hypothesis doesn't violate the positive definiteness of R.

PatriziaTravella↗

DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗

NASA NTRS1985

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗

NASA NTRS2000

Effect of Medium Symmetries in Limiting the Number of Parameters Estimated with Polarimetric Interferometry

The addition of interferometric backscattering pairs to the conventional polarimetric SAR data over forests and other vegetated areas increases the dimensionality of the data space, in principle enabling the estimation of a larger number of vegetation parameters. Without regard to the sensitivity of these data to vegetation scattering parameters, this paper poses the question: Will increasing the data channels as such result in a one-to-one increase in the number of parameters that can be estimated, or do vegetation and data properties inherently limit that number otherwise? In this paper, the complete polarimetric interferometric covariance matrix is considered and various symmetry properties of the scattering medium are used to study whether any of the correlation pairs can be eliminated. The number of independent pairs has direct consequences in their utility in parameter estimation schemes, since the maximum number of parameters that can be estimated cannot exceed the number of unique measurements. The independent components of the polarimetric interferometric SAR (POL/INSAR) data are derived for media with reflection, rotation, and azimuth symmetries, which are often encountered in vegetated surfaces. Similar derivations have been carried out before for simple polarimetry, i.e., zero baseline. This paper extends those to the interferometric case of general nonzero baselines. It is shown that depending on the type of symmetries present, the number of independent available measurements that can be used to estimate medium parameters will vary. In particular, whereas in the general case there are 27 mathematically independent measurements possible from a polarimetric interferometer, this number can be reduced to 15, 9, and 6 if the medium has reflection, rotation, or azimuthal symmetries, respectively. The results can be used in several ways in the interpretation of SAR data and the development of parameter estimation schemes, which will be discussed at the presentation. Recent POL/INSAR data from the JPL AIRSAR over a forested area will be used to demonstrate the results of this derivation. This work was performed by the Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, under contract from the National Aeronautics and Space Administration.

Moghaddam, Mahta↗

NASA NTRS2019

Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy

Solar induced chlorophyll fluorescence has been shown to be increasingly an useful proxy for the estimation of gross primary productivity (GPP), at a range of spatial scales. Here, we explore the seasonality in a continuous time series of canopy solar induced fluorescence (hereafter SiF) and its relation to canopy gross primary production (GPP), canopy light use efficiency (LUE), and direct estimates of leaf level photochemical efficiency in an evergreen canopy. SiF was calculated using infilling in two bands from the incoming and reflected radiance using a pair of Ocean Optics USB2000+ spectrometers operated in a dual field of view mode, sampling at a 30 min time step using custom written automated software, from early spring through until autumn in 2011. The optical system was mounted on a tower of 18 m height adjacent to an eddy covariance system, to observe a boreal forest ecosystem dominated by Scots pine. (Pinus sylvestris) A Walz MONITORING-PAM, multi fluorimeter system, was simultaneously mounted within the canopy adjacent to the footprint sampled by the optical system. Following correction of the SiF data for O2 and structural effects, SiF, SiF yield, LUE, the photochemicsl reflectance index (PRI), and the normalized difference vegetation index (NDVI) exhibited a seasonal pattern that followed GPP sampled by the eddy covariance system. Due to the complexities of solar azimuth and zenith angle (SZA) over the season on the SiF signal, correlations between SiF, SiF yield, GPP, and LUE were assessed on SZA <50° and under strictly clear sky conditions. Correlations found, even under these screened scenarios, resulted around ~r2 = 0.3. The diurnal responses of SiF, SiF yield, PAM estimates of effective quantum yield (ΔF/Fm′), and meteorological parameters demonstrated some agreement over the diurnal cycle. The challenges inherent in SiF retrievals in boreal evergreen ecosystems are discussed.

photosynthetic↗

DOE OSTI2022

Stronger fertilization effects on aboveground versus belowground plant properties across nine U.S. grasslands

Abstract Increased nutrient inputs due to anthropogenic activity are expected to increase primary productivity across terrestrial ecosystems, but changes in allocation aboveground versus belowground with nutrient addition have different implications for soil carbon (C) storage. Thus, given that roots are major contributors to soil C storage, understanding belowground net primary productivity (BNPP) and biomass responses to changes in nutrient availability is essential to predicting carbon–climate feedbacks in the context of interacting global environmental changes. To address this knowledge gap, we tested whether a decade of nitrogen (N) and phosphorus (P) fertilization consistently influenced aboveground and belowground biomass and productivity at nine grassland sites spanning a wide range of climatic and edaphic conditions in the continental United States. Fertilization effects were strong aboveground, with both N and P addition stimulating aboveground biomass at nearly all sites (by 30% and 36%, respectively, on average). P addition consistently increased root production (by 15% on average), whereas other belowground responses to fertilization were more variable, ranging from positive to negative across sites. Site‐specific responses to P were not predicted by the measured covariates. Atmospheric N deposition mediated the effect of N fertilization on root biomass and turnover. Specifically, atmospheric N deposition was positively correlated with root turnover rates, and this relationship was amplified with N addition. Nitrogen addition increased root biomass at sites with low N deposition but decreased it at sites with high N deposition. Overall, these results suggest that the effects of nutrient supply on belowground plant properties are context dependent, particularly with regard to background N supply rates, demonstrating that site conditions must be considered when predicting how grassland ecosystems will respond to increased nutrient loading from anthropogenic activity.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1986

Precise orbit determination of a geosynchronous satellite with differential very-long-baseline interferometry

The use of differential VLBI (DVLBI) to determine the orbit of a U.S. geosynchronous satellite with a position accuracy of a few meters is examined. DVLBI was employed to track the geosynchronous satellite DSCS-II; two observation channels with 2 MHz bandwidth each in a 2 GHz band were used to receive the telemetry signal of the satellite and quasar radio signals. The random observation, delay, station location, and quasar position errors are discussed. The correlation processing of the quasar and satellite signals is described. The accuracy of the estimated satellite position is evaluated using a covariance analysis. It is determined that the satellite position was estimated with an accuracy of about 4 m at the center of the observation arc. It is noted that the modeling error of the solar radiation pressure factor, ionospheric correction error, and observation error need to be corrected in order to provide more accurate orbit determinations.

Shiomi, Tadashi↗

NASA NTRS1992

UV variability and the nature of the continuum source in NGC5548

The short-wave UV continuum of NGC5548 varied quasi-periodically during the IUE monitoring from December 1988 to August 1989. The spectral index of this SWP continuum varied in tight correlation with it, inviting comparison with models to search for clues to the physical nature of the continuum source. The observed covariation of spectral index and continuum are consistent with the model of a black body of constant area (that of a sphere the size of Neptune's orbit) but variable temperature ( 18,000K to 25,000 K). A thin accretion disk around a black hole may also be consistent with this data if the accretion rate can be variable on time scales of 100 days while remaining uniform across the disk. There is marginal evidence for hysteresis in which spectral index changes lead changes in the SWP flux, although the CCF peaks at zero time lag.

Stoner, R. E.↗