Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5

Mesoscale Raised Rim Depressions (MRRDs) on Earth: A Review of the Characteristics, Processes, and Spatial Distributions of Analogs for Mars

Fields of mesoscale raised rim depressions (MRRDs) of various origins are found on Earth and Mars. Examples include rootless cones, mud volcanoes, collapsed pingos, rimmed kettle holes, and basaltic ring structures. Correct identification of MRRDs on Mars is valuable because different MRRD types have different geologic and/or climatic implications and are often associated with volcanism and/or water, which may provide locales for biotic or prebiotic activity. In order to facilitate correct identification of fields of MRRDs on Mars and their implications, this work provides a review of common terrestrial MRRD types that occur in fields. In this review, MRRDs by formation mechanism, including hydrovolcanic (phreatomagmatic cones, basaltic ring structures), sedimentological (mud volcanoes), and ice-related (pingos, volatile ice-block forms) mechanisms. For each broad mechanism, we present a comparative synopsis of (i) morphology and observations, (ii) physical formation processes, and (iii) published hypothesized locations on Mars. Because the morphology for MRRDs may be ambiguous, an additional tool is provided for distinguishing fields of MRRDs by origin on Mars, namely, spatial distribution analyses for MRRDs within fields on Earth. We find that MRRDs have both distinguishing and similar characteristics, and observation that applies both to their mesoscale morphology and to their spatial distribution statistics. Thus, this review provides tools for distinguishing between various MRRDs, while highlighting the utility of the multiple working hypotheses approach.

Burr, Devon M.↗

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗

Magnetic Helicity Condensation and the Solar Cycle

Solar filaments exhibit a global chirality pattern where dextral/sinistral filaments, corresponding to negative/ positive magnetic helicity, are dominant in the northern/southern hemisphere. This pattern is opposite to the sign of magnetic helicity injected by differential rotation along east–west oriented polarity inversion lines, posing a major conundrum for solar physics. A resolution of this problem is offered by the magnetic helicity-condensation model of Antiochos. To investigate the global consequences of helicity condensation for the hemispheric chirality pattern, we apply a temporally and spatially averaged statistical approximation of helicity condensation. Realistic magnetic field configurations in both the rising and declining phases of the solar cycle are simulated. For the helicity-condensation process, we assume convective cells consisting of positive/negative vorticities in the northern/southern hemisphere that inject negative/positive helicity. The magnitude of the vorticity is varied as a free parameter, corresponding to different rates of helicity injection. To reproduce the observed percentages of dominant and minority filament chiralities, we find that a vorticity of magnitude 2.5 x 10(exp −6)/s is required. This rate, however, is insufficient to produce the observed unimodal profile of chirality with latitude. To achieve this, a vorticity of at least 5 x 10(exp −6)/s is needed. Our results place a lower limit on the small-scale helicity injection required to dominate differential rotation and reproduce the observed hemispheric pattern. Future studies should aim to establish whether the helicity injection rate due to convective flows and/or flux emergence across all latitudes of the Sun is consistent with our results.

Mackay, Duncan H.↗

Geostationary Lightning Mapper Flash Characteristics of Electrified Snowfall Events

This study examines characteristics of lightning in snowfall events (i.e., thundersnow, 20 TSSN) from the perspective of the Geostationary Lightning Mapper (GLM) and the National Environmental Satellite Data and Information Service (NESDIS) merged Snowfall Rate (mSFR) product. A thundersnow detection algorithm (TDA) was derived from the GLM and mSFR which resulted in a probability of detection (POD) of 66.7% when compared to the Meteorological Terminal Air Report (METAR) reports of TSSN. However, using the TDA an additional 2,175 lightning flashes within detected snowfall were identified that were not observed by the METAR reports, indicating that TSSN has been under reported in previous literature. TSSN flashes observed by GLM have mean flash areas, durations, and total optical energy outputs of 754 km2, 402 ms, and 1,342 fJ, which are between the 50th and 99th percentile values for all flashes within the GLM field of view. A comparison with data from the National Lightning Detection Network (NLDN) indicated that the NLDN had at least one cloud or ground flash detection in 1,709 of the 2,214 flashes observed by GLM in snowfall. An average of 5.85 NLDN flashes were assigned to a single GLM flash when the NLDN flash data were constrained by the GLM flash duration and spatial footprint. Statistically significant (p<0.01) differences in flash area and flash energy were found between flashes that were observed by the NLDN and those that were not. Additionally, when GLM was combined with the NLDN, at least 11.1% of flashes involved a tall human-made object like an antenna or wind turbine.

