Search NASA⌕ Search

SEARCH · Search NASA

Results for “Standard deviation modeling”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗
DOE OSTI2025

Search for a new scalar resonance decaying to a Higgs boson and another new scalar particle in the final state with two bottom quarks and two photons in proton-proton collisions at $\sqrt{s} = 13$ TeV

A search is presented for a new scalar resonance, X, decaying to a standard model Higgs boson and another new scalar particle, Y, in the final state where the Higgs boson decays to a $\text{b}\overline{\text{b} }$ pair, while the Y particle decays to a pair of photons. The search is performed in the mass range 240–1000 GeV for the resonance X, and in the mass range 70–800 GeV for the particle Y, using proton-proton collision data collected by the CMS experiment at $\sqrt{s}=13$ TeV, corresponding to an integrated luminosity of 132 fb −1 . In general, the data are found to be compatible with the standard model expectation. Observed (expected) upper limits at 95% confidence level on the product of the production cross section and the relevant branching fraction are extracted for the X → YH process, and are found to be within the range of 0.05–2.69 (0.08–1.94) fb, depending on m X and m Y . The most significant deviation from the background-only hypothesis is observed for X and Y masses of 300 and 77 GeV, respectively, with a local (global) significance of 3.33 (0.65) standard deviations.

Beyond Standard Model↗
DOE OSTI2025

A versatile machine learning workflow for high-throughput analysis of supported metal catalyst particles

Accurate and efficient characterization of nanoparticles (NPs), particularly regarding particle size distribution, is essential for advancing our understanding of their structure-property relationship and facilitating their design for various applications. In this study, we introduce a novel two-stage artificial intelligence (AI)-driven workflow for NP analysis that leverages prompt engineering techniques from state-of-the-art single-stage object detection and large-scale vision transformer (ViT) architectures. This methodology is applied to transmission electron microscopy (TEM) and scanning TEM (STEM) images of heterogeneous catalysts, enabling high-resolution, high-throughput analysis of particle size distributions for supported metal catalyst NPs. The model's performance in detecting and segmenting NPs is validated across diverse heterogeneous catalyst systems, including various metals (Ru, Cu, PtCo, and Pt), supports (silica (SiO 2 ), γ-alumina (γ-Al 2 O 3 ), and carbon black), and particle diameter size distributions with mean and standard deviations ranging from 1.6 ± 0.2 nm to 9.7 ± 4.6 nm. The proposed machine learning (ML) methodology achieved an average F1 overlap score of 0.91 ± 0.01 and demonstrated the ability to disentangle overlapping NPs anchored on catalytic support materials. The segmentation accuracy is further validated using the Hausdorff distance and robust Hausdorff distance metrics, with the 90th percent of the robust Hausdorff distance showing errors within 0.4 ± 0.1 nm to 1.4 ± 0.6 nm. In conclusion, our AI-assisted NP analysis workflow demonstrates robust generalization across diverse datasets and can be readily applied to similar NP segmentation tasks without requiring costly model retraining.

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

Measurement of 𝑊⁢𝑊⁢𝑍 and 𝑍⁢𝐻 Production Cross Sections at $\sqrt{𝑠}$ = 13 and 13.6 TeV

A measurement is presented of the cross section in proton-proton collisions for the production of two 𝑊 bosons and one 𝑍 boson. It is based on data recorded by the CMS experiment at the CERN LHC at center-of-mass energies $\sqrt{𝑠}$ = 13 and 13.6 TeV, corresponding to an integrated luminosity of 200 fb −1 . Events with four charged leptons (electrons or muons) in the final state are selected. Both nonresonant 𝑊⁢𝑊⁢𝑍 production and 𝑍⁢𝐻 production, with the Higgs boson decaying into two 𝑊 bosons, are reported. For the first time, the two processes are measured separately in a simultaneous fit. Combining the two modes, signal strengths relative to the standard model (SM) predictions of 0.75$^{+0.34}_{−0.29}$ and 1.74$^{+0.71}_{−0.60}$ are measured for $\sqrt{𝑠}$ = 13 and 13.6 TeV, respectively. The observed (expected) significance for the triboson signal is 3.8 (2.5) standard deviations for $\sqrt{𝑠}$ = 13.6 TeV, thus providing the first evidence for triboson production at this center-of-mass energy. Combining the two modes and the two center-of-mass energies, the inclusive signal strength relative to the SM prediction is measured to be 1.03$^{+0.31}_{−0.28}$, with an observed (expected) significance of 4.5 (5.0) standard deviations.

