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 757 records · Page 42

Statistical study of ULF wave occurrence in the dayside magnetosphere

Ultralow-frequency (ULF) waves are observed almost everywhere in the dayside magnetosphere. The mechanism by which these waves are generated and transformed in the dayside magnetosphere is still not understood. Here we report a statistical study of these waves based on magnetic field data from the International Sun-Earth Explorer 1 (ISEE 1) spacecraft. Data from the first traversal of the spacecraft through the entire dayside magnetosphere have been examined to determine the spatial distribution of wave occurrence. Successive 20-min segments of data were transformed to a field-aligned coordinate system. The parallel component was detrended and all three components of the field spectrally analyzed. Wave occurrence was defined by the presence of significant peaks in the power spectra. Wave events were categorized by three wave frequency bands: Pc 3 with T approximately 10-45 s; Pc 4 with T approximately 45-150 s; the short-period part of the Pc 5 wave band with T approximately 150-324 s. Properties of the spectral peaks were then entered into a data base. The data base was next sorted to determine the spatial occurrence pattern for the waves. Our results show that Pc 3 waves most frequently occur just outside synchronous orbit and are approximately centered on local noon. Pc 4 waves have a similar distribution with its peak further out. Pc 5 waves have high occurrence rate at the two flanks of the magnetosphere. Peaks in spectra obtained near the magnetopause are less clearly defined than those deeper in the magnetosphere.

Cao, M.↗

The Ultra-Luminous X-ray Source Population from the Chandra Archive of Galaxies

One hundred fifty-four discrete non-nuclear Ultra-Luminous X-ray (ULX) sources, with spectroscopically-determined intrinsic X-ray luminosities greater than 1 e39 ergs/s, are identified in 82 galaxies observed with Chandra's Advanced CCD Imaging Spectrometer. Source positions, X-ray luminosities, and spectral and timing characteristics are tabulated. Statistical comparisons between these X-ray properties and those of the weaker discrete sources in the same fields (mainly neutron star and stellar-mass black hole binaries) are made. Sources above approximately le38 ergs per second display similar spatial, spectral, color, and variability distributions. In particular, there is no compelling evidence in the sample for a new and distinct class of X-ray object such as the intermediate-mass black holes. 83% of ULX candidates have spectra that can be described as absorbed power laws with index = 1.74 and column density = 2.24e21 l per square centimeter, or approximately 5 times the average Galactic column. About 20% of the ULX's have much steeper indices indicative of a soft, and likely thermal, spectrum. The locations of ULXs in their host galaxies are strongly peaked towards their galaxy centers. The deprojected radial distribution of the ULX candidates is somewhat steeper than an exponential disk, indistinguishable from that of the weaker sources. About 5--15% of ULX candidates are variable during the Chandra observations (which average 39.5 ks). Comparison of the cumulative X-ray luminosity functions of the ULXs to Chandra Deep Field results suggests approximately 25% of the sources may be background objects including 14% of the ULX candidates in the sample of spiral galaxies and 44% of those in elliptical galaxies implying the elliptical galaxy ULX population is severely compromised by background active galactic nuclei. Correlations with host galaxy properties confirm the number and total X-ray luminosity of the ULXs are associated with recent star formation and with galaxy merging and interactions. The preponderance of ULXs in star-forming galaxies as well as their similarities to less-luminous sources suggest they originate in a young but short-lived population such as the high-mass X-ray binaries with a smaller contribution (based on spectral slope) from recent supernovae. The number of ULXs in elliptical galaxies scales with host galaxy mass and can be explained most simply as the high-luminosity end of the low-mass X-ray binary population.

Swartz, Douglas A.↗

Prediction of Acoustic Environments from Horizontal Rocket Firings

