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 1,099 records · Page 61

Spectral characteristics and the extent of paleosols of the Palouse formation

Three spectral models defining the spatial distribution of soil areas by levels of amorphous iron, organic carbon, and the ratio of amorphous iron to organic carbon were developed and field verification studies were conducted. The models used particular Thematic Mapper band ratios selected by statistical correlation with soil chemical data. The ability of the models to indicate erosion severity and to differentiate between iron enriched and carbonate paleosols is discussed. In addition, the effect of vegetation cover on paleosols is addressed.

Frazier, B. E.↗

Structural analysis of airborne flux estimates over a region

Aircraft-based observations of turbulence fields of velocity, moisture, and temperature are used to study coherent turbulent structures that dominate turbulent transfer of moisture and heat above three different eco-systems. Flux traces are defragmented, to reconstruct the presumed full size (along the sampled transect) of these structures, and flux traces are simplified by elimination of those that contribute negligibly to the flux estimate. Structures are analyzed in terms of size, spatial distribution, and contribution to the flux, in the four 'quadrant' modes of eddy-covariance transfer (excess up/down and deficit up/down). The effect of nonlinear detrending of moisture and temperature data on this 'structural analysis,' over surfaces with heterogeneous surface wetness, is also examined. Results over grassland, wetland, and moist and dry agricultural land, show that nonlinear detrending may provide a more physically realistic description of structures. Significant differences are observed between structure size and associated relative flux contribution, between moist and dry areas, with smaller structures playing a more important role over the moist areas. Structure size generally increases with height, as spatial reorganization from smaller structures into larger ones takes place. This coincides with a gradual loss of surface 'signature' (position and clustering of plumes above localized source areas). The data are expected to provide a basis for an eventual statistical description of boundary-layer transfer events , and help to interpret the link between boundary-layer transfer and hydrological surface conditions.

Caramori, Paulo↗

Vertical Scales of Turbulence at the Mount Wilson Observatory

The vertical scales of turbulence at the Mount Wilson Observatory are inferred from data from the University of California at Berkeley Infrared Spatial Interferometer (ISI), by modeling path length fluctuations observed in the interferometric paths to celestial objects and those in instrumental ground-based paths. The correlations between the stellar and ground-based path length fluctuations and the temporal statistics of those fluctuations are modeled on various timescales to constrain the vertical scales. A Kolmogorov-Taylor turbulence model with a finite outer scale was used to simulate ISI data. The simulation also included the white instrumental noise of the interferometer, aperture-filtering effects, and the data analysis algorithms. The simulations suggest that the path delay fluctuations observed in the 1992-1993 ISI data are largely consistent with being generated by refractivity fluctuations at two characteristic vertical scales: one extending to a height of 45 m above the ground, with a wind speed of about 1 m/ s, and another at a much higher altitude, with a wind speed of about 10 m/ s. The height of the lower layer is of the order of the dimensions of trees and other structures near the interferometer, which suggests that these objects, including elements of the interferometer, may play a role in generating the lower layer of turbulence. The modeling indicates that the high- attitude component contributes primarily to short-period (less than 10 s) fluctuations, while the lower component dominates the long-period (up to a few minutes) fluctuations. The lower component turbulent height, along with outer scales of the order of 10 m, suggest that the baseline dependence of long-term interferometric, atmospheric fluctuations should weaken for baselines greater than a few tens of meters. Simulations further show that there is the potential for improving the seeing or astrometric accuracy by about 30%-50% on average, if the path length fluctuations in the lower component are directly calibrated. Statistical and systematic effects induce an error of about 15 m in the estimate of the lower component turbulent altitude.

Treuhaft, Robert N.↗

Mathematical Inversion of Lightning Data: Techniques and Applications

A survey of some interesting mathematical inversion studies dealing with radio, optical, and electrostatic measurements of lightning are presented. A discussion of why NASA is interested in lightning, what specific physical properties of lightning are retrieved, and what mathematical techniques are used to perform the retrievals are discussed. In particular, a relatively new multi-station VHF time-of-arrival (TOA) antenna network is now on-line in Northern Alabama and will be discussed. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The LMA supports on-going ground-validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. The LMA also provides detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and offers interesting comparisons with other meteorological/geophysical datasets. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. A new channel mapping retrieval algorithm is introduced for this purpose. To characterize the spatial distribution of retrieval errors, the algorithm has been applied to analyze literally tens of millions of computer-simulated lightning VHF point sources that have been placed at various ranges, azimuths, and altitudes relative to the LMA network. Statistical results are conveniently summarized in high-resolution, color-coded, error maps.

