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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Monthly Oceanic Rainfall Derived from Space-Borne Microwave Data

A technique for estimating monthly oceanic rainfall rate using multi-channel microwave measurements has been developed. This technique is applied to the Special Sensor Microwave Imager (SSM/I) data since 1987 to infer monthly rainfall for the Global Precipitation Climatology Project (GPCP). A modified version of this algorithm is now being applied to the TRMM Microwave Imager (TMI) data. There are three prominent features of this algorithm. First, the knowledge of the form of the rainfall intensity probability density function used to augment the measurements. Second, utilizing a linear combination of the 19.35 and 22.235 GHz channels to deemphasize the effect of water vapor. Third, an objective technique has been developed to estimate the rain layer thickness from the 19.35 and 22.23 5 GHz brightness temperature histograms, TMI data with better spatial resolution and 24 hour sampling (vs. sun-synchronized sampling, which is limited to two narrow intervals of local solar time for DMSP satellites) prompt us to study the similarity and difference between these two rainfall estimates. One year of rainfall data (January to December 1998) are used in this study. The mean differences of these two rain estimates are typically less than 15%. Statistical tests are administrated to evaluate the differences between rainfall estimates from SSM/I and TMI data. Differences between TMI and SSM/I rain estimates are discussed in the context of radiative transfer calculations and rainfall estimate algorithm.

Chang, Alfred T. C.↗
NASA NTRS1986

M.I.T./Canadian vestibular experiments on the Spacelab-1 mission: 4. Space motion sickness: symptoms, stimuli, and predictability

Space sickness symptoms were observed by 4 specially trained observers on Spacelab-1. Three reported persistent symptoms, and vomited repeatedly during the first and/or second day of flight. Head movements on all axes were provocative, particularly in pitch and roll. Head acceleration data recorded from 2 symptomatic crewmen showed that after several hours of physical activity in orbit, symptoms appeared, and thereafter both crewmen were compelled to limit head movements. Firm body contact with motionless surfaces helped alleviate symptoms. When crewmembers floated into unfamiliar body orientations in the cabin, inherent ambiguities in static visual orientation cues sometimes produced spatial reorientation episodes which were also provocative. Symptoms largely resembled those of other forms of prolonged motion sickness, superimposed upon other symptoms attributable to fluid shift. All 4 eventually used anti-motion sickness drugs. When they did, vomiting frequency was reduced. By the 4th day, symptoms subsided, and head accelerations again increased in magnitude and variability. Sickness intensity in orbit was not predicted by statistically concordant results of 6 acute preflight susceptibility tests. However, results from a longer duration preflight prism goggles test showed an apparent correlation. All subjects were asymptomatic making head movements in parabolic flight 4 days after the mission, but not 1 year later. Overall, results support the view that space sickness is a motion sickness.

Non-NASA Center↗
NASA NTRS2010

Downscaling Soil Moisture in the Southern Great Plains Through a Calibrated Multifractal Model for Land Surface Modeling Applications

Accounting for small-scale spatial heterogeneity of soil moisture (theta) is required to enhance the predictive skill of land surface models. In this paper, we present the results of the development, calibration, and performance evaluation of a downscaling model based on multifractal theory using aircraft!based (800 m) theta estimates collected during the southern Great Plains experiment in 1997 (SGP97).We first demonstrate the presence of scale invariance and multifractality in theta fields of nine square domains of size 25.6 x 25.6 sq km, approximately a satellite footprint. Then, we estimate the downscaling model parameters and evaluate the model performance using a set of different calibration approaches. Results reveal that small-scale theta distributions are adequately reproduced across the entire region when coarse predictors include a dynamic component (i.e., the spatial mean soil moisture ) and a stationary contribution accounting for static features (i.e., topography, soil texture, vegetation). For wet conditions, we found similar multifractal properties of soil moisture across all domains, which we ascribe to the signature of rainfall spatial variability. For drier states, the theta fields in the northern domains are more intermittent than in southern domains, likely because of differences in the distribution of vegetation coverage. Through our analyses, we propose a regional downscaling relation for coarse, satellite-based soil moisture estimates, based on ancillary information (static and dynamic landscape features), which can be used in the study area to characterize statistical properties of small-scale theta distribution required by land surface models and data assimilation systems.

Mascaro, Giuseppe↗
NASA NTRS2018

Comparisons of Two Spatial Implementations of a Crop Model Using Remotely Sensed Observations over Southeastern United States

