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 451 records · Page 25

Molecules and chains in a strong magnetic field - Statistical treatment

A Thomas-Fermi-Dirac-Weizsaecker statistical model is developed and employed to investigate diatomic molecules and infinite molecular chains in strong magnetic fields. The standard magnetic Thomas-Fermi-Dirac kinetic, potential, and exchange energy functionals are supplemented by a gradient correction to the kinetic energy. The numerical method used for solving this system in two spatial dimensions is detailed. Numerical solutions for a wide range of magnetic strengths and elements are presented to demonstrate the robustness, as well as the limitations, of the statistical approach. These calculations qualitatively reproduce many of the results of detailed quantum mechanical treatments. For example, the fractional binding energy is greatest for low atomic numbers and for strong magnetic fields.

Abrahams, Andrew M.↗

Statistical Estimation of Orbital Debris Populations with a Spectrum of Object Size

Orbital debris is a real concern for the safe operations of satellites. In general, the hazard of debris impact is a function of the size and spatial distributions of the debris populations. To describe and characterize the debris environment as reliably as possible, the current NASA Orbital Debris Engineering Model (ORDEM2000) is being upgraded to a new version based on new and better quality data. The data-driven ORDEM model covers a wide range of object sizes from 10 microns to greater than 1 meter. This paper reviews the statistical process for the estimation of the debris populations in the new ORDEM upgrade, and discusses the representation of large-size (greater than or equal to 1 m and greater than or equal to 10 cm) populations by SSN catalog objects and the validation of the statistical approach. Also, it presents results for the populations with sizes of greater than or equal to 3.3 cm, greater than or equal to 1 cm, greater than or equal to 100 micrometers, and greater than or equal to 10 micrometers. The orbital debris populations used in the new version of ORDEM are inferred from data based upon appropriate reference (or benchmark) populations instead of the binning of the multi-dimensional orbital-element space. This paper describes all of the major steps used in the population-inference procedure for each size-range. Detailed discussions on data analysis, parameter definition, the correlation between parameters and data, and uncertainty assessment are included.

Xu, Y. -l↗

Spaceborne Lidar Retrievals of PM2.5 for Air Quality Studies and Applications

Fine particulate matter (PM2.5) substantially contributes to air pollution and negatively affects human health. While many studies have investigated the use of passive column-integrated aerosol optical depth to infer surface PM2.5, the use of lidar observations for air quality characterization is not nearly as extensive. Lidar measurements are critical, however, due to the vertical aerosol information they provide, including near the surface. In this presentation, we first provide an overview of various lidar-based approaches for estimating PM2.5 concentrations and then discuss how lidar measurements can assist other air quality applications. For example, estimates of PM2.5 have been obtained in a physics-based approach through CALIOP near-surface aerosol extinction retrievals, assumptions on the mass extinction efficiency, and incorporating other parameters (an aerosol hygroscopic growth factor and PM2.5/PM10 ratio). Application of this algorithm over the contiguous United States (CONUS) from 2006 to 2018 yielded larger PM2.5 values over the eastern and western CONUS (~10-15 μg/m³) and lower PM2.5 levels in the central CONUS (~5 μg/m³). These spatial patterns were similar to those from gridded PM2.5 concentrations obtained through in situ measurements at ground stations operated by the US Environmental Protection Agency. In another approach, the Cloud Aerosol Transport System (CATS) lidar was used with the Goddard Earth Observing System (GEOS) model in a 1D ensemble-based variational technique to obtain PM2.5 over the US and Europe, and the spatial patterns of the CATS/GEOS based PM2.5 concentrations generally captured those from surface stations (with corresponding hourly EPA PM2.5 vs CATS PM2.5 statistics of R=0.4 and bias=1.5 μg/m³). In our recent work, as part of the Models, In situ, and Remote sensing of Aerosols (MIRA) Working Group, we have applied both the CALIOP and CATS/GEOS based approaches over the highly polluted country of India during the post-monsoon season (September-October 2016). We derived elevated levels of two-month mean PM2.5 (~100 μg/m³) in northern India, especially near New Delhi. These high PM2.5 concentrations in the Indo-Gangetic plain are driven in large part from the seasonal burning of crop residue and meteorological conditions typical at this time of the year, such as low wind speeds and a shallow boundary layer. While the satellite-derived PM2.5 moderately replicates (R = ~0.7-0.9) the spatial variability in the two-month mean of surface in situ PM2.5 from monitoring sites operated by the Central and State Pollution Control Boards, we show results from specific scenes for which there are large deviations between the satellite-derived PM2.5 and in situ measurements. Other current work on this topic focuses on developing PM2.5 estimates using airborne high spectral resolution lidar measurements through machine learning regression algorithms and involves several parameters (e.g., aerosol extinction, color ratio, lidar ratio). Application of this method over major metropolitan areas in the US and Asia have resulted in high correlations (R = 0.93) with surface measurements. This airborne lidar approach can be adapted to spaceborne lidar measurements, and all three of these approaches can be applied to ESA’s EarthCARE Atmospheric Lidar instrument, setting the stage for the future Cloud Aerosol Lidar for Global Scale Observations of the Ocean-Land Atmosphere System (CALIGOLA) mission. Ultimately, beyond estimates of PM2.5, the aerosol vertical distribution from lidars can benefit studies involving passive sensor approaches for PM2.5 proxies (including from geostationary satellites), wildfire smoke plume injection heights, volcanic emissions (e.g., ash height retrievals), and aerosol/air quality model assimilation, evaluation, and forecasts.

