Search NASA⌕ Search

SEARCH · Search NASA

Results for “Total variance”

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 217 records · Page 12

The Influence of Observation Errors on Analysis Error and Forecast Skill Investigated with an Observing System Simulation Experiment

The Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework is used to explore the response of analysis error and forecast skill to observation quality. In an OSSE, synthetic observations may be created that have much smaller error than real observations, and precisely quantified error may be applied to these synthetic observations. Three experiments are performed in which synthetic observations with magnitudes of applied observation error that vary from zero to twice the estimated realistic error are ingested into the Goddard Earth Observing System Model (GEOS-5) with Gridpoint Statistical Interpolation (GSI) data assimilation for a one-month period representing July. The analysis increment and observation innovation are strongly impacted by observation error, with much larger variances for increased observation error. The analysis quality is degraded by increased observation error, but the change in root-mean-square error of the analysis state is small relative to the total analysis error. Surprisingly, in the 120 hour forecast increased observation error only yields a slight decline in forecast skill in the extratropics, and no discernable degradation of forecast skill in the tropics.

Observation Error↗

Review of a Draft - Avoiding the Impossible: Re-focusing a Non-Feasible Mission 2-Hrs into a 3-Day Engineering Session

Concurrent engineering offers a great many benefits to engineers and mission designers throughout the world of aerospace. The only downside of concurrent engineering, and this is somewhat unavoidable, is that you don’t know the results of a design session until the end when it is completed. Usually, this is not a problem – you wouldn’t start building a spacecraft before the design is finished. However, within mass and cost constrained systems, you may end up with a final design that although technically sound – is not feasible due to mass or cost limits. Employing in-session mass and cost models with flexible inputs that refine their estimates and variance as more detailed information comes in throughout a design session allows major design changes to be made when the probability of breaching a mass or cost cap exceeds a threshold level. This enables mission designers to re-focus the study, and avoid spending 3-days with 15 engineers designing a non-feasible mission. By understanding key correlations and nested relationships within mass or cost, and specifically mass or cost allocations per mission element by mission type, it’s possible to get flexible-input, statistically based mass and cost estimates very early in the design process. Baseline models are seeded using mission characteristics and general parameters (outer planetary orbiter-probe mission, $500M cost cap for example) to provide a rough estimate of the expected mass or cost. As information gets solidified during the session, it gets added to the model and the estimates are updated. Continuing the orbiter-probe mission example, modeling probe heat shield cost as a percent of total probe cost, and probe cost as a percent of total flight system cost, and total flight system cost as a percent of total mission cost allows a design team to roll-up solidified information to estimate the probability of fitting within a mass or cost constraint early in a concurrent design session. When only the heat shield cost is known, the variance of the final estimate is higher, whereas when the full probe gets defined, naturally, the variance of the estimate decreases. A methodology, model, verification and demo implementation for cost limit breach are presented.

Hihn, Jairus M.↗

Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data

Forest stand structure and biomass data were collected using conventional forest inventory techniques in tropical, subtropical, and warm temperate forest biomes. The feasibility of detecting tropical forest successional age class and total biomass differences using Landsat-Thematic mapper (TM) data, was evaluated. The Normalized Difference Vegetation Index (NDVI) calculated from Landsat-TM data were not significantly correlated with forest regeneration age classes in the mountain terrain of the Luquillo Experimental Forest, Puerto Rico. The low sun angle and shadows cast on steep north and west facing slopes reduced spectral reflectance values recorded by TM orbital altitude. The NDVI, calculated from low altitude aircraft scanner data, was significatly correlated with forest age classes. However, analysis of variance suggested that NDVI differences were not detectable for successional forests older than approximately 15-20 years. Also, biomass differences in young successional tropical forest were not detectable using the NDVI. The vegetation index does not appear to be a good predictor of stand structure variables (e.g., height, diameter of main stem) or total biomass in uneven age, mixed broadleaf forest. Good correlation between the vegetation index and low biomass in even age pine plantations were achieved for a warm temperate study site. The implications of the study for the use of NDVI for forest structure and biomass estimation are discussed.

