Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

New Approaches to Quantifying Transport Model Error in Atmospheric CO2 Simulations

In recent years, much progress has been made in observing CO2 distributions from space. However, the use of these observations to infer source/sink distributions in inversion studies continues to be complicated by difficulty in quantifying atmospheric transport model errors. We will present results from several different experiments designed to quantify different aspects of transport error using the Goddard Earth Observing System, Version 5 (GEOS-5) Atmospheric General Circulation Model (AGCM). In the first set of experiments, an ensemble of simulations is constructed using perturbations to parameters in the model s moist physics and turbulence parameterizations that control sub-grid scale transport of trace gases. Analysis of the ensemble spread and scales of temporal and spatial variability among the simulations allows insight into how parameterized, small-scale transport processes influence simulated CO2 distributions. In the second set of experiments, atmospheric tracers representing model error are constructed using observation minus analysis statistics from NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA). The goal of these simulations is to understand how errors in large scale dynamics are distributed, and how they propagate in space and time, affecting trace gas distributions. These simulations will also be compared to results from NASA's Carbon Monitoring System Flux Pilot Project that quantified the impact of uncertainty in satellite constrained CO2 flux estimates on atmospheric mixing ratios to assess the major factors governing uncertainty in global and regional trace gas distributions.

Ott, L.↗
NASA NTRS2020

Turbulent mixing in supercritical jets: effect of compressibility factor and inflow condition

Fuel injection and turbulent mixing at supercritical pressures determines ignition and combustion in numerous engineering applications. Flow evolution under such conditions is characterized by strong non-linear coupling between dynamics, transport coefficients, and thermodynamics. Experimental studies observe that the jets injected at supercritical pressures exhibit significantly different dynamics from the jets at subcritical conditions, owing to the lack of distinct liquid and gas phases in supercritical state. Thus, the averaged flow quantities such as the potential core length, jet spatial growth rate and velocity decay profiles differ in the two conditions, resulting in different mixed-fluid distributions. In this study, turbulent jet direct numerical simulations (DNS) are performed to examine the variations in flow statistics between injection of Nitrogen (N₂) in Nitrogen (N₂) at both subcritical (perfect-gas) and supercritical conditions. In all cases, isothermal round jets at Reynolds number (Re_{D}), based on jet diameter (D) and jet orifice velocity (U₀), of 5000 are considered. For mixing analyses, a passive scalar transported with the flow is examined.

Sharan, Nek↗
NASA NTRS2021

Assessment of the Advanced Very High Resolution Radiometer (AVHRR) for Snowfall Retrieval in High Latitudes Using CloudSat and Machine Learning

Precipitation retrieval is a challenging topic, especially in high latitudes (HL), and current precipitation products face ample challenges over these regions. This study investigates the potential of the Advanced Very High Resolution Radiometer (AVHRR) for snowfall retrieval in HL using CloudSat radar information and machine learning (ML). With all the known limitations, AVHRR observations should be considered for HL snowfall retrieval because 1) AVHRR data have been continuously collected for about four decades on multiple platforms with global coverage, and similar observations will likely continue in the future; 2) current passive microwave satellite precipitation products have several issues over snow and ice surfaces; and 3) good coincident observations between AVHRR and CloudSat are available for training ML algorithms. Using ML, snowfall rate was retrieved from AVHRR’s brightness temperature and cloud probability, as well as auxiliary information provided by numerical reanalysis. The results indicate that the ML-based retrieval algorithm is capable of detection and estimation of snowfall with comparable or better statistical scores than those obtained from the Atmospheric Infrared Sounder (AIRS) and two passive microwave sensors contributing to the Global Precipitation Measurement (GPM) mission constellation. The outcomes also suggest that AVHRR-based snowfall retrievals are spatially and temporally reasonable and can be considered as a quantitatively useful input to the merged precipitation products that require frequent sampling or long-term records.

Mohammad Reza↗
NASA NTRSDate not supplied

A NASA GISTEMP Observational Uncertainty Ensemble: Regional and Monthly Uncertainty

