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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Review of Bayesian Networks for Spatial Data

We report Bayesian networks are a popular class of multivariate probabilistic models as they allow for the translation of prior beliefs about conditional dependencies between variables to be easily encoded into their model structure. Due to their widespread usage, they are often applied to spatial data for inferring properties of the systems under study and also generating predictions for how these systems may behave in the future. We review published research on methodologies for representing spatial data with Bayesian networks and also summarize the application areas for which Bayesian networks are employed in the modeling of spatial data. We find that a wide variety of perspectives are taken, including a GIS-centric focus on efficiently generating geospatial predictions, a statistical focus on rigorously constructing graphical models controlling for spatial correlation, as well as a range of problem-specific heuristics for mitigating the effects of spatial correlation and dependency arising in spatial data analysis. Special attention is also paid to potential future directions for integration of Bayesian networks with spatial processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A species’ response to spatial climatic variation does not predict its response to climate change

The dominant paradigm for assessing ecological responses to climate change assumes that future states of individuals and populations can be predicted by current, species-wide performance variation across spatial climatic gradients. However, if the fates of ecological systems are better predicted by past responses to in situ climatic variation through time, this current analytical paradigm may be severely misleading. Empirically testing whether spatial or temporal climate responses better predict how species respond to climate change has been elusive, largely due to restrictive data requirements. Here, we leverage a newly collected network of ponderosa pine tree-ring time series to test whether statistically inferred responses to spatial versus temporal climatic variation better predict how trees have responded to recent climate change. When compared to observed tree growth responses to climate change since 1980, predictions derived from spatial climatic variation were wrong in both magnitude and direction. This was not the case for predictions derived from climatic variation through time, which were able to replicate observed responses well. Future climate scenarios through the end of the 21st century exacerbated these disparities. These results suggest that the currently dominant paradigm of forecasting the ecological impacts of climate change based on spatial climatic variation may be severely misleading over decadal to centennial timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An evaluation of air quality in major urban areas of India

Rapid economic growth and burgeoning population have contributed to enhanced levels of PM 2.5 concentrations in urban regions of India. Evaluation of ambient air quality facilitates the assessment of effectiveness of emission control measures and early identification of new sources. This study provides a comprehensive statistical analysis of PM 2.5 concentrations in key urban areas across India, including Delhi, Kolkata, Mumbai, Chennai, Hyderabad, and several regional centers. Data from 2017 to 2023 was analyzed using trend analysis, cluster analysis, principal component analysis, and geostatistical interpolation to understand spatiotemporal variations and sources. The analysis reveals significant differences in spatial distribution of PM 2.5 concentrations with high annual averages in urban regions in Indo-Gangetic plain (82–123 μg m −3 ) and relatively lower concentrations (29–46 μg m −3 ) in southern urban areas of Kerala, Tamil Nadu and Andhra Pradesh. Delhi state had the highest 24-averaged PM 2.5 concentrations (112 μg m −3 ) followed by urban regions in Uttar Pradesh, Bihar and West Bengal (94 μg m −3 ). Trend analysis from 2017 to 2023 revealed an overall 2.5% decline in site-wide PM2.5 concentrations, with the exception of Ludhiana, which exhibited a consistent annual increase of 10%. Principal component analysis (PCA) attributes 30% of the variance to wintertime emissions, 13% to biomass burning, and 18% to the regional haze in the northern Indo-Gangetic Plain. Different analyses clearly demonstrates the contribution of biomass burning to pollution in Delhi and surrounding cities. Transboundary pollution to Kolkata is likely from the highly polluted region in Indo-Gangetic Plain. Coastal cities of Mumbai and Chennai has relatively lower pollution attributed to the influence of sea breeze dilution, with mostly local contribution and some potential transport from upwind industry clusters. Hyderabad also has local contribution due to high density of vehicular traffic and local small industries. This study shows that mitigation efforts targeting clusters of regions should be undertaken to curb the high PM2.5 pollution. Policy measures should be implemented both at local and the intra-state level to address shared sources and transport of pollution.

Hysplitbacktrajectories↗
NASA NTRS2002

On the Choice of Variable for Atmospheric Moisture Analysis