Sader, Steven A.↗

Evaluating Retention Index Score Assumptions to Refine GC–MS Metabolite Identification

As metabolomics grows into a high-throughput and high demand research field, current metrics for the identification of small molecules in gas chromatography–mass spectrometry (GC–MS) still require manual verification. Though steps have been taken to improve scoring metrics by combining spectral similarity (SS) and retention index (RI), the problem persists. A large body of literature has analyzed and refined SS scores, but few studies have explicitly studied improvements to RI scores. Here, in this study, we examined whether uninvestigated assumptions of the RI score are valid and propose ways to improve them. Query RIs were matched to library RI with a generous window of ±35 to avoid unintentional removal of valid compound identifications. Each match was manually verified as a true positive (TP), true negative, or unknown. Metabolites with at least 30 TP identifications were included in downstream analyses, resulting in a total of 87 metabolites from samples of varying complexity and type (e.g., amino acid mixtures, human urine, fungal species, and so on.). Our results showed that the RI score assumptions of normality, consistent variance across metabolites, and a mean error centered at 0 are often violated. We demonstrated through a cross-validation analysis that modifying these underlying assumptions according to empirical metabolite-specific distributions improved the TP and negative rankings. Further, we statistically determined the minimum number of samples required to estimate distributional parameters for scoring metrics. Overall, this work proposes a robust statistical pipeline to reduce the time bottleneck of metabolite identification by improving RI scores and thus minimize the effort to complete manual verification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Including massive neutrinos in thermal Sunyaev Zeldovich power spectrum and cluster counts analyses

We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the total neutrino mass and data combination. In particular, we find that using the transfer functions for cold dark matter (CDM) + baryons in the computation of the halo mass function, instead of the transfer functions including neutrino perturbations, as prescribed in recent work, yields an ≈0.25 percent downward shift of the σ 8 constraint from tSZ power spectrum data, with a fiducial neutrino mass Σm ν = 0.06 eV. In ΛCDM, with an X-ray mass bias corresponding to the expected hydrostatic mass bias, i.e. (1 - b) ≃ 0.8, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1, and KiDS+VIKING-450. In νΛCDM, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as (1 - b) = 0.64 ± 0.04 (68 percent CL). For forecasts, we find that competitive neutrino mass measurements using cosmic variance limited SZ power spectrum require masking the heaviest clusters and probing the small-scale SZ power spectrum up to ℓ max ≈ 10 4 . Although this is challenging, we find that SZ power spectrum can realistically be used to tightly constrain intracluster medium properties: we forecast a 2 percent determination of the X-ray mass bias by combining CMB-S4 and our mock SZ power spectrum with ℓ max = 10 3 .

Astronomy & Astrophysics↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Crop classification using airborne radar and Landsat data

NASA 13.3 GHz airborne radar data from a soil moisture measurement analysis is used to investigate the statistical nature of the radar backscattering coefficient for bare ground and three different crop types, and to evaluate the crop classification rates using Landsat data alone or combined with the airborne survey. The scatterometer was a fan-beam Doppler system, VV polarized, and is considered only for 50 deg angles of incidence. A total of 36 fields were covered a week apart by the aircraft and Landsat, and Rayleigh statistics were used in the frequency averaging to eliminate fluctuations due to random fluctuations. Within-field variances were calculated for the Landsat and the radar data and used to design optimum crop classification procedures. The Landsat Band 4 readings were 67% accurate, and an increase in accuracy of 10% was achieved by the addition of the radar data.

Ulaby, F. T.↗

A study of satellite-derived moisture with emphasis on the Gulf of Mexico