The historical global temperature record is an essential data product for quantifying the variability and change of the Earth system. In recent years, better characterization of observational uncertainty in global and hemispheric trends has become available, but the methodologies are not necessarily applicable to analyses at smaller regional areas, or monthly means, where station sparsity and other systematic issues contribute to greater uncertainty. This work details a gridded uncertainty ensemble of historical temperature anomalies from the Goddard Institute for Space Studies (GISS) Surface Temperature product (GISTEMP) product. This ensemble characterizes the complex spatial and temporal correlation structure of uncertainty in gridded historical temperature, enabling proper uncertainty propagation for climate and social science at regional and monthly scales. This work details the methodology for generating the uncertainty ensemble, key statistics of the uncertainty evolution over space and time, and provides best practices for using the uncertainty ensemble in future studies. Summary statistics from the uncertainty ensemble are in good agreement with production GISTEMP. Two applications of the uncertainty ensemble are also presented. First, the warmest year on record is shown to most likely be 2016 with a 53.2% chance and 2020 as the second most likely with a 44.4% chance. Second, it is shown that the arctic is warming 2.5 - 5 times faster than the globe, significantly faster than the regularly quoted twice as fast.

GISTEMP↗
NASA NTRS2017

Statistical Approaches for the Definition of Landslide Rainfall Thresholds and their Uncertainty Using Rain Gauge and Satellite Data

Models for forecasting rainfall-induced landslides are mostly based on the identification of empirical rainfall thresholds obtained exploiting rain gauge data. Despite their increased availability, satellite rainfall estimates are scarcely used for this purpose. Satellite data should be useful in ungauged and remote areas, or should provide a significant spatial and temporal reference in gauged areas. In this paper, the analysis of the reliability of rainfall thresholds based on rainfall remote sensed and rain gauge data for the prediction of landslide occurrence is carried out. To date, the estimation of the uncertainty associated with the empirical rainfall thresholds is mostly based on a bootstrap resampling of the rainfall duration and the cumulated event rainfall pairs (D,E) characterizing rainfall events responsible for past failures. This estimation does not consider the measurement uncertainty associated with D and E. In the paper, we propose (i) a new automated procedure to reconstruct ED conditions responsible for the landslide triggering and their uncertainties, and (ii) three new methods to identify rainfall threshold for the possible landslide occurrence, exploiting rain gauge and satellite data. In particular, the proposed methods are based on Least Square (LS), Quantile Regression (QR) and Nonlinear Least Square (NLS) statistical approaches. We applied the new procedure and methods to define empirical rainfall thresholds and their associated uncertainties in the Umbria region (central Italy) using both rain-gauge measurements and satellite estimates. We finally validated the thresholds and tested the effectiveness of the different threshold definition methods with independent landslide information. The NLS method among the others performed better in calculating thresholds in the full range of rainfall durations. We found that the thresholds obtained from satellite data are lower than those obtained from rain gauge measurements. This is in agreement with the literature, where satellite rainfall data underestimate the 'ground' rainfall registered by rain gauges.

landslide prediction↗
NASA NTRS2000

The Impact of ENSO on Extratropical Low Frequency Noise in Seasonal Forecasts

