Empirical fitting of forward backscattering models for multitemporal retrieval of soil moisture from radar data at L-band

A multitemporal algorithm, originally conceived for the C-band radar aboard the Sentinel-1 satellite, has been updated to retrieve soil moisture from L-band radar data, such as those provided by the National Aeronautics and Space Administration Soil Moisture Active/Passive (SMAP) mission. This type of algorithm may deliver more accurate soil moisture maps that mitigate the effect of roughness and vegetation changes. Within the multitemporal inversion scheme based on the Bayesian maximum a posteriori probability (MAP) criterion, a dense time series of radar measurements is integrated to invert a forward backscattering model. The model calibration and validation tasks have been accomplished using the data collected during the SMAP validation experiment 12 spanning several soil conditions (pasture, wheat, corn, and soybean). The data have been used to update the forward model for bare soil scattering at L-band and to tune a simple vegetation scattering model considering two different classes of vegetation: those producing mainly single scattering effects and those characterized by a significant multiple scattering involving terrain surface and vegetation elements interaction. The algorithm retrievals showed a root mean square difference (RMSD) around 5% over bare soil, soybean, and cornfields. As for wheat, a bias was observed; when removed, the RMSD went down from 7.7% to 5%

Development of SMAP Mission Cal/Val Activities

The Soil Moisture Active Passive (SMAP) mission is a NASA directed mission to map global land surface soil moisture and freeze-thaw state. Instrument and mission details are shown. The key SMAP soil moisture product is provided at 10 km resolution with 0.04cubic cm/cubic cm accuracy. The freeze/thaw product is provided at 3 km resolution and 80% frozen-thawed classification accuracy. The full list of SMAP data products is shown

Retrieving landscape freeze/thaw state fromSoil Moisture Active Passive (SMAP) radar and radiometer measurements

Over one-third of the global land area undergoes a seasonal transition between predominantly frozen and non-frozen conditions each year, with the land surface freeze/thaw (FT) state a significant control on hydrological and biospheric processes over northern land areas and at high elevations. The NASA Soil Moisture Active Passive (SMAP) mission produced a daily landscape FT product at 3-km spatial resolution derived from ascending and descending orbits of SMAP high-resolution L-band (1.4 GHz) radar measurements. Following the failure of the SMAP radar in July 2015, coarser (36-km) footprint SMAP radiometer inputs were used to develop an alternative daily passive microwave freeze/thaw product. In this study, in situ observations are used to examine differences in the sensitivity of the 3-km radar versus the 36-km radiometer measurements to the landscape freeze/thaw state during the period of overlapping instrument operation. Assessment of the retrievals at high-latitude SMAP core validation sites showed excellent agreement with in situ flags, exceeding the 80% SMAP mission accuracy requirement. Similar performance was found for the radar and radiometer products using both air temperature and soil temperature derived FT reference flags. There was a tendency for SMAP thaw retrievals to lead the surface flags due to the influence of wet snow cover conditions on both the radar and radiometer signal. Comparison with other satellite derived FT products showed those derived from passive measurements (SMAP radiometer; Aquarius radiometer; Advanced Microwave Scanning Radiometer - 2) retrieved less frozen area than the active products (SMAP radar; Aquarius radar)

Global evaluation of SMAP/Sentinel-1 soil moisture products