Visible-Infrared Spin Scan Radiometer (VISSR) Atmospheric Sounder (VAS) moisture retrievals are compared to the National Meteorological Center Regional Analysis and Forecast System (RAFS) 12-h forecast and to 1200 UTC rawinsondes over the U.S. and the Gulf of Mexico on a daily basis for nearly 1.5 years. The principal objective is to determine what information the current moisture retrievals add to that available from the RAFS and surface data. The data are examined from the climatological perspective, that is, total precipitable water over the seasons for three geographical regions, and also for synoptic applications, that is, vertical and horizontal resolution. VAS retrievals are found to be systematically too moist at higher values. The variance of the VAS soundings more closely agrees with the rawinsonde at locations around the Gulf of Mexico than the RAFS. An examination of a case (6 June 1989) over the Gulf of Mexico region comparing three layers of VAS-derived moisture to the RAFS forecast shows the former capable of outperforming the latter in both the horizontal and, to some extent, the vertical frame of reference.

Schreiner, Anthony J.↗

Total-sky and clear-sky variability of outgoing longwave radiation and shortwave radiation and its relation to cloud fluctuation statistics

Correlation statistics are used to investigate some of the characteristic temporal and spatial scales of the day-to-day changes in cloudiness which tend to dominate the variance of outgoing LW radiation and SW radiation fields. The temporal decorrelation time of the radiation field as measured by the lag time required for the autocorrelation of a region to fall to 0.4 is about 12 hours in the midlatitude storm track. Decorrelation time approaches one to two days when the averaging area of the region is increased significantly. Decorrelation times are also longer in the equatorial and tropical regions. The short decorrelation times are a manifestation of the creation, breakup, and movement of cloud systems. Cross-correlation analysis of the radiation field in the midlatitude storm track region and the subtropical dry zone shows advection speeds that range from about 25 m/sec in the storm tracks to about 8 m/sec in the subtropical dry zone.

Bess, T. D.↗

Can Remotely Sensed Snow Disappearance Explain Seasonal Water Supply?

Understanding the relationship between remotely sensed snow disappearance and seasonal water supply may become vital in coming years to supplement limited ground based, in situ measurements of snow in a changing climate. For the period 2001–2019, we investigated the relationship between satellite derived Day of Snow Disappearance (DSD)—the date at which snow has completely disappeared—and the seasonal water supply, i.e., the April—July total streamflow volume, for 15 snow dominated basins across the western U.S. A Monte Carlo framework was applied, using linear regression models to evaluate the predictive skill—defined here as a model’s ability to accurately predict seasonal flow volumes—of varied predictors, including DSD and in situ snow water equivalent (SWE), across a range of spring forecast dates. In all basins there is a statistically significant relationship between mean DSD and seasonal water supply (p ≤ 0.05), with mean DSD explaining roughly half of the variance. Satellite-based model skill improves later in the forecast season, surpassing the skill of in-situ-based (SWE) models in skill in 10 of the 15 basins by the latest forecast date. We found little to no correlation between model error and basin characteristics such as elevation and the ratio of snow water equivalent to total precipitation. Despite a relatively short data record, this exploratory analysis shows promise for improving seasonal water supply prediction, in particular for snow dominated basins lacking in situ observations.

snow remote sensing↗

Allometric and Mobile Terrestrial LiDAR Modeling of Aboveground Woody Biomass of Populus in Coppice Production

Poplars ( Populus spp.) and their hybrids are increasingly being grown in coppice production to generate bioenergy feedstocks at frequent intervals. Allometric equations are re-quired to predict aboveground biomass (AGB) of coppiced individuals with minimal field measurements. Likewise, remote sensing tools like LiDAR (light detection and ranging) can be used if models are available to predict AGB from point cloud data. Therefore, this study sought to develop equations to predict dry woody AGB from field measurements and LiDAR data from coppiced poplar field trials containing eastern cottonwood ( P. del-toides ) and hybrid poplar taxa. We found that taxa-specific allometric models containing the summed basal area of the three largest stems in the coppice provided the best predictive model, with stem height and stem count failing to provide additional explanatory power. The best predictive LiDAR-based model was independent of taxa but had slightly lower adjusted R 2 and higher RMSE than the allometric model. It contained four parameters including crown volume, leaf area index, variance of height returns, and the top point density (i.e., density metric 9 or the proportion of points in the highest point interval when the point cloud is evenly divided into ten vertical intervals). In total, these models can be used to quickly and efficiently estimate dry woody AGB of Populus coppice systems for bioenergy feedstock production.

