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 487 records · Page 27

Imaging the 3.4-micron feature in Comet Levy (1990c)

The IRTF's ProtoCAM has been used to observe the organic emission feature peaking near 3.36 microns in Comet Levy. 2D images of the spatial distribution of the 3.4-micron organic emission exhibit brightness profiles statistically indistinguishable from that of the dust continuum, within 1225 of the cometary nucleus; this suggests a progenitor for the emission feature which is not a photodissociation fragment of some other species.

Klavetter, James J.↗

Factor analysis as a tool for spectral line component separation 21cm emission in the direction of L1780

The spectra of the 21cm HI radiation from the direction of L1780, a small high-galactic latitude dark/molecular cloud, were analyzed by multivariate methods. Factor analysis was performed on HI (21cm) spectra in order to separate the different components responsible for the spectral features. The rotated, orthogonal factors explain the spectra as a sum of radiation from the background (an extended HI emission layer), and from the L1780 dark cloud. The coefficients of the cloud-indicator factors were used to locate the HI 'halo' of the molecular cloud. Our statistically derived 'background' and 'cloud' spectral profiles, as well as the spatial distribution of the HI halo emission distribution were compared to the results of a previous study which used conventional methods analyzing nearly the same data set.

Toth, L. V.↗

Comparison of Surface Ground Temperature from Satellite Observations and the Off-Line Land Surface GEOS Assimilation System

The surface ground temperature (Tg) is an important meteorological variable, because it represents an integrated thermal state of the land surface determined by a complex surface energy budget. Furthermore, Tg affects both the surface sensible and latent heat fluxes. Through these fluxes. the surface budget is coupled with the atmosphere above. Accurate Tg data are useful for estimating the surface radiation budget and fluxes, as well as soil moisture. Tg is not included in conventional synoptical weather station reports. Currently, satellites provide Tg estimates globally. It is necessary to carefully consider appropriate methods of using these satellite data in a data assimilation system. Recently, an Off-line Land surface GEOS Assimilation (OLGA) system was implemented at the Data Assimilation Office at NASA-GSFC. One of the goals of OLGA is to assimilate satellite-derived Tg data. Prior to the Tg assimilation, a thorough investigation of satellite- and model-derived Tg, including error estimates, is required. In this study we examine the Tg from the n Project (ISCCP DI) data and the OLGA simulations. The ISCCP data used here are 3-hourly DI data (2.5x2.5 degree resolution) for 1992 summer months (June, July, and August) and winter months (January and February). The model Tg for the same periods were generated by OLGA. The forcing data for this OLGA 1992 simulation were generated from the GEOS-1 Data Assimilation System (DAS) at Data Assimilation Office NASA-GSFC. We examine the discrepancies between ISCCP and OLGA Tg with a focus on its spatial and temporal characteristics, particularly on the diurnal cycle. The error statistics in both data sets, including bias, will be estimated. The impact of surface properties, including vegetation cover and type, topography, etc, on the discrepancies will be addressed.

Yang, R.↗

Stratospheric Analysis and Forecasting in the Northern Winter of 1999/2000: The NASA DAO's GEOS-3 System

An evaluation is presented of the performance in the northern winter 1999/2000 of the GEOS-3 troposphere-stratosphere data assimilation system (DAS). The impacts of the two main input data types are assessed: upper-air soundings (sondes) provide wind and temperature information and satellite-based (Tiros Operational Vertical Sounders: TOVS) give estimates of the thermal structure. It is shown that in the low stratosphere (300-70hPa) the analyses are generally slightly warmer than the sonde data, but colder than the TOVS data; this relationship reverses between 70 and 10 hPa. There are geographical biases, related to the spatial and temporal coverage of the observation types and to the statistical weights assigned to them in the DAS. Forecasts show a tendency to reduce zonal asymmetries in the atmospheric flow and to suppress stratospheric temperature minima. In the DAS, the analysis increments compensate for this, but it leads to important biases in the multi-day forecasts. The analysis increments are as large as the diabatic forcing in the lower polar stratosphere, indicating a substantial model bias. The results provide important insights into the roles of different data types and the circulation model in producing accurate analyses for studies of polar chemistry and physical processes.

Pawson, Steven↗

Data Assimilation and Data Fusion for Planetary Atmospheres