MAP/Sentinel-1 soil moisture is the latest SMAP (Soil Moisture Active Passive) product derived from synergistic utilization of the radiometry observations of SMAP and radar backscattering data of Sentinel-1. This product is the first and only global soil moisture (SM) map at 1 km and 3 km spatial resolutions. In this paper, we evaluated the SMAP/Sentinel-1 SM product from different viewpoints to better understand its quality, advantages, and likely limitations. A comparative analysis of this product and in situ measurements, for the time period March 2015 to January 2022, from 35 dense and sparse SM networks and 561 stations distributed around the world was carried out. We examined the effects of land cover, vegetation fraction, water bodies, urban areas, soil characteristics, and seasonal climatic conditions on the performance of active–passive SMAP/Sentinel-1 in estimating the SM. We also compared the performance metrics of enhanced SMAP (9 km) and SMAP/Sentinel-1 products (3 km) to analyze the effects of the active–passive disaggregation algorithm on various features of the SMAP SM maps. Results showed satisfactory agreement between SMAP/Sentinel-1 and in situ SM measurements for most sites (r values between 0.19 and 0.95 and ub-RMSE between 0.03 and 0.17), especially for dense sites without representativeness errors. Thanks to the vegetation effect correction applied in the active–passive algorithm, the SMAP/Sentinel-1 product had the highest correlation with the reference data in grasslands and croplands. Results also showed that the accuracy of the SMAP/Sentinel-1 SM product in different networks is independent of the presence of water bodies, urban areas, and soil types.Peer ReviewedPostprint (published version

Soil Moisture ActivePassive (SMAP) L-Band Microwave Radiometer Post-Launch Calibration

The SMAP microwave radiometer is a fully-polarimetric L-band radiometer flown on the SMAP satellite in a 6 AM/ 6 PM sun-synchronous orbit at 685 km altitude. Since April, 2015, the radiometer is under calibration and validation to assess the quality of the radiometer L1B data product. Calibration methods including the SMAP L1B TA2TB (from Antenna Temperature (TA) to the Earths surface Brightness Temperature (TB)) algorithm and TA forward models are outlined, and validation approaches to calibration stability/quality are described in this paper including future work. Results show that the current radiometer L1B data satisfies its requirements

Parameter Optimization of a Discrete Scattering Model by Integration of Global Sensitivity Analysis Using SMAP Active and Passive Observations

Active and passive microwave signatures respond differently to the land surface and provide complementary information on the characteristics of the observed scenes. The objective of this paper is to explore the synergy of active radar and passive radiometer observations at the same spatial scale to constrain a discrete radiative transfer model, the Tor Vergata (TVG) model, to gain insights into the microwave scattering and emission mechanisms over grasslands. The TVG model can simultaneously simulate the backscattering coefficient and emissivity with a set of input parameters. To calibrate this model, in situ soil moisture and temperature data collected from the Maqu area in the northeastern region of the Tibetan Plateau, interpolated leaf area index (LAI) data from the Moderate Resolution Imaging Spectroradiometer LAI eight-day products, and concurrent and coincident Soil Moisture Active Passive (SMAP) radar and radiometer observations are used. Because this model needs numerous input parameters to be driven, the extended Fourier amplitude sensitivity test is first applied to conduct global sensitivity analysis (GSA) to select the sensitive and insensitive parameters. Only the most sensitive parameters are defined as free variables, to separately calibrate the active-only model (TVG-A), the passive-only model (TVG-P), and the active and passive combined model (TVG-AP). The accuracy of the calibrated models is evaluated by comparing the SMAP observations and the model simulations. The results show that TVG-AP can well reproduce the backscattering coefficient and brightness temperature, with correlation coefficients of 0.87, 0.89, 0.78, and 0.43 and root-mean-square errors of 0.49 dB, 0.52 dB, 7.20 K, and 10.47 K for σ HH⁰ , σ VV⁰ , TBH, and TBV, respectively. In contrast, TVG-A and TVG-P can only accurately model the backscattering coefficient and brightness temperature, respectively. Without any modifications of the calibrated parameters, the error metrics computed from the validation data are slightly worse than those of the calibration data. These results demonstrate the feasibility of the synergistic use of SMAP active radar and passive radiometer observations under the unified framework of a physical model. In addition, the results demonstrate the necessity and effectiveness of applying GSA in model optimization. It is expected that these findings can contribute to the development of model-based soil moisture retrieval methods using active and passive microwave remote sensing data

Evaluation of SMAP Freeze/Thaw Retrieval Accuracy at Core Validation Sites in the Contiguous United States

Seasonal freeze-thaw (FT) impacts much of the northern hemisphere and is an important control on its water, energy, and carbon cycle. Although FT in natural environments extends south of 45°N, FT studies using the L-band have so far been restricted to boreal or greater latitudes. This study addresses this gap by applying a seasonal threshold algorithm to Soil Moisture Active Passive (SMAP) data (L3_SM_P) to obtain a FT product south of 45°N (‘SMAP FT’), which is then evaluated at SMAP core validation sites (CVS) located in the contiguous United States (CONUS). SMAP landscape FT retrievals are usually in good agreement with 0–5 cm soil temperature at SMAP grids containing CVS stations (\u3e70%). The accuracy could be further improved by taking into account specific overpass time (PM), the grid-specific seasonal scaling factor, the data aggregation method, and the sampling error. Annual SMAP FT extent maps compared to modeled soil temperatures derived from the Goddard Earth Observing System Model Version 5 (GEOS-5) show that seasonal FT in CONUS extends to latitudes of about 35–40°N, and that FT varies substantially in space and by year. In general, spatial and temporal trends between SMAP and modeled FT were similar

Evaluation of Sentinel-1, SMAP and SMOS surface soil moisture products for distributed eco-hydrological modelling in Mediterranean forest basins

[EN] Reliable distributed hydrological modeling, especially in semi-arid areas, must consider the inclusion of surface soil moisture (SSM) spatial information during the calibration process. This variable plays a key role in the evapotranspiration processes that determine the hydrological cycle. The coarse resolution of the SSM estimates by satellite remote sensing has restricted the application of this approach to only large basins, focusing most of the studies in the consideration of simply the temporal dynamics of this variable. The growing efforts in providing higher spatial resolution through disaggregating methodologies or new sensor estimates facilitates the application of this spatial approach to small basins. This paper explores the applicability of the currently available satellite surface soil moisture estimates for distributed eco-hydrological modelling in Mediterranean forest basins. On one hand, this study contributes to fill the existing research gap on the use of remote sensing SSM spatial patterns within the distributed hydrological modelling framework in small basins. On the other hand, it serves as an indirect validation method for the spatial performance of satellite SSM products. To achieve this goal, we implemented the eco-hydrological model TETIS in three case studies named: Hozgarganta (southern Spain), Ceira (western Portugal) and Carraixet (eastern Spain). The SSM estimates selected for comparison were Sentinel-1 SSM provided by the Copernicus Global Land Services (CGLS), SMAP SSM disaggregated using Sentinel-1 (SPL2SMAP_S) provided by the National Aeronautics and Space Administration (NASA), SMOS SSM provided by the Barcelona Expert Center (BEC), and SMOS and SMAP SSM disaggregated using the DISPATCH algorithm provided by Lobelia Earth. The methodology employed involved a multi-objective and multi-variable calibration in terms of remote sensing SSM spatial patterns and in-situ streamflow, using the Spatial Efficiency Metric (SPAEF) and the Nash-Sutcliffe efficiency index (NSE) respectively. Before model calibration a sensitivity analysis of the most influent variables was performed. The temporal and spatial comparison of the reference SSM products revealed inconsistencies amongst products. The disaggregating methodology determined the spatial agreement to a greater degree than the sensor itself (i.e. SMAP, SMOS). In spite of the differences amongst products, the multi-objective calibration approach proposed increased the robustness of the hydrological modelling.This study was founded by the Spanish AEI within the program WaterJPI through the project iAqueduct (PCI2019-103729) , by the EC Life project ResilientForests (LIFE17 CCA/ES/000063) , and by the project Water4Cast funded by Generalitat Valenciana (PROMETEO/2021/074) . We also acknowledge the following hydrometeorological data providers institutions: SiAR, SAIH-HIDROSUR, SAIH Jucar and SNIRH.Gomis-Cebolla, J.; Garcia-Arias, A.; Perpinyà-Vallès, M.; Francés, F. (2022). Evaluation of Sentinel-1, SMAP and SMOS surface soil moisture products for distributed eco-hydrological modelling in Mediterranean forest basins. Journal of Hydrology. 608:1-19. https://doi.org/10.1016/j.jhydrol.2022.12756911960

Mapping Soil Moisture from Remotely Sensed and In-situ Data with Statistical Methods

Soil moisture is an important factor for accurate prediction of agricultural productivity and rainfall runoff with hydrological models. Remote sensing satellites such as Soil Moisture Active Passive (SMAP) offer synoptic views of soil moisture distribution at a regional-to-global scale. To use the soil moisture product from these satellites, however, requires a downscaling of the data from an usually large instantaneous field of view (i.e. 36 km) to the watershed analysis scales ranging from 30 m to 1 km. In addition, validation of the soil moisture products using the ground station observations without an upscaling treatment would lead to cross-level fallacy. In the literature of geographical analysis, scale is one of the top research concens because of the needs for multi-source geospatial data fusion. This dissertation research introduced a multi-level soil moisture data assimilation and processing methodology framework based on spatial information theories. The research contains three sections: downscaling using machine learning and geographically weighted regression, upscaling ground network observation to calibrate satellite data, and spatial and temporal multi-scale data assimilation using spatio-temporal interpolation. (1) Soil moisture downscaling In the first section, a downscaling method is designed using 1-km geospatial data to obtain subpixel soil moisture from the 9-km soil moisture product of the SMAP satellite. The geospatial data includes normalized difference vegetation index (NDVI), land surface temperature (LST), gross primary productivity (GPP), topographical moisture index (TMI), with all resampled to 1-km resolution. The machine learning algorithm – random forest was used to create a prediction model of the soil moisture at a 1-km resolution. The 1-km soil moisture product was compared with the ground samples from the West Texas Mesonet (WTM) station data. The residual was then interpolated to compensate the unpredicted variability of the model. The entire process was based on the concept of regression kriging- where the regression was done by the random forest model. Results show that the downscaling approach was able to achieve better accuracy than the current statistical downscaling methods. (2) Station network data upscaling The Texas Soil Observation Network (TxSON) network was designed to test the feasibility of upscaling the in-situ data to match the scale of the SMAP data. I advanced the upscaling method by using the Voronoi polygons and block kriging with a Gaussian kernel aggregation. The upscaling algorithm was calibrated using different spatial aggregation parameters, such as the fishnet cell size and Gaussian kernel standard deviation. The use of the kriging can significantly reduce the spatial autocorrelation among the TxSON stations because of its declustering ability. The result proved the new upscaling method was better than the traditional ones. (3) Multi-scale data fusion in a spatio-temporal framework None of the current works for soil moisture statistical downscaling honors time and space equally. It is important, however, that the soil moisture products are consistent in both domains. In this section, the space-time kriging model for soil moisture downscaling and upscaling computation framework designed in the last two sections is implemented to create a spatio-temporal integrated solution to soil moisture multi-scale mapping. The present work has its novelty in using spatial statistics to reconcile the scale difference from satellite data and ground observations, and therefore proposes new theories and solutions for dealing with the modifiable areal unit problem (MAUP) incurred in soil moisture mapping from satellite and ground stations