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 973 records · Page 54

Advances in Landslide Hazard Forecasting: Evaluation of Global and Regional Modeling Approach

A prototype global satellite-based landslide hazard algorithm has been developed to identify areas that exhibit a high potential for landslide activity by combining a calculation of landslide susceptibility with satellite-derived rainfall estimates. A recent evaluation of this algorithm framework found that while this tool represents an important first step in larger-scale landslide forecasting efforts, it requires several modifications before it can be fully realized as an operational tool. The evaluation finds that the landslide forecasting may be more feasible at a regional scale. This study draws upon a prior work's recommendations to develop a new approach for considering landslide susceptibility and forecasting at the regional scale. This case study uses a database of landslides triggered by Hurricane Mitch in 1998 over four countries in Central America: Guatemala, Honduras, EI Salvador and Nicaragua. A regional susceptibility map is calculated from satellite and surface datasets using a statistical methodology. The susceptibility map is tested with a regional rainfall intensity-duration triggering relationship and results are compared to global algorithm framework for the Hurricane Mitch event. The statistical results suggest that this regional investigation provides one plausible way to approach some of the data and resolution issues identified in the global assessment, providing more realistic landslide forecasts for this case study. Evaluation of landslide hazards for this extreme event helps to identify several potential improvements of the algorithm framework, but also highlights several remaining challenges for the algorithm assessment, transferability and performance accuracy. Evaluation challenges include representation errors from comparing susceptibility maps of different spatial resolutions, biases in event-based landslide inventory data, and limited nonlandslide event data for more comprehensive evaluation. Additional factors that may improve algorithm performance accuracy include incorporating additional triggering factors such as tectonic activity, anthropogenic impacts and soil moisture into the algorithm calculation. Despite these limitations, the methodology presented in this regional evaluation is both straightforward to calculate and easy to interpret, making results transferable between regions and allowing findings to be placed within an inter-comparison framework. The regional algorithm scenario represents an important step in advancing regional and global-scale landslide hazard assessment and forecasting.

Kirschbaum, Dalia B.↗

Numerical Model Sensitivity to Heterogeneous Satellite Derived Vegetation Roughness

The sensitivity of a mesoscale weather prediction model to a 1 km satellite-based vegetation roughness initialization is investigated for a domain within the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database, specific to the domain, that was previously derived using a physically based procedure and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and iii) a MODIS climatologic roughness database that like (i) is a function only of land cover type, but possesses domain specific mean values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations within a domain including Oklahoma, Eastern Arkansas, and Northwest Louisiana during a 4-day period within IHOP 2002. Sensitivity analysis compares the impact the three roughness initializations on time-series temperature, precipitation probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE). Overall, the results indicate that, for the current investigation, replacement of the standard look-up table values with the satellite-derived values statistically improves model performance for most observed variables. Such natural roughness heterogeneity enhances the surface wind speed, PBL height and TKE production up to 10 percent, with a lesser effect over grassland, and greater effect over mixed land cover domains.

Jasinski, Michael↗

An AeroCom–AeroSat study: intercomparison of satellite AOD datasets for aerosol model evaluation