The implications of using different control variables for the analysis of moisture observations in a global atmospheric data assimilation system are investigated. A moisture analysis based on either mixing ratio or specific humidity is prone to large extrapolation errors, due to the high variability in space and time of these parameters and to the difficulties in modeling their error covariances. Using the logarithm of specific humidity does not alleviate these problems, and has the further disadvantage that very dry background estimates cannot be effectively corrected by observations. Relative humidity is a better choice from a statistical point of view, because this field is spatially and temporally more coherent and error statistics are therefore easier to obtain. If, however, the analysis is designed to preserve relative humidity in the absence of moisture observations, then the analyzed specific humidity field depends entirely on analyzed temperature changes. If the model has a cool bias in the stratosphere this will lead to an unstable accumulation of excess moisture there. A pseudo-relative humidity can be defined by scaling the mixing ratio by the background saturation mixing ratio. A univariate pseudo-relative humidity analysis will preserve the specific humidity field in the absence of moisture observations. A pseudorelative humidity analysis is shown to be equivalent to a mixing ratio analysis with flow-dependent covariances. In the presence of multivariate (temperature-moisture) observations it produces analyzed relative humidity values that are nearly identical to those produced by a relative humidity analysis. Based on a time series analysis of radiosonde observed-minus-background differences it appears to be more justifiable to neglect specific humidity-temperature correlations (in a univariate pseudo-relative humidity analysis) than to neglect relative humidity-temperature correlations (in a univariate relative humidity analysis). A pseudo-relative humidity analysis is easily implemented in an existing moisture analysis system, by simply scaling observed-minus background moisture residuals prior to solving the analysis equation, and rescaling the analyzed increments afterward.

Dee, Dick P.↗
NASA NTRS1993

Threshold effects on V/V(max) for gamma-ray bursts

The interpretation of the observed gamma-ray burst V/V(max) statistic in terms of spatial distributions is model-dependent. Detection of gamma-ray bursts requires the counting rate in one or more detectors to exceed a threshold C(lim) determined from a time-dependent background rate B(t). The sampling depth of the burst detector is thus time-dependent, and, if burst sources are nonuniform in space, the observed V/V(max) distribution will be affected by B(t). We demonstrate this effect with a simple geometric distribution of standard candles and argue that V/V(max) statistic without information on threshold variations is insufficient for rigorous data analysis. Peak count rates and threshold values must be given separately for all events in order to facilitate a meaningful comparison of observations with theoretical distribution models.

Hartmann, D. H.↗
NASA NTRS2022

Canadian and Alaskan Wildfire Smoke Particle Properties, Their Evolution and Controlling Factors, From Satellite Observations

The optical and chemical properties of biomass burning (BB) smoke particles greatly affect the impact that wildfires have on climate and air quality. Previous work has demonstrated some links between smoke properties and factors such as fuel type and meteorology. However, the factors controlling BB particle speciation at emission are not adequately understood nor are the factors driving particle aging during atmospheric transport. As such, modeling wildfire smoke impacts on climate and air quality remains challenging. The potential to provide robust, statistical characterizations of BB particles based on ecosystem type and ambient environmental conditions with remote sensing data is investigated here. Space-based Multi-angle Imaging SpectroRadiometer (MISR) observations, combined with the MISR Research Aerosol (RA) algorithm and the MISR Interactive Explorer (MINX) tool, are used to retrieve smoke plume aerosol optical depth (AOD) and to provide constraints on plume vertical extent; smoke age; and particle size, shape, light-absorption properties, and absorption spectral dependence. These tools are applied to numerous wildfire plumes in Canada and Alaska, across a range of conditions, to create a regional inventory of BB particle-type temporal and spatial distribution. We then statistically compare these results with satellite measurements of fire radiative power (FRP) and land cover characteristics, as well as short-term climate, meteorological, and drought information from the Modern-Era Retrospective analysis for Research and Applications (MERRA-2) reanalysis and the North American Drought Monitor. We find statistically significant differences in the retrieved smoke properties based on land cover type, with fires in forests producing the thickest plumes containing the largest, brightest particles and fires in savannas and grasslands exhibiting the opposite. Additionally, the inferred dominant aging mechanisms and the timescales over which they occur vary systematically between land types. This work demonstrates the potential of remote sensing to constrain BB particle properties and the mechanisms governing their evolution over entire ecosystems. It also begins to realize this potential, as a means of improving regional and global climate and air quality modeling in a rapidly changing world.

Katherine T. Junghenn Noyes↗
NASA NTRS1980

Interfaces between statistical analysis packages and the ESRI geographic information system

Interfaces between ESRI's geographic information system (GIS) data files and real valued data files written to facilitate statistical analysis and display of spatially referenced multivariable data are described. An example of data analysis which utilized the GIS and the statistical analysis system is presented to illustrate the utility of combining the analytic capability of a statistical package with the data management and display features of the GIS.

