Results for “Spatial statistics” — Page 15

DOE OSTI2026

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Downscaling↗

DOE OSTI2026

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Environmental sciences↗

NASA NTRS1976

Urban land use monitoring from computer-implemented processing of airborne multispectral data

Machine processing techniques were applied to multispectral data obtained from airborne scanners at an elevation of 600 meters over central Indianapolis in August, 1972. Computer analysis of these spectral data indicate that roads (two types), roof tops (three types), dense grass (two types), sparse grass (two types), trees, bare soil, and water (two types) can be accurately identified. Using computers, it is possible to determine land uses from analysis of type, size, shape, and spatial associations of earth surface images identified from multispectral data. Land use data developed through machine processing techniques can be programmed to monitor land use changes, simulate land use conditions, and provide impact statistics that are required to analyze stresses placed on spatial systems.

Todd, W. J.↗

NASA NTRS1986

Statistics of the geomagnetic secular variation for the past 5Ma

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

NASA NTRS1988

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

DOE OSTI2025

Atacama Cosmology Telescope: DR6 gravitational lensing and SDSS BOSS cross-correlation measurement and constraints on gravity with the 𝐸 𝐺 statistic

We derive new constraints on the 𝐸 𝐺 statistic as a test of gravity, combining the cosmic microwave background (CMB) lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS BOSS) CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the 𝐸 𝐺 statistic can distinguish Λ⁢ CDM from alternative models of gravity. We find 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢1$^{+0.06}_{−0.05}$ for Atacama Cosmology Telescope (ACT) and CMASS data at 68.28% confidence level, and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢9$^{+0.14}_{−0.11}$ for the ACT and LOWZ. Systematic errors are estimated to be 3% and 4%, respectively. Including CMB lensing information from Planck PR4 results in 𝐸 𝐺 ⁡(𝑧 eff = 0.555) = 0.3⁢4$^{+0.05}_{−0.05}$ with CMASS and 𝐸 𝐺 ⁡(𝑧 eff = 0.316) = 0.4⁢3$^{+0.11}_{−0.09}$ with LOWZ. These are consistent with predictions for the Λ⁢ CDM model that best fits the Planck CMB anisotropy and SDSS BOSS baryon acoustic oscillations (BAO), where 𝐸$^{GR}_{𝐺⁡}$(𝑧 eff =0.555) =0.401 ± 0.005 for CMB lensing combined with CMASS and 𝐸$^{GR}_{𝐺}$⁡(𝑧 eff = 0.316) = 0.452 ± 0.005 combined with LOWZ. We also find 𝐸 𝐺 to be scale independent, with probability to exceed >5%, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Improving seasonal precipitation forecasts in the Western United States through statistical downscaling

Abstract Seasonal precipitation forecasts in the western United States are critical resources for water resource management, especially during winter. While current seasonal forecasting systems provide monthly precipitation forecasts operationally, their coarse resolution limits their effectiveness in capturing the localized precipitation patterns and snowpack conditions essential for water resource managers in the mountainous regions. Here, analog statistical downscaling is demonstrated as an effective approach to enhance the spatial resolution of operational seasonal forecasts provided by the North American Multi-Model Ensemble. Downscaling was performed by building an analog ‘library’, in which corresponding model forecasts and observed values during the training period were stored. In the testing period, unseen model forecasts referenced the closest historical forecast from the analog library and applied the corresponding observational value for each point. This analysis indicates that downscaled products can capture localized features more accurately than the original coarse resolution forecasts, reducing forecast error across the western United States. Moreover, downscaling individual ensemble members—rather than downscaling the ensemble mean—further reduces forecasting error for their multi-model ensemble mean products. The greatest error reductions in the downscaled product, measured by root mean squared error (RMSE), were observed at low to mid-elevations (500–2000 meters), with 50%–70% improvement relative to the original forecasts. In the higher elevations (2000 meters and above), changes in RMSE relative to the original forecast were limited to 10%–30% improvements. The improvement is more substantial for forecast systems with 10 ensemble members compared to that with 4 members, but this relationship does not hold for the system with 24 ensemble members. These findings show that analog statistical downscaling can effectively address the spatial limitations of seasonal precipitation forecasts with minimal computational cost, providing a valuable framework for enhancing coarse resolution forecasting products while providing insights into the timing of ensemble mean calculations during the downscaling process.

Vernon, B. (ORCID:0009000891670689)↗

NASA NTRS2010

Optical and Microphysical Retrievals of Marine Stratocumulus Clouds off the Coast of Namibia from Satellite and Aircraft