To better understand and characterize current uncertainties in the important observational constraint of climate models of aerosol optical depth (AOD), we evaluate and intercompare 14 satellite products, representing nine different retrieval algorithm families using observations from five different sensors on six different platforms. The satellite products (super-observations consisting of 1°×1° daily aggregated retrievals drawn from the years 2006, 2008 and 2010) are evaluated with AErosol RObotic NETwork (AERONET) and Maritime Aerosol Network (MAN) data. Results show that different products exhibit different regionally varying biases (both under- and overestimates) that may reach ±50 %, although a typical bias would be 15 %–25 % (depending on the product). In addition to these biases, the products exhibit random errors that can be 1.6 to 3 times as large. Most products show similar performance, although there are a few exceptions with either larger biases or larger random errors. The intercomparison of satellite products extends this analysis and provides spatial context to it. In particular, we show that aggregated satellite AOD agrees much better than the spatial coverage (often driven by cloud masks) within the 1°×1° grid cells. Up to ∼50 % of the difference between satellite AOD is attributed to cloud contamination. The diversity in AOD products shows clear spatial patterns and varies from 10 % (parts of the ocean) to 100 % (central Asia and Australia). More importantly, we show that the diversity may be used as an indication of AOD uncertainty, at least for the better performing products. This provides modellers with a global map of expected AOD uncertainty in satellite products, allows assessment of products away from AERONET sites, can provide guidance for future AERONET locations and offers suggestions for product improvements. We account for statistical and sampling noise in our analyses. Sampling noise, variations due to the evaluation of different subsets of the data, causes important changes in error metrics. The consequences of this noise term for product evaluation are discussed.

Aerosol↗

Advances in gamma-ray burst astronomy

The paper reviews the current status of observational research into the highly puzzling problem of cosmic gamma-ray bursts. Up to now there is no convincing temporal or spatial correlation with any known celestial processes or objects. The typical burst source strength is somewhere between 10 to the 26th power ergs and 10 to the 52nd power ergs. A list of 39 confirmed and 9 unconfirmed cosmic gamma-ray bursts observed by satellite is given, showing times of occurrence and in some cases, the size in erg per square centimeter. Several approaches to the problem of source object identification are discussed: (1) accumulation of observation statistics with their present poor resolution and research for correlative phenomena, (2) mapping out celestial source distributions with moderate resolution in order to search for galactic or other anisotropies in direction or to look for repeating source regions, and (3) very high-accuracy localization of the source directions of one or several bursts in order to pinpoint a tiny fraction of the celestial sphere for correlative radio, optical and X-ray studies. Planned future instrumentation for deep space probes and multiple-balloon studies is described.

Cline, T. L.↗

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed, 2D surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a narrow band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary-layer flow. Random irregularities which are spatially homogeneous in nature, and also possess a continuous spectrum in the streamwise direction, lead to instability amplitudes that are intermediate to those caused by the periodic and isolated irregularities, respectively, being, in fact, of the same order as the geometrical mean of the amplitudes in the latter two cases. A physical explanation for these asymptotic scalings is given, in addition to providing an analytical expression for the expected value of the instability amplitude for an ensemble of statistically irregular surfaces with random phase distributions.

Choudhari, Meelan↗

Assessment of Global Cloud Datasets from Satellites: Project and Database Initiated by the GEWEX Radiation Panel

Clouds cover about 70% of the Earth's surface and play a dominant role in the energy and water cycle of our planet. Only satellite observations provide a continuous survey of the state of the atmosphere over the whole globe and across the wide range of spatial and temporal scales that comprise weather and climate variability. Satellite cloud data records now exceed more than 25 years in length. However, climatologies compiled from different satellite datasets can exhibit systematic biases. Questions therefore arise as to the accuracy and limitations of the various sensors. The Global Energy and Water cycle Experiment (GEWEX) Cloud Assessment, initiated in 2005 by the GEWEX Radiation Panel, provided the first coordinated intercomparison of publically available, standard global cloud products (gridded, monthly statistics) retrieved from measurements of multi-spectral imagers (some with multiangle view and polarization capabilities), IR sounders and lidar. Cloud properties under study include cloud amount, cloud height (in terms of pressure, temperature or altitude), cloud radiative properties (optical depth or emissivity), cloud thermodynamic phase and bulk microphysical properties (effective particle size and water path). Differences in average cloud properties, especially in the amount of high-level clouds, are mostly explained by the inherent instrument measurement capability for detecting and/or identifying optically thin cirrus, especially when overlying low-level clouds. The study of long-term variations with these datasets requires consideration of many factors. A monthly, gridded database, in common format, facilitates further assessments, climate studies and the evaluation of climate models.

Stubenrauch, C. J.↗

