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 793 records · Page 44

Self-organization of cosmic radiation pressure instability

Under some circumstances the absorption of radiation momentum by an absorbing medium opens the possibility of a dynamical instability, sometimes called 'mock gravity'. Here, a simplified abstract model is studied in which the radiation source is assumed to remain spatially uniform, there is no reabsorption or reradiated light, and no forces other than radiative pressure act on the absorbing medium. It is shown that this model displays the unique feature of being not only unstable, but also self-organizing. The structure approaches a statistical dynamical steady state which is almost independent of initial conditions. In this saturated state the absorbers are concentrated in thin walls around empty bubbles; as the instability develops the big bubbles get bigger and the small ones get crushed and disappear. A linear analysis shows that to first order the thin walls are indeed stable structures. It is speculated that this instability may play a role in forming cosmic large-scale structure.

Hogan, Craig J.↗

Validation of Satellite-Based Precipitation Products from TRMM to GPM

The global precipitation measurement mission (GPM) has been in operation for seven years and continues to provide a vast quantity of global precipitation data at finer temporospatial resolutions with improved accuracy and coverage. GPM’s signature algorithm, the integrated multisatellite retrievals for GPM (IMERG) is a next-generation of precipitation product expected for wide variety of research and operational applications. This study evaluates the latest version (V06B) of IMERG and its predecessor, the tropical rainfall measuring mission (TRMM) multisatellite precipitation (TMPA) 3B42 (V7) using ground-based and gauge-corrected multiradar multisensor system (MRMS) precipitation products over the conterminous United States (CONUS). The spatial distributions of all products are analyzed. The error characteristics are further examined for 3B42 and IMERG in winter and summer by an error decomposition approach, which partitions total bias into hit bias, biases due to missed precipitation and false precipitation. The volumetric and categorical statistical metrics are used to quantitatively evaluate the performance of the two satellite-based products. All products show a similar precipitation climatology with some regional differences. The two satellite-based products perform better in the eastern CONUS than in the mountainous Western CONUS. The evaluation demonstrates the clear improvement in IMERG precipitation product in comparison with its predecessor 3B42, especially in reducing missed precipitation in winter and summer, and hit bias in winter, resulting in better performance in capturing lighter and heavier precipitation.

MRMS↗

Estimating pixel-level uncertainty in ocean color retrievals from MODIS

The spectral distribution of marine remote sensing reflectance, R(rs), is the fundamental measurement of ocean color science, from which a host of bio-optical and biogeochemical properties of the water column can be derived. Estimation of uncertainty in these derived properties is thus dependent on knowledge of the uncertainty in satellite-retrieved R(rs) (u(c)(R(rs))) at each pixel. Uncertainty in R(rs), in turn, is dependent on the propagation of various uncertainty sources through the R(rs) retrieval process, namely the atmospheric correction (AC). A derivative-based method for uncertainty propagation is established here to calculate the pixel-level uncertainty in R(rs), as retrieved using NASA’s multiple-scattering epsilon (MSEPS) AC algorithm and verified using Monte Carlo (MC) analysis. The approach is then applied to measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite, with uncertainty sources including instrument random noise, instrument systematic uncertainty, and forward model uncertainty. The uc(Rrs) is verified by comparison with statistical analysis of coincident retrievals from MODIS and in situ Rrs measurements, and our approach performs well in most cases. Based on analysis of an example 8-day global products, we also show that relative uncertainty in R(rs) at blue bands has a similar spatial pattern to the derived concentration of the phytoplankton pigment chlorophyll-a (chl-a), and around 7.3%, 17.0%, and 35.2% of all clear water pixels (chl-a ≤ 0.1 mg/cu.m) with valid u(c)(R(rs)) have a relative uncertainty ≤ 5% at bands 412 nm, 443 nm, and 488 nm respectively, which is a common goal of ocean color retrievals for clear waters. While the analysis shows that u(c)(R(rs)) calculated from our derivative-based method is reasonable, some issues need further investigation, including improved knowledge of forward model uncertainty and systematic uncertainty in instrument calibration.

Pixel-level uncertainty↗

Use of a diffuser for producing texture masks for image processing

The use of a diffuser in synthesizing binary texture-variance masks for applications in image processing is investigated. The binary texture-variance masks are used to separate various objects in an image on the basis of their spatial texture variance. The diffuser is required to reduce the speckling present in the mask. The reduction in speckling is described quantitatively as a reduction of the misclassification probability associated with the separation of the objects into several texture classes. General results are furnished for the statistics of the field scattered by the diffuser, and these are used to obtain expressions for the misclassification probabilities in the presence of a diffuser.

Stark, H.↗