Masuoka, E.↗
DOE OSTI2023

A Statistical Framework for Evaluating Rain Microphysics in Model Simulations and Disdrometer Observations

Abstract Statistical analyses of a large disdrometer data set and a diverse set of model simulations for convection using the Regional Atmospheric Modeling System were conducted, with the mutual goal of providing insights into precipitation formation and microphysical processes. We demonstrate that a two‐moment bulk microphysical model successfully captures the dominant observed modes of variability in rainfall related to rainfall intensity and raindrop size distributions. The model reproduced the general distribution of observed precipitation groups (PGs) derived from Principal Component Analysis. The multi‐variable analysis also uncovered some shortcomings in the model as well as limitations of the disdrometer data. The model solutions were constrained in their predicted drop size distributions (DSDs) due to the fixed DSD parameters assumed in a two‐moment microphysics scheme. A case study from the Mid‐latitude Continental Clouds and Convection Experiment field project demonstrated how model results can be used to contextualize the disdrometer observations which are limited in sample size, spatial coherence, and detection of small drops and low drop concentrations. The case study also showed that the spatial patterns of the statistically derived PGs revealed by the model are consistent with the hypothesized microphysical processes that determine surface rain DSDs. This work demonstrates how leveraging the strengths of observations and models together can improve our understanding and representation of rain microphysical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Parametric, Frequency-Domain Approach for Clutter Analysis & Rejection in Remote Sensing

A novel approach is presented for parametric analysis of remotely-sensed ground and cloud clutter. A spatial-frequency-domain clutter model is generated from an extensive, one-year database of weather imagery and statistics are given for each spatial frequency. This approach is useful for the analysis and design of spatial and temporal clutter-rejection filters, which can also be analyzed in this domain.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Estimation of context for statistical classification of multispectral image data

Recent investigations have demonstrated the effectiveness of a contextual classifier that combines spatial and spectral information employing a general statistical approach. This statistical classification algorithm exploits the tendency of certain ground cover classes to occur more frequently in some spatial contexts than in others. Indeed, a key input to this algorithm is a statistical characterization of the context: the context function. An unbiased estimator of the context function is discussed which, besides having the advantage of statistical unbiasedness, has the additional advantage over other estimation techniques of being amenable to an adaptive implementation in which the context-function estimate varies according to local contextual information. Results from applying the unbiased estimator to the contextual classification of three real Landsat data sets are presented and contrasted with results from noncontextual classifications and from contextual classifications utilizing other context-function estimation techniques.

Tilton, J. C.↗
DOE OSTI2022

Special Issue: Geostatistics and Machine Learning

Abstract Recent years have seen a steady growth in the number of papers that apply machine learning methods to problems in the earth sciences. Although they have different origins, machine learning and geostatistics share concepts and methods. For example, the kriging formalism can be cast in the machine learning framework of Gaussian process regression. Machine learning, with its focus on algorithms and ability to seek, identify, and exploit hidden structures in big data sets, is providing new tools for exploration and prediction in the earth sciences. Geostatistics, on the other hand, offers interpretable models of spatial (and spatiotemporal) dependence. This special issue on Geostatistics and Machine Learning aims to investigate applications of machine learning methods as well as hybrid approaches combining machine learning and geostatistics which advance our understanding and predictive ability of spatial processes.

58 GEOSCIENCES↗
DOE OSTI2023

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Benchmarking structural evolution methods for training of machine learned interatomic potentials

When creating training data for machine-learned interatomic potentials (MLIPs), it is common to create initial structures and evolve them using molecular dynamics (MD) to sample a larger configuration space. Here, we benchmark two other modalities of evolving structures, contour exploration (CE) and dimer-method (DM) searches against MD for their ability to produce diverse and robust density functional theory training data sets for MLIPs. We also discuss the generation of initial structures which are either from known structures or from random structures in detail to further formalize the structure-sourcing processes in the future. The polymorph-rich zirconium-oxygen composition space is used as a rigorous benchmark system for comparing the performance of MLIPs trained on structures generated from these structural evolution methods. Using Behler–Parrinello neural networks as our MLIP models, we find that CE and the DM searches are generally superior to MD in terms of spatial descriptor diversity and statistical accuracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unraveling Cu Chemical Signature in CdTe by Spectral Fluorescence Mapping

