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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comment on “Five Decades of Observed Daily Precipitation Reveal Longer and More Variable Drought Events Across Much of the Western United States”

Abstract Changes in precipitation patterns with climate change could have important impacts on human and natural systems. Zhang et al. (2021, https://doi.org/10.1029/2020gl092293 ) report trends in daily precipitation patterns over the last five decades in the western United States, focusing on meteorological drought. They report that dry intervals (calculated at the annual or seasonal level) have increased across much of the southwestern U.S., with statistical assessment suggesting the results are statistically robust. However, Zhang et al. (2021, https://doi.org/10.1029/2020gl092293 ) preprocess their annual (or seasonal) averages to compute 5‐year moving window averages before using established statistical techniques for trend analysis that assume independence about some fixed trend. Here we show that the moving window preprocessing violates that independence assumption and inflates the statistical significance of their trend estimates. This raises questions about the robustness of their results. We conclude by discussing the difficulty of adjusting for spatial structure when assessing time trends in a regional context.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

At Risk Population Estimates for Belarus, Poland and Slovakia with Machine Learning

High-resolution gridded population modeling is crucial for various applications, including disaster response planning, infectious disease spread modeling, climate change impact estimation, policy development, and more. Multiple gridded population datasets have been developed, each tailored to meet specific objectives. Among them, LandScan Global dataset is designed to represent ambient and unwarned population distributions. However, this dataset relies on a statistical approach that requires manual adjustments, making it time consuming and labour intensive. Existing machine learning (ML) methods often train and test at different spatial resolutions, potentially leading to inflated results, and they rely on Census population totals for disaggregation. To address these limitations, in this study we developed population estimates using ML models trained and tested at a consistent 30 arc-second resolution (≈1 square kilometer), specifically using Random Forest (RF) and XGBoost. These models were trained on 2020 datum to predict for 2021 for three countries: Belarus, Poland, and Slovakia. Our findings show that both RF (MAE varies from 5.75 to 13.25) and XGBoost (MAE varies from 8.15 to 23.44) model performance is close to LandScan Global estimates. Furthermore, neither of the models performed the best across all grid cells: the RF model was more effective in areas with lower populations, while XGBoost excelled in more densely populated regions. The proposed approach can be used for countries where the Census data is not available.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗
DOE OSTI2025

Environmental Conditions Affecting Global Mesoscale Convective System Occurrence

Abstract The ERA5 environments of mesoscale convective systems (MCSs), tracked from satellite observations, are assessed over a 20-yr period. The use of a large set of MCS tracks allows us to robustly test the sensitivity of the results to factors such as region, latitude, and diurnal cycle. We aim to provide novel information on environments of observed MCSs for assessments of global atmospheric models and to improve their ability to simulate MCSs. Statistical analysis of all tracked MCSs is performed in two complementary ways. First, we investigate the environments when an MCS has occurred at different spatial scales before and after MCS formation. Several environmental variables are found to show marked changes before MCS initiation, particularly over land. The vertically integrated moisture flux convergence shows a robust signal across different regions and when considering MCS initiation diurnal cycle. We also found spatial scale dependence of the environments between 200 and 500 km, providing new evidence of a natural length scale for use with MCS parameterization. In the second analysis, the likelihood of MCS occurrence for given environmental conditions is evaluated, by considering all environments and determining the probability of being in an MCS core or shield region. These are compared to analogous non-MCS environments, allowing discrimination between conditions suitable for MCS and non-MCS occurrence. Three environmental variables are found to be useful predictors of MCS occurrence: total column water vapor, midlevel relative humidity, and total column moisture flux convergence. Such relations could be used as trigger conditions for the parameterization of MCSs, thereby strengthening the dependence of the MCS scheme on the environment. Significance Statement Large storm systems called mesoscale convective systems form across Earth. These are collections of thunderstorms, with associated high-level clouds that produce substantial, lighter rainfall and modulate Earth’s energy balance. They produce hazardous weather conditions, such as floods and high winds, and are responsible for a high percentage of rainfall in many regions globally. We investigate the environmental conditions under which they form, so that we can understand the spatial extent of the environment which is important for their formation, and also where and when the effects of these storms might be felt. The novel information generated here should help improve the representation of these storms in weather and climate models, improving the prediction of rainfall, thunderclouds, and high-level clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

A laser velocimeter system for large-scale aerodynamic testing

A unique laser velocimeter was developed that is capable of sensing two orthogonal velocity components from a variable remote distance of 2.6 to 10 m for use in the 40- by 80-Foot and 80- by 120-Foot Wind Tunnels and the Outdoor Aerodynamic Research Facility at Ames Research Center. The system hardware, positioning instrumentation, and data acquisition equipment are described in detail; system capabilities and limitations are discussed; and expressions for systematic and statistical accuracy are developed. Direct and coupled laboratory measurements taken with the system are compared with measurements taken with a laser velocimeter of higher spatial resolution, and sample data taken in the open circuit exhaust flow of a 1/50-scale model of the 80- by 120-Foot Wind Tunnel are presented.

Reinath, M. S.↗
NASA NTRS1987

Evaluation of several classification schemes for mapping forest cover types in Michigan

Landsat MSS data were evaluated for mapping forest cover types in the northern Lower Peninsula of Michigan. The study examined seasonal variations, interpretation procedures and vegetation composition/distribution and their effect on overall classification accuracy and ability to identify individual pine species. Photographic images were used for visual interpretations while digital analysis was performed using a common (ERDAS) microcomputer image processing system. The classification schemes were evaluated using contingency tables and were ranked using the KAPPA statistic. The various classification schemes were ranked differentially according to study site location. Visual interpretation procedures ranked best, or least accurate, depending on the spatial distribution and complexity of the forest cover. Supervised classification techniques were more accurate than unsupervised clustering over all sites and seasons. Maximum likelihood classification of June data was superior to any digital classification technique of February data. The study indicates that classification accuracy is more dependent on the composition and distribution of forests in the northern lower Peninsula of Michigan than on the selection of a particular classification scheme.

Hudson, W. D.↗
NASA NTRS2012

Geospatial Method for Computing Supplemental Multi-Decadal U.S. Coastal Land-Use and Land-Cover Classification Products, Using Landsat Data and C-CAP Products

This paper discusses the development and implementation of a geospatial data processing method and multi-decadal Landsat time series for computing general coastal U.S. land-use and land-cover (LULC) classifications and change products consisting of seven classes (water, barren, upland herbaceous, non-woody wetland, woody upland, woody wetland, and urban). Use of this approach extends the observational period of the NOAA-generated Coastal Change and Analysis Program (C-CAP) products by almost two decades, assuming the availability of one cloud free Landsat scene from any season for each targeted year. The Mobile Bay region in Alabama was used as a study area to develop, demonstrate, and validate the method that was applied to derive LULC products for nine dates at approximate five year intervals across a 34-year time span, using single dates of data for each classification in which forests were either leaf-on, leaf-off, or mixed senescent conditions. Classifications were computed and refined using decision rules in conjunction with unsupervised classification of Landsat data and C-CAP value-added products. Each classification's overall accuracy was assessed by comparing stratified random locations to available reference data, including higher spatial resolution satellite and aerial imagery, field survey data, and raw Landsat RGBs. Overall classification accuracies ranged from 83 to 91% with overall Kappa statistics ranging from 0.78 to 0.89. The accuracies are comparable to those from similar, generalized LULC products derived from C-CAP data. The Landsat MSS-based LULC product accuracies are similar to those from Landsat TM or ETM+ data. Accurate classifications were computed for all nine dates, yielding effective results regardless of season. This classification method yielded products that were used to compute LULC change products via additive GIS overlay techniques.

Spruce, J. P.↗
NASA NTRS2023

An Unbiased CO Survey Toward the Northern Region of the Small Magellanic Cloud with the Atacama Compact Array. II. CO Cloud Catalog

The nature of molecular clouds and their statistical behavior in subsolar metallicity environments are not fully explored yet. We analyzed data from an unbiased CO (J = 2–1) survey at the spatial resolution of ∼2 pc in the northern region of the Small Magellanic Cloud with the Atacama Compact Array to characterize the CO cloud properties. A cloud-decomposition analysis identified 426 spatially/velocity-independent CO clouds and their substructures. Based on the cross-matching with known infrared catalogs by Spitzer and Herschel, more than 90% CO clouds show spatial correlations with point sources. We investigated the basic properties of the CO clouds and found that the radius–velocity linewidth (R–σ v ) relation follows the Milky Way-like power-law exponent, but the intercept is ∼1.5 times lower than that in the Milky Way. The mass functions (dN/dM) of the CO luminosity and virial mass are characterized by an exponent of ∼1.7, which is consistent with previously reported values in the Large Magellanic Cloud and in the Milky Way.

Local Group↗
DOE OSTI2022

A physical method for downscaling land surface temperatures using surface energy balance theory

Fine-resolution land surface temperature (LST) derived from thermal infrared remote sensing images is a good indicator of surface water status and plays an essential role in the exchange of energy and water between land and atmosphere. A physical surface energy balance (SEB)-based LST downscaling method (DTsEB) is developed to downscale coarse remotely sensed thermal infrared LST products with fine-resolution visible and near-infrared data. Here, the DTsEB method is advantageous for its ability to mechanically interrelate surface variables contributing to the spatial variation of LST, to quantitatively weigh the contributions of each related variable within a physical framework, and to efficaciously avoid the subjective selection of scaling factors and the establishment of statistical regression relationships. The applicability of the DTsEB method was tested by downscaling 12 scenes of 990 m Moderate Resolution Imaging Spectroradiometer (MODIS) and aggregated Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) LST products to 90 m resolution at six overpass times between 2005 and 2015 over three 9.9 km by 9.9 km cropland (mixed by grass, tree, and built-up land) study areas. Three typical LST downscaling methods, namely the widely applied TsHARP, the later developed least median square regression downscaling (LMS) and the geographically weighted regression (GWR), were introduced for intercomparison. The results showed that the DTsEB method could more effectively reconstruct the subpixel spatial variations in LST within the coarse-resolution pixels and achieve a better downscaling accuracy than the TsHARP, LMS and GWR methods. The DTsEB method yielded, on average, root mean square errors (RMSEs) of 2.01 K and 1.42 K when applied to the MODIS datasets and aggregated ASTER datasets, respectively, which were lower than those obtained with the TsHARP method, with average RMSEs of 2.41 K and 1.71 K, the LMS method, with average RMSEs of 2.35 K and 1.63 K, and the GWR method, with average RMSEs of 2.38 K and 1.64 K, respectively. The contributions of the related surface variables to the subpixel spatial variation in the LST varied both spatially and temporally and were different from each other. In summary, the DTsEB method was demonstrated to outperform the TsHARP, LMS, and GWR methods and could be used as a good alternative for downscaling LST products from coarse to fine resolution with high robustness and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Harnessing the Population Statistics of Subhalos to Search for Annihilating Dark Matter