Global food security is one of the most pressing issues of the current century, particularly for developing nations. Agricultural simulation models can be a key component in testing new technologies, seeds and cultivars etc. However, inaccurate input information, model related errors and the mode of implementation can also add to model uncertainties. In this study, the crop model is implemented in two separate fashions: a)gridded (GriDSSAT model) and b) using random spatial ensembles (RHEAS model). This is done in the Southeastern US to evaluate and understand the modelperformance over a region data availabilities. Once the model performance is assessed, multiple satellite based earth observation parameters such as soil moisture, vegetation index etc. can be assimilated into crop models to reduce input and model related uncertainties particularly in data limited regions. In this study, the National Agricultural Statistical Services (NASS) reported yield data at county levels are used for comparison andvalidation purposes. The GriDSSAT model estimation of corn yields in comparison with the reported NASS yields showed an overall RMSD of nearly 3720 (kg/ha) whereas RMSD for the RHEAS model implementation was 3550 (kg/ha). Overall the GriDSSAT model had negative bias of nearly 2400 kg/ha (except for 2013) while RHEAS had a slight positive bias of 400 kg/ha (approx.).

SERVIR↗
NASA NTRS2021

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗
DOE OSTI2022

Computing Thermodynamic Properties of Fluids Augmented by Nanoconfinement: Application to Pressurized Methane

Nanoconfined fluids exhibit remarkably different thermodynamic behavior compared to the bulk phase. These confinement effects render predictions of thermodynamic quantities of nanoconfined fluids challenging. In particular, confinement creates a spatially varying density profile near the wall that is primarily responsible for adsorption and capillary condensation behavior. Significant fluctuations in thermodynamic quantities, inherent in such nanoscale systems, coupled to strong fluid–wall interactions give rise to this near-wall density profile. Empirical models have been proposed to explain and model these effects, yet no first-principles based formulation has been developed. We present a statistical mechanics framework that embeds such a coupling to describe the effect of the fluid–wall interaction in amplifying the near-wall density behavior for compressible gases at elevated pressures such as pressurized methane in confinement. We show that the proposed theory predicts accurately the adsorbed layer thickness as obtained with small-angle neutron scattering measurements. Furthermore, the predictions of density under confinement from the proposed theory are shown to be in excellent agreement with available experimental and atomistic simulations data for a range of temperatures for nanoconfined methane. While the framework is presented for evaluating the near-wall density, owing to its rigorous foundation in statistical mechanics, the proposed theory can also be generalized for predicting phase-transition and nonequilibrium transport of nanoconfined fluids.

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

Using Temporal Deep Learning Models to Estimate Daily Snow Water Equivalent Over the Rocky Mountains

Abstract In this study we construct and compare three different deep learning (DL) models for estimating daily snow water equivalent (SWE) from high‐resolution gridded meteorological fields over the Rocky Mountain region. To train the DL models, Snow Telemetry (SNOTEL) station‐based SWE observations are used as the prediction target. All DL models produce higher median Nash‐Sutcliffe Efficiency (NSE) values than a conceptual SWE model and interpolated gridded data sets, although mean squared errors also tend to be higher. Sensitivity of the SWE prediction to the model's input variables is analyzed using an explainable artificial intelligence (XAI) method, yielding insight into the physical relationships learned by the models. This method reveals the dominant role precipitation and temperature play in snowpack dynamics. In applying our models to estimate SWE throughout the Rocky Mountains, an extrapolation problem arises since the statistical properties of SWE (e.g., annual maximum) and geographical properties of individual grid points (e.g., elevation) differ from the training data. This problem is solved by normalizing the SWE with its historical maximum value to alleviate extrapolation for all tested DL models. Our work shows that the DL models are promising tools for estimating SWE, and sufficiently capture relevant physical relationships to make them useful for spatial and temporal extrapolation of SWE values.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1991

Venus ionospheric tail rays - Spatial distributions and interplanetary magnetic field control

The overall properties of Venus ionospheric tail rays (such as density, spatial extent, and distribution) and their relationship to the draped magnetic field configuration behind the planet were investigated using measurements obtained by the Pioneer Venus Orbiter Langmuir probe, a magnetometer, and a plasma-wave detector. The results suggest that tail rays are a normal feature of the steady solar wind interaction with Venus and are not generally associated with a central tail plasma sheet. The statistics of the tail rays occurrence point toward the existence of a distributed terminator ionosphere source, consistent with findings of Brace et al. (1990).

Ong, M.↗
NASA NTRS1992

Detection of x ray sources in PROS