AGB↗

Variability in daily, zonal mean lower-stratospheric temperatures

Satellite data from the microwave sounding unit (MSU) channel 4, when carefully merged, provide daily zonal anomalies of lower-stratosphere temperature with a level of precision between 0.01 and 0.08 C per 2.5 deg latitude band. Global averages of these daily zonal anomalies reveal the prominent warming events due to volcanic aerosol in 1982 (El Chichon) and 1991 (Mt. Pinatubo), which are on the order of 1 C. The quasibiennial oscillation (QBO) may be extracted from these zonal data by applying a spatial filter between 15 deg N and 15 deg S latitude, which resembles the meridional curvature. Previously published relationships between the QBO and the north polar stratospheric temperatures during northern winter are examined but were not found to be reproduced in the MSU4 data. Sudden stratospheric warmings in the north polar region are represented in the MSU4 data for latitudes poleward of 70 deg N. In the Southern Hemisphere, there appears to be a moderate relationship between total ozone concentration and MSU4 temperatures, though it has been less apparent in 1991 and 1992. In terms of empirical modes of variability, the authors find a strong tendency in EOF 1 (39.2% of the variance) for anomalies in the Northern Hemisphere polar regions to be counterbalanced by anomalies equatorward of 40 deg N and 40 deg S latitudes. In addition, most of the modes revealed significant power in the 15-20 day period band.

Christy, John R.↗

Interpretation of TOMS Observations of Tropical Tropospheric Ozone with a Global Model and In Situ Observations

We interpret the distribution of tropical tropospheric ozone columns (TTOCs) from the Total Ozone Mapping Spectrometer (TOMS) by using a global three-dimensional model of tropospheric chemistry (GEOS-CHEM) and additional information from in situ observations. The GEOS-CHEM TTOCs capture 44% of the variance of monthly mean TOMS TTOCs from the convective cloud differential method (CCD) with no global bias. Major discrepancies are found over northern Africa and south Asia where the TOMS TTOCs do not capture the seasonal enhancements from biomass burning found in the model and in aircraft observations. A characteristic feature of these northern topical enhancements, in contrast to southern tropical enhancements, is that they are driven by the lower troposphere where the sensitivity of TOMS is poor due to Rayleigh scattering. We develop an efficiency correction to the TOMS retrieval algorithm that accounts for the variability of ozone in the lower troposphere. This efficiency correction increases TTOC's over biomass burning regions by 3-5 Dobson units (DU) and decreases them by 2-5 DU over oceanic regions, improving the agreement between CCD TTOCs and in situ observations. Applying the correction to CCD TTOCs reduces by approximately DU the magnitude of the "tropical Atlantic paradox" [Thompson et al, 2000], i.e. the presence of a TTOC enhancement over the southern tropical Atlantic during the northern African biomass burning season in December-February. We reproduce the remainder of the paradox in the model and explain it by the combination of upper tropospheric ozone production from lightning NOx, peristent subsidence over the southern tropical Atlantic as part of the Walker circulation, and cross-equatorial transport of upper tropospheric ozone from northern midlatitudes in the African "westerly duct." These processes in the model can also account for the observed 13-17 DU persistent wave-1 pattern in TTOCs with a maximum above the tropical Atlantic and a minimum over the tropical Pacific during all seasons. The photochemical effects of mineral dust have only a minor role on the modeled distribution of TTOCs, including over northern Africa, due to multiple competing effects. The photochemical effects of mineral dust globally decease annual mean OH concentrations by 9%. A global lightning NOx source of 6 Tg N yr(sup -1) in the model produces a simulation that is most consistent with TOMS and in situ observations.

Martin, Randall V.↗

Resolved gravity waves in the tropical stratosphere: Impact of horizontal resolution and deep convection parametrization