This study examines the uncertainty in forecasts of the January-February-March (JFM) mean extratropical circulation, and how that uncertainty is modulated by the El Nino/Southern Oscillation (ENSO). The analysis is based on ensembles of hindcasts made with an Atmospheric General Circulation Model (AGCM) forced with sea surface temperatures observed during; the 1983 El Nino and 1989 La Nina events. The AGCM produces pronounced interannual differences in the magnitude of the extratropical seasonal mean noise (intra-ensemble variability). The North Pacific, in particular, shows extensive regions where the 1989 seasonal mean noise kinetic energy (SKE), which is dominated by a "PNA-like" spatial structure, is more than twice that of the 1983 forecasts. The larger SKE in 1989 is associated with a larger than normal barotropic conversion of kinetic energy from the mean Pacific jet to the seasonal mean noise. The generation of SKE due to sub-monthly transients also shows substantial interannual differences, though these are much smaller than the differences in the mean flow conversions. An analysis of the Generation of monthly mean noise kinetic energy (NIKE) and its variability suggests that the seasonal mean noise is predominantly a statistical residue of variability resulting from dynamical processes operating on monthly and shorter times scales. A stochastically-forced barotropic model (linearized about the AGCM's 1983 and 1989 base states) is used to further assess the role of the basic state, submonthly transients, and tropical forcing, in modulating the uncertainties in the seasonal AGCM forecasts. When forced globally with spatially-white noise, the linear model generates much larger variance for the 1989 base state, consistent with the AGCM results. The extratropical variability for the 1989 base state is dominanted by a single eigenmode, and is strongly coupled with forcing over tropical western Pacific and the Indian Ocean, again consistent with the AGCM results. Linear calculations that include forcing from the AGCM variance of the tropical forcing and submonthly transients show a small impact on the variability over the Pacific/North American region compared with that of the base state differences.

Schubert, Siegfried D.↗
NASA NTRS1997

A Statistical Study of the Global Structure of the Ring Current

In this paper we derive the average configuration of the ring current as a function of the state of the magnetosphere as indicated by the Dst index. We sort magnetic field data from the Combined Release and Radiation Effects Satellite (CRRES) by spatial location and by the Dst index in order to produce magnetic field maps. From these maps we calculate local current systems by taking the curl of the magnetic field. NN7e find both the westward (outer) and the eastward (inner) components of the ring current. We find that the ring current intensity varies linearly with D.St as expected, and that the ring current is asymmetric for all Dst values. The azimuthal peak of the ring current is located in the afternoon sector for quiet conditions, and near midnight for disturbed conditions. The ring current also moves closer to the Earth during disturbed conditions. We are able to recreate the Dst index by integrating the magnetic perturbations caused by the ring current. We find that we needed to apply a 20 nT offset to Dst, and assume a perfectly conducting Earth to obtain an optimal agreement between the computed and the observed Dst. We interpret the 20 nT offset as the magnetic field generated by the quiet time ring current used as baseline in computing Dst.

Jorgensen, A. M.↗
NASA NTRSDate not supplied

Australian Bushfire in 2020: Accessing MERRA-2 Data in GES DISC Remotely through OPeNDAP & Calculating Statistics with Python3

The evolution and transport of thick haze from the 2020 Australian bushfire is tracked using the NASA Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data archived at the NASA GES DISC data center. The MERRA-2 provides global data at 0.5 x 0.625 spatial resolution since the year 1980. In this use case, with xarray, a python3 library, we remotely accessed the hourly aerosol optical depth (AOD) and PM2.5 data through the OpeNDAP service provided by GES DISC. We also derived weekly data from hourly ones.

Xiaohua Pan↗
NASA NTRS1972

A review of some aspects of inferring the ozone profile by inversion of ultraviolet radiance measurements

The mathematical inversion of light-scattering observations to obtain the atmospheric ozone profile is discussed in terms of the filtering properties of the physical and mathematical processes for different spatial scales. Within this context, it is shown that the physical process of scattering acts as a low pass filter; which transfers large-scale profile information efficiently to radiance observations but very strongly attenuates small- or fine-scale profile information. To avoid domination of the mathematical inversion by the random error of radiance measurements, an equivalent spatial-scale filtering in the inversion procedure is essential. The available evidence suggests that mathematical inversion to obtain the low level ozone profile below 25 to 30 km is either inferior to or no better than the statistical estimation of the profile using total ozone as predictor. However, inversion profiles for high level ozone above 25 to 30 km appear to have moderately good accuracy.

Mateer, C. L.↗
NASA NTRS1972

An experimental investigation of the power spectrum of phase modulation induced on a satellite radio signal by the ionosphere

The power spectrum of phase modulation imposed upon satellite radio signals by the inhomogeneous F-region of the ionosphere (100 - 500 km) was studied. Tapes of the S-66 Beacon B Satellite recorded during the period 1964 - 1966 were processed to yield or record the frequency of modulation induced on the signals by ionospheric dispersion. This modulation is produced from the sweeping across the receiving station as the satellite transits of the two dimensional spatial phase pattern are produced on the ground. From this a power spectrum of structure sizes comprising the diffracting mechanism was determined using digital techniques. Fresnel oscillations were observed and analyzed along with some comments on the statistical stationarity of the shape of the power spectrum observed.

Moser, D. T.↗
NASA NTRS1973

Comparative evaluation of spatial features in automatic land use classification from photographic imagery

Five spatial features have been evaluated for their applicability in automatic land use classification of photographic imagery. Data arrays of approximately 10,000 square meters were classified on the basis of an 8 by 8 point data array. Statistical features, information features, sequency features and texture features were evaluated using a distance to prototype classifier and an adaptive classifier. The results indicate approximately 70% accuracy in the classification.

Herzog, J. H.↗
NASA NTRS1985

Influence of spatial variability of hydraulic characteristics of soils on surface parameters obtained from remote sensing data in infrared and microwaves

The correct interpretation of thermal and hydraulic soil parameters infrared from remotely sensed data (thermal infrared, microwaves) implies a good understanding of the causes of their temporal and spatial variability. Given this necessity, the sensitivity of the surface variables (temperature, moisture) to the spatial variability of hydraulic soil properties is tested with a numerical model of heat and mass transfer between bare soil and atmosphere. The spatial variability of hydraulic soil properties is taken into account in terms of the scaling factor. For a given soil, the knowledge of its frequency distribution allows a stochastic use of the model. The results are treated statistically, and the part of the variability of soil surface parameters due to that of soil hydraulic properties is evaluated quantitatively.

Brunet, Y.↗
NASA NTRS1986

Utilizing remote sensing of Thematic Mapper data to improve our understanding of estuarine processes and their influence on the productivity of estuarine-dependent fisheries

LANDSAT thematic mapper (TM) data are being used to refine and validate a stochastic spatial computer model to be applied to coastal resource management problems in Louisiana. Two major aspects of the research are: (1) the measurement of area of land (or emergent vegetation) and water and the length of the interface between land and water in TM imagery of selected coastal wetlands (sample marshes); and (2) the comparison of spatial patterns of land and water in the sample marshes of the imagery to that in marshes simulated by a computer model. In addition to activities in these two areas, the potential use of a published autocorrelation statistic is analyzed.

Browder, J. A.↗
NASA NTRS1993

Cloud properties from the analysis of AVHRR observations for FIRE 2

Preliminary results are presented for cloud properties from the analysis of AVHRR observations for FIRE 2. The properties were obtained from a combination of the spatial coherence method and a multispectral retrieval scheme. Geographically gritted fields for the number of cloud layers were produced. For single layered cloud systems, fractional cloud cover, cloud emission temperature, cloud emissivity, and particle size were retrieved. Statistics on the properties of upper-level clouds and the Coffeeville cloud conditions are presented.

Lin, Xijian↗
NASA NTRS1993

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗
NASA NTRS2007

Diurnal Lightning Distributions as Observed by the Optical Transient Detector (OTD) and the Lightning Imaging Sensor (LIS)

Data obtained from the Optical Transient Detector (April 1995 to March 2000) and the Lightning Imaging Sensor (December 1997 to December 2005) satellites (70 and 35 inclination low earth orbits, respectively) are used to statistically determine the number of flashes in the annual and seasonal diurnal cycle as a function of local and universal time. The data are further subdivided by season, land versus ocean, northern versus southern hemisphere, and other spatial (e.g., continents) and temporal (e.g., time of peak diurnal amplitude) categories. The data include corrections for detection efficiency and instrument view time. Continental results display strong diurnal variation, with a lightning peak in the late afternoon and a minimum in late morning. In regions of the world dominated by large mesoscale convective systems the peak in the diurnal curve shifts toward late evening or early morning hours. The maximum diurnal flash rate occurs in June-August, corresponding to the Northern Hemisphere summer, while the minimum occurs in December-February. Summer lightning dominates over winter activity and springtime lightning dominates over autumn activity at most continental locations. This latter behavior occurs especially strongly over the Amazon region in South America in September-November. Oceanic lightning activity in winter and autumn tends to exceed that in summer and spring. Global lightning is well correlated in phase but not in amplitude with the Carnegie curve. The diurnal flash rate varies about 4-35 percent about the mean, while the Carnegie curve varies around 4-15 percent.

Bailey, Jeff C.↗
NASA NTRS2026

CH3OH and HCN in Interstellar Comet 3I/ATLAS Mapped with the ALMA Atacama Compact Array: Distinct Outgassing Behaviors and a Remarkably High CH3OH/HCN Production Rate Ratio

We report the detection of methanol (CH3OH) toward interstellar comet 3I/ATLAS using the At-acama Compact Array of the Atacama Large Millimeter/Submillimeter Array (ALMA) on UT 2025 August 28, September 18 and 22, and October 1, and of hydrogen cyanide (HCN) on September 12 and 15. These observations spanned pre-perihelion heliocentric distances (r H ) of 2.6 – 1.7 au. The molecules showed outgassing patterns distinct from one another, with HCN production being depleted in the sunward hemisphere of the coma, whereas CH3OH was enhanced in that direction. Statistical analysis of molecular scale lengths in 3I/ATLAS indicated that CH3OH included production from coma sources at Lp > 258 km at 99% confidence. However low signal-to-noise on long baselines, which sample emission on small spatial scales closest to the nucleus, prevented definitively ruling out CH3OH as purely a parent species. In contrast, HCN production was indistinguishable from direct nucleus sublimation. The CH3OH production rate increased sharply from August through October, including an uptick near the inner edge of the H2O sublimation zone at rH = 2 au. Compared to comets measured to date at radio wavelengths, the derived CH3OH/HCN ratios in 3I/ATLAS of 124+30−34 and 79+11−14on September 12 and 15, respectively, are among the most enriched values measured in any comet, surpassed only by anomalous solar system comet C/2016 R2 (PanSTARRS)

Comae↗
NASA NTRS2017

On the Impact of Granularity of Space-Based Urban CO2 Emissions in Urban Atmospheric Inversions: A Case Study for Indianapolis, IN

Quantifying greenhouse gas (GHG) emissions from cities is a key challenge towards effective emissions management. An inversion analysis from the INdianapolis FLUX experiment (INFLUX) project, as the first of its kind, has achieved a top-down emission estimate for a single city using CO2 data collected by the dense tower network deployed across the city. However, city-level emission data, used as a priori emissions, are also a key component in the atmospheric inversion framework. Currently, fine-grained emission inventories (EIs) able to resolve GHG city emissions at high spatial resolution, are only available for few major cities across the globe. Following the INFLUX inversion case with a global 1x1 km ODIAC fossil fuel CO2 emission dataset, we further improved the ODIAC emission field and examined its utility as a prior for the city scale inversion. We disaggregated the 1x1 km ODIAC non-point source emissions using geospatial datasets such as the global road network data and satellite-data driven surface imperviousness data to a 3030 m resolution. We assessed the impact of the improved emission field on the inversion result, relative to priors in previous studies (Hestia and ODIAC). The posterior total emission estimate (5.1 MtC/yr) remains statistically similar to the previous estimate with ODIAC (5.3 MtC/yr). However, the distribution of the flux corrections was very close to those of Hestia inversion and the model-observation mismatches were significantly reduced both in forward and inverse runs, even without hourly temporal changes in emissions. EIs reported by cities often do not have estimates of spatial extents. Thus, emission disaggregation is a required step when verifying those reported emissions using atmospheric models. Our approach offers gridded emission estimates for global cities that could serves as a prior for inversion, even without locally reported EIs in a systematic way to support city-level Measuring, Reporting and Verification (MRV) practice implementation.

EIs↗