Koshak, William↗

Influence of Satellite-Based Heterogeneous Vegetation Momentum Roughness on Mesoscale Model Dynamics During IHOP 2002

The sensitivity of mesoscale weather prediction model to a vegetation roughness initialization is investigated for the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database previously derived using a physically based procedure and MODIS imagery, and iii) a MODIS climatologic roughness database that possesses the same spatial heterogeneity as (i) but with mean land class values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations from meteorological stations within the Oklahoma mesonet and surrounding region during IHOP20O2. A sensitivity analysis on the impact the MODIS-based roughness fields is also made through a time-series intercomparison of temperature bias, probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE) the results that, for the current replacement of the standard land-cover type based roughness values with the satellite-derived fields statistically improves model performance for most of the observed variables. Further, the satellite-based roughness enhances the surface wind speed, PBL height and TKE production on the order of 3 to l0 percent, with a lesser effect over grassland and cropland domains, and the greater effect over mixed land cover domains

Jasinski, Michael↗

Gridded National Inventory of U.S. Methane Emissions

We present a gridded inventory of US anthropogenic methane emissions with 0.1 deg x 0.1 deg spatial resolution, monthly temporal resolution, and detailed scale dependent error characterization. The inventory is designed to be onsistent with the 2016 US Environmental Protection Agency (EPA) Inventory of US Greenhouse Gas Emissionsand Sinks (GHGI) for 2012. The EPA inventory is available only as national totals for different source types. We use a widerange of databases at the state, county, local, and point source level to disaggregate the inventory and allocate the spatial and temporal distribution of emissions for individual source types. Results show large differences with the EDGAR v4.2 global gridded inventory commonly used as a priori estimate in inversions of atmospheric methane observations. We derive grid-dependent error statistics for individual source types from comparison with the Environmental Defense Fund (EDF) regional inventory for Northeast Texas. These error statistics are independently verified by comparison with the California Greenhouse Gas Emissions Measurement (CALGEM) grid-resolved emission inventory. Our gridded, time-resolved inventory provides an improved basis for inversion of atmospheric methane observations to estimate US methane emissions and interpret the results in terms of the underlying processes.

Methane↗

Introduction to Data Assimilation

Atmospheric data assimilation is a class of techniques used for producing descriptions of fields of air temperature, pressure, humidity, wind, etc. on a spatial grid or in terms of a finite functional representation. These are then used to initialize numerical weather forecasts or to analyze the atmosphere for other purposes. The techniques combine past, present, and even future observations in an approximate statistically optimal way. Various types of statistical or physically-based models and their corresponding adjoints are employed to relate diverse fields in both time and space and to relate what is observed to what is being analyzed. Computationally, the problem is very demanding and onstraining on the techniques that can be employed on a routine basis.

Atmospheric data assimilation↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

IRA Energy Community Data Layers