Convectively generated gravity waves (CGGWs) are important for numerical weather prediction due to their effect on the quasi-biennial oscillation (QBO) in the stratosphere. Here, using global ECMWF IFS simulations at TCo7999 (or 1.25 km), TCo2559 (or 3.9 km) and TCo1279 (or 7.8 km) horizontal resolutions, sensitivity of resolved CGGWs to the horizontal resolution and to the explicit versus parametrized representation of deep convection is elucidated during the westerly shear phase of the QBO. Parametrized deep convection is found to inhibit CGGWs, resulting in a twofold reduction in CGGW forcing. When deep convection is explicitly resolved, the total CGGW forcing is almost unchanged across the horizontal resolutions. However, the contribution of long and mesoscale CGGWs (with horizontal wavelengths 100 km ≤ $λ_h$ < $1,900$ km) to the total CGGW forcing decreases and the contribution of smaller-scale CGGWs (with $λ_h$ < $100$ km) increases as the horizontal resolution increases. At the maximum CGGW forcing altitude, at TCo7999 resolution 43% of the total CGGW forcing is due to long and mesoscale waves, whereas at TCo2559 and TCo1279 resolutions their contribution is 70% and 90%, respectively. While CGGW forcing from long and mesoscale waves is similar at TCo7999 resolution with explicit deep convection and at TCo1279 resolution with parametrized deep convection, CGGW forcing from these waves is artificially enhanced at TCo1279 and TCo2559 resolutions with explicit deep convection. This is due to the explicit deep convection being too strong and having too much variance for 100 km ≤ $λ_h$ < $1,900$ km. Therefore, parametrizations of deep convection and CGGWs (to account for forcing from waves with $λ_h$ < $100$ km) are required even at TCo2559 resolution. Additionally, resolved CGGW forcing at TCo7999 resolution is examined for the easterly shear phase of the QBO; similar to the westerly shear phase, the smaller-scale waves contribute 55% to the total CGGW forcing at the maximum CGGW forcing altitude.

54 ENVIRONMENTAL SCIENCES↗

Pulmonary-impedance power spectral analysis: A facile means of detecting radiation-induced gastrointestinal distress and performance decrement in man

Changes in respiratory variance revealed by power spectral analysis of the pulmonary impedance pneumogram can be used to detect and measure stresses directly or indirectly affecting human respiratory function. When gastrointestinal distress occurred during a series of 5 total-body exposures of 30 R at a rate of 1.5 R/min, it was accompanied by typical shifts in pulmonary impedance power spectra. These changes did not occur after protracted exposure of 250 R (30 R daily) at 1.5 R/hr that failed to cause radiation sickness. This system for quantitating respiratory effort can also be used to detect alterations in one's ability to perform under controlled exercise conditions.

Rick, R. C.↗

The geoid spectrum from altimetry

Satellite altimetry information from the world's major oceans was analyzed to arrive at a geoid power spectrum. Using the equivalent of about 7 revolutions of data (mostly from GEOS-3) the power spectrum of the sea surface generally follows the expected values from Kaula's rule applied to the geoid. Analysis of overlapping altimetry arcs (and oceanographic data) shows that the surface spectrum is dominated by the geoid to about 500 cycles (40 km half wavelength) but that sea state departures are significant starting at about 250 cycles (80 km). Estimates of geopotential variances from a derived (smooth) geoid spectrum show significantly less power than Kaula's rule to about 60 cycles, but somewhat more from there to about 400 cycles. At less than 40 km half wavelength, the total power in the marine geoid may be negligible.

Wagner, C. A.↗

Vertical-structure effects on planetary microwave brightness temperature measurements - Applications to the lunar regolith

An attempt to constrain the effects of vertical variations in dielectric properties on lunar microwave observations is presented. A numerical approach for deriving the reflectivity and microwave weighting function of a vertically varying half-space is used, assuming variance in the dielectric properties with depth only, and negligible magnetic effects. The cases of continuous and stratified models of vertical structures are discussed, and a concentration of emitted energy in upper layers is found. The total emitted energy oscillates, varying with the thickness of the upper soil layer, but averaging out interference effects due to random variations in the substrate depth. Consideration is also given to the vertical structure effects on the lunation-mean disk-center brightness temperature, its variations, and the regolith electrical loss, and predicted reflectivity effects by feasible models of the lunar regolith dielectric profile.

Keihm, S. J.↗