The Milky Way's dark matter halo is expected to host numerous low-mass subhalos with no detectable associated stellar component. Such subhalos are invisible unless their dark matter annihilates to visible states such as photons. One of the established methods for identifying candidate subhalos is to search for individual unassociated gamma-ray sources with properties consistent with the dark matter expectation. However, robustly ruling out an astrophysical origin for any such candidate is challenging. In this work, we present a complementary approach that harnesses information about the entire population of subhalos—such as their spatial and mass distribution in the Galaxy—to search for a signal of annihilating dark matter. Using simulated data, we show that the collective emission from subhalos can imprint itself in a unique way on the statistics of observed photons, even when individual subhalos may be too dim to be resolved on their own. Additionally, we demonstrate that, for the models we consider, the signal can be identified even in the face of unresolved astrophysical point-source emission of extragalactic and Galactic origin. Furthermore, this establishes a new search technique for subhalos that is complementary to established methods, and that could have important ramifications for gamma-ray dark matter searches using observatories such as the Fermi Large Area Telescope and the Cerenkov Telescope Array.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

The spectral correlation of functional clouds: a statistical test for theoretical models

We compute the spectral correlation function (SCF) of CO J=1-0 maps of moleculary cloud complexes. The SCF is a power law over approximately an order of magnitude in spatial separation in every map. The power law slope of the SCF, its normalization, and the spectral line width averaged over the whole map, are computed for all the observational maps.