lightning↗

Evaluating Rock Pool Hydroperiod Fluctuation using Climate Variables to Inform Habitat Monitoring and Protection in the Western Sonoran Desert

As warming and drying trends continue to impact the greater American Southwest, effective ecosystem management increasingly relies upon understanding the relationships between climate and water resources. Ephemeral freshwater rock pools, known as tinajas, have great ecological and cultural importance as some of the only sources of surface water in the western Sonoran Desert (WSD). Tinaja flooding and drying cycles, known as hydroperiods, vary based on meteorologic and climatologic conditions; however, a lack of extensive research relating climatic impacts to tinajas puts these critical ecosystems further at risk. Tinajas throughout the WSD are monitored by the National Park Service (NPS) using resource intensive strategies including in situ trail camera observation and direct measurement. To aid NPS monitoring efforts, this research used remotely sensed climate data to analyze spatiotemporal climate trends and relationships between climate variables and tinaja hydroperiods in the WSD between 1979–2022. Using Aqua and Terra Moderate Resolution Imaging Spectroradiometers (MODIS), University of Idaho Gridded Surface Meteorological Dataset (gridMET), and OpenET data, the project analyzed land surface temperature, evapotranspiration, precipitation, wind velocity, and solar radiation. The team generated climate anomaly time series and climate normal maps for the WSD, identifying statistically significant spatial and temporal trends. These data were then compared to daily, qualitative in situ hydroperiod observations taken between 2019–2022. This work will be used by the NPS to inform the monitoring and protection of tinajas in the WSD. Findings contribute to a limited body of research concerning climate and tinajas, which are often overlooked despite their disproportionate ecological importance.

Annie Britton↗

Algorithm for Identifying Erroneous Rain-Gauge Readings

An algorithm analyzes rain-gauge data to identify statistical outliers that could be deemed to be erroneous readings. Heretofore, analyses of this type have been performed in burdensome manual procedures that have involved subjective judgements. Sometimes, the analyses have included computational assistance for detecting values falling outside of arbitrary limits. The analyses have been performed without statistically valid knowledge of the spatial and temporal variations of precipitation within rain events. In contrast, the present algorithm makes it possible to automate such an analysis, makes the analysis objective, takes account of the spatial distribution of rain gauges in conjunction with the statistical nature of spatial variations in rainfall readings, and minimizes the use of arbitrary criteria. The algorithm implements an iterative process that involves nonparametric statistics.

Rickman, Doug↗

A new technique for enlargement and reconstruction of digital sensor imagery

A new image enlargement and reconstruction routine, the digital image enlarging balanced reconstruction algorithm (DIEBRA), is described. Using a highly modified form of 'balanced' two-dimensional polynomial interpolation, this program enlarges digital imagery, creating synthetic high spatial resolution images. Statistical analysis shows the DIEBRA-generated imagery to be significantly closer to true high spatial resolution imagery at all frequencies than imagery generated by a cubic convolutional reconstruction filter.

Rhea, W. Joseph↗

Classifying Urban Land Covers Using Local Indices of Spatial Complexity

