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 937 records · Page 52

Unsupervised Classification of Global Radar Units on Venus

Characterization of the Venusian surface in terms of its radar properties was accomplished by application of an unsupervised, linear discriminant algorithm to two Pioneer-Venus (PV) Orbiter radar data sets: the RMS-slope (surface roughness) and reflectivity. Both databases were spatially filtered to the same effective resolution of 100 km prior to classification. A recent supervised classification study using these data was based on presupposed morphologic significance of selected data ranges. The knowledge of both Venusian geology and the geologic significance of the radar data is so limited that the data warrant a more unsupervised approach; for this study a linear discriminant classifier was chosen. This approach is purely statistical, thereby removing any observer bias. Statistical significance of the resulting clusters was evaluated by an ancillary program in which an F test utilizing the Mahalanobis' distance.

Kozak, R. C.↗

A scanning defect mapping system for semiconductor characterization

We have developed an optical scanning system that generates maps of the spatial distributions of defects in single and polycrystalline silicon wafers. This instrument, called Scanning Defect Mapping System, utilizes differences in the scattering characteristics of dislocation etch pits and grain boundaries from a defect-etched sample to identify and count them. This system simultaneously operates in the dislocation mode and the grain boundary (GB) mode. In the 'dislocation mode,' the optical scattering from the etch pits is used to statistically count dislocations, while ignoring the GB's. Likewise, in the 'grain boundary mode' the system only recognizes the local scattering from the GB's to generate grain boundary distributions. The information generated by this instrument is valuable for material quality control, identifying mechanisms of defect generation and the nature of thermal stresses during the crystal growth, and the solar cell process design.

Sopori, Bushnan L.↗

On the Lighthill relationship and sound generation from isotropic turbulence

In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.

Zhou, YE↗

Strong gravitational lensing statistics as a test of cosmogonic scenarios

Gravitational lensing statistics can provide a direct and powerful test of cosmic structure formation theories. Since lensing tests, directly, the magnitude of the nonlinear mass density fluctuations on lines of sight to distant objects, no issues of 'bias' (of mass fluctuations with respect to galaxy density fluctuations) exist here, although lensing observations provide their own ambiguities of interpretation. We develop numerical techniques for generating model density distributions with the very large spatial dynamic range required by lensing considerations and for identifying regions of the simulations capable of multiple image lensing in a conservative and computationally efficient way that should be accurate for splittings significantly larger than 3 seconds. Applying these techniques to existing standard Cold dark matter (CDM) (Omega = 1) and Primeval Baryon Isocurvature (PBI) (Omega = 0.2) simulations (normalized to the Cosmic Background Explorer Satellite (COBE) amplitude), we find that the CDM model predicts large splitting (greater than 8 seconds) lensing events roughly an order-of-magnitude more frequently than the PBI model. Under the reasonable but idealized assumption that lensing structrues can be modeled as singular isothermal spheres (SIS), the predictions can be directly compared to observations of lensing events in quasar samples. Several large splitting (Delta Theta is greater than 8 seconds) cases are predicted in the standard CDM model (the exact number being dependent on the treatment of amplification bias), whereas none is observed. In a formal sense, the comparison excludes the CDM model at high confidence (essentially for the same reason that CDM predicts excessive small-scale cosmic velocity dispersions.) A very rough assessment of low-density but flat CDM model (Omega = 0.3, Lambda/3H(sup 2 sub 0) = 0.7) indicates a far lower and probably acceptable level of lensing. The PBI model is consistent with, but not strongly tested by, the available lensing data, and other open models would presumably do as well as PBI. These preliminary conclusions and the assumptions on which they are based can be tested and the analysis can be applied to other cosmogonic models by straightforward extension of the work presented here.

Cen, Renyue↗

The Effect of Strain Rate on the Evolution of Plane Wakes Subjected to Irrotational Strains

Direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of irrotational plane strain applied at three different strain rates have been generated. The strain geometry is such that the flow is compressed in the streamwise direction and expanded in the cross-stream direction with the spanwise direction being unstrained. This geometry is the temporally evolving analogue of a spatially evolving wake in an adverse pressure gradient. A pseudospectral numerical method with up to 16 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2,000. Although the evolutions of many statistics are nearly collapsed when plotted against total strain, there are some differences owing to the different strain rate histories. The impact of strain-rate on the wake spreading rate, the peak velocity deficit, the Reynolds stress profiles, and the flow structure is examined.