Travis D Toth↗

The effects of longitudinal roughness elements upon the turbulent boundary layer

Velocity profiles and bursting statistics were measured in the presence of cylindrical longitudinal roughness elements (LREs) and compared to similar work by Johansen and Smith (1986). An algorithm was devised to detect low speed streaks (LSS) from the hot-wire rake data, allowing an estimation of the LREs' ability to reduce spatial randomness of the sublayer structure. The velocity profiles agree with Johansen and Smith's result that the effects of the LREs are felt only for y(+) less than about 15. The analysis of the bursting statistics and LSS algorithm output shows that the flow relaxes back to a flat plate boundary layer between the LREs. A crude first approximation is that the LREs simply impose the no-slip boundary condition at an elevation equivalent to their diameter. While the LREs might serve as a nucleation site for the LSS, spatial randomness of the structure is affected only very close to the LREs.

Blackwelder, R. F.↗

A stochastic-dynamic model for global atmospheric mass-field statistics

Global atmospheric mass field error correlations based on satellite observations and on numerical forecasts show strong and systematic latitude dependence. A model for the latitude dependent spatial correlation structure of mass field forecast errors is derived from dynamical considerations. Three methods of solution were tested. In the first method, the equation was solved by expansion in spherical harmonics, and the correlation function was computed analytically using the expansion coefficients. In the second method, the finite difference equivalent of the equation was solved using a fast poisson solver. The correlation function was computed using stratified sampling of the individual realizations. In the third method, a higher order equation was derived, and solved directly in finite differences by two successive applications of the fast poisson solver. The three methods were compared for accuracy and efficiency, and the third method was chosen as clearly superior.

Ghil, M.↗

Recognition and characterization of hierarchical interstellar structure. II - Structure tree statistics

A new method of image analysis is described, in which images partitioned into 'clouds' are represented by simplified skeleton images, called structure trees, that preserve the spatial relations of the component clouds while disregarding information concerning their sizes and shapes. The method can be used to discriminate between images of projected hierarchical (multiply nested) and random three-dimensional simulated collections of clouds constructed on the basis of observed interstellar properties, and even intermediate systems formed by combining random and hierarchical simulations. For a given structure type, the method can distinguish between different subclasses of models with different parameters and reliably estimate their hierarchical parameters: average number of children per parent, scale reduction factor per level of hierarchy, density contrast, and number of resolved levels. An application to a column density image of the Taurus complex constructed from IRAS data is given. Moderately strong evidence for a hierarchical structural component is found, and parameters of the hierarchy, as well as the average volume filling factor and mass efficiency of fragmentation per level of hierarchy, are estimated. The existence of nested structure contradicts models in which large molecular clouds are supposed to fragment, in a single stage, into roughly stellar-mass cores.

Houlahan, Padraig↗

Estimating noise levels of remotely sensed measurements from satellites using spatial structure analysis

A technique is presented whereby the noise level of satellite measurements of the atmosphere and earth can be estimated. The technique analyzes a spatial array of data measured by a satellite instrument. A minimum of about 200 satellite measurements is required, preferably in a regular pattern. Statistical structure analysis is used to describe a combination of the mean gradient and noise in the data. The noise level is then estimated by separating out the gradient information and leaving only the noise. Results are presented for four satellite sounding instruments, and effective blackbody or brightness temperature noise levels were compared to prelaunch specifications or inflight calibrations for each instrument. Comparisons showed that in the absence of cloud-contaminated measurements (in the case of infrared data) and away from the highly variable ground surface, the noise level of various satellite instruments can be obtained without the need for calibration data. The noise levels imply how much spatial averaging is possible, without smearing the detected geophysical gradient, and how much is necessary, to meet the absolute signal accuracy requirements for the intended use of the satellite measurements.

Hillger, Donald W.↗

On the structure of turbulent wakes and merging shear layers of multielement airfoils