turbulence kinematics↗
DOE OSTI2024

Autocalibration of the E3SM Version 2 Atmosphere Model Using a PCA-Based Surrogate for Spatial Fields

Global Climate Model tuning (calibration) is a tedious and time-consuming process, with high-dimensional input and output fields. Experts typically tune by iteratively running climate simulations with hand-picked values of tuning parameters. Many, in both the statistical and climate literature, have proposed alternative calibration methods, but most are impractical or difficult to implement. We present a practical, robust, and rigorous calibration approach on the atmosphere-only model of the Department of Energy's Energy Exascale Earth System Model (E3SM) version 2. Our approach can be summarized into two main parts: (a) the training of a surrogate that predicts E3SM output in a fraction of the time compared to running E3SM, and (b) gradient-based parameter optimization. To train the surrogate, we generate a set of designed ensemble runs that span our input parameter space and use polynomial chaos expansions on a reduced output space to fit the E3SM output. We use this surrogate in an optimization scheme to identify values of the input parameters for which our model best matches gridded spatial fields of climate observations. To validate our choice of parameters, we run E3SMv2 with the optimal parameter values and compare prediction results to expertly-tuned simulations across 45 different output fields. This flexible, robust, and automated approach is straightforward to implement, and we demonstrate that the resulting model output matches present day climate observations as well or better than the corresponding output from expert tuned parameter values, while considering high-dimensional output and operating in a fraction of the time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coherent Random Lasing in Subwavelength Quasi‐2D Perovskites