The problem of detecting discrete sources in x-ray images has much in common with the problem of automatic source detection at other wavelengths. In all cases, one searches for positive brightness enhancements exceeding a certain threshold, which appear consistent with what one expects for a point source, in the presence of a (possibly) spatially variable background. Multidimensional point spread functions (e.g., dependent on detector position and photon energy) are also common. At the same time, the problem in x-ray astronomy has some unique aspects. For example, for typical x-ray exposures in current or recent observatories, the number of available pixels far exceeds the number of actual x-ray events, so Poisson, rather than Gaussian statistics apply. Further, extended cosmic x-ray sources are common, and one often desires to detect point sources in the vicinity or even within bright, diffuse x-ray emission. Finally, support structures in x-ray detectors often cast sharp shadows in x-ray images making it necessary to detect sources in a region of rapidly varying exposure. We have developed a source detection package within the IRAF/PROS environment which attempts to deal with some of the problems of x-ray source detection. We have patterned our package after the successful Einstein Observatory x-ray source detection programs. However, we have attempted to improve the flexibility and accessibility of the functions and to provide a graphical front-end for the user. Our philosophy has been to use standard IRAF tasks whenever possible for image manipulation and to separate general functions from mission-specific ones. We will report on the current status of the package and discuss future developments, including simulation tasks, to allow the user to assess detection efficiency and source significance, tasks to determine source intensity, and alternative detection algorithms.

Deponte, J.↗
NASA NTRS2005

Small-Scale Drop-Size Variability: Empirical Models for Drop-Size-Dependent Clustering in Clouds

By analyzing aircraft measurements of individual drop sizes in clouds, it has been shown in a companion paper that the probability of finding a drop of radius r at a linear scale l decreases as l(sup D(r)), where 0 less than or equals D(r) less than or equals 1. This paper shows striking examples of the spatial distribution of large cloud drops using models that simulate the observed power laws. In contrast to currently used models that assume homogeneity and a Poisson distribution of cloud drops, these models illustrate strong drop clustering, especially with larger drops. The degree of clustering is determined by the observed exponents D(r). The strong clustering of large drops arises naturally from the observed power-law statistics. This clustering has vital consequences for rain physics, including how fast rain can form. For radiative transfer theory, clustering of large drops enhances their impact on the cloud optical path. The clustering phenomenon also helps explain why remotely sensed cloud drop size is generally larger than that measured in situ.

Marshak, Alexander↗
NASA NTRS2016

Automated Storm Tracking and the Lightning Jump Algorithm Using GOES-R Geostationary Lightning Mapper (GLM) Proxy Data

This study develops a fully automated lightning jump system encompassing objective storm tracking, Geostationary Lightning Mapper proxy data, and the lightning jump algorithm (LJA), which are important elements in the transition of the LJA concept from a research to an operational based algorithm. Storm cluster tracking is based on a product created from the combination of a radar parameter (vertically integrated liquid, VIL), and lightning information (flash rate density). Evaluations showed that the spatial scale of tracked features or storm clusters had a large impact on the lightning jump system performance, where increasing spatial scale size resulted in decreased dynamic range of the system's performance. This framework will also serve as a means to refine the LJA itself to enhance its operational applicability. Parameters within the system are isolated and the system's performance is evaluated with adjustments to parameter sensitivity. The system's performance is evaluated using the probability of detection (POD) and false alarm ratio (FAR) statistics. Of the algorithm parameters tested, sigma-level (metric of lightning jump strength) and flash rate threshold influenced the system's performance the most. Finally, verification methodologies are investigated. It is discovered that minor changes in verification methodology can dramatically impact the evaluation of the lightning jump system.

lightning jump↗
DOE OSTI2025

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗
DOE OSTI2024

A Stereo Camera Simulator for Large-Eddy Simulations of Continental Shallow Cumulus Clouds Based on Three-Dimensional Path-Tracing

The complex spatial and temporal structure of cumulus clouds complicates their representation in weather and climate models. Classic meteorological instrumentation struggles to fully capture these features. Networks of multiple high-resolution hemispheric cameras are increasingly used to fill this data gap, and provide information on this missing multi-dimensional spatial information. In this study, a path-tracing algorithm is used to generate virtual camera images of resolved clouds in large-eddy simulations (LES). These images are then used as a camera network simulator, allowing reconstructions of three-dimensional cloud edges from the model output. Because the actual LES cloud field is fully known, the combined path-tracing and reconstruction method can be statistically analyzed. The method is applied to LES realizations of summertime shallow cumulus at the Jülich Observatory for Cloud Evolution (JOYCE), Germany, which also routinely operates a camera network. We find that the path-tracing method allows accurate reconstruction of up to 70% of the visible cloud edges. Additional sensitivity tests show that the method is robust for changes in its hyperparameters. The sensitivity to cloud optical thickness is also investigated, finding a cloud boundary placement error of approximately 182 m. This error can be considered typical for cloud boundary reconstruction using real stereo camera imagery. The results provide proof of principle for future use of the method for evaluating LES clouds against camera network imagery, and for further optimizing the configuration of such camera networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modes and Mechanisms of Pacific Decadal-Scale Variability