Rogers, Michael M.↗

[Overview of the GSFC Effort to Build Microwave Detectors for CMB Polarization]

Filled arrays of bolometers are currently being employed for use in astronomy from the far-infrared through millimeter portion of the electromagnetic spectrum. Because of the large range of wavelengths for which such detectors are applicable, the size of a pixel (p) relative to the wavelength (lambda) will vary according to the specific application of a given available technology. As the pixel size becomes comparable to a wavelength of the absorbed radiation, correlations in the radiation field at the pixel can affect the imaging properties of the system independent of the correlations introduced by the front-end optics. We study the dependence of image fidelity and induced polarization on the size of the pixel by employing a formalism in which diffraction due to the pixel boundary is treated by propagating the second-order statistical correlations of the radiation field through a model optical system. We construct simulated polarized images of square pixels for various ratios of p lambda. For the limit in which few modes are supported by the pixel (p lambda lesssiml), we find that the diffraction due to the pixel edges is non-negligible and hence must be considered along with the telescope diffraction pattern in modeling the ultimate spatial resolution of an imaging system. For the case in which the pixel is over-moded (p lambda ggl), the geometric limit is approached as expected. This technique provides a quantitative approach to understand and optimize the imaging properties of planar detectors in the few-mode limit.

Chuss, David↗

An Initial Assessment of a SMAP Soil Moisture Disaggregation Scheme Using TIR Surface Evaporation Data over the Continental United States

The Soil Moisture Active Passive (SMAP) mission is dedicated toward global soil moisture mapping. Typically, an L-band microwave radiometer has spatial resolution on the order of 36-40 km, which is too coarse for many specific hydro-meteorological and agricultural applications. With the failure of the SMAP active radar within three months of becoming operational, an intermediate (9-km) and finer (3-km) scale soil moisture product solely from the SMAP mission is no longer possible. Therefore, the focus of this study is a disaggregation of the 36-km resolution SMAP passive-only surface soil moisture (SSM) using the Soil Evaporative Efficiency (SEE) approach to spatial scales of 3-km and 9-km. The SEE was computed using thermal-infrared (TIR) estimation of surface evaporation over Continental U.S. (CONUS). The disaggregation results were compared with the 3 months of SMAP-Active (SMAP-A) and Active/Passive (AP) products, while comparisons with SMAP-Enhanced (SMAP-E), SMAP-Passive (SMAP-P), as well as with more than 180 Soil Climate Analysis Network (SCAN) stations across CONUS were performed for a 19 month period. At the 9-km spatial scale, the TIR-Downscaled data correlated strongly with the SMAP-E SSM both spatially (r = 0.90) and temporally (r = 0.87). In comparison with SCAN observations, overall correlations of 0.49 and 0.47; bias of −0.022 and −0.019 and unbiased RMSD of 0.105 and 0.100 were found for SMAP-E and TIR-Downscaled SSM across the Continental U.S., respectively. At 3-km scale, TIR-Downscaled and SMAP-A had a mean temporal correlation of only 0.27. In terms of gain statistics, the highest percentage of SCAN sites with positive gains (>55%) was observed with the TIR-Downscaled SSM at 9-km. Overall, the TIR-based downscaled SSM showed strong correspondence with SMAP-E; compared to SCAN, and overall both SMAP-E and TIR-Downscaled performed similarly, however, gain statistics show that TIR-Downscaled SSM slightly outperformed SMAP-E.

Disaggregation↗

Simulation of ocean SAR images via phase history generation