Abstract Quasi‐2D lead halide perovskites have garnered increasing interest as lasing gain media. Relatively simple fabrication, high refractive index, and unique quantum well structure encourage their use in traditional cavity lasers and cavity‐free systems called random lasers (RLs). Despite tremendous advances reported thus far, coherent random lasing in quasi‐2D perovskite subwavelength films has not been reported. Consequently, coherent optical feedback mechanisms in quasi‐2D perovskite systems are still unexplored. Here, this work reports the observation of coherent random lasing in subwavelength quasi‐2D perovskite films. Statistical analysis of spectral measurements reveals Lévy‐like intensity fluctuations, replica symmetry breaking confirms random lasing, and the coherent modes are studied with spectral and spatial correlation techniques. The observed coherent lasing modes are found to be extended states that arise from the random crystal grain structure during fabrication and span the entire pump volume. These modes out‐compete diffusive lasing due to their coherence.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Polishing the Gold Standard: The Role of Orbital Choice in CCSD(T) Vibrational Frequency Prediction

While CCSD(T) with spin-restricted Hartree-Fock (RHF) orbitals has long been lauded for its ability to accurately describe closed-shell interactions, the performance of CCSD(T) on open-shell species is much more erratic, especially when using a spin-unrestricted HF (UHF) reference. Previous studies have shown improved treatment of open-shell systems when a non-HF set of molecular orbitals, like Brueckner or Kohn-Sham density functional theory (DFT) orbitals, is used as a reference. Inspired by the success of regularized orbital-optimized second-order Møller-Plesset perturbation theory (κ-OOMP2) orbitals as reference orbitals for MP3, we investigate the use of κ-OOMP2 orbitals and various DFT orbitals as reference orbitals for CCSD(T) calculations of the corrected ground-state harmonic vibrational frequencies of a set of 36 closed-shell (29 neutrals, 6 cations, 1 anion) and 59 open-shell diatomic species (38 neutrals, 15 cations, 6 anions). The aug-cc-pwCVTZ basis set is used for all calculations. The use of κ-OOMP2 orbitals in this context alleviates difficult cases observed for both UHF orbitals and OOMP2 orbitals. Removing two multireference systems and 12 systems with ambiguous experimental data leaves a pruned data set. Overall performance on the pruned data set highlights CCSD(T) with a B97 orbital reference (CCSD(T):B97), CCSD(T) with a κ-OOMP2 orbital reference (CCSD(T):κ-OOMP2), and CCSD(T) with a B97M-rV orbital reference (CCSD(T):B97M-rV) with RMSDs of 8.48 cm -1 , and 8.50 cm -1 , and 8.75 cm -1 respectively, outperforming CCSD(T):UHF by nearly a factor of 5. Moreover, the performance on the closed- and open-shell subsets shows these methods are able to treat open-shell and closed-shell systems with comparable accuracy and robustness. CCSD(T) with RHF orbitals is seen to improve upon UHF for the closed-shell species, while spatial symmetry breaking in a number of restricted open-shell HF (ROHF) references leads CCSD(T) with ROHF reference orbitals to exhibit the poorest statistical performance of all methods surveyed for open-shell species. The use of κ-OOMP2 orbitals has also proven useful in diagnosing multireference character that can hinder the reliability of CCSD(T).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Stochastic relativistic viscous hydrodynamics from the Metropolis algorithm