The modes of Pacific decadal-scale variability (PDV), traditionally defined as statistical patterns of variance, reflect to first order the ocean's integration (i.e., reddening) of atmospheric forcing that arises from both a shift and a change in strength of the climatological (time-mean) atmospheric circulation. While these patterns concisely describe PDV, they do not distinguish among the key dynamical processes driving the evolution of PDV anomalies, including atmospheric and ocean teleconnections and coupled feedbacks with similar spatial structures that operate on different timescales. In this review, we synthesize past analysis using an empirical dynamical model constructed from monthly ocean surface anomalies drawn from several reanalysis products, showing that the PDV modes of variance result from two fundamental low-frequency dynamical eigenmodes: the North Pacific–central Pacific (NP-CP) and Kuroshio–Oyashio Extension (KOE) modes. Both eigenmodes highlight how two-way tropical–extratropical teleconnection dynamics are the primary mechanisms energizing and synchronizing the basin-scale footprint of PDV. While the NP-CP mode captures interannual- to decadal-scale variability, the KOE mode is linked to the basin-scale expression of PDV on decadal to multidecadal timescales, including contributions from the South Pacific.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Effects of spatial variability and scale on areal -average evapotranspiration

This paper explores the effect of spatial variability and scale on areally-averaged evapotranspiration. A spatially-distributed water and energy balance model is employed to determine the effect of explicit patterns of model parameters and atmospheric forcing on modeled areally-averaged evapotranspiration over a range of increasing spatial scales. The analysis is performed from the local scale to the catchment scale. The study area is King's Creek catchment, an 11.7 sq km watershed located on the native tallgrass prairie of Kansas. The dominant controls on the scaling behavior of catchment-average evapotranspiration are investigated by simulation, as is the existence of a threshold scale for evapotranspiration modeling, with implications for explicit versus statistical representation of important process controls. It appears that some of our findings are fairly general, and will therefore provide a framework for understanding the scaling behavior of areally-averaged evapotranspiration at the catchment and larger scales.

Famiglietti, J. S.↗
NASA NTRS1994

Time scales of pattern evolution from cross-spectrum analysis of advanced very high resolution radiometer and coastal zone color scanner imagery

We have selected square subareas (110 km on a side) from coastal zone color scanner (CZCS) and advanced very high resolution radiometer (AVHRR) images for 1981 in the California Current region off northern California for which we could identify sequences of cloud-free data over periods of days to weeks. We applied a two-dimensional fast Fourier transformation to images after median filtering, (x, y) plane removal, and cosine tapering. We formed autospectra and coherence spectra as functions of a scalar wavenumber. Coherence estimates between pairs of images were plotted against time separation between images for several wide wavenumber bands to provide a temporal lagged coherence function. The temporal rate of loss of correlation (decorrelation time scale) in surface patterns provides a measure of the rate of pattern change or evolution as a function of spatial dimension. We found that patterns evolved (or lost correlation) approximately twice as rapidly in upwelling jets as in the 'quieter' regions between jets. The rapid evolution of pigment patterns (lifetime of about 1 week or less for scales of 50-100 km) ought to hinder biomass transfer to zooplankton predators compared with phytoplankton patches that persist for longer times. We found no significant differences between the statistics of CZCS and AVHRR images (spectral shape or rate of decorrelation). In addition, in two of the three areas studied, the peak correlation between AVHRR and CZCS images from the same area occurred at zero lag, indicating that the patterns evolved simutaneously. In the third area, maximum coherence between thermal and pigment patterns occurred when pigment images lagged thermal images by 1-2 days, mirroring the expected lag of high pigment behind low temperatures (and high nutrients) in recently upwelled water. We conclude that in dynamic areas such as coastal upwelling systems, the phytoplankton cells (identified by pigment color patterns) behave largely as passive scalars at the mesoscale and that growth, death, and sinking of phytoplankton collectively play at most a mariginal role in determining the spectral statistics of the pigment patterns.

Denman, Kenneth L.↗
NASA NTRS2016

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