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 703 records · Page 39

Sensitivity of Lower Stratospheric Assimilated Ozone on Error Covariance Modeling and Data Selection

Assimilated ozone is produced at the NASA/Goddard Data Assimilation Office by blending ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line transport model. The current system tends to overestimate the amount of lower stratospheric ozone. This is a region where ozone plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. Since these are all important applications of assimilated ozone products, effort is being directed towards reducing this bias. The SBUV ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone in the lower stratosphere, and its vertical distribution in particular, are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Accordingly, the sensitivity of the assimilated lower stratospheric ozone fields to changes in the TOMS error-covariance modeling and the SBUV data selection has been investigated. The use of a spatially correlated TOMS error covariance model led to improvements in the product. However, withholding the SBUV/2 data for the layer between 63 and 126 hPa typically degraded the product, a result which vindicates the use of this layer ozone product, despite its known errors. These efforts to improve the lower stratospheric distribution will be extended to include a more advanced forecast error covariance model, and by assimilating ozone products from new instruments on Envisat and EOS Aura.

Stajner, Ivanka↗

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗

Large Contribution of Meteorological Factors to Inter-Decadal Changes in Regional Aerosol Optical Depth

Aerosol optical depth (AOD) has become a crucial metric for assessing global climate change. Although global and regional AOD trends have been studied extensively, it remains unclear what factors are driving the inter-decadal variations in regional AOD and how to quantify the relative contribution of each dominant factor. This study used a long-term (1980–2016) aerosol dataset from the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) reanalysis, along with two satellite-based AOD datasets (MODIS/Terra and MISR) from 2001 to 2016, to investigate the long-term trends in global and regional aerosol loading. Statistical models based on emission factors and meteorological parameters were developed to identify the main factors driving the inter-decadal changes of regional AOD and to quantify their contribution. Evaluation of the MERRA-2 AOD with the ground-based measurements of AERONET indicated significant spatial agreement on the global scale (r= 0.85, root-mean-square error = 0.12, mean fractional error = 38.7 %, fractional gross error = 9.86 % and index of agreement = 0.94). However, when AOD observations from the China Aerosol Remote Sensing Network (CARSNET) were employed for independent verification, the results showed that MERRA-2 AODs generally underestimated CARSNET AODs in China (relative mean bias = 0.72 and fractional gross error =−34.3 %). In general, MERRA-2 was able to quantitatively reproduce the annual and seasonal AOD trends on both regional and global scales, as observed by MODIS/Terra, although some differences were found when compared to MISR. Over the 37-year period in this study, significant decreasing trends were observed over Europe and the eastern United States. In contrast, eastern China and southern Asia showed AOD increases, but the increasing trend of the former reversed sharply in the most recent decade. The statistical analyses suggested that the meteorological parameters explained a larger proportion of the AOD variability (20.4 %–72.8 %) over almost all regions of interest (ROIs) during 1980–2014 when compared with emission factors (0 %–56 %). Further analysis also showed that SO2 was the dominant emission factor, explaining 12.7 %–32.6 % of the variation in AOD over anthropogenic-aerosol-dominant regions, while black carbon or organic carbon was the leading factor over the biomass-burning-dominant (BBD) regions, contributing 24.0 %–27.7 % of the variation. Additionally, wind speed was found to be the leading meteorological parameter, explaining 11.8 %–30.3 % of the variance over the mineral-dust-dominant regions, while ambient humidity (including soil moisture and relative humidity) was the top meteorological parameter over the BBD regions, accounting for 11.7 %–35.5 % of the variation. The results of this study indicate that the variation in meteorological parameters is a key factor in determining the inter-decadal change in regional AOD.

Huizheng Che↗

Multi-Scale Fractal Analysis of Image Texture and Pattern