We propose an algorithm for simulating stochastic relativistic fluid dynamics based on Metropolis updates. Each step of the algorithm begins with an update based on ideal hydrodynamics. This is followed by proposing random (spatial) momentum transfers between fluid cells, keeping the total energy fixed. These proposals are then accepted or rejected using the change in entropy as a statistical weight. The algorithm reproduces relativistic viscous hydrodynamics in the “density frame,” which is a formulation of viscous hydrodynamics we review and clarify here. This formulation is first order in time and requires no auxiliary dynamical fields such as Π 𝜇⁢𝜈 . The only parameters are the shear and bulk viscosities and the equation of state. Here, by adopting the 3+1 split of general relativity, we extend the Metropolis algorithm to general space-time coordinates, such as Bjorken coordinates, which are commonly used to simulate heavy-ion collisions.

Hydrodynamic noise↗
NASA NTRS1972

Ground pattern analysis in the Great Plains

The author has identified the following significant results. During preliminary analysis of the ERTS-1 RBV band-2 images it was possible to detect many geologic features and to delineate several geologic formation boundaries. Also, the spatial frequency characteristics of some different geologic formations are distinguishable from each other. At the present time there is not enough data to provide any statistical evaluation of these differences, but it is hoped that more quantitative results will be produced within the next few months. Multispectral analysis will begin when the retrospective order for the other spectral bands of the first images is received. Since many geologic formation boundaries are visible on the single band that has been analyzed so far, even greater visibility is expected with the spectral information available in the additional bands.

Egbert, D. D.↗
NASA NTRS1984

A laser velocimeter system for the Ames 40- by 80-foot and 80by 120-foot wind tunnels

A unique laser velocimeter has been developed that is capable of sensing two orthogonal velocity components from a variable remote distance of 2.6 to 10 m for use in the Ames 40- by 80-Foot and 80- by 120-Foot wind tunnels, as well as in the Ames Outdoor Aerodynamic Research Facility. The system hardware, positioning instrumentation, and data acquisition equipment are described in detail. System capabilities and limitations are discussed, and expressions for systematic and statistical accuracy are developed. Direct and coupled laboratory measurements of a turbulent swirling flow are compared with measurements taken with a laser velocimeter of higher spatial resolution. Sample measurements taken in the open-circuit exhaust flow of the 80- by 120-Foot Wind Tunnel 1/50-scale model are presented.

Reinath, M. S.↗