The overarching goal of this Cooperative Agreement was to develop a model and procedures for the data assimilation of planetary spacecraft atmospheric observations. Data assimilation - in its application to weather analysis and prediction - is the process of finding an initial state of the meteorological variables (winds, temperatures, pressures, etc.) of an atmosphere, which, when propagated forward in time using a deterministic general circulation model, reproduces all of the available observations over that time to within the measurement and computational errors. With this definition, data assimilation is seen to be a natural extension of well-known least-squares minimization techniques. The primary complication arises from the scale of the problem: For the Martian atmosphere with the available nadir-viewing Thermal Emission Spectrometer data from Mars Global Surveyor, approximately 1,000,000 individual measurements of channel radiances (in the 15-micrometer region, where these radiances relate directly to the surface and atmospheric temperature) were made per day. A suitable general circulation model for dealing with this data set has on the order of 20,000 independent variables. After some spatial and temporal averaging of the data - which provides a necessary statistical estimate of the representativeness of the measurements, a crucial issue in data assimilation - the problem reduces in scale to the solution of approximately 50,000 equations for the 20,000 variables.

Houben, Howard↗

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗

Statistical model of the stimulated forward Brillouin scattering driven by a randomized laser beam in plasma

The modeling of a spatially incoherent laser beam remains a central problem of the parametric instabilities in the context of inertial confinement fusion. This letter gives a simplified and comprehensive overview of the recent analytical developments regarding the modeling of these laser beams and a comparison with a dedicated experiment. Our model accounts for the first time for the statistical standard deviation of the gain and accurately captures the entanglement between wave mixing processes and the speckle correlations thus resolving the longstanding contradictions between the random phase approximation and the model of independent speckles. It is successfully compared to a recent laser beam spray experiment and the associated paraxial simulations, demonstrating that backscattering predictions require accounting for the beam spray. Furthermore, our framework thus provides a way to evaluate and guide the analysis of parametric instabilities in high laser energy experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗

Hybrid chemical characterization of latent images in EUV resist with 12 nm half-pitch features

With the advancement of high numerical aperture extreme ultraviolet (EUV) lithography, the new platform will enable chipmakers to achieve critical dimensions of 8 nm. However, resist materials face significant challenges in delivering increased sensitivity while managing rising stochastic variations. We aim to develop comprehensive techniques to characterize the chemical profile of latent images, stored in EUV resists after exposure and postexposure baking, which is essential for understanding the origin of stochastic effects. Infrared photo-induced force microscopy (IR PiFM) is a bimodal atomic force microscopy technique combined with an infrared light source, allowing for simultaneous sub-5 nm topographic and chemical characterization within a localized environment. Critical-dimension resonant soft X-ray scatterometry (CD-RSoXS) provides statistical data that reveal structural and chemical information for comparative analysis. For the first time, IR PiFM has been used to chemically map the latent images of EUV resists (after exposure and postexposure baking) at a record high resolution of 12 nm half-pitch, enabling nondestructive analysis of patterns prior to development. Furthermore, CD-RSoXS offers direct experimental observation and comparison of exposed, postexposure baked, and developed patterns, which align with the IR PiFM results. We demonstrate that the IR PiFM technique offers valuable insights into both high spatial resolution and local chemical information simultaneously. In addition, CD-RSoXS provides statistical results that support our main findings. This hybrid metrology approach leverages a multifaceted dataset by combining the most reliable information from each source, which is essential for a comprehensive understanding of the stochastic effects in EUV lithography processes.

O’Reilly, Padraic↗

Long-term Average Spectral and Spatial Distributions of Plasmaspheric Hiss Observed by the Akebono and IMAGE Satellites