Fractals embody important ideas of self-similarity, in which the spatial behavior or appearance of a system is largely independent of scale. Self-similarity is defined as a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. An ideal fractal (or monofractal) curve or surface has a constant dimension over all scales, although it may not be an integer value. This is in contrast to Euclidean or topological dimensions, where discrete one, two, and three dimensions describe curves, planes, and volumes. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution. However, most geographical phenomena are not strictly self-similar at all scales, but they can often be modeled by a stochastic fractal in which the scaling and self-similarity properties of the fractal have inexact patterns that can be described by statistics. Stochastic fractal sets relax the monofractal self-similarity assumption and measure many scales and resolutions in order to represent the varying form of a phenomenon as a function of local variables across space. In image interpretation, pattern is defined as the overall spatial form of related features, and the repetition of certain forms is a characteristic pattern found in many cultural objects and some natural features. Texture is the visual impression of coarseness or smoothness caused by the variability or uniformity of image tone or color. A potential use of fractals concerns the analysis of image texture. In these situations it is commonly observed that the degree of roughness or inexactness in an image or surface is a function of scale and not of experimental technique. The fractal dimension of remote sensing data could yield quantitative insight on the spatial complexity and information content contained within these data. A software package known as the Image Characterization and Modeling System (ICAMS) was used to explore how fractal dimension is related to surface texture and pattern. The ICAMS software was verified using simulated images of ideal fractal surfaces with specified dimensions. The fractal dimension for areas of homogeneous land cover in the vicinity of Huntsville, Alabama was measured to investigate the relationship between texture and resolution for different land covers.

Emerson, Charles W.↗

Satellite disintegration dynamics

The subject of satellite disintegration is examined in detail. Elements of the orbits of individual fragments, determined by DOD space surveillance systems, are used to accurately predict the time and place of fragmentation. Dual time independent and time dependent analyses are performed for simulated and real breakups. Methods of statistical mechanics are used to study the evolution of the fragment clouds. The fragments are treated as an ensemble of non-interacting particles. A solution of Liouville's equation is obtained which enables the spatial density to be calculated as a function of position, time and initial velocity distribution.

Dasenbrock, R. R.↗

The velocity structure and turbulence at the center of the Orion Nebula

Using the Coude Feed System of the Kitt Peak National Observatory at its highest possible dispersion and a CCD detector with wide dynamic range it was possible to examine in detail the radial and turbulent velocities of the Orion Nebula at the seeing limit. For the first time, multiple components of the spectral lines were resolved and the Gaussian profiles fit to them. Using four different observing runs, two circular areas 2 minutes in diameter were observed, the first centered on the theta sup 1 Orionis and the second centered on theta sup 2 Orionis. The velocity resolution of 1 km/s combined with the angular resolution of 2 seconds allowed observation in detail at line splitting as well as the location and measurement of high internal motions, about -50 km/s with respect to the bulk of ionized matter. Because over 900 different positions were measured, modern statistical techniques were used for the analysis of turbulence, so that the dispersion relation and the structure function could be calculated and compared with predictions of the Kolmogorov theory. The random motions showed correlation with spatial separation, as expected for turbulence; but comparison with predictions of the classical theory for nebular turbulence showed poor agreement. If Kolmogorov theory applies, then the energy must be input at many scales.

Castaneda, Hector O.↗

The systematic radial downflow in the transition region of the quiet sun from limb-to-limb observations of the C IV resonance lines

This paper presents absolute velocities of C IV 1548, 1550 nm measured as a function of position along the solar equator, which was free of both active regions and coronal holes, and uniformly representative of the quiet sun. These observations were made with moderate spatial resolution (18 arcsec) using an EUV spectrometer dedicated to measuring absolute wavelengths (velocities) by direct comparison with a platinum spectrum generated on board the sounding rocket. On the assumption that systematic horizontal motions cancel statistically so that the line-of-sight velocities approach zero at the limb, a net radial downflow of 7.5 + or - 1.0 km/sec was found. The assumption was tested using the wavelength reference and found to be valid within the absolute accuracy of the rest wavelengths of the C IV lines.

Rottman, Gary J.↗

Modeling of the reactant conversion rate in a turbulent shear flow

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Frankel, S. H.↗

Numerical simulation of two-dimensional spatially-developing mixing layers