A method for simulating an ocean synthetic-aperture radar (SAR) image is illustrated for a simple internal wave current pattern. The method calculates both the amplitude image and the radar signal history. The simulation model used consists of six stages. In the first stage, the full wave spectrum in two spatial coordinates is calculated from the wind speed and direction using the action spectral density equation. In the second stage, the pixel size is selected and used as the basis for dividing the spectrum into large and small scale motions. Realizations for the large-scale ocean surface height and velocity are then calculated. In the third stage, the sensor wavelength and geometry are used to calculate the small scale statistics (radar cross section, coherence time, root-mean square (RMS), radial velocity and RMS slope). In the fourth stage, an autoregressive method is used to generate a realization of the surface reflectivity history that is consistent with the radial velocity and the radial velocity variance. In the fifth stage, the signal history is generated by summing the reflectivities at the proper times with the appropriate antenna weighting function for the SAR. As a consequence of this process, speckle is automatically included in the signal. In the final stage, the SAR image can be created by using any of the traditional ways to process the signal history, including variable focusing and multilook processing.

Bennett, John R.↗

Supernova remnant rich fields in the Carina spiral arm

The analysis of the ROSAT PSPC data on five fields containing supernova remnants several months ago was completed. Dr. Una Hwang, prepared a paper describing our results which was published last August ('An X-ray Study of Five Supernova Remnants in the Carina Spiral Arm', by Hwang and Markert, 1994, Ap. J., 431, p. 819). Hwang's earlier analysis of this data became part of her PhD thesis ( 'X-ray Studies of Supernova Remnants', February 1994, MIT). A copy of the Hwang and Markert paper is appended. The results of the study are well-summarized in the Hwang and Markert paper: the best spatial-spectral X-ray study yet made of the intriguing SNR G296.1-0.7 was obtained. This study showed interesting spectral variations over the surface of the object.It was also determined the gross physical properties of G296 based on its X-ray emission. Four other fields were also examined. For three of these the upper limits to the radio object were determined, and in one case found a weak, but statistically significant X-ray object coincident with the peak of the radio flux.

Markert, Thomas↗

Spatial Distribution of Large Cloud Drops

By analyzing aircraft measurements of individual drop sizes in clouds, we have shown in a companion paper (Knyazikhin et al., 2004) that the probability of finding a drop of radius r at a linear scale l decreases as l(sup D(r)) where 0 less than or equal to D(r) less than or equal to 1. This paper shows striking examples of the spatial distribution of large cloud drops using models that simulate the observed power laws. In contrast to currently used models that assume homogeneity and therefore a Poisson distribution of cloud drops, these models show strong drop clustering, the more so the larger the drops. The degree of clustering is determined by the observed exponents D(r). The strong clustering of large drops arises naturally from the observed power-law statistics. This clustering has vital consequences for rain physics explaining how rain can form so fast. It also helps explain why remotely sensed cloud drop size is generally biased and why clouds absorb more sunlight than conventional radiative transfer models predict.

Marshak, A.↗

Passive/active Microwave Soil Moisture Change Disaggregation Using Smapvex12 Data

The SMAPVEX12 (Soil Moisture Active Passive (SMAP) Validation Experiment 2012) experiment was conducted during June-July 2012 in Manitoba, Canada with the goal of collecting remote sensing data and ground measurements for the development and testing of soil moisture retrieval algorithms under varying vegetation and soil conditions for the SMAP satellite. The aircraft based soil moisture data provided by the passive/active microwave sensor PALS (Passive and Active L-band System) has a nominal spatial resolution of 1600 m. However, this resolution is not compatible with agricultural, meteorological and hydrological studies that require high spatial resolutions and this issue can be solved by soil moisture disaggregation. The soil moisture disaggregation algorithm integrates radiometer soil moisture retrievals and high-resolution radar observations and it can provide soil moisture estimates at a finer scale than the radiometer data alone. In this study, a change detection algorithm was used for disaggregation of coarse resolution passive microwave soil moisture retrievals with radar backscatter coefficients obtained from the higher spatial resolution UAVSAR (Unmanned Air Vehicle Synthetic Aperture Radar) at crop field scale. The accuracy of the disaggregated change in soil moisture was evaluated using ground based soil moisture measurements collected during SMAPVEX12 campaign. The results showed that soil moisture spatial variabilities were better characterized by the disaggregated change in soil moisture estimates at 5 m / 800 m resolution as well as good agreement with in situ measurements. It also showed that VWC (Vegetation Water Content) did not have a big impact on disaggregation algorithm performance, with R2 of the disaggregated results ranging 0.628-0.794. The 5 m and 800m resolution disaggregated soil moisture did no show significant difference in statistical performance variables.