Drought Indicators Based on Model Assimilated GRACE Terrestrial Water Storage Observations

The Gravity Recovery and Climate Experiment (GRACE) twin satellites observe time variations in Earth's gravity field which yield valuable information about changes in terrestrial water storage (TWS). GRACE is characterized by low spatial (greater than 150,000 square kilometers) and temporal (greater than 10 day) resolution but has the unique ability to sense water stored at all levels (including groundwater) systematically and continuously. The GRACE Data Assimilation System (GRACE-DAS), based on the Catchment Land Surface Model (CLSM) enhances the value of the GRACE water storage data by enabling spatial and temporal downscaling and vertical decomposition into moisture 39 components (i.e. groundwater, soil moisture, snow), which individually are more useful for scientific applications. In this study, GRACE-DAS was applied to North America and GRACE-based drought indicators were developed as part of a larger effort that investigates the possibility of more comprehensive and objective identification of drought conditions by integrating spatially, temporally and vertically disaggregated GRACE data into the U.S. and North American Drought Monitors. Previously, the Drought Monitors lacked objective information on deep soil moisture and groundwater conditions, which are useful indicators of drought. Extensive datasets of groundwater storage from USGS monitoring wells and soil moisture from the Soil Climate Analysis Network (SCAN) were used to assess improvements in the hydrological modeling skill resulting from the assimilation of GRACE TWS data. The results point toward modest, but statistically significant, improvements in the hydrological modeling skill across major parts of the United States, highlighting the potential value of GRACE assimilated water storage field for improving drought detection.

Houborg, Rasmus↗

Heliocentric Distance and Solar Activity Dependence of Sustained Quasi-radial Interplanetary Magnetic Field Occurrence

Planets close to their stars experience an interplanetary magnetic field (IMF) that is dominantly quasi-radial. Our solar system serves as a laboratory to study how the occurrence of quasi-radial IMF varies away from the star and under different stellar activities. Furthermore, on time and spatial scales relevant to magnetospheric physics, solar wind variability prevails in the form of structures generated both at the Sun and locally in the interplanetary space. The stationary Parker spiral model only approximates the large-scale structure of the IMF. Deviations from the Parker spiral often result in strongly radial magnetic fields that give rise to kinetic foreshock turbulence, which in turn can impact planetary magnetospheres. The relative significance of this type of interaction can be estimated statistically based on the occurrence rate of cases where the IMF is directed along the radial direction, leading to the entire day-side magnetosphere being downstream of the ion foreshock. We use observations covering radial distances from 0.1 to 10 au and more than 2 solar cycles to quantify the prevalence of radial IMFs throughout the heliosphere. Near Earth's orbit, it is found that the occurrence rates of quasi-radial and southward IMF orientations are similar, and that the Pearson correlation coefficient is ${{ \mathcal R }}_{{xy}}\sim -0.7$ calculated between quasi-radial IMF occurrence rate and solar activity. A negative correlation is demonstrated for radial distances extending to at least Mars but not to Saturn.

Brandon L Burkholder↗

Statistical Distribution of Inflation on Lava Flows: Analysis of Flow Surfaces on Earth and Mars