Spatial power spectra of the crustal geomagnetic field and core geomagnetic field

Equations providing numerical values of the geomagnetic field spherical harmonic spatial power spectrum as defined by Lowes (1966, 1974) are obtained and this power spectrum is related to various other power spectra. Equations relating the spherical harmonic spatial power spectrum to average great circle power spectra for components of the vector magnetic field in the radial direction, along the great circle track and perpendicular to the first two directions are derived under the assumption that the sources of the field are internal. A statistical model for the crustal and core geomagnetic fields is proposed and used to derive equations for the expected main and crustal spherical harmonic power spectra. The model equations are then compared with observations to determine a scale factor which is then used to obtain an estimate for the core radius and a great circle power spectrum for the field component perpendicular to the great circle and radial directions which are in good agreement with observations. The predicted spherical harmonic power spectrum for the crustal field is found to be consistent with POGO satellite and aircraft data. Other possible models for the crustal and core geomagnetic fields are also briefly considered.

Mcleod, M. G.↗

Star formation in active galaxies and quasars

The observational evidence for a causal or statistical link between star formation and active galactic nuclei are reviewed. The chief difficulty is in quantitatively ascertaining the star formation rate in active galaxies: most of the readily observable manifestations of star formation superficially resemble those of an active nucleus. Careful multi-wavelength spatially-resolved observations demonstrate that many Seyfert galaxies are undergoing star formation. A survey of CO emission from Seyferts (interpreted in conjunction with IRAS data) suggests that type 2 Seyferts have unusually high rates of star formation, but type 1 Seyferts do not. Recent work also suggests that many powerful radio galaxies may be actively forming stars: radio galaxies with strong emission-lines often have blue colors and strong far-infrared emission. Determining the star formation rate in the host galaxies of quasars is especially difficult. Multicolor imaging and long-slit spectroscopy suggests that many of the host galaxies of radio-loud quasars are blue, and a cold interstellar medium has been detected in some quasar hosts.

Heckman, Timothy M.↗

A New Fitting Procedure for the Blended He I 1083 nm Multiplet

The He I 1083 nm multiplet is a powerful tool for observing the outer solar atmosphere but is difficult to analyze because the lines are weak, highly variable, and spectrally blended, both internally and with other neighboring solar and telluric lines. After separation from nearby spectral features, two components of the He I multiplet are resolved. Fitting these lines with two unconstrained Gaussian profiles always gives a ratio of major to minor component of less than half the value which would be expected for optically thin lines. One possibility for explaining the discrepancy between the weakness of the lines and the ratio of the spectral components is to assume that the line formation regions are concentrated in laterally unresolved, optically thick structures with small filling factor. However, we present here a least-squares fitting technique using cubic splines with fixed breakpoints with the constraint that the blend is the sum of three identically shaped profiles shifted in wavelength according to the atomic structure of the blend and weighted by the corresponding statistics weights, in agreement with optically thin line formation. The basis functions for the fitting procedure have no built-in spectral symmetry or shape. The resulting underlying profiles tend to be asymmetric with excess absorption to the blue, consistent with formation by "hot" and "cool" spatial elements within the observed volume, with the hotter regions having differential motion toward the observer. The results agree well with NASA/XSO Spectromagnetograph observations in quiet sun and coronal holes. Partial funding of this research was provided through the NASA Sun-Earth Connection SR&T program.

Jones, H. P.↗

Small-Scale Variability of Large Cloud Drops

Cloud droplet size distribution is one of the most fundamental subjects in cloud physics. Understanding of spatial distribution and small-scale fluctuations of cloud droplets is essential for both cloud physics and atmospheric radiation. For cloud physics, it relates to the coalescence growth of raindrops while for radiation, it has a strong impact on a cloud's radiative properties. Most of the existing cloud radiation and precipitation formation models assume that the mean number of drops with a given radius varies proportionally to volume. The analysis of microphysical data on liquid water drop sizes shows that, for sufficiently small volumes, the number is proportional to the drop size dependent power of the volume. For abundant small drops present, the exponent is 1 as assumed in the conventional approach. However, for rarer large drops, the exponents fall below unity. At small scales, therefore, the mean number of large drops decreases with volume at a slower rate than the conventional approach assumes, suggesting more large drops at these scales than conventional models account for; their impact is consequently underestimated. Size dependent models of spatial distribution of cloud drops that simulate the observed power laws show strong drop clustering, the more so the larger the drops. The degree of clustering is determined by the observed exponents. The strong clustering of large drops arises naturally from the observed power-law statistics. Current theories of photon-cloud interaction and warm rain formation will need radical revision in order to produce these statistics; their underlying equations are unable to yield the observed power law.