Gauge bosons↗
NASA NTRS2007

Mars Global Reference Atmospheric Model (Mars-GRAM 2005) Applications for Mars Science Laboratory Mission Site Selection Processes

The new Mars-GRAM auxiliary profile capability, using data from TES observations, mesoscale model output, or other sources, allows a potentially higher fidelity representation of the atmosphere, and a more accurate way of estimating inherent uncertainty in atmospheric density and winds. Figure 3 indicates that, with nominal value rpscale=1, Mars-GRAM perturbations would tend to overestimate observed or mesoscale-modeled variability. To better represent TES and mesoscale model density perturbations, rpscale values as low as about 0.4 could be used. Some trajectory model implementations of Mars-GRAM allow the user to dynamically change rpscale and rwscale values with altitude. Figure 4 shows that an mscale value of about 1.2 would better replicate wind standard deviations from MRAMS or MMM5 simulations at the Gale, Terby, or Melas sites. By adjusting the rpscale and rwscale values in Mars-GRAM based on figures such as Figure 3 and 4, we can provide more accurate end-to-end simulations for EDL at the candidate MSL landing sites.

Justh, H. L.↗
NASA NTRS2023

Characterization of Errors in Satellite-Based HCHO/NO 2 Tropospheric Column Ratios With Respect to Chemistry, Column-to-PBL Translation, Spatial Representation, and Retrieval Uncertainties

The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.

Amir H. Souri↗
DOE OSTI2025

Observation of W + W – γ production in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector and constraints on anomalous quartic gauge-boson couplings

This Letter reports the observation of W + W – γ triboson production in 140 fb –1 of data collected by the ATLAS detector from proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV at the LHC. Events with an opposite-charge eμ pair, a high transverse-momentum photon, and significant missing transverse momentum are considered. The observed (expected) significance of the signal is 5.9 (6.0) standard deviations. The measured fiducial cross-section, defined for the W + W – γ → e ± μ ∓ $ν\bar{ν}$γ final state is 6.2 ± 0.8 (stat.) ± 0.6 (sys.) fb, in good agreement with the Standard Model prediction of 6.1$^{+1.0}_{–0.7}$ fb. Constraints on the Wilson coefficients of 13 dimension-8 operators describing physics beyond the Standard Model through anomalous quartic gauge-boson couplings are derived using the effective field theory framework.

Anomalous quartic gauge couplings↗
NASA NTRS1973

Lunar polarization studies at 3.1 mm wavelength

Observations of the distribution of linearly polarized lunar thermal emission were made at a wavelength of 3.1 mm with The University of Texas 4.88 m parabolic reflector (0.042 deg HPBW). A shadow corrected, rough surface, thermal emission model for a homogeneous moon was least-squares-fitted to the polarization data. Results indicate an effective lunar dielectric constant of 1.34 plus or minus 0.04 with surface roughness characterized by a standard deviation of 17 (plus or minus 5) deg for surface slopes with a normal probability density, independent of lunar phase. A comparison of these results with published values at other wavelengths suggests that the effective lunar dielectric constant, as obtained by lunar emission measurements, decreases with decreasing wavelength of observation. This wavelength dependence may be interpreted in terms of an inhomogeneous surface and/or a surface that possesses intermediate scale surface roughness.

White, T. L.↗
NASA NTRS1996

Squeezed States, Uncertainty Relations and the Pauli Principle in Composite and Cosmological Models

The importance of not only uncertainty relations but also the Pauli exclusion principle is emphasized in discussing various 'squeezed states' existing in the universe. The contents of this paper include: (1) Introduction; (2) Nuclear Physics in the Quark-Shell Model; (3) Hadron Physics in the Standard Quark-Gluon Model; (4) Quark-Lepton-Gauge-Boson Physics in Composite Models; (5) Astrophysics and Space-Time Physics in Cosmological Models; and (6) Conclusion. Also, not only the possible breakdown of (or deviation from) uncertainty relations but also the superficial violation of the Pauli principle at short distances (or high energies) in composite (and string) models is discussed in some detail.

Terazawa, Hidezumi↗
NASA NTRS2012

Vandenberg Air Force Base Upper Level Wind Launch Weather Constraints