In recent years, advances in research and engineering have led to more powerful launch vehicles which can reach areas of space not yet explored. These more powerful vehicles yield acoustic environments potentially destructive to the vehicle or surrounding structures. Therefore, it has become increasingly important to be able to predict the acoustic environments created by these vehicles in order to avoid structural and/or competent failure. The current industry standard technique for predicting launch-induced acoustic environments was developed by Eldred in the early 1970's and is published in NASA SP-80721. Recent work2 has shown Eldred's technique to be inaccurate for current state-of-the-art launch vehicles. Due to the high cost of full-scale and even sub-scale rocket experiments, very little rocket noise data is available. Furthermore, much of the work thought to be applicable to rocket noise has been done with heated jets. Tam3,4 has done an extensive amount of research on jets of different nozzle exit shape, diameter, velocity, and temperature. Though the values of these parameters, especially exit velocity and temperature, are often very low compared to these values in rockets, a lot can be learned about rocket noise from jet noise literature. The turbulent nature of jet and rocket exhausts is quite similar. Both exhausts contain turbulent structures of varying scale-termed the fine and large scale turbulence by Tam. The finescale turbulence is due to small eddies from the jet plume interacting with the ambient atmosphere. According to Tam et al., the noise radiated by this envelope of small-scale turbulence is statistically isotropic. Hence, one would expect the noise from the small scale turbulence of the jet to be nearly omni-directional. The coherent nature of the large-scale turbulence results in interference of the noise radiated from different spatial locations within the jet. This interference-whether it is constructive or destructive-results in highly directional noise radiation. Tam3 has proposed a model to predict the acoustic environment due to jets and while it works extremely well for jets, it was found to be inappropriate for rockets8. A model to predict the acoustic environment due to a launch vehicle in the far-field which incorporates concepts from both Eldred and Tam was created. This was done using five sets of horizontally fired rocket data, obtained between 2008 and 2012. Three of these rockets use solid propellant and two use liquid propellant. Through scaling analysis, it is shown that liquid and solid rocket motors exhibit similar spectra at similar amplitudes. This model is accurate for these five data sets within 5 dB of the measured data for receiver angles of 30deg to 160deg (with respect to the downstream exhaust centerline). The model uses the following vehicle parameters: nozzle exit diameter and velocity, radial distance from source to receiver, receiver angle, mass flow rate, and acoustic efficiency.

Giacomoni, Clothilde↗

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗

A statistical study of plasma sheet dynamics using ISEE 1 and 2 energetic particle flux data

Plasma sheet dynamics during substorms are studied by analyzing 461 cases of transient dropout events of the 1.5 and 6-keV particle fluxes detected by ISEE 1 and 2 satellites. The instruments for detecting low- and high-energy particles are described. The spatial distribution of flux dropout events, and the events' relationship to magnetospheric activity level are examined. Substorm events without observed flux dropout events are investigated. The data reveal that the flux dropout distribution is isotropic, between 12-23 earth radii, and is present in the entire nightside plasma sheet; and the substorms without flux dropout are more frequent near earth and magnetospheric flanks. It is observed that tailward of 12 earth radii the flux dropout events and substorms without flux dropout are similar. The Chao et al. (1977) MHD rarefaction wave propagation model and the Hones (1973, 1980) near-tail, X-type magnetic neutral line formation model are discussed and compared to the experimental data. It is noted that neither model explains the plasma sheet dynamics observed.

Dandouras, J.↗

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. A diffraction pattern analysis of MSS images led to the development of spatial signatures for farm land, urban areas and mountains. Four spatial features are employed to describe the spatial characteristics of image cells in the digital data. Three spectral features are combined with the spatial features to form a seven dimensional vector describing each cell. Then, the classification of the feature vectors is accomplished by using the maximum likelihood criterion. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. It was also determined that for a given geographic area the statistics of the classes remain invariable for a period of a month, but vary substantially between seasons. Three ERTS-1 images from the Phoenix, Arizona area were processed, and recognition rates between 85% and 100% were obtained for the terrain classes of desert, farms, mountains, and urban areas. To eliminate the need for training data, a new clustering algorithm has been developed. Seven ERTS-1 images from four test sites have been processed through the clustering algorithm, and high recognition rates have been achieved for all terrain classes.

Gramenopoulos, N.↗

High-beta turbulence in two-dimensional magnetohydrodynamics

Equations of ideal magnetohydrodynamics are used to study incompressible turbulent flows in a specified geometry where all the field quantities vary with only two spatial dimensions. The procedures adopted are basically those of Kraichnan (1967), in which classical equilibrium ensembles are built around constants of the motion identified from the Fourier-transformed equations of motion. Once the constants of the motion are identified, the statistical formulation of the problem is presented in a phase space whose coordinates are the real and imaginary parts of the Fourier coefficients. Canonical ensembles are constructed in this phase space by classical arguments. The theory developed permits isolation of some qualitatively new gross physical effects which have so far not been calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long-wavelength, self-generated, magnetic fields for a wide range of initial parameters.