Marshak, Alexander↗

Global Patterns of Lightning Properties Derived by OTD and LIS

The satellite instruments Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) provide unique empirical data about the frequency of lightning flashes around the globe (OTD), and the tropics (LIS), which 5 has been used before to compile a well received global climatology of flash rate densities. Here we present a statistical analysis of various additional lightning properties derived from OTD/LIS, i.e. the number of so-called "events" and "groups" per flash, as well as 10 the mean flash duration, footprint and radiance. These normalized quantities, which can be associated with the flash "strength", show consistent spatial patterns; most strikingly, oceanic flashes show higher values than continental flashes for all properties. Over land, regions with high (Eastern US) 15 and low (India) flash strength can be clearly identified. We discuss possible causes and implications of the observed regional differences. Although a direct quantitative interpretation of the investigated flash properties is difficult, the observed spatial patterns provide valuable information for the 20 interpretation and application of climatological flash rates. Due to the systematic regional variations of physical flash characteristics, viewing conditions, and/or measurement sensitivities, parametrisations of lightning NOx based on total flash rate densities alone are probably affected by regional biases.

Beirle, Steffen↗

Daily Kilometer-Scale MODIS Satellite Maps of PM2.5 Describe Wintertime Episodes

The San Joaquin Valley (SJV) suffers from severe health-endangering episodes of PM2.5 aerosol loadings in wintertime; episodes last approximately 5 days and differ in geographical distribution and composition. PM2.5 stations are scattered; consequently the use of remote sensing to map variable regional patterns of these varying respirable aerosol concentrations is desirable. High-precision AOT retrievals can capture column particulate loading. However,PM2.5 mapping is challenging due to several reasons: particularly thin mixed layers (ML) and thus relatively low aerosol optical thickness (AOT) close to current measurement limits, variable and a typical composition of the aerosols, and complex surface bidirectional reflectance. However, the West does present some advantages in analysis. Air basins are isolated from long-distance transport, and experience predominant strong meteorological subsidence. Thus these Western basin regions have fewer problematic cases of overriding aerosol layers detached from the surface. To counter such local overriding, Chu et al. have described an approach for the Eastern US, and He et al have described a synoptic classification approach useful in Shanghai. The Bay Area Air Quality Management District (BAAQMD) expands our experience with the use of AOT, with lower PM2.5 and several isolated sub-basins. We have prepared daily maps of episodes in each region. We present also a sequence of increasingly detailed statistical models, AOT initially appears to contribute little information; however, inclusion of weather information reveals its utility. Lyapustin and Wang's MultiAngle Implementation of Atmospheric Correction (MAIAC) retrieval for AOT provided the most useful operational remote sensing information for these regions. It provides high (1-km) spatial resolution maps and a high percentage of availability. Empirical regression methods have found that random effects regression models (aka mixed effects models, ME) employing AOT provide good estimates of ground PM2.5 concentrations.Here, we attempt to extend these methods and evaluate the usefulness of AOT with greater physical analysis, based on DISCOVER-AQ4 experience.

Chatfield, Robert B.↗

A New Parametrization of the Rain Drop Size Distribution

This paper revisits the problem of finding a parametric form for the rain drop size distribution (DSD) which 1) is an appropriate model for tropical rainfall, and 2) involves statistically independent parameters. Using TOGA/COARE data, we derive a parametrization which meets these criteria. This new parametrization is an improvement on the one that was derived in [3] using TRMM ground truth data from Darwin, Australia. The new COARE data allows us to verify that the spatial variability fo the two

rain↗

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗

Visualizing Spatially Varying Distribution Data

Box plot is a compact representation that encodes the minimum, maximum, mean, median, and quarters information of a distribution. In practice, a single box plot is drawn for each variable of interest. With the advent of more accessible computing power, we are now facing the problem of visual icing data where there is a distribution at each 2D spatial location. Simply extending the box plot technique to distributions over 2D domain is not straightforward. One challenge is reducing the visual clutter if a box plot is drawn over each grid location in the 2D domain. This paper presents and discusses two general approaches, using parametric statistics and shape descriptors, to present 2D distribution data sets. Both approaches provide additional insights compared to the traditional box plot technique

Kao, David↗

Field Penetration in a Rectangular Box Using Numerical Techniques: An Effort to Obtain Statistical Shielding Effectiveness

This paper emphasizes the application of numerical methods to explore the ideas related to shielding effectiveness from a statistical view. An empty rectangular box is examined using a hybrid modal/moment method. The basic computational method is presented followed by the results for single- and multiple observation points within the over-moded empty structure. The statistics of the field are obtained by using frequency stirring, borrowed from the ideas connected with reverberation chamber techniques, and extends the ideas of shielding effectiveness well into the multiple resonance regions. The study presented in this paper will address the average shielding effectiveness over a broad spatial sample within the enclosure as the frequency is varied.