Radar Backscatter↗

Retrieval of Aerosol Optical Depth Under Thin Cirrus from MODIS: Application to an Ocean Algorithm

A strategy for retrieving aerosol optical depth (AOD) under conditions of thin cirrus coverage from the Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. We adopt an empirical method that derives the cirrus contribution to measured reflectance in seven bands from the visible to shortwave infrared (0.47, 0.55, 0.65, 0.86, 1.24, 1.63, and 2.12 μm, commonly used for AOD retrievals) by using the correlations between the top-of-atmosphere (TOA) reflectance at 1.38 micron and these bands. The 1.38 micron band is used due to its strong absorption by water vapor and allows us to extract the contribution of cirrus clouds to TOA reflectance and create cirrus-corrected TOA reflectances in the seven bands of interest. These cirrus-corrected TOA reflectances are then used in the aerosol retrieval algorithm to determine cirrus-corrected AOD. The cirrus correction algorithm reduces the cirrus contamination in the AOD data as shown by a decrease in both magnitude and spatial variability of AOD over areas contaminated by thin cirrus. Comparisons of retrieved AOD against Aerosol Robotic Network observations at Nauru in the equatorial Pacific reveal that the cirrus correction procedure improves the data quality: the percentage of data within the expected error +/-(0.03 + 0.05 ×AOD) increases from 40% to 80% for cirrus-corrected points only and from 80% to 86% for all points (i.e., both corrected and uncorrected retrievals). Statistical comparisons with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) retrievals are also carried out. A high correlation (R = 0.89) between the CALIOP cirrus optical depth and AOD correction magnitude suggests potential applicability of the cirrus correction procedure to other MODIS-like sensors.

cirrus correction↗

Satellite-Based Drought Reporting on the Navajo Nation

The Navajo Nation (NN) is the largest reservation in the US, and faces challenges related to water management during long-term and widespread drought episodes. The Navajo Nation is a federally recognized tribe, which has boundaries within Arizona, New Mexico, and Utah. The Navajo Nation has a land area of over 70,000 square kilometers. The Navajo Nation Department of Water Resources (NNDWR) reports on drought and climatic conditions through the use of regional Standardized Precipitation Index (SPI) values and a network of in-situ rainfall, streamflow, and climate data. However, these data sources lack the spatial detail and consistent measurements needed to provide a coherent understanding of the drought regime within the Nation's regional boundaries. This project, as part of NASA's Western Water Applications Office (WWAO), improves upon the recently developed Drought Severity Assessment Tool (DSAT) to ingest satellite-based precipitation data to generate SPI values for specific administrative boundaries within the reservation. The tool aims to: (1) generate SPI values and summary statistics for regions of interest on various timescales, (2) to visualize SPI values within a web-map application, and (3) produce maps and comparative statistical outputs in the format required for annual drought reporting. The co-development of the DSAT with NN partners is integral to increasing the sustained use of Earth Observations for water management applications. This tool will provide data to support the NN in allocation of drought contingency dollars to the regions most adversely impacted by declines in water availability.

Satellite-based↗

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

The land-water interface of coastal marshes may influence the production of estuarine-dependent fisheries more than the area of these marshes. To test this hypothesis, a spatial model was created to explore the dynamic relationship between marshland-water interface and level of disintegration in the decaying coastal marshes of Louisiana's Barataria, Terrebonne, and Timbalier basins. Calibrating the model with Landsat Thematic Mapper satellite imagery, a parabolic relationship was found between land-water interface and marsh disintegration. Aggregated simulation data suggest that interface in the study area will soon reach its maximum and then decline. A statistically significant positive linear relationship was found between brown shrimp catch and total interface length over the past 28 years. This relationship suggests that shrimp yields will decline when interface declines, possibly beginning about 1995.

Browder, Joan A.↗

Methods of Attenuation Correction for the TRMM Precipitation Radar Data