The radiation belt slot region is known to result from losses of energetic electrons by enhanced pitch-angle scattering by whistler mode waves associated with plasmaspheric hiss emission. The distributions of whistler mode waves in the slot L range are therefore important for understanding the electron radiation belt. The sources and distributions of the waves are, however, still controversial. In the present study, using the Akebono/MCA data [1989-20051 and the IMAGE/RPI data [2000-20051, we have constructed the average plasmaspheric hiss spectral distributions over a broad frequency range. In addition, we have investigated the spatial distributions of plasmaspheric hiss with the wave map technique [Green et a1.(2005)]. Our study shows that the broadband plasmaspheric hiss are distributed in the frequency range of 100Hz to several kHz, and exhibit a broad intensive peak. The frequency of the intensity peak tends to increase with magnetic latitude. The frequencies of the most intense waves in the nominal slot L range (2<3) during quiet times (Dst>-50nT) are found to be between 300Hz and 600Hz on average. During high storm activity (Dst <-150nT), however, the peak frequencies become slightly lower. The intensity of plasmaspheric hiss clearly depends on substorm activity as measured by the AE index, consistent with Meredith et a1.(2004). The hiss wave intensity maps also show a strong local time asymmetry. The large amplitude waves are observed at 6:OO-19:OO MLT. From our extensive analysis, we have also found an L dependence of hiss activity, with the larger amplitude waves being observed at lower L during substorm active conditions. The same tendency can be found for solar activity. The average intensities of the waves during 1989-1991 and 2000-2001 are a few dB larger than those during 1992-1997 and 2005. The most intense waves are observed at lower L during high solar activity. The statistical study on spectrum features of the plasmaspheric hiss together with the spatial distribution show clear dependences of storm, substorm and solar activities.

Fung, Shing↗

Myriad World Baseline: Global Geodemographic Estimates

The LandScan Myriad World Baseline (MWB) method produces global, residential (nighttime/home-location) gridded geodemographic estimates based on 5-year age/gender cohorts—at 30-arcsecond (≈1 km) resolution. MWB is designed to fill gaps where detailed, georeferenced survey data (e.g., Demographic and Health Surveys (DHS)) are missing or outdated, and to provide a baseline that can support human security analysis, including consequence assessment, “patterns of life” modeling, and scenario-based population futures. MWB’s workflow spatializes household-level age/gender characteristics from the GLOPOP-S dataset by conflating household and gridded expected relative wealth adapted from Global Gridded Relative Deprivation Index (GRDI), then adjusts them to a target year of interest. Age/gender estimates are then applied to harmonize lowest-administrative-level statistics with LandScan residential counts, yielding final geodemographic estimates. Two validation case studies are presented: Ghana (2021) and Tokyo/Kanagawa, Japan (2020), illustrating spatial variability in demographic cohorts and comparing MWB outputs to official gridded statistics. Results show close overall alignment relative to validation criteria including population pyramids and age-dependency ratios.

Tuccillo, Joe [ORNL] (ORCID:0000000259300943)↗

Nanoscale Spatial Distribution of Supported Nanoparticles Controls Activity and Stability in Powder Catalysts for CO Oxidation and Photocatalytic H 2 Evolution

Here we discuss, supported metal nanoparticles are essential components of high-performing catalysts, and their structures are intensely researched. In comparison, nanoparticle spatial distribution in powder catalysts is conventionally not quantified, and the influence of this collective property on catalyst performance remains poorly investigated. Here, we demonstrate a general colloidal self-assembly method to control uniformity of nanoparticle spatial distribution on common industrial powder supports. We quantify distributions on the nanoscale using image statistics and show that the type of nanospatial distribution determines not only the stability, but also the activity of heterogeneous catalysts. Widely investigated systems (Au–TiO 2 for CO oxidation thermocatalysis and Pd–TiO 2 for H 2 evolution photocatalysis) were used to showcase the universal importance of nanoparticle spatial organization. Spatially and temporally resolved microkinetic modeling revealed that nonuniformly distributed Au nanoparticles suffer from local depletion of surface oxygen, and therefore lower CO oxidation activity, as compared to uniformly distributed nanoparticles. Nanoparticle spatial distribution also determines the stability of Pd–TiO 2 photocatalysts, because nonuniformly distributed nanoparticles sinter while uniformly distributed nanoparticles do not. This work introduces new tools to evaluate and understand catalyst collective (ensemble) properties in powder catalysts, which thereby pave the way to more active and stable heterogeneous catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scale Dependence of Multiplier Distributions for Particle Concentration, Enstrophy, and Dissipation in the Inertial Range of Homogeneous Turbulence

Turbulent flows preferentially concentrate inertial particles depending on their stopping time or Stokes number, which can lead to significant spatial variations in the particle concentration. Cascade models are one way to describe this process in statistical terms. Here, we use a direct numerical simulation (DNS) dataset of homogeneous, isotropic turbulence to determine probability distribution functions (PDFs) for cascade multipliers, which determine the ratio by which a property is partitioned into subvolumes as an eddy is envisioned to decay into smaller eddies. We present a technique for correcting effects of small particle numbers in the statistics. We determine multiplier PDFs for particle number, flow dissipation, and enstrophy, all of which are shown to be scale dependent. However, the particle multiplier PDFs collapse when scaled with an appropriately defined local Stokes number. As anticipated from earlier works, dissipation and enstrophy multiplier PDFs reach an asymptote for sufficiently small spatial scales. From the DNS measurements, we derive a cascade model that is used it to make predictions for the radial distribution function (RDF) for arbitrarily high Reynolds numbers, Re, finding good agreement with the asymptotic, infinite Re inertial range theory of Zaichik and Alipchenkov [New J. Phys. 11, 103018 (2009)]. We discuss implications of these results for the statistical modeling of the turbulent clustering process in the inertial range for high Reynolds numbers inaccessible to numerical simulations.

turbulent flows↗

High-Resolution South American Wind Resource Data Downscaled with Generative Machine Learning Conditioned on Near-Surface Observations

High-resolution historical wind data was developed for the entirety of South America using the innovative Super-Resolution for Renewable Resource Data (sup3r) machine learning framework. The publicly available Sup3rWind South America dataset represents a significant advancement in wind resource data generation, leveraging generative machine learning conditioned on near-surface observations from the Meteorological Assimilation Data Ingest System (MADIS) to efficiently and accurately downscale coarse reanalysis data from the European Centre for Medium-Range Weather Forecasts (ERA5). This approach produces fine-scale, spatially and temporally coherent wind and meteorological fields hundreds of times more computationally efficient than traditional numerical weather modeling methods, enabling access to high-fidelity wind information across both continental and offshore regions. Sup3rWind South America builds on the earlier Sup3rWind Ukraine dataset through improvements in model architecture and outputs conditioned on near-surface observation inputs. As with the Ukraine data release, this dataset includes wind speed, wind direction, temperature, relative humidity, and pressure at a horizontal resolution of ~2 km, representing a 15x spatial enhancement relative to the 31 km ERA5 grid. Wind speed and direction are provided at 5-minute resolution, a 12x temporal refinement compared to the hourly ERA5 data, while temperature, relative humidity, and pressure remain at hourly resolution. The data covers all years from 2005 to 2024. Before downscaling, ERA5 inputs were bias-corrected using long-term monthly means and a limited number of quality-controlled observations to align large-scale statistics with regional conditions. The resulting dataset is the first publicly available high-resolution timeseries wind record that provides full spatial coverage of South America. Model validation demonstrates strong agreement with observations across several statistical metrics, consistent with other state-of-the-art high-resolution wind resource datasets. The potential applications of Sup3rWind South America span renewable energy resource assessment, energy system modeling, and grid resilience analysis. The 20-year record and high spatial and temporal resolution support accurate estimation of long-term energy yield and the economic feasibility of potential wind development sites. Continuous coverage across both continental and offshore regions enables comprehensive site prospecting within exclusive economic zones. The 2 km, 5-minute resolution data provide the spatial and temporal variability required for power system simulation, operational planning, and regional risk assessments.

17 WIND ENERGY↗

Effective Assimilation of SMAP Observations Using Statistical Techniques

Statistical techniques permit the retrieval of soil moisture estimates in a model climatology while retaining the spatial and temporal signatures of the satellite observations. As a consequence, the need for bias correction prior to an assimilation of these estimates is reduced, which could result in a more effective use of the independent information provided by the satellite observations. In this study, a statistical neural network (NN) retrieval algorithm is calibrated using SMAP brightness temperature observations and modeled soil moisture estimates (similar to those used to calibrate the SMAP Level 4 DA system). Daily values of surface soil moisture are estimated using the NN and then assimilated into the NASA Catchment model. The skill of the assimilation estimates is assessed based on a comprehensive comparison to in situ measurements from the SMAP core and sparse network sites as well as the International Soil Moisture Network. The NN retrieval assimilation is found to significantly improve the model skill, particularly in areas where the model does not represent processes related to agricultural practices. Additionally, the NN method is compared to assimilation experiments using traditional bias correction techniques. The NN retrieval assimilation is found to more effectively use the independent information provided by SMAP resulting in larger model skill improvements than assimilation experiments using traditional bias correction techniques.

Kolassa, J.↗

The ERTS-1 investigation (ER-600). Volume 6: ERTS-1 signature extension analysis, July 1972 - June 1973

Feature classification, spatially and temporally, was extended over the Houston test site area. The Earth Resources Technology Satellite (ERTS-1) multispectral scanner data from August, September, and October 1972, of five widely separated lakes were used as statistical training fields and test sites. Short term temporal (same day to 36 days) and moderately long term spatial (within and between three ERTS multispectral scanner frames) signature extensions have been verified with respect to large relatively homogeneous features. The most significant feature dependent variable affecting spatial and short term extension was water turbidity. Long term signature extension will require a model to compensate or modify the ERTS-1 multispectral scanner data for significant sun angle changes. The presence of atmospheric haze changed the absolute signature but always by approximately the same amount so that the measured water signature was always the same. The normally occurring variations in atmospheric haze conditions had no major effect on the water signatures in this study.

Erb, R. B.↗