Bunting, Charles F.↗

Backscatter Characteristics of the Winter Sea Ice Cover in the Beaufort

The microwave backscatter statistics of the sea ice cover in the winter Beaufort Sea were characterized using C-band Synthetic Aperture Radar (SAR) data collected by ERS-1. Sea ice backscatter signatures were sampled from SAR image data collected during the winter of 1991/92 and winter of 1992/93. The spatial and temporal variabilities of the backscatter signatures of different ice types are discussed. The results show considerable seasonal stability of the backscatter signature of multi-year ice as well as first-year ice. Small amplitude regional variations of the multi-year ice backscatter can be observed. Consistent contrast between multiyear ice and first-year ice is maintained throughout the season.

Kwok, R.↗

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

Twentieth Century Redistribution in Climatic Drivers of Global Tree Growth

Energy and water limitations of tree growth remain insufficiently understood at large spatiotemporal scales, hindering model representation of interannual or longer-term ecosystem processes. By assessing and statistically scaling the climatic drivers from 2710 tree-ring sites, we identified the boreal and temperate land areas where tree growth during 1930–1960 CE responded positively to temperature (20.8 ± 3.7 Mio km2; 25.9 ± 4.6%), precipitation (77.5 ± 3.3 Mio km2; 96.4 ± 4.1%), and other parameters. The spatial manifestation of this climate response is determined by latitudinal and altitudinal temperature gradients, indicating that warming leads to geographic shifts in growth limitations. We observed a significant (P < 0.001) decrease in temperature response at cold-dry sites between 1930–1960 and 1960–1990 CE, and the total temperature-limited area shrunk by −8.7 ± 0.6 Mio km2. Simultaneously, trees became more limited by atmospheric water demand almost worldwide. These changes occurred under mild warming, and we expect that continued climate change will trigger a major redistribution in growth responses to climate.

Babst, Flurin↗

Synergistic Use of Hyperspectral UV-Visible OMI and Broadband Meteorological Imager MODIS Data for a Merged Aerosol Product

The retrieval of optimal aerosol datasets by the synergistic use of hyperspectral ultraviolet(UV)–visible and broadband meteorological imager (MI) techniques was investigated. The Aura Ozone Monitoring Instrument (OMI) Level 1B (L1B) was used as a proxy for hyperspectral UV–visible instrument data to which the Geostationary Environment Monitoring Spectrometer (GEMS) aerosol algorithm was applied. Moderate-Resolution Imaging Spectroradiometer (MODIS) L1B and dark target aerosol Level 2 (L2) data were used with a broadband MI to take advantage of the consistent time gap between the MODIS and the OMI. First, the use of cloud mask information from the MI infrared (IR) channel was tested for synergy. High-spatial-resolution and IR channels of the MI helped mask cirrus and sub-pixel cloud contamination of GEMS aerosol, as clearly seen in aerosol optical depth (AOD) validation with Aerosol Robotic Network (AERONET) data. Second, dust aerosols were distinguished in the GEMS aerosol-type classification algorithm by calculating the total dust confidence index (TDCI) from MODIS L1B IR channels. Statistical analysis indicates that the Probability of Correct Detection (POCD) between the forward and inversion aerosol dust models (DS) was increased from 72% to 94% by use of the TDCI for GEMS aerosol-type classification, and updated aerosol types were then applied to the GEMS algorithm. Use of the TDCI for DS type classification in the GEMS retrieval procedure gave improved single-scattering albedo (SSA) values for absorbing fine pollution particles (BC) and DS aerosols. Aerosol layer height (ALH) retrieved from GEMS was compared with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data, which provides high-resolution vertical aerosol profile information. The CALIOP ALH was calculated from total attenuated backscatter data at 1064 nm, which is identical to the definition of GEMS ALH. Application of the TDCI value reduced the median bias of GEMS ALH data slightly. The GEMS ALH bias approximates zero, especially for GEMS AOD values of>~0.4 and GEMS SSA values of<~0.95.Finally, the AOD products from the GEMS algorithm and MI were used in aerosol merging with the maximum-likelihood estimation method, based on a weighting factor derived from the standard deviation of the original AOD products. With the advantage of the UV–visible channel in retrieving aerosol properties over bright surfaces, the combined AOD products demonstrated better spatial data availability than the original AOD products, with comparable accuracy. Furthermore, pixel-level error analysis of GEMS AOD data indicates improvement through MI synergy.

aerosol↗