The structure of the attached turbulent flow in the vicinity of a NACA 4412 airfoil equipped with a single-slotted flap was studied. The airfoil/flap configuration was tested at a Mach number of 0.06 and a Reynolds number of 1.3 x 10 to the 6th in a 7- by 10-Foot Wind Tunnel. Surface-pressure measurements were made on the main airfoil and on the flap. Detailed measurements, obtained using a high-spatial-resolution laser Doppler anemometer, were made of the mean velocity flow field and of the second-order statistical quantities (Reynolds stresses) in the boundary layers, wakes, and merging shear layers. The experimental observations are compared with theoretical predictions of pressure, mean velocity, and Reynolds stress.

Olson, L. E.↗

Use of MAGSAT anomaly data for crustal structure and mineral resources in the US Midcontinent

The analysis and preliminary interpretation of investigator-B MAGSAT data are addressed. The data processing included: (1) removal of spurious data points; (2) statistical smoothing along individual data tracks, to reduce the effect of geomagnetic transient disturbances; (3) comparison of data profiles spatially coincident in track location but acquired at different times; (4) reduction of data by weighted averaging to a grid with 1 deg xl deg latitude/longitude spacing, and with elevations interpolated and weighted to a common datum of 400 km; (5) wavelength filtering; and (6) reduction of the anomaly map to the magnetic pole. Agreement was found between a magnitude data anomaly map and a reduce-to-the-pole map supporting the general assumption that, on a large scale (long wavelength), it is induced crustal magnetization which is responsible for major anamalies. Anomalous features are identified and explanations are suggested with regard to crustal structure, petrologic characteristics, and Curie temperature isotherms.

Carmichael, R. S.↗

Galactic neutron stars and gamma-ray bursts

The association of the gamma-ray burst phenomenon with Galactic neutron stars is investigated statistically, and the distance to presently observable gamma-ray bursts is constrained. This is done by calculating the spatial distribution and kinematic properties of a sample of Population I neutron stars and comparing their properties to those of observed gamma-ray bursts. Both brightness distribution and angular distribution on the celestial sphere are employed. Current observations suggest negligible correlation on small angular scales, large-scale isotropy, and a radial distribution that is approximately uniform to the distances currently observed. It is concluded that the distribution of Population I neutron stars is consistent with current observations if gamma-ray bursts occur continuously throughout the neutron star lifetime and if the current sampling depth is at least about 150 pc but not more than roughly 2 kpc.

Hartmann, Dieter↗

Structural properties of compact groups

We report the results of a systematic study of galaxies in the regions of Hickson compact groups. Our sample is composed of the 22 Hickson groups which are located in the southern hemisphere and have cz less than 9000 km/s. Making use of digitized images of IIIa-J plates that cover an area of 0.5 x 0.5 deg around each group, we were able to detect and classify images down to a magnitude limit of 19.5 in the B band. This limit is typically three magnitudes fainter than previous studies. Most groups show a statistically significant excess of fainter galaxies compared to the background. These fainter galaxies typically have a somewhat more extended spatial distribution than the brighter galaxies originally classified by Hickson. Our data suggest that Hickson groups have a wide range in density and radius, ranging from very compact structures with overdensities of the order of 10(exp 2) and crossing times of roughly 0.01 H(sub 0 sup -1), to much more diffuse structures, similar to loose groups, with overdensities of about 3 and crossing times of roughly 0.5 H(sub 0 sup -1).

De Carvalho, R. R.↗

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗

Evaluation of Terms in the Water Vapor Budget Using Airborne Dial and In Situ Measurements from the Southern Great Plans 1997 Experiment

The Southern Great Plains (SGP97) field experiment was conducted in Oklahoma during June and July 1997 primarily to validate soil moisture retrieval algorithms using microwave radiometer measurements from aircraft as well as in situ surface measurements. One important objective of the SGP97 experiment plan was to examine the effect of soil moisture on the evolution of the atmospheric boundary layer (ABL) and clouds over the Southern Great Plains during the warm season. To support boundary layer studies during SGP97. the NASA Langley Research Center's Lidar Atmospheric Sensing Experiment (LASE) was flown on a NASA-P3 aircraft in conjunction with the Electronically Scanned Thinned Array Radiometer (ESTAR). The LASE instrument is an airborne, downward-looking differential absorption lidar (DIAL) system capable of measuring water vapor concentration as well as aerosol backscatter with high horizontal and vertical resolution in the ABL. Here, we will demonstrate how the LASE data can be used to determine water vapor statistics and most of the water vapor budget terms in the ABL. This information can then be related to spatial variations in soil moisture and the surface energy budget. The extensive surface and aircraft in situ measurements conducted during SGP97 provide information on the ABL that cannot be retrieved from the LASE data alone and also offer an excellent opportunity to validate the remote water vapor budget measurements with LASE.

Senff, Christoph J.↗