The surface morphology of a lava flow results from processes that take place during the emplacement of the flow. Certain types of features, such as tumuli, lava rises and lava rise pits, are indicators of flow inflation or endogenous growth of a lava flow. Tumuli in particular have been identified as possible indicators of tube location, indicating that their distribution on the surface of a lava flow is a junction of the internal pathways of lava present during flow emplacement. However, the distribution of tumuli on lava flows has not been examined in a statistically thorough manner. In order to more rigorously examine the distribution of tumuli on a lava flow, we examined a discrete flow lobe with numerous lava rises and tumuli on the 1969 - 1974 Mauna Ulu flow at Kilauea, Hawaii. The lobe is located in the distal portion of the flow below Holei Pali, which is characterized by hummocky pahoehoe flows emplaced from tubes. We chose this flow due to its discrete nature allowing complete mapping of surface morphologies, well-defined boundaries, well-constrained emplacement parameters, and known flow thicknesses. In addition, tube locations for this Mauna Ulu flow were mapped by Holcomb (1976) during flow emplacement. We also examine the distribution of tumuli on the distal portion of the hummocky Thrainsskjoldur flow field provided by Rossi and Gudmundsson (1996). Analysis of the Mauna Ulu and Thrainsskjoldur flow lobes and the availability of high-resolution MOC images motivated us to look for possible tumuli-dominated flow lobes on the surface of Mars. We identified a MOC image of a lava flow south of Elysium Mons with features morphologically similar to tumuli. The flow is characterized by raised elliptical to circular mounds, some with axial cracks, that are similar in size to the tumuli measured on Earth. One potential avenue of determining whether they are tumuli is to look at the spatial distribution to see if any patterns similar to those of tumuli-dominated terrestrial flows can be identified. Since tumuli form by the injection of lava beneath a crust, the distribution of tumuli on a flow should represent the distribution of thermally preferred pathways beneath the surface of the crust. That distribution of thermally preferred pathways may be a function of the evolution of a basaltic lava flow. As a longer-lived flow evolves, initially broad thermally preferred pathways would evolve to narrower, more well-defined tube-like pathways. The final flow morphology clearly preserves the growth of the flow over time, with inflation features indicating pathways that were not necessarily contemporaneously active. Here, we test using statistical analysis whether this final flow morphology produces distinct distributions that can be used to readily determine the distribution of thermally preferred pathways beneath the surface of the crust.

Glazel, L. S.↗

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗

The reduction, verification and interpretation of MAGSAT magnetic data over Canada

Analysis of Investigator B tapes was extended to use a statistical screening of individual data points in each grid cell over the map area. The mean value of data from each cell was computed and all data values more than 4nT from the mean were rejected. The revised mean and standard deviation was computed. The average standard deviation for all map grid cells was 1.4nT for the scalar anomaly map and about 2nT for the vector component anomaly maps. At this stage, the grid cell is larger and spatial resolution lower for the vector maps than for the scalar maps, because the vector components are considerably more disturbed and more data points were rejected. An approach for altitude adjustment is in progress so that corrections can be made and data reduction procede. To evaluate the use of scalar and vector MAGSAT data for the production of regional magnetic charts, charts of the north, east, vertical, and scalar field derived from MAGSAT and from Canadian aeromagnetic data were compared. It appears possible to retrieve crustal signals from the horizontal vector components despite external disturbance fields at high latitude.

Coles, R. L.↗

A comparison of measured and calculated optical properties of atmospheric aerosols at infrared wavelengths

Measurements of 10.6-micron lidar backscatter were compared with calculated backscatter based on nearly simultaneous observations of stratospheric and tropospheric aerosol size distributions. It was found that there is better agreement in the troposphere, even though the uncertainties of the calculation are greater for this region due to the variables in both the spatial concentration and the physical makeup of the aerosol. A second comparison study was made to test the consistency of the mean tropospheric extinction values at 1.02 micron (as reported by the SAGE satellite) with the values calculated from an ensemble of 400 measured size distributions thought to be representative of midcontinental tropospheric aerosol. The two methods produce consistent results within the expected degree of uncertainty. The ensemble of 400 'proven' size distributions is then used to calculate a statistical relationship between the 1.02-micron extinction and the 10.6-micron backscatter.

Rosen, James M.↗

Thermal and Driven Stochastic Growth of Langmuir Waves in the Solar Wind and Earth's Foreshock

