Synthetic Aperture Radar has shown its large potential for retrieving soil moisture maps at regional scales. However, since the backscattered signal is determined by several surface characteristics, the retrieval of soil moisture is an ill-posed problem when using single configuration imagery. Unless accurate surface roughness parameter values are available, retrieving soil moisture from radar backscatter usually provides inaccurate estimates. The characterization of soil roughness is not fully understood, and a large range of roughness parameter values can be obtained for the same surface when different measurement methodologies are used. In this paper, a literature review is made that summarizes the problems encountered when parameterizing soil roughness as well as the reported impact of the errors made on the retrieved soil moisture. A number of suggestions were made for resolving issues in roughness parameterization and studying the impact of these roughness problems on the soil moisture retrieval accuracy and scale.

The objective of this paper is to propose a new semi-empirical radar backscattering model for bare soil surfaces based on the Dubois model. A wide dataset of backscattering coefficients extracted from synthetic aperture radar (SAR) images and in situ soil surface parameter measurements (moisture content and roughness) is used. The retrieval of soil parameters from SAR images remains challenging because the available backscattering models have limited performances. Existing models, physical, semi-empirical, or empirical, do not allow for a reliable estimate of soil surface geophysical parameters for all surface conditions. The proposed model, developed in HH, HV, and VV polarizations, uses a formulation of radar signals based on physical principles that are validated in numerous studies. Never before has a backscattering model been built and validated on such an important dataset as the one proposed in this study. It contains a wide range of incidence angles (18°–57°) and radar wavelengths (L, C, X), well distributed, geographically, for regions with different climate conditions (humid, semi-arid, and arid sites), and involving many SAR sensors. The results show that the new model shows a very good performance for different radar wavelengths (L, C, X), incidence angles, and polarizations (RMSE of about 2 dB). This model is easy to invert and could provide a way to improve the retrieval of soil parameters.

The aim of this paper is to evaluate the most used radar backscattering models (Integral Equation Model “IEM”, Oh, Dubois, and Advanced Integral Equation Model “AIEM”) using a wide dataset of SAR (Synthetic Aperture Radar) data and experimental soil measurements. These forward models reproduce the radar backscattering coefficients ( ${\mathsf{\sigma }}^0$ ) from soil surface characteristics (dielectric constant, roughness) and SAR sensor parameters (radar wavelength, incidence angle, polarization). The analysis dataset is composed of AIRSAR, SIR-C, JERS-1, PALSAR-1, ESAR, ERS, RADARSAT, ASAR and TerraSAR-X data and in situ measurements (soil moisture and surface roughness). Results show that Oh model version developed in 1992 gives the best fitting of the backscattering coefficients in HH and VV polarizations with RMSE values of 2.6 dB and 2.4 dB, respectively. Simulations performed with the Dubois model show a poor correlation between real data and model simulations in HH polarization (RMSE = 4.0 dB) and better correlation with real data in VV polarization (RMSE = 2.9 dB). The IEM and the AIEM simulate the backscattering coefficient with high RMSE when using a Gaussian correlation function. However, better simulations are performed with IEM and AIEM by using an exponential correlation function (slightly better fitting with AIEM than IEM). Good agreement was found between the radar data and the simulations using the calibrated version of the IEM modified by Baghdadi (IEM_B) with bias less than 1.0 dB and RMSE less than 2.0 dB. These results confirm that, up to date, the IEM modified by Baghdadi (IEM_B) is the most adequate to estimate soil moisture and roughness from SAR data.

The objective of this paper is to extend the semi-empirical calibration of the backscattering Integral Equation Model (IEM) initially proposed for Synthetic Aperture Radar (SAR) data at C- and X-bands to SAR data at L-band. A large dataset of radar signal and in situ measurements (soil moisture and surface roughness) over bare soil surfaces were used. This dataset was collected over numerous agricultural study sites in France, Luxembourg, Belgium, Germany and Italy using various SAR sensors (AIRSAR, SIR-C, JERS-1, PALSAR-1, ESAR). Results showed slightly better simulations with exponential autocorrelation function than with Gaussian function and with HH than with VV. Using the exponential autocorrelation function, the mean difference between experimental data and Integral Equation Model (IEM) simulations is +0.4 dB in HH and −1.2 dB in VV with a Root Mean Square Error (RMSE) about 3.5 dB. In order to improve the modeling results of the IEM for a better use in the inversion of SAR data, a semi-empirical calibration of the IEM was performed at L-band in replacing the correlation length derived from field experiments by a fitting parameter. Better agreement was observed between the backscattering coefficient provided by the SAR and that simulated by the calibrated version of the IEM (RMSE about 2.2 dB).

Approximately, 30% of the Sentinel-1 (S-1) swath over land is imaged with incidence angles higher than 40°. Still, the interplay among the scattering mechanisms taking place at such a high incidence and their implications on the backscatter information content is often disregarded. This article investigates, through an experimental and numerical study, the S-1 sensitivity to the surface soil moisture (SSM) over agricultural fields observed at low (~33°) and high (~43°) incidence angles and quantifies the impact of the incidence angle on the SSM retrieval accuracy. The study sites are the Apulian Tavoliere (Italy) and REd de MEDición de la HUmedad del Suelo (REMEDHUS) (Spain), which are both instrumented with a hydrologic network continuously measuring SSM. At low incidence angles, results confirm that for crops such as wheat and barley, dominated in C-band by surface scattering, there exists a good sensitivity of S-1 VV to SSM. At high incidence angles, the sensitivity to SSM holds through the combination of the soil attenuated and double bounce scattering. Conversely, over crops dominated by volume scattering, such as sugar beet, the S-1 VV signal is not correlated with the in situ SSM observations, neither at low nor at high incidence. For all the crops, the sensitivity of S-1 to SSM in VH is found significantly lower than in VV. The impact of the incidence angle on the SSM retrieval has been studied with a recursive algorithm based on a short-term change detection approach. An upper and lower bounds for the worsening of the S-1 VV retrieval performance at far versus near range observations have been estimated. In the worst-case scenario, the root mean square error (RMSE) increases from ~0.056 m³/m³, at low incidence, to ~0.071 m³/m³, at high incidence. The mechanism that lowers the retrieval accuracy at high incidence angles is further investigated in the synthetic experiment and its impact on the RMSE is estimated in terms of the volume scattering contribution.

The retrieval of soil moisture information with spatially and temporally high resolution from Synthetic Aperture Radar (SAR) observations is still a challenge. By using multi-orbit Sentinel-1 C-band time series, we present a novel approach for estimating volumetric soil moisture content for agricultural areas with a temporal resolution of one to two days, based on a short-term change detection method. By applying an incidence angle normalization and a Fourier Series transformation, the effect of varying incidence angles on the backscattering signal could be reduced. As the C-band co-polarized backscattering signal is prone to vegetational changes, it is used in this study for the vegetational correction of its related backscatter ratios. The retrieving algorithm was implemented in a cloud-processing environment, enabling a potential global and scalable application. Validated against eight in-situ cosmic ray neutron probe stations across the Rur catchment (Germany) as well as six capacitance stations at the Apulian Tavoliere (Italy) site for the years 2018 to 2020, the method achieves a correlation coefficient of R of 0.63 with an unbiased Root Mean Square Error of 0.063 m³/m³.