X-ray absorption spectroscopy (XAS) has been shown to be a powerful tool to unravel the chemical environment of a given atom within a matrix. When used in correlative X-ray microscopy approaches, XAS allows one to probe with nanoscale precision regions of particular interest in an absorber. Herein, we use X-ray absorption near edge structure (XANES) to evaluate the chemical environment of Cu atoms within a CdTe solar cell. The reconstruction of XANES spectra from XRF maps have unfolded 2D maps of Cu chemical structures. In this work, we found that most Cu atoms exist in Cu 2 Te and Cu 1.4 Te phase. Moreover, we found traces of CuTe, Cu 2 O, CuO, Cu 2 S, CuS, and metallic Cu phase. Investigating Cu chemical structures at different performing areas, we found no observable correlation between Cu chemical structures and electrical performance. This approach allows tracking of Cu chemical structures along with electrical performance and elemental distribution simultaneously, with high spatial resolution in a statistically practical way.

CdTe↗
DOE OSTI2023

Nano-infrared imaging of metal insulator transition in few-layer 1T-TaS 2

Abstract Among the family of transition metal dichalcogenides, 1T-TaS 2 stands out for several peculiar physical properties including a rich charge density wave phase diagram, quantum spin liquid candidacy and low temperature Mott insulator phase. As 1T-TaS 2 is thinned down to the few-layer limit, interesting physics emerges in this quasi 2D material. Here, using scanning near-field optical microscopy, we perform a spatial- and temperature-dependent study on the phase transitions of a few-layer thick microcrystal of 1T-TaS 2 . We investigate encapsulated air-sensitive 1T-TaS 2 prepared under inert conditions down to cryogenic temperatures. We find an abrupt metal-to-insulator transition in this few-layer limit. Our results provide new insight in contrast to previous transport studies on thin 1T-TaS 2 where the resistivity jump became undetectable, and to spatially resolved studies on non-encapsulated samples which found a gradual, spatially inhomogeneous transition. A statistical analysis suggests bimodal high and low temperature phases, and that the characteristic phase transition hysteresis is preserved down to a few-layer limit.

42 ENGINEERING↗
DOE OSTI2023

Feasibility of Correlated Extensive Air Shower Detection with a Distributed Cosmic-Ray Network

We explore the sensitivity offered by a global network of cosmic-ray detectors to a novel, unobserved phenomenon: widely separated simultaneous extended air showers. Existing localized observatories work independently to observe individual showers, offering insight into the source and nature of ultrahigh-energy cosmic rays. However no current observatory is large enough to provide sensitivity to anticipated processes such as the Gerasimova–Zatsepin effect or potential new physics that generate simultaneous air showers separated by hundreds to thousands of kilometers. A global network of consumer electronics (the Cosmic Rays Found In Smartphones (CRAYFIS) experiment), may provide a novel opportunity for observation of such phenomena. Two user scenarios are explored. In the first, with maximal user adoption, we find that statistically significant discoveries of spatially separated but coincident showers are possible within a couple years. In the second, more practical adoption model with 10 6 active devices, we find a worldwide CRAYFIS to be sensitive to novel "burst" phenomena where many simultaneous extensive air showers (EAS) occur at once.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The Earth Model Column Collaboratory (EMC 2 ) v1.1: an open-source ground-based lidar and radar instrument simulator and subcolumn generator for large-scale models

Abstract. Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

58 GEOSCIENCES↗
NASA NTRS1973

The detection of intermediate size magnetic anomalies in Cosmos-49 and OGO-2, 4, and 6 data

Benkova, Dolginov, and Simonenko have recently reported the presence of intermediate size magnetic anomalies in the data from COSMOS-49 and hypothesized a crustal and/or upper mantle origin. The spherical harmonic models of the internal potential function were examined, based on the OGO-2, 4, and 6 data (POGO (10/68) and later models), and verified the locations and amplitudes of those anomalies whose wavelengths approximate 4000 km. The comparison was made by subtracting a field model developed with a truncated series of n* = 9 from one computed with n* = 11 and generating a residual map equivalent to the COSMOS-49 data. The patterns of delta F so computed from POGO were then compared with the IZMIRAN maps and also were analyzed statistically, in both the spatial and frequency domains, using residuals computed from the raw COSMOS-49 data with the n* = 9 COSMOS-49 field model as reference. The two sets of data were thus derived from completely independent sets of observations and field references. The two patterns are shown to agree very well over the whole earth surface up to the 50 deg latitude limit of COSMOS-49.

Regan, R. D.↗