Assessing Satellite-Based Fire Data for use in the National Emissions Inventory

Biomass burning is significant to emission estimates because: (1) it can be a major contributor of particulate matter and other pollutants; (2) it is one of the most poorly documented of all sources; (3) it can adversely affect human health; and (4) it has been identified as a significant contributor to climate change through feedbacks with the radiation budget. Additionally, biomass burning can be a significant contributor to a regions inability to achieve the National Ambient Air Quality Standards for PM 2.5 and ozone, particularly on the top 20% worst air quality days. The United States does not have a standard methodology to track fire occurrence or area burned, which are essential components to estimating fire emissions. Satellite imagery is available almost instantaneously and has great potential to enhance emission estimates and their timeliness. This investigation compares satellite-derived fire data to ground-based data to assign statistical error and helps provide confidence in these data. The largest fires are identified by all satellites and their spatial domain is accurately sensed. MODIS provides enhanced spatial and temporal information, and GOES ABBA data are able to capture more small agricultural fires. A methodology is presented that combines these satellite data in Near-Real-Time to produce a product that captures 81 to 92% of the total area burned by wildfire, prescribed, agricultural and rangeland burning. Each satellite possesses distinct temporal and spatial capabilities that permit the detection of unique fires that could be omitted if using data from only one satellite.

Soja, Amber J.↗

The Atacama Cosmology Telescope: Cosmology from Galaxy Clusters Detected Via the Sunyaev-Zel'dovich Effect

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sample used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cosmology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives (sigma)8 = 0.851 +/- 0.115 and w = -1.14 +/- 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation obtained from numerical simulations and calibrated by X-ray observations, we find (sigma)8 + 0.821 +/- 0.044 and w = -1.05 +/- 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations plus type Ia supernova which give (sigma)8 = 0.802 +/- 0.038 and w = -0.98 +/- 0.053. A stacking analysis of the clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and the cosmological parameters derived from them are broadly consistent with current models.

Sehgal, Neelima↗

Estimating Uncertainties in the Multi-Instrument SBUV Profile Ozone Merged Data Set

The MOD data set is uniquely qualified for use in long-term ozone analysis because of its long record, high spatial coverage, and consistent instrument design and algorithm. The estimated MOD uncertainty term significantly increases the uncertainty over the statistical error alone. Trends in the post-2000 period are generally positive in the upper stratosphere, but only significant at 1-1.6 hPa. Remaining uncertainties not yet included in the Monte Carlo model are Smoothing Error ( 1 from 10 to 1 hPa) Relative calibration uncertainty between N11 and N17Seasonal cycle differences between SBUV records.

Uncertainties↗

Unifying Radiative Transfer Models in Computer Graphics and Remote Sensing, Part II: A Differentiable, Polarimetric Forward Model and Validation

The constellation of Earth-observing satellites continuously collects measurements of scattered radiance, which must be transformed into geophysical parameters in order to answer fundamental scientific questions about the Earth. Retrieval of these parameters requires highly flexible, accurate, and fast forward and inverse radiative transfer models. Existing forward models used by the remote sensing community are typically accurate and fast, but sacrifice flexibility by assuming the atmosphere or ocean is composed of plane-parallel layers. Monte Carlo forward models can handle more complex scenarios such as 3D spatial heterogeneity, but are relatively slower. We propose looking to the computer graphics community for inspiration to improve the statistical efficiency of Monte Carlo forward models and explore new approaches to inverse models for remote sensing. In Part 2 of this work, we demonstrate that Monte Carlo forward models in computer graphics are capable of sufficient accuracy for remote sensing by extending Mitsuba 3, a forward and inverse modeling framework recently developed in the computer graphics community, to simulate simple atmosphere-ocean systems and show that our framework is capable of achieving error on par with codes currently used by the remote sensing community on benchmark results.

Radiative transfer↗

Unifying radiative transfer models in computer graphics and remote sensing, Part I: A Survey

The constellation of Earth-observing satellites continuously collects measurements of scattered radiance, which must be transformed into geophysical parameters in order to answer fundamental scientific questions about the Earth. Retrieval of these parameters requires highly flexible, accurate, and fast forward and inverse radiative transfer models. Existing forward models used by the remote sensing community are typically accurate and fast, but sacrifice flexibility by assuming the atmosphere or ocean is composed of plane-parallel layers. Monte Carlo forward models can handle more complex scenarios such as 3D spatial heterogeneity, but are relatively slower. We propose looking to the computer graphics community for inspiration to improve the statistical efficiency of Monte Carlo forward models and explore new approaches to inverse models for remote sensing. In Part 1 of this work, we examine the evolution of radiative transfer models in computer graphics and highlight recent advancements that have the potential to push forward models in remote sensing beyond their current periphery of realism.

Radiative transfer↗