The surface reference technique (SRT) has been studied extensively both theoretically and experimentally over the last decade. It is only with the launch of the Tropical Rain Measuring Mission (TRMM) Precipitation Radar, however, that we can begin to test directly whether the technique provides a reliable means of attenuation correction for spaceborne weather radars. Preliminary results indicate that the method yields results that are qualitatively reasonable when the rain rate is moderate or high and when the surface provides a stable reference value. The structure of the normalized radar cross section of the surface (NRCS), however, is highly complex in the sense that the statistics change with background type, incidence angle, location, and season. As a consequence of this, the reliability of the path attenuation estimate is not proportional simply to the amount of rain attenuation but to the rain attenuation relative to the inherent variability of the reference target. To illustrate the behavior of the SRT, several overpasses of the Hurricane Bonnie are shown. We also show that the proper choice of reference data set (spatial, temporal and global) can be critical to the success of the method.

Meneghini, Robert↗

Multi-Group Maximum Entropy Model for Translational Non-Equilibrium

The aim of the current work is to describe a new model for flows in translational non- equilibrium. Starting from the statistical description of a gas proposed by Boltzmann, the model relies on a domain decomposition technique in velocity space. Using the maximum entropy principle, the logarithm of the distribution function in each velocity sub-domain (group) is expressed with a power series in molecular velocity. New governing equations are obtained using the method of weighted residuals by taking the velocity moments of the Boltzmann equation. The model is applied to a spatially homogeneous Boltzmann equation with a Bhatnagar-Gross-Krook1(BGK) model collision operator and the relaxation of an initial non-equilibrium distribution to a Maxwellian is studied using the model. In addition, numerical results obtained using the model for a 1D shock tube problem are also reported.

Entropy↗

Analysis techniques for two-dimensional infrared data

In order to evaluate infrared detection and remote sensing systems, it is necessary to know the characteristics of the observational environment. For both scanning and staring sensors, the spatial characteristics of the background may be more of a limitation to the performance of a remote sensor than system noise. This limitation is the so-called spatial clutter limit and may be important for systems design of many earth application and surveillance sensors. The data used in this study is two dimensional radiometric data obtained as part of the continuing NASA remote sensing programs. Typical data sources are the Landsat multi-spectral scanner (1.1 micrometers), the airborne heat capacity mapping radiometer (10.5 - 12.5 micrometers) and various infrared data sets acquired by low altitude aircraft. Techniques used for the statistical analysis of one dimensional infrared data, such as power spectral density (PSD), exceedance statistics, etc. are investigated for two dimensional applicability. Also treated are two dimensional extensions of these techniques (2D PSD, etc.), and special techniques developed for the analysis of 2D data.

Winter, E. M.↗

Estimation of Supersonic Stage Separation Aerodynamics of Winged-Body Launch Vehicles Using Response Surface Methods

A collection of statistical and mathematical techniques referred to as response surface methodology was used to estimate the longitudinal stage separation aerodynamic characteristics of a generic, bimese, winged multi-stage launch vehicle configuration using data obtained on small-scale models at supersonic speeds in the NASA Langley Research Center Unitary Plan Wind Tunnel. The simulated Mach 3 staging was dominated by multiple shock wave interactions between the orbiter and booster vehicles throughout the relative spatial locations of interest. This motivated a partitioning of the overall inference space into several contiguous regions within which the separation aerodynamics were presumed to be well-behaved and estimable using cuboidal and spherical central composite designs capable of fitting full second-order response functions. The primary goal was to approximate the underlying overall aerodynamic response surfaces of the booster vehicle in belly-to-belly proximity to the orbiter vehicle using relatively simple, lower-order polynomial functions that were piecewise-continuous across the full independent variable ranges of interest. The quality of fit and prediction capabilities of the empirical models were assessed in detail, and the issue of subspace boundary discontinuities was addressed. The potential benefits of augmenting the central composite designs to full third order using computer-generated D-optimality criteria were also evaluated. The usefulness of central composite designs, the subspace sizing, and the practicality of fitting low-order response functions over a partitioned inference space dominated by highly nonlinear and possibly discontinuous shock-induced aerodynamics are discussed.

Erickson, Gary E.↗