Data, geospatial data resources, and the linked mapping tool and web services reflect data for two types of potentially qualifying energy communities: 1) Census tracts and directly adjoining tracts that have had coal mine closures since 1999 or coal-fired electric generating unit retirements since 2009. These census tracts qualify as energy communities. 2) Metropolitan statistical areas (MSAs) and non-metropolitan statistical areas (non-MSAs) that are energy communities for 2023 and 2024, along with their fossil fuel employment (FFE) status. Additional information on energy communities and related tax credits can be accessed on the Interagency Working Group on Coal & Power Plant Communities & Economic Revitalization Energy Communities website (https://energycommunities.gov/energy-community-tax-credit-bonus/). Use limitations: these spatial data and mapping tool may not be relied upon by taxpayers to substantiate a tax return position or for determining whether certain penalties apply and will not be used by the IRS for examination purposes. The mapping tool does not reflect the application of the law to a specific taxpayer’s situation, and the applicable Internal Revenue Code provisions ultimately control.

Census Tract↗

Wildfire Smoke Particle Properties and Evolution, From Space-Based Multi-Angle Imaging II: The Williams Flats Fire during the FIREX-AQ Campaign

Although the characteristics of biomass burning events and the ambient ecosystem determine emitted smoke composition, the conditions that modulate the partitioning of black carbon (BC) and brown carbon (BrC) formation are not well understood, nor are the spatial or temporal frequency of factors driving smoke particle evolution, such as hydration, coagulation, and oxidation, all of which impact smoke radiative forcing. In situ data from surface observation sites and aircraft field campaigns offer deep insight into the optical, chemical, and microphysical traits of biomass burning (BB) smoke aerosols, such as single scattering albedo (SSA) and size distribution, but cannot by themselves provide robust statistical characterization of both emitted and evolved particles. Data from the NASA Earth Observing System’s Multi-Angle Imaging SpectroRadiometer (MISR) instrument can provide at least a partial picture of BB particle properties and their evolution downwind, once properly validated. Here we use in situ data from the joint NOAA/NASA 2019 Fire Influence on Regional to Global Environments Experiment-Air Quality (FIREX-AQ) field campaign to assess the strengths and limitations of MISR-derived constraints on particle size, shape, light-absorption, and its spectral slope, as well as plume height and associated wind vectors. Based on the satellite observations, we also offer inferences about aging mechanisms effecting downwind particle evolution, such as gravitational settling, oxidation, secondary particle formation, and the combination of particle aggregation and condensational growth. This work builds upon our previous study, adding confidence to our interpretation of the remote-sensing data based on an expanded suite of in situ measurements for validation. The satellite and in situ measurements offer similar characterizations of particle property evolution as a function of smoke age for the 06 August Williams Flats Fire, and most of the key differences in particle size and absorption can be attributed to differences in sampling and changes in the plume geometry between sampling times. Whereas the aircraft data provide validation for the MISR retrievals, the satellite data offer a spatially continuous mapping of particle properties over the plume, which helps identify trends in particle property downwind evolution that are ambiguous in the sparsely sampled aircraft transects. The MISR data record is more than two decades long, offering future opportunities to study regional wildfire plume behavior statistically, where aircraft data are limited or entirely lacking.

biomass burning↗

Mapping Electrical Crosstalk in Pixelated Sensor Arrays

Electronic coupling effects such as Inter-Pixel Capacitance (IPC) affect the quantitative interpretation of image data from CMOS, hybrid visible and infrared imagers alike. Existing methods of characterizing IPC do not provide a map of the spatial variation of IPC over all pixels. We demonstrate a deterministic method that provides a direct quantitative map of the crosstalk across an imager. The approach requires only the ability to reset single pixels to an arbitrary voltage, different from the rest of the imager. No illumination source is required. Mapping IPC independently for each pixel is also made practical by the greater S/N ratio achievable for an electrical stimulus than for an optical stimulus, which is subject to both Poisson statistics and diffusion effects of photo-generated charge. The data we present illustrates a more complex picture of IPC in Teledyne HgCdTe and HyViSi focal plane arrays than is presently understood, including the presence of a newly discovered, long range IPC in the HyViSi FPA that extends tens of pixels in distance, likely stemming from extended field effects in the fully depleted substrate. The sensitivity of the measurement approach has been shown to be good enough to distinguish spatial structure in IPC of the order of 0.1%.

focal plane array↗

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantifying resolution in cosmological N- body simulations using self-similarity

ABSTRACT We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using abacus , we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 per cent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times, the lower cut-off scale on convergence evolves in comoving units as a−1/2 (a being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Nuclear rotation at the fission limit in Rf 254

A ground-state rotational band in the fissile nucleus Rf254 was observed for the first time. Levels up to spin 14 and excitation energy of 1.56 MeV were observed. The Rf254 nuclei were produced using the Pb206(Ti50,2n) fusion-evaporation reaction. It is the weakest reaction channel ever studied using in-beam γ-ray spectroscopic methods. The reaction products were separated from the beam in the Argonne gas-filled analyzer (AGFA). The Rf254 nuclei were implanted into a double-sided Si strip detector at the AGFA focal plane and tagged with subsequent ground-state spontaneous fission decays using temporal and spatial correlations. Prompt γ rays in coincidence with the Rf254 recoils were detected in the Gammasphere array of Ge detectors. In order to identify the ground-state rotational band in Rf254, a method for identifying rotational bands in low statistics γ-ray spectra was developed. The deduced Rf254 kinematic moment of inertia is smaller compared to neighboring even-even nuclei. This is most likely associated with a slightly lower quadrupole deformation and stronger pairing correlations in Rf254. The behavior of the moment of inertia as a function of rotational frequency is similar to that of the lighter N=150 isotones Fm250 and No252.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerical Evidence for Many-Body Localization in Two and Three Dimensions

Disorder and interactions can lead to the breakdown of statistical mechanics in certain quantum systems, a phenomenon known as many-body localization (MBL). Much of the phenomenology of MBL emerges from the existence of localized-bits, or l-bits, a set of conserved quantities that are spatially localized and binary (i.e., possess only ±1 eigenvalues). While MBL and l-bits are known to exist in one-dimensional systems, their existence in dimensions greater than one is a key open question. To tackle this question, we develop an algorithm that can find approximate binary l-bits in arbitrary dimensions by adaptively generating a basis of operators in which to represent the l-bit. We use the algorithm to study four models: the one-, two-, and three-dimensional disordered Heisenberg models and the two-dimensional disordered hard-core Bose-Hubbard model. For all four of the models studied, our algorithm finds high-quality l-bits at large disorder strength and rapid qualitative changes in the distributions of l-bits in particular ranges of disorder strengths, suggesting the existence of MBL transitions. Furthermore, these transitions in the one-dimensional Heisenberg model and two-dimensional Bose-Hubbard model coincide well with past estimates of the critical disorder strengths in these models which further validates the evidence of MBL-like behavior in the other two and three-dimensional models we examine. In addition to finding MBL-like behavior in higher dimensions, our algorithm can be used to probe MBL in various geometries and dimensionality.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MIC-DP: A Scalable Correlation-Aware Differential Privacy Framework for High-Dimensional Data

Conventional differential privacy (DP) assumes record independence, limiting effectiveness on real-world datasets with temporal, spatial, or structural correlations. These dependencies undermine privacy guarantees and degrade utility in domains like healthcare, IoT, and smart city analytics. We propose Maximum Information Correlated Differential Privacy (MIC-DP), a novel framework that dynamically calibrates noise based on statistical dependencies. MIC-DP uses the Maximum Information Coefficient (MIC) to capture both linear and nonlinear correlations without explicit modeling, enabling adaptive sensitivity adjustment and improved privacy–utility trade-offs. Evaluations on healthcare (MIMIC), demographic (ACI), and synthetic datasets show that MIC-DP reduces mean absolute error (MAE) by up to 5.2% under strict privacy budgets (ϵ≤1), with aggregate utility improvements reaching 18% across datasets and evaluation metrics. MIC-DP provides formal (ϵ,δ)-privacy guarantees, scales efficiently with feature count, and supports deployment in moderate-scale, privacy-sensitive applications. Its tunable performance and runtime efficiency make MIC-DP suitable for privacy-sensitive applications where low-latency analytics and strong privacy guarantees must coexist. These results demonstrate MIC-DP’s effectiveness as a correlation-aware solution for practical DP.

Yang, Wenjun [Univ. of Washington, Tacoma, WA (Uni↗

Assessing sensitivities of climate model weighting to multiple methods, variables, and domains in the south-central United States

Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.

54 ENVIRONMENTAL SCIENCES↗

The LANDSAT-1 multispectral scanner as a tool in the classification of inland lakes

Relationships between LANDSAT-1 multispectral scanner (MSS) data and the trophic status of a group of lakes in the north-northeastern part of the United States were studied by predicting the magnitudes of two trophic state indicators, estimating lake position on a multivariate trophic scale, and automatically classifying lakes according to their trophic state. Initially, the principal component ordination was employed with 100 lakes. MSS data for some 20 lakes was then extracted from computer-compatible tapes (CCT) using a binary marking technique. The output was in the form of descriptive statistics and photographic concatenations. Color ratios were incorporated into regression models for the prediction of Secchi disc transparency, chlorophyll a, and lake position on the tropic scale. Results indicate that the LANDSAT-1 system, although handicapped by low spectral and spatial resolutions as well as excessive cloud cover, can be used as a supplemental data source in lake survey programs.

Boland, D. H. P.↗