The 30th Operational Support Squadron Weather Flight (30 OSSWF) provides comprehensive weather services to the space program at Vandenberg Air Force Base (VAFB) in California. One of their responsibilities is to monitor upper-level winds to ensure safe launch operations of the Minuteman III ballistic missile. The 30 OSSWF tasked the Applied Meteorology Unit (AMU) to analyze VAFB sounding data with the goal of determining the probability of violating (PoV) their upper-level thresholds for wind speed and shear constraints specific to this launch vehicle, and to develop a tool that will calculate the PoV of each constraint on the day of launch. In order to calculate the probability of exceeding each constraint, the AMU collected and analyzed historical data from VAFB. The historical sounding data were retrieved from the National Oceanic and Atmospheric Administration Earth System Research Laboratory archive for the years 1994-2011 and then stratified into four sub-seasons: January-March, April-June, July-September, and October-December. The maximum wind speed and 1000-ft shear values for each sounding in each subseason were determined. To accurately calculate the PoV, the AMU determined the theoretical distributions that best fit the maximum wind speed and maximum shear datasets. Ultimately it was discovered that the maximum wind speeds follow a Gaussian distribution while the maximum shear values follow a lognormal distribution. These results were applied when calculating the averages and standard deviations needed for the historical and real-time PoV calculations. In addition to the requirements outlined in the original task plan, the AMU also included forecast sounding data from the Rapid Refresh model. This information provides further insight for the launch weather officers (LWOs) when determining if a wind constraint violation will occur over the next few hours on day of launch. The interactive graphical user interface (GUI) for this project was developed in Microsoft Excel using Visual Basic for Applications. The GUI displays the critical sounding data easily and quickly for the LWOs on day of launch. This tool will replace the existing one used by the 30 OSSWF, assist the LWOs in determining the probability of exceeding specific wind threshold values, and help to improve the overall upper winds forecast for the launch customer.

Shafer, Jaclyn A.↗
NASA NTRS2012

OH PLIF Visualization of the UVa Supersonic Combustion Experiment: Configuration A

Hydroxyl radical (OH) planar laser-induced fluorescence (PLIF) measurements were performed in the University of Virginia s dual-mode scramjet experiment. The test section was set up in configuration A, which includes a Mach 2 nozzle, combustor, and extender section. Hydrogen fuel was injected through an unswept compression ramp at two different equivalence ratios. Through the translation of the optical system and the use of two separate camera views, the entire optical range of the combustor was accessed. Single-shot, average, and standard deviation images of the OH PLIF signal are presented at several streamwise locations. The results show the development of a highly turbulent flame structure and provide an experimental database to be used for numerical model assessment.

Johansen, Craig T.↗
NASA NTRS2013

OH PLIF Visualization of the UVa Supersonic Combustion Experiment: Configuration A

Hydroxyl radical (OH) planar laser-induced fluorescence (PLIF) measurements were performed in the University of Virginia supersonic combustion experiment. The test section was set up in configuration A, which includes a Mach 2 nozzle, combustor, and extender section. Hydrogen fuel was injected through an unswept compression ramp at two different equivalence ratios. Through the translation of the optical system and the use of two separate camera views, the entire optically accessible range of the combustor was imaged. Single-shot, average, and standard deviation images of the OH PLIF signal are presented at several streamwise locations. The results show the development of a highly turbulent flame structure and provide an experimental database to be used for numerical model assessment.

Johansen, Craig T.↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2026

Combination of searches for nonresonant Higgs boson pair production in proton–proton collisions at $\sqrt{s}$ = 13 TeV

This paper presents a combination of searches for the nonresonant production of Higgs boson pairs (HH) in proton–proton collisions at a centre-of-mass energy of 13 TeV. The dataset was collected by the CMS experiment at the LHC from 2016 to 2018 and corresponds to a total integrated luminosity of 138 fb –1 . The observed (expected) upper limit on the inclusive HH production cross-section relative to the standard model (SM) prediction is found to be 3.5 (2.5). Assuming all other Higgs boson couplings are equal to their SM values, the Higgs boson trilinear self-coupling modifier κ λ = λ 3 /λ$^{SM}_{3}$ is constrained in the range –1.35 ≤ κ λ ≤ 6.37 at 95% confidence level (CL). Similarly, for the coupling modifier κ 2V , which governs the interactions between two vector bosons and two Higgs bosons, we have excluded κ 2V = 0 at more than five standard deviations for all values of κ λ . At 95% CL, assuming that the other couplings are equal to their SM values, κ 2V is constrained in the range 0.64 ≤ κ 2V ≤ 1.40. This work also investigates HH production in several new physics scenarios, using the Higgs effective field theory (HEFT) framework. The HEFT framework is further exploited to study various ultraviolet complete models with an extended Higgs sector and set constraints on specific parameters. An extrapolation of the results in the expected integrated luminosity following the high-luminosity upgrade of the LHC is also reported.