Two-dimensional, incompressible, spatially developing mixing layer simulations are performed at Re = 10(exp 2) and 10(exp 4) with two classes of perturbations applied at the inlet boundary; combinations of discrete modes from linear stability theory, and a broad spectrum of modes derived from experimentally measured velocity spectra. The effect of the type and strength of inlet perturbations on vortex dynamics and time-averaged properties are explored. Two-point spatial velocity and autocorrelations are used to estimate the size and lifetime of the resulting coherent structures and to explore possible feedback effects. The computed time-averaged properties such as mean velocity profiles, turbulent statistics, and spread rates show good agreement with experimentally measured values. It is shown that by forcing with a broad spectrum of modes derived from an experimental energy spectrum many experimentally observed phenomena can be reproduced by a 2-D simulation. The strength of the forcing merely affected the length required for the dominant coherent structures to become fully-developed. Thus intensities comparable to those of the background turbulence in many wind tunnel experiments produced the same results, given sufficient simulation length.

Wilson, R. V.↗

Observational Review of Gamma-Ray Bursts

The observed gamma-ray temporal, spectral, intensity and spatial distribution characteristics of GRBs (Gamma Ray Bursts), primarily from data obtained from the Compton Observatory, will be described. The talk will concentrate on recent studies of burst properties, correlations of GRB parameters and other statistical studies that have only recently come to light with the unprecedented sample of over two thousand GRBs, along with some mention of studies in progress by members of the BATSE (Burst and Transient Source Experiment) team. Recent studies of possible observational biases, un-triggered GRBs and threshold calculations for BATSE will also be described.

Fishman, Gerald J.↗

Evaluation of the Global Multi-Resolution Terrain Elevation Data 2010 (GMTED2010) Using ICESat Geodetic Control

Supported by NASA's Earth Surface and Interior (ESI) Program, we are producing a global set of Ground Control Points (GCPs) derived from the Ice, Cloud and land Elevation Satellite (ICESat) altimetry data. From February of 2003, to October of 2009, ICESat obtained nearly global measurements of land topography (+/- 86deg latitudes) with unprecedented accuracy, sampling the Earth's surface at discrete approx.50 m diameter laser footprints spaced 170 m along the altimetry profiles. We apply stringent editing to select the highest quality elevations, and use these GCPs to characterize and quantify spatially varying elevation biases in Digital Elevation Models (DEMs). In this paper, we present an evaluation of the soon to be released Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010). Elevation biases and error statistics have been analyzed as a function of land cover and relief. The GMTED2010 products are a large improvement over previous sources of elevation data at comparable resolutions. RMSEs for all products and terrain conditions are below 7 m and typically are about 4 m. The GMTED2010 products are biased upward with respect to the ICESat GCPs on average by approximately 3 m.

Carabajal, Claudia C.↗

Air Quality Over the Canadian Oil Sands: A First Assessment Using Satellite Observations

Results from the first assessment of air quality over the Canadian oil sands -- one ofthe largest industrial undertakings in human history -- using satellite remote sensing observations of two pollutants, nitrogen dioxide (N0O) and sulfur dioxide (SO2), are presented. High-resolution maps were created that revealed distinct enhancements in both species over an area (roughly 30 km x 50 km) of intensive surface mining at scales of a few kilometers. The magnitude of these enhancements, quantified in terms of total mass, are comparable to the largest seen in Canada from individual sources. The rate of increase in NO2 between 2005 and 2010 was assessed at 10.4 +/- 3.5%/year and resulted from increases both in local values as well as the spatial extent of the enhancement. This is broadly consistent with both surface-measurement trends and increases in annual bitumen production. An increase in SO2 was also found, but given larger uncertainties, it is not statistically significant.

McLinden, C. A.↗

Statistical Perspectives on Stratospheric Transport

Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Sparling, L. C.↗

Statistical analysis of atmospheric turbulence about a simulated block building

An array of towers instrumented to measure the three components of wind speed was used to study atmospheric flow about a simulated block building. Two-point spacetime correlations of the longitudinal velocity component were computed along with two-point spatial correlations. These correlations are in good agreement with fundamental concepts of fluid mechanics. The two-point spatial correlations computed directly were compared with correlations predicted by Taylor's hypothesis and excellent agreement was obtained at the higher levels which were out of the building influence. The correlations fall off significantly in the building wake but recover beyond the wake to essentially the same values in the undisturbed, higher regions.