Though the emphasis of the Southern Africa Regional Science Initiative 2000 (SAFARI-2000) dry season campaign was largely on emission sources and transport, the assemblage of aircraft (including the high altitude NASA ER-2 remote sensing platform and the University of Washington CV-580, UK MRF C-130, and South African Weather Bureau JRA in situ aircrafts) provided a unique opportunity for cloud studies. Therefore, as part of the SAFARI initiative, investigations were undertaken to assess regional aerosol-cloud interactions and cloud remote sensing algorithms. In particular, the latter part of the experiment concentrated on marine boundary layer stratocumulus clouds off the southwest coast of Africa. Associated with cold water upwelling along the Benguela current, the Namibian stratocumulus regime has received limited attention but appears to be unique for several reasons. During the dry season, outflow of continental fires and industrial pollution over this area can be extreme. From below, upwelling provides a rich nutrient source for phytoplankton (a source of atmospheric sulfur through DMS production as well as from decay processes). The impact of these natural and anthropogenic sources on the microphysical and optical properties of the stratocumulus is unknown. Continental and Indian Ocean cloud systems of opportunity were also studied during the campaign. SAFARI 2000 aircraft flights off the coast of Namibia were coordinated with NASA Terra Satellite overpasses for synergy with the Moderate Resolution Imaging Spectroradiometer (MODIS) and other Terra instruments. MODIS was developed by NASA and launched onboard the Terra spacecraft on December 18, 1999 (and Aqua spacecraft on May 4, 2002). Among the remote sensing algorithms developed and applied to this sensor are cloud optical and microphysical properties that include cloud thermodynamic phase, optical thickness, and effective particle radius of both liquid water and ice clouds. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, and fundamental atmospheric research. The archived MODIS Collection 5 cloud products processing stream will be used to analyze low water cloud scenes off the Namibian and Angolan coasts during SAFARI 2000 time period, as well as other years. Pixel-level Terra and Aqua MODIS retrievals (l. km spatial resolution at nadir) and gridded (1' uniform grid) statistics of cloud optical thickness and effective particle radius will be presented, including joint probability distributions between the two quantities. In addition, perspectives from the MODIS Airborne Simulator, which flew on the ER-2 during SAFARI 2000 providing high spatial resolution retrievals (50 m at nadir), will be presented as appropriate. The H-SAF Program requires an experimental operational European-centric Satellite Precipitation Algorithm System (E-SPAS) that produces medium spatial resolution and high temporal resolution surface rainfall and snowfall estimates over the Greater European Region including the Greater Mediterranean Basin. Currently, there are various types of experimental operational algorithm methods of differing spatiotemporal resolutions that generate global precipitation estimates. This address will first assess the current status of these methods and then recommend a methodology for the H-SAF Program that deviates somewhat from the current approach under development but one that takes advantage of existing techniques and existing software developed for the TRMM Project and available through the public domain.

Platnick, Steven E.↗

NASA NTRS1973

Far field effects of weak spatially periodic inhomogeneities.

Small variations in the properties of an acoustic medium produce cumulative effects that significantly alter wave forms propagating over long distances. To investigate some aspects of these effects, without statistical models, the one-dimensional wave equation with a weak spatially periodic variation in the sound speed is treated. An approximation is developed valid for such weak inhomogeneities over large distances based upon the approximate treatment of a linear differential equation with a periodic coefficient. It is applied to the propagation of an initial step wave and identifies the dispersive nature of the medium which accounts for the cumulative effects.

Friedman, M. B.↗

NASA NTRS2004

Spatially Complete Global Surface Albedos Derived from Terra/MODIS Data

Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. Recent production of land surface anisotropy, diffuse bihemispherical (white-sky) albedo and direct beam directional hemispherical (black-sky) albedo from observations acquired by the MODIS instruments aboard NASA s Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface's radiative characteristics. Cloud cover, which cutails retrievals, and the presence of ephemeral and seasonal snow limit the snow-free data to approximately half the global land surfaces on an annual equal-angle basis. This precludes the MOD43B3 albedo products from being used in some remote sensing and ground-based applications, climate models, and global change research projects. An ecosystem-dependent temporal interpolation technique is described that has been developed to fill missing or seasonally snow-covered data in the official MOD43B3 albedo product. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data. The resulting snow-free value-added products provide the scientific community with spatially and temporally complete global white- and black-sky surface albedo maps and statistics. These products are stored on 1'(approximately 10 km) and coarser resolution equal-angle grids, and are computed for the first seven MODIS wavelengths, ranging from 0.47 through 2.1 microns, and for three broadband wavelengths, 0.3-0.7,0.3-5.0 and 0.7-5.0 microns.

King, Michael D.↗

NASA NTRS1992

Multi-resolution processing for fractal analysis of airborne remotely sensed data

Fractal geometry is increasingly becoming a useful tool for modeling natural phenomenon. As an alternative to Euclidean concepts, fractals allow for a more accurate representation of the nature of complexity in natural boundaries and surfaces. Since they are characterized by self-similarity, an ideal fractal surface is scale-independent; i.e. at different scales a fractal surface looks the same. This is not exactly true for natural surfaces. When viewed at different spatial resolutions parts of natural surfaces look alike in a statistical manner and only for a limited range of scales. Images acquired by NASA's Thermal Infrared Multispectral Scanner are used to compute the fractal dimension as a function of spatial resolution. Three methods are used to determine the fractal dimension - Schelberg's line-divider method, the variogram method, and the triangular prism method. A description of these methods and the results of applying these methods to a remotely-sensed image is also presented. Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km (high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 1 Feb. 1992. Radiosonde data was collected for that duration to profile the characteristics of the atmosphere. This corresponds to 3 different pixel sizes - 5m, 15m, and 30m. After, simulating different spatial sampling intervals within the same image for each of the 3 image sets, the results are cross-correlated to compare the extent of detail and complexity that is obtained when data is taken at lower spatial intervals.

Jaggi, S.↗

NASA NTRS1994

An X-ray study of five supernova remnants in the Carina spiral arm

The ROSAT Position Sensitive Proportional Counter (PSPC) is used to perform an exploratory study of four fields in the Carina spiral arm containing five radio supernova remnants, only one of which has previously been studied in X-rays. We present upper limits for the detection in X-rays of G298.5 - 0.3, G298.6 - 0.0, and G299.0+0.2, and report a 4 sigma detection of G296.8-0.3. In addition, we present detailed spatial and spectral analysis of the bright X-ray remnant G296.1-0.7, which has previously been studied by both the Einstein IPC and EXOSAT LE/CMA. We detect relatively slight, but statistically significant, variations in the spectrum across the remnant via spatially resolved spectral fits and a study of the spatial variation of hardness ratios. In general, the spectrum is characteristic of a thermal plasma with kT about 0.2 keV and N(sub H) about 1.5 x 10(exp 21/sq. cm). The total X-ray emitting mass is estimated to be about 250 solar mass for an optically estimated distance of 4 kpc to the remnant. At this distance, the linear dimensions of the remnant are roughly 35 - 50 pc, implying an age on the order of 20,000 yr. Assuming that X-ray and radio brightnesses are related by SIGMA(sub R) proportional to SIGMA(exp 0.69)(sub X) and that the four radio remnants have X-ray spectral characteristics similar to G296.1-0.7, we find that the column densities to these sources must be several times 10(exp 22)/sq cm in order to explain their low X-ray count rates. This column density is considerably in excess of the X-ray fitted column density to G296.1-0.7, but is comparable to the total column densities in H I measured via the 21 cm line in the directions to all five remnants. This implies that G296.1 - 0.7 is at a significantly smaller distance than the other remnants.

Hwang, Una↗

DOE OSTI2026

Short-term electricity load forecasting: Application-driven evaluation of machine learning models across spatial and temporal scales

As we transition towards a decarbonized economy, the integration of variable renewable energy resources and new demands (e.g., electric vehicles, heat pumps) into the electricity grid places unprecedented pressure on grid operators to effectively anticipate and manage peak load. In this context, machine learning algorithms are proving to be indispensable for accurate short-term load forecasting, a crucial task to address these challenges. This study benchmarks 6 machine learning algorithms, including three neural networks and three tree-based algorithms, across various levels of spatial aggregation and time horizons (1, 4, 8, 24, and 48 h). The central contribution of this work is the comparison and analysis of load forecasting models not only based on statistical metrics, but also based on a novel error metric, which evaluates the cost implications of forecast errors for power system stakeholders. Results show that tree-based models outperform neural networks, based on statistical metrics, and yield less skewed error distributions for most spatial scales. However, through the lens of the novel error metric, neural networks are the more competitive choice, especially for forecast horizons that exceed 8 h. The study concludes with actionable recommendations to grid operators and highlights the need for the development of error metrics that link forecasting accuracy to operational costs. To promote transparency and open science, the datasets and Python code are open-sourced via a supplementary repository.

Houben, Nikolaus↗

NASA NTRS1976

Wheat signature modeling and analysis for improved training statistics

The author has identified the following significant results. The spectral, spatial, and temporal characteristics of wheat and other signatures in LANDSAT multispectral scanner data were examined through empirical analysis and simulation. Irrigation patterns varied widely within Kansas; 88 percent of wheat acreage in Finney was irrigated and 24 percent in Morton, as opposed to less than 3 percent for western 2/3's of the State. The irrigation practice was definitely correlated with the observed spectral response; wheat variety differences produced observable spectral differences due to leaf coloration and different dates of maturation. Between-field differences were generally greater than within-field differences, and boundary pixels produced spectral features distinct from those within field centers. Multiclass boundary pixels contributed much of the observed bias in proportion estimates. The variability between signatures obtained by different draws of training data decreased as the sample size became larger; also, the resulting signatures became more robust and the particular decision threshold value became less important.

Nalepka, R. F.↗

NASA NTRS2004

Application of Geostatistical Simulation to Enhance Satellite Image Products

With the deployment of Earth Observing System (EOS) satellites that provide daily, global imagery, there is increasing interest in defining the limitations of the data and derived products due to its coarse spatial resolution. Much of the detail, i.e. small fragments and notches in boundaries, is lost with coarse resolution imagery such as the EOS MODerate-Resolution Imaging Spectroradiometer (MODIS) data. Higher spatial resolution data such as the EOS Advanced Spaceborn Thermal Emission and Reflection Radiometer (ASTER), Landsat and airborne sensor imagery provide more detailed information but are less frequently available. There are, however, both theoretical and analytical evidence that burn scars and other fragmented types of land covers form self-similar or self-affine patterns, that is, patterns that look similar when viewed at widely differing spatial scales. Therefore small features of the patterns should be predictable, at least in a statistical sense, with knowledge about the large features. Recent developments in fractal modeling for characterizing the spatial distribution of undiscovered petroleum deposits are thus applicable to generating simulations of finer resolution satellite image products. We will present example EOS products, analysis to investigate self-similarity, and simulation results.

Hlavka, Christine A.↗

DOE OSTI2026

Structured illumination for surface-resolved grazing-incidence X-ray scattering

Grazing-incidence (GI) scattering techniques are widely used to characterize thin films, offering high surface sensitivity and insight into morphology and structure. However, these approaches typically provide statistical averaged information due to elongated footprint or limited spatial resolution due to beam size. Here we introduce a method that combines structured illumination with GI X-ray scattering and leverages our computational imaging approach to resolve local structural details. We demonstrate that our method captures local features of an organic semiconductor thin film without the need for sample rotation as in tomography. The method expands GI techniques from statistical averaging to high-resolution imaging, thereby providing the capability for detailed analysis of local material properties, such as domain shape, orientation and polymorphism, which are critical for advancing material design towards more efficient and tailored materials.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1983

The spatial correlation function of rich clusters of galaxies

A series of objective statistical estimators is applied to directly study the three-dimensional distribution of rich clusters, using a recently completed redshift sample of 104 Abell clusters out to distance class four and galactic latitude of 30 deg or more. The sample and the relevant positional and redshift selection effects are discussed, and the two-point angular correlation function and nearest-neighbor distribution of clusters in the sample are determined, as well as in the larger D = 5+6 sample and in several subsamples of different richness classes and latitude zones. The redshift correlation function of clusters for various angular separations on the sky is examined as a test for the reality of the observed angular correlations. The two-point spatial correlation function of clusters in the sample is determined, and the consistency among the three observed correlation functions is tested. A model that describes all the observed correlations is presented.

Bahcall, N. A.↗

NASA NTRS2016

Tropospheric Ozone Change from 1980 to 2010 Dominated by Equatorward Redistribution of Emissions

Ozone is an important air pollutant at the surface, and the third most important anthropogenic greenhouse gas in the troposphere. Since 1980, anthropogenic emissions of ozone precursors methane, non-methane volatile organic compounds, carbon monoxide and nitrogen oxides (NOx) have shifted from developed to developing regions. Emissions have thereby been redistributed equatorwards, where they are expected to have a stronger effect on the tropospheric ozone burden due to greater convection, reaction rates and NOx sensitivity. Here we use a global chemical transport model to simulate changes in tropospheric ozone concentrations from 1980 to 2010, and to separate the influences of changes in the spatial distribution of global anthropogenic emissions of short-lived pollutants, the magnitude of these emissions, and the global atmospheric methane concentration. We estimate that the increase in ozone burden due to the spatial distribution change slightly exceeds the combined influences of the increased emission magnitude and global methane. Emission increases in Southeast, East and South Asia may be most important for the ozone change, supported by an analysis of statistically significant increases in observed ozone above these regions. The spatial distribution of emissions dominates global tropospheric ozone, suggesting that the future ozone burden will be determined mainly by emissions from low latitudes.

Zhang, Yuqiang↗