Statistical distributions of Langmuir wave fields in the solar wind and the edge of Earth's foreshock are analyzed and compared with predictions for stochastic growth theory (SGT). SGT quantitatively explains the solar wind, edge, and deep foreshock data as pure thermal waves, driven thermal waves subject to net linear growth and stochastic effects, and as waves in a pure SGT state, respectively, plus radiation near the plasma frequency f(sub p). These changes are interpreted in terms of spatial variations in the beam instability's growth rate and evolution toward a pure SGT state. SGT analyses of field distributions are shown to provide a viable alternative to thermal noise spectroscopy for wave instruments with coarse frequency resolution, and to separate f(sub p) radiation from Langmuir waves.

Cairns, Iver H.↗

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

Radar Cross Section of Orbital Debris Objects

This discussion is concerned with the radar-data analysis and usage involved in the building of model orbital debris (OD) populations in the near-Earth environment, focusing on radar cross section (RCS). While varying with radar wavelength, physical dimension, material composition, overall shape and structure, the RCS of an irregular object is also strongly dependent on its spatial orientation. The historical records of observed RCSs for cataloged OD objects in the Space Surveillance Network are usually distributed over an RCS range, forming respective characteristic patterns. The National Aeronautics and Space Administration (NASA) Size Estimation Model provides an empirical probability-density function of RCS as a function of “effective diameter” (or “characteristic length”), which makes it feasible to predict possible RCS distributions for a given model OD population and to link data with model from a statistical perspective. The discussion also includes application of the widely used method of moments (MoM) and the Generalized Multi-particle Mie-solution (GMM) in the prediction of the RCS of arbitrarily shaped objects. Theoretical calculation results for an aluminum cube are compared with corresponding experimental measurements.

Xu, Yu-Lin↗

Statistical Properties of the Population of the Galactic Center Filaments II: The Spacing between Filaments

We carry out a population study of magnetized radio filaments in the Galactic centre using MeerKAT data by focusing on the spacing between the filaments that are grouped. The morphology of a sample of 43 groupings containing 174 magnetized radio filaments are presented. Many grouped filaments show harp-like, fragmented cometary tail-like, or loop-like structures in contrast to many straight filaments running mainly perpendicular to the Galactic plane. There are many striking examples of a single filament splitting into two prongs at a junction, suggestive of a flow of plasma along the filaments. Spatial variations in spectral index, brightness, bending, and sharpening along the filaments indicate that they are evolving on a 105−6-yr time-scale. The mean spacings between parallel filaments in a given grouping peaks at ∼16 arcsec. We argue by modeling that the filaments in a grouping all lie on the same plane and that the groupings are isotropically oriented in 3D space. One candidate for the origin of filamentation is interaction with an obstacle, which could be a compact radio source, before a filament splits and bends into multiple filaments. In this picture, the obstacle or sets the length scale of the separation between the filaments. Another possibility is synchrotron cooling instability occurring in cometary tails formed as a result of the interaction of cosmic ray driven Galactic centre outflow with obstacles such as stellar winds. In this picture, the mean spacing and the mean width of the filaments are expected to be a fraction of a parsec, consistent with observed spacing.

plasmas↗

The spatial distribution of gamma-ray bursts from Pioneer Venus Orbiter

The homogeneity of the burst source population in space with the log N(greater than P sub c)-log P sub c and log (dN/dP sub c)-log P sub c relations is studied on the basis of a sample of 225 observed gamma-ray bursts between September 1978 and July 1988 from the Pioneer Venus Orbiter. The relations N(greater than P sub c) is approximately equal to P sub c exp -1.5 and dN/dP sub c is approximately equal to P sub c exp -2.5 are expected if the gamma bursters are uniformly distributed in space. The Kolmogorov-Smirnov statistic gives a probability of 0.93 that the observed integral distribution N(greater than P sub c) is drawn from the expected uniform distribution P sub c exp -1.5. A least-squares 1 standard deviation fit to the observed differential distribution yields dN/dP sub c of about P sub c -2.42 +/- 0.07. Both are consistent with a uniform distribution in space of the parent population of burst sources.

Chuang, K. W.↗

Sampling errors in the estimation of empirical orthogonal functions

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

North, G. R.↗