BSM↗
DOE OSTI2026

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗
NASA NTRS1992

Estimating structural attributes of Douglas-fir/western hemlock forest stands from Landsat and SPOT imagery

Relationships between spectral and texture variables derived from SPOT HRV 10 m panchromatic and Landsat TM 30 m multispectral data and 16 forest stand structural attributes is evaluated to determine the utility of satellite data for analysis of hemlock forests west of the Cascade Mountains crest in Oregon and Washington, USA. Texture of the HRV data was found to be strongly related to many of the stand attributes evaluated, whereas TM texture was weakly related to all attributes. Data analysis based on regression models indicates that both TM and HRV imagery should yield equally accurate estimates of forest age class and stand structure. It is concluded that the satellite data are a valuable source for estimation of the standard deviation of tree sizes, mean size and density of trees in the upper canopy layers, a structural complexity index, and stand age.

Cohen, Warren B.↗
NASA NTRS1998

A Simple Model of Cirrus Horizontal Inhomogeneity and Cloud Fraction

A simple model of horizontal inhomogeneity and cloud fraction in cirrus clouds has been formulated on the basis that all internal horizontal inhomogeneity in the ice mixing ratio is due to variations in the cloud depth, which are assumed to be Gaussian. The use of such a model was justified by the observed relationship between the normalized variability of the ice water mixing ratio (and extinction) and the normalized variability of cloud depth. Using radar cloud depth data as input, the model reproduced well the in-cloud ice water mixing ratio histograms obtained from horizontal runs during the FIRE2 cirrus campaign. For totally overcast cases the histograms were almost Gaussian, but changed as cloud fraction decreased to exponential distributions which peaked at the lowest nonzero ice value for cloud fractions below 90%. Cloud fractions predicted by the model were always within 28% of the observed value. The predicted average ice water mixing ratios were within 34% of the observed values. This model could be used in a GCM to produce the ice mixing ratio probability distribution function and to estimate cloud fraction. It only requires basic meteorological parameters, the depth of the saturated layer and the standard deviation of cloud depth as input.

Smith, Samantha A.↗
NASA NTRS2008

Field Observations of Soil Moisture Variability across Scales

In this study, over 36,000 ground-based soil moisture measurements collected during the SGP97, SGP99, SMEX02, and SMEX03 field campaigns were analyzed to characterize the behavior of soil moisture variability across scales. The field campaigns were conducted in Oklahoma and Iowa in the central USA. The Oklahoma study region is sub-humid with moderately rolling topography, while the Iowa study region is humid with low-relief topography. The relationship of soil moisture standard deviation, skewness and the coefficient of variation versus mean moisture content was explored at six distinct extent scales, ranging from 2.5 m to 50 km. Results showed that variability generally increases with extent scale. The standard deviation increased from 0.036 cm3/cm3 at the 2.5-m scale to 0.071 cm3/cm3 at the 50-km scale. The log standard deviation of soil moisture increased linearly with the log extent scale, from 16 m to 1.6 km, indicative of fractal scaling. The soil moisture standard deviation versus mean moisture content exhibited a convex upward relationship at the 800-m and 50-km scales, with maximum values at mean moisture contents of roughly 0.17 cm3/cm3 and 0.19 cm3/cm3, respectively. An empirical model derived from the observed behavior of soil moisture variability was used to estimate uncertainty in the mean moisture content for a fixed number of samples at the 800-m and 50-km scales, as well as the number of ground-truth samples needed to achieve 0.05 cm3/cm3 and 0.03 cm3/cm3 accuracies. The empirical relationships can also be used to parameterize surface soil moisture variations in land surface and hydrological models across a range of scales. To our knowledge, this is the first study to document the behavior of soil moisture variability over this range of extent scales using ground-based measurements. Our results will contribute not only to efficient and reliable satellite validation, but also to better utilization of remotely sensed soil moisture products for enhanced modeling and prediction.

Famiglietti, James S.↗