The skewed statistical distributions of land cover types in complex, heterogeneous urban areas limits the effectiveness of traditional spectrally based maximum-likelihood classifiers. This work examines the utility of fractal dimension and Moran's I index of spatial autocorrelation in segmenting high-resolution panchromatic and lower-resolution multispectral imagery. Tools available in the Image Characterization and Modeling System (ICAMS) were used to analyze multi-temporal and multi-platform imagery of Atlanta, Georgia. In this example, land cover change trajectories from forest or grassland to built up land covers lead to decreased spatial autocorrelation. In lower resolution imagery such as Landsat MSS, the complex details of forested land covers and urbanized areas are smoothed, and texture-based change detection is less effective. Although segmentation of panchromatic images is possible using fractal dimension or Moran's I, widely differing land covers often yield similar values of these indices. Better results are obtained when a surface of local fractal dimension or spatial autocorrelation is combined as an additional layer in a supervised maximum-likelihood multispectral classification.

Arumugam, Mahesh↗

Comment on “Advanced Testing of Low, Medium, and High ECS CMIP6 GCM Simulations Versus ERA5-T2m” by N. Scafetta (2022)

Scafetta (2022, https://doi.org/10.1029/2022gl097716) purports to test Coupled Model Intercomparison Project Phase 6 (CMIP6) climate models through a comparison of temperature changes over three decades. Unfortunately, the paper contains numerous conceptual and statistical errors that undermine all of the conclusions. First, no uncertainty is given for the observational temperature difference, making it impossible to assess compatibility with any model result. Second, the CMIP6 data are the ensemble means for each model, but the metric being tested is sensitive to the internal variability and so the full ensemble for each model must be used. When this is corrected, the conclusion that “all models with ECS > 3.0°C overestimate the observed global surface warming” is not sustained. Third, the statistical test in Section 2 would reject all models even in a perfect model setup given sufficient ensemble members, thus the second conclusion “that spatial t-statistics rejects the data-model agreement” is also not sustainable.

CMIP6↗

Downscaling Satellite Precipitation with Emphasis on Extremes: A Variational 1-Norm Regularization in the Derivative Domain

The increasing availability of precipitation observations from space, e.g., from the Tropical Rainfall Measuring Mission (TRMM) and the forthcoming Global Precipitation Measuring (GPM) Mission, has fueled renewed interest in developing frameworks for downscaling and multi-sensor data fusion that can handle large data sets in computationally efficient ways while optimally reproducing desired properties of the underlying rainfall fields. Of special interest is the reproduction of extreme precipitation intensities and gradients, as these are directly relevant to hazard prediction. In this paper, we present a new formalism for downscaling satellite precipitation observations, which explicitly allows for the preservation of some key geometrical and statistical properties of spatial precipitation. These include sharp intensity gradients (due to high-intensity regions embedded within lower-intensity areas), coherent spatial structures (due to regions of slowly varying rainfall),and thicker-than-Gaussian tails of precipitation gradients and intensities. Specifically, we pose the downscaling problem as a discrete inverse problem and solve it via a regularized variational approach (variational downscaling) where the regularization term is selected to impose the desired smoothness in the solution while allowing for some steep gradients(called 1-norm or total variation regularization). We demonstrate the duality between this geometrically inspired solution and its Bayesian statistical interpretation, which is equivalent to assuming a Laplace prior distribution for the precipitation intensities in the derivative (wavelet) space. When the observation operator is not known, we discuss the effect of its misspecification and explore a previously proposed dictionary-based sparse inverse downscaling methodology to indirectly learn the observation operator from a database of coincidental high- and low-resolution observations. The proposed method and ideas are illustrated in case studies featuring the downscaling of a hurricane precipitation field.

Hurricanes↗

Development of a “Nature Run” for Observing System Simulation Experiments (OSSEs) for Snow Mission Development

Snow is a fundamental component of global and regional water budgets, particularly in mountainous areas and regions downstream that rely on snowmelt for water resources. Land surface models (LSMs) are commonly used to develop spatially distributed estimates of snow water equivalent (SWE) and runoff. However, LSMs are limited by uncertainties in model physics and parameters, among other factors. In this study, we describe the use of model calibration tools to improve snow simulations within the Noah-MP LSM as the first step in an observing system simulation experiment (OSSE). Noah-MP is calibrated against the University of Arizona (UA) SWE product over a western Colorado domain. With spatially varying calibrated parameters, we run calibrated and default Noah-MP simulations for water years 2010–20. By evaluating both simulations against the UA dataset, we show that calibration decreases domain averaged temporal RMSE and bias for snow depth from 0.15 to 0.13 m and from −0.036 to −0.0023 m, respectively, and improves the timing of snow ablation. Increased snow simulation performance also improves estimates of model-simulated runoff in four of six study basins, though only one has statistically significant improvement. Spatially distributed Noah-MP snow parameters perform better than default uniform values. We demonstrate that calibrating variables related to snow albedo calculations and rain–snow partitioning, among other processes, is a necessary step for creating a nature run that reasonably approximates true snow conditions for the OSSEs. Additionally, the inclusion of a snowfall scaling term can address biases in precipitation from meteorological forcing datasets, further improving the utility of LSMs for generating reliable spatiotemporal estimates of snow.

Melissa L. Wrzesien↗

Development and Evaluation of High-Resolution Climate Simulations Over the Mountainous Northeastern United States

The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.

climatology↗

The Use of Gridded Fossil Fuel CO2 Emissions (FFCO2) Inventory for Climate Mitigation Applications: Errors, Uncertainties, and Current and Future Challenges

Emission Inventory (EI) is a fundamental tool to monitor global compliance of greenhouse gases (GHGs) emissions reduction actions. Inventory guidelines provide a best practice to help EI compilers to make comparable national emission estimates, in spite of the differences in data availability across countries and regions. There are a variety of sources of errors and uncertainties, however, that originate beyond what the inventory guidelines define. For example, spatially-explicit EIs, which are a key product for atmospheric modeling applications, are often developed for research purposes, and there are no specific guidelines to disaggregate emission estimates from country scale. On top of that, EIs are fundamentally prone to systematic biases due to the simple calculation methodology and thus an objective evaluation (e.g. atmospheric top-down estimates) is needed to assure the accuracy of the estimates. ODIAC is a global high-resolution (1x1 km) fossil fuel carbon dioxide (CO2) gridded EI that is now often used in atmospheric CO2 modeling. ODIAC is based on disaggregation of national emission estimates made by CDIAC, which is the well accepted standard in the community. The ODIAC emission data product is updated on an annual basis using best available statistical data. Subnational spatial emission patterns are estimated using power plant profiles and satellite-observations of nighttime lights. In addition to the conventional CDIAC gridded data product, ODIAC carries international bunker emissions (shipping and aviation), which allows flux inversion modelers to accurately impose the global total fossil fuel emissions and their horizontal and vertical distribution. We have extensively evaluated ODIAC emissions using fine-grained EIs as well as a high-resolution atmospheric model simulation across different scales (national, subnational/regional, and urban policy relevant) with a focus on the uncertainties associated with the emission disaggregation. We have examined the use of NASA's Black Marble Suomi-NPP/VIIRS nightlight data.

Oda, Tomohiro↗

Scales of North Atlantic wind stress curl determined from the comprehensive ocean-atmosphere data set

Nineteen years of wind data over the North Atlantic are used to calculate a field of wind stress curl. An empirical orthogonal function (EOF) analysis is performed on this field, resulting in spatial patterns of wind stress curl and associated time series. A Monte Carlo technique is used to establish the statistical significance of each spatial pattern, and the associated time series are spectrally analyzed. The first four statistically significant EOF modes represent more than 50 percent of the curl variance, and the spatial patterns of curl associated with these modes exhibit the major elements of North Atlantic climatology. Most of the time series spectral variance is contained in annual and semiannual frequencies. The features observed include the individual annual variation of the subtropical high and the subpolar low, the annual oscillation of intensity between pressure centers, the influence of localized strong SST gradients and associated cyclogenesis regions, and the constant nature of the trades.

Ehret, Laura L.↗

Atmospheric pressure, density, temperature and wind variations between 50 and 200 km

Data on atmospheric pressure, density, temperature and winds between 50 and 200 km were collected from sources including Meteorological Rocket Network data, ROBIN falling sphere data, grenade release and pitot tube data, meteor winds, chemical release winds, satellite data, and others. These data were analyzed by a daily difference method and results on the distribution statistics, magnitude, and spatial structure of the irregular atmospheric variations are presented. Time structures of the irregular variations were determined by the analysis of residuals from harmonic analysis of time series data. The observed height variations of irregular winds and densities are found to be in accord with a theoretical relation between these two quantities. The latitude variations (at 50 - 60 km height) show an increasing trend with latitude. A possible explanation of the unusually large irregular wind magnitudes of the White Sands MRN data is given in terms of mountain wave generation by the Sierra Nevada range about 1000 km west of White Sands. An analytical method is developed which, based on an analogy of the irregular motion field with axisymmetric turbulence, allows measured or model correlation or structure functions to be used to evaluate the effective frequency spectra of scalar and vector quantities of a spacecraft moving at any speed and at any trajectory elevation angle.

Justus, C. G.↗

Short and long periodic atmospheric variations between 25 and 200 km

Previously collected data on atmospheric pressure, density, temperature and winds between 25 and 200 km from sources including Meteorological Rocket Network data, ROBIN falling sphere data, grenade release and pitot tube data, meteor winds, chemical release winds, satellite data, and others were analyzed by a daily difference method and results on the distribution statistics, magnitude, and spatial structure of gravity wave and planetary wave atmospheric variations are presented. Time structure of the gravity wave variations were determined by the analysis of residuals from harmonic analysis of time series data. Planetary wave contributions in the 25-85 km range were discovered and found to have significant height and latitudinal variation. Long period planetary waves, and seasonal variations were also computed by harmonic analysis. Revised height variations of the gravity wave contributions in the 25 to 85 km height range were computed. An engineering method and design values for gravity wave magnitudes and wave lengths are given to be used for such tasks as evaluating the effects on the dynamical heating, stability and control of spacecraft such as the space shuttle vehicle in launch or reentry trajectories.

Justus, C. G.↗

User oriented ERTS-1 images

Photographic reproduction of ERTS-1 images are capable of displaying only a portion of the total information available from the multispectral scanner. Methods are being developed to generate ERTS-1 images oriented towards special users such as agriculturists, foresters, and hydrologists by applying image enhancement techniques and interactive statistical classification schemes. Spatial boundaries and linear features can be emphasized and delineated using simple filters. Linear and nonlinear transformations can be applied to the spectral data to emphasize certain ground information. An automatic classification scheme was developed to identify particular ground cover classes such as fallow, grain, rape seed or various vegetation covers. The scheme applies the maximum likelihood decision rule to the spectral information and classifies the ERTS-1 image on a pixel by pixel basis. Preliminary results indicate that the classifier has limited success in distinguishing crops, but is well adapted for identifying different types of vegetation.

Shlien, S.↗

The local complex of O and B stars. I - Distribution of stars and interstellar dust

The O-B5 stars, supergiants, young clusters, and associations within 1 kpc of the sun populate two flat systems inclined to each other by 19 to 22 deg. The historical background, statistical significance, composition, spatial arrangement of the contents, and interstellar extinction in the two belts are discussed. A more or less random distribution in space and in age characterizes the O-B5 stars of the 'galactic belt', which is aligned nearly along the Milky Way. The 'Gould belt' is inclined to the Milky Way (north in Sco-Oph and south in Orion), and exhibits a projected distribution of O-B5 stars in its mean plane that resembles a 'dragonfly', with five major features defining it. A crude 'diameter' of the system is 750 to 1000 pc, and the sun's position is eccentric, lying toward Ophiuchus. The nuclear age of the system, while not unique, may be characterized as 30 m.y. from the spectral type of the broad main-sequence turnup near B2-5. Most of the O-B2 stars and youngest stellar groups near the sun belong to the Gould belt, but both belts have approximately equal space densities of B3-B5 stars and similar average values of interstellar extinction.-

Stothers, R.↗