Steely, S. L., Jr.↗

Photon statistics of a two-mode squeezed vacuum

We investigate the general case of the photon distribution of a two-mode squeezed vacuum and show that the distribution of photons among the two modes depends on four parameters: two squeezing parameters, the relative phase between the two oscillators and their spatial orientation. The distribution of the total number of photons depends only on the two squeezing parameters. We derive analytical expressions and present pictures for both distributions.

Schrade, Guenter↗

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. For the recognition of terrain types, spatial signatures are developed from the diffraction patterns of small areas of ERTS-1 images. This knowledge is exploited for the measurements of a small number of meaningful spatial features from the digital Fourier transforms of ERTS-1 image cells containing 32 x 32 picture elements. Using these spatial features and a heuristic algorithm, the terrain types in the vicinity of Phoenix, Arizona were recognized by the computer with a high accuracy. Then, the spatial features were combined with spectral features and using the maximum likelihood criterion the recognition accuracy of terrain types increased substantially. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. Nonlinear transformations of the feature vectors are required so that the terrain class statistics become approximately Gaussian. 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.

Gramenopoulos, N.↗

Steady induction effects in geomagnetism. Part 1B: Geomagnetic estimation of steady surficial core motions: A non-linear inverse problem

The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation (SV) indicated by models of the observed geomagnetic field is examined in the source-free mantle/frozen-flux core (SFI/VFFC) approximation. This inverse problem is non-linear because solutions of the forward problem are deterministically chaotic. The SFM/FFC approximation is inexact, and neither the models nor the observations they represent are either complete or perfect. A method is developed for solving the non-linear inverse motional induction problem posed by the hypothesis of (piecewise, statistically) steady core surface flow and the supposition of a complete initial geomagnetic condition. The method features iterative solution of the weighted, linearized least-squares problem and admits optional biases favoring surficially geostrophic flow and/or spatially simple flow. Two types of weights are advanced radial field weights for fitting the evolution of the broad-scale portion of the radial field component near Earth's surface implied by the models, and generalized weights for fitting the evolution of the broad-scale portion of the scalar potential specified by the models.

Voorhies, Coerte V.↗

A Historic View of Solar Coronal Mass Ejections (CMEs)

We present a historic overview of CME observations, ending with concepts for future measurement capabilities. One of the first detections of what we now call a CME was provided by instrumentation on OSO-7 and reported by Tousey (1973); but the phrase "corona) mass ejection" was coined after the Skylab/ATM coronagraph detected dozens of the transients over its nine month observing run (e.g., Munro et al., 1979). Pre-discovery identification of likely CMEs were then reported in historic eclipse photographs and drawings (e.g., Eddy, 1974; Cliver, 1989). Multi-year observations followed with groundbased MLSO MK3/4 coronagraph (1980-present), and spacebased missions: Solwind (1979-1985), SMM (1980-1989), SOHO LASCO/EIT (1996-present), SMEI (2003-present), and STEREO SECCHI (2006-present). The Spartan 201 coronagraph flew in space multiple times. The influential Gosling (1993) "solar flare myth" manuscript identified CMEs as the cause of the most severe geomagnetic storms, thus cementing their importance in Sun-Earth connection studies. A new window into CMEs was opened with the launch of SOHO in 1995 when the UVCS spectrometer began returning plasma diagnostics of a significant number of events (e.g., Ciaravella et al., 2006). What about the future for CME research? Statistical properties of the UVCS CME observations are forthcoming; the STEREO mission should continue to return views of CMEs from unique vantage points; and the recent launch of SDO should provide new insights into the small spatial scale dynamics of activity associated with CMEs. Several new observing techniques have been demonstrated at total eclipses, and inclusion on spacebased platforms in the future could also prove valuable for understanding CMEs. A common element of several recent proposals is to image the white-light corona with extremely high spatial resolution. The momentary glimpses of the corona during total solar eclipses have shown fine structure that is not captured by global models, and dynamics of these structured elements may be important to resolve the question of CME initiation.

SaintCyr, Orville C.↗