Fyfe, D.↗

The possible modifications of the Hisse model for pure LANDSAT agricultural data

An idea, due to A. Feiveson, is presented for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field). Theoretical arguments are given which show that any significant refinement of the model beyond Feiveson's proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Guassian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗

The possible modifications of the HISSE model for pure LANDSAT agricultural data

A method for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field) is discussed. Theoretical arguments are given which show that any significant refinement of the model beyond this proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Gaussian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.↗

Recent Advances on INSAR Temporal Decorrelation: Theory and Observations Using UAVSAR

We review our recent advances in understanding the role of temporal decorrelation in SAR interferometry and polarimetric SAR interferometry. We developed a physical model of temporal decorrelation based on Gaussian-statistic motion that varies along the vertical direction in forest canopies. Temporal decorrelation depends on structural parameters such as forest height, is sensitive to polarization and affects coherence amplitude and phase. A model of temporal-volume decorrelation valid for arbitrary spatial baseline is discussed. We tested the inversion of this model to estimate forest height from model simulations supported by JPL/UAVSAR data and lidar LVIS data. We found a general good agreement between forest height estimated from radar data and forest height estimated from lidar data.

polarimetric SAR interferometry↗

Correlation techniques and measurements of wave-height statistics

Statistical measurements of wave height fluctuations have been made in a wind wave tank. The power spectral density function of temporal wave height fluctuations evidenced second-harmonic components and an f to the minus 5th power law decay beyond the second harmonic. The observations of second harmonic effects agreed very well with a theoretical prediction. From the wave statistics, surface drift currents were inferred and compared to experimental measurements with satisfactory agreement. Measurements were made of the two dimensional correlation coefficient at 15 deg increments in angle with respect to the wind vector. An estimate of the two-dimensional spatial power spectral density function was also made.

Guthart, H.↗

The morphology of energetic O/+/ ions during two magnetic storms - Temporal variations

Results are presented for a statistical analysis of the temporal variations of precipitating O(+) and H(+) ions in the energy range 0.7-12 keV during the magnetic storms on Dec. 16-18, 1971. Emphasis is on the temporal variations of parameters describing the intensity, average energy, and spatial location of the precipitation zones of the two ionic species. It is shown that the intensity of the precipitating O(+) ions correlate well with the geomagnetic indices which measure the strength of magnetospheric substorm activity and the strength of the storm time ring current, indicating that a previously unknown strong coupling mechanism existed between the magnetosphere and the ionosphere during the storm period. Correlations are found between the locations of the O(+) and H(+) precipitation zones and between the average energies of the two ionic species. Precipitation is shown to be an important loss mechanism for ring current ions with energy less than 12 keV during the magnetic storm period studied.

Sharp, R. D.↗

Edge detection - Image-plane versus digital processing

To optimize edge detection with the familiar Laplacian-of-Gaussian operator, it has become common to implement this operator with a large digital convolution mask followed by some interpolation of the processed data to determine the zero crossings that locate edges. It is generally recognized that this large mask causes substantial blurring of fine detail. It is shown that the spatial detail can be improved by a factor of about four with either the Wiener-Laplacian-of-Gaussian filter or an image-plane processor. The Wiener-Laplacian-of-Gaussian filter minimizes the image-gathering degradations if the scene statistics are at least approximately known and also serves as an interpolator to determine the desired zero crossings directly. The image-plane processor forms the Laplacian-of-Gaussian response by properly combining the optical design of the image-gathering system with a minimal three-by-three lateral-inhibitory processing mask. This approach, which is suggested by Marr's model of early processing in human vision, also reduces data processing by about two orders of magnitude and data transmission by up to an order of magnitude.

Huck, Friedrich O.↗

Calibration of the DSCOVR EPIC Visible and NIR Channels using Multiple LEO Radiometers

The first five years of operation of the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) at the Lagrange one point have produced results that uniquely complement the data from currently operating low orbit Earth-observing instruments. In this paper we describe an updated unified approach to EPIC calibration. In this approach, calibration coefficients and their trends were obtained by comparing EPIC observations to the measurements from polar orbiting radiometers. In this study L1B reflectances from Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua and Terra satellites, Multi-angle Imaging Spectroradiometer (MISR) onboard Terra and Visible Infrared Imaging Radiometer (VIIRS) onboard the Suomi National Polar-orbiting Partnership (Suomi NPP) spacecraft were used to infer calibration coefficients for four EPIC visible and near-infrared channels: 443 nm, 551 nm, 680 nm, and 780 nm. EPIC Version three measurements made between June 2015 and August 2020 were used for comparison. The calibration procedure identifies the most homogeneous low Earth orbit radiometer scenes matching scattering angles that are temporarily and spatially collocated with EPIC observations. These scenes are used to determine reflectance to count (R/C) ratios in spectrally analogous channels. Seasonal average R/C ratios were analyzed to obtain EPIC calibration gains and trends. The trends for the full dataset period are not statistically significant except in the 443 nm channel. No significant changes in calibration were found after the instrument’s exit from safe hold in March 2020. The R/C ratios were also used to determine the differences in EPIC gains resulting from separate calibrations: against MODIS Aqua or Terra, as well as against forward or aftward MISR cameras. Statistical tests indicate that the differences between the two datasets are not significant except in the 780 nm channels where Aqua-derived coefficients may be around 2% lower compared to Terra. The dependence of EPIC calibration gains on the instrument scattering angle and on DSCOVR-Earth distance were investigated. Lastly, model Low Earth Orbit (LEO) reflectances calculated to match the EPIC viewing geometry were employed to study how EPIC calibration coefficients depend on EPIC-LEO viewing geometry differences. The effect of LEO and EPIC angular mismatch on calibration was shown to be small.

deep space climate observatory/earth polychromatic↗

Data Science and the Knowledge Discovery Adventure

This talk will cover the important steps involved in the data science and knowledge discovery process: • Initial fact gathering (interview domain experts, review reports, articles, state-of-the-art) • Identify the problem (prediction, classification, statistical analysis, etc.) • Survey supporting data sources • Understand the data (numerical, categorical, text, sampling rate, data quality issues, etc.) • Selecting relevant features and sources • Acquire the data (set up agreements with the data stewards, APIs to download, etc.) • Merge data sources (temporal, spatial, common key, other ontologies...) • Feature Engineering (non linear domain knowledge or physics-based relationships) • Build data processing pipeline (may need to tap into data stream, develop parallel processing algorithm, federated learning etc.) • Build model and test (tune hyper-parameters, cross validation.) • Analyze/Validate results (do the results make sense. Does it answer the original question). • Deploy/Publish (Monitor and assess benefits)

Data science↗

Modeling of topographic effects on Antarctic sea ice using multivariate adaptive regression splines

The role of seafloor topography in the spatial variations of the southern ocean sea ice cover as observed (every other day) by the Nimbus 7 scanning multichannel microwave radiometer satellite in the years 1980, 1983, and 1984 is studied. Bottom bathymetry can affect sea ice surface characteristics because of the basically barotropic circulation of the ocean south of the Antarctic Circumpolar current. The main statistical tool used to quantify this effect is a local nonparametric regression model of sea ice concentration as a function of the depth and its first two derivatives in both meridional and zonal directions. First, we model the relationship of bathymetry to sea ice concentration in two sudy areas, one over the Maud Rise and the other over the Ross Sea shelf region. The multiple correlation coefficient is found to average 44% in the Maud Rise study area and 62% in the Ross Sea study area over the years 1980, 1983, and 1984. Second, a strategy of dividing the entire Antarctic region into an overlapping mosaic of small areas, or windows is considered. Keeping the windows small reduces the correlation of bathymetry with other factors such as wind, sea temperature, and distance to the continent. We find that although the form of the model varies from window to window due to the changing role of other relevant environmental variables, we are left with a spatially consistent ordering of the relative importance of the topographic predictors. For a set of three representative days in the Austral winter of 1980, the analysis shows that an average of 54% of the spatial variation in sea ice concentration over the entire ice cover can be attributed to topographic variables. The results thus support the hypothesis that there is a sea ice to bottom bathymetry link. However this should not undermine the considerable influence of wind, current, and temperature which affect the ice distribution directly and are partly responsible for the observed bathymetric effects.

De Veaux, Richard D.↗