In a context of high stress on water resources and agricultural production at the global level, together with climate change marked by an increase in the frequency of these events, drought is ...considered to be a strong threat both socially and economically. The Mediterranean region is a hot spot of climate change; it is also characterized by a scarcity of water resources that places intense pressure on agricultural productivity. This article analyzes the potential for using multiple remote sensing tools in the quantification and predictability of drought in Northwest Africa. Three satellite products are considered: the Normalized Difference Vegetation Index (NDVI), Soil Moisture Index (SWI), and Land Surface Temperature (LST). A discussion of the variability of these products and their inter-correlation is presented, illustrating a generally high consistency between them. Statistical anomaly indices are then computed and a drought severity mapping is presented. The results illustrate in particular a high percentage of dry conditions in the region studied during the last ten years (2007-2017). Finally, we propose the use of the analog statistical approach to identify similar evolutions of the three variables in the past. Although this technique is not a forecast, it provides a strong indication of the plausible future trajectory of a given hydrological season.
The recent deployment of ESA's Sentinel operational satellites has established a new paradigm for remote sensing applications. In this context, Sentinel-1 radar images have made it possible to ...retrieve surface soil moisture with a high spatial and temporal resolution. This paper presents two methodologies for the retrieval of soil moisture from remotely-sensed SAR images, with a spatial resolution of 100 m. These algorithms are based on the interpretation of Sentinel-1 data recorded in the VV polarization, which is combined with Sentinel-2 optical data for the analysis of vegetation effects over a site in Urgell (Catalunya, Spain). The first algorithm has already been applied to observations in West Africa by Zribi et al., 2008, using low spatial resolution ERS scatterometer data, and is based on change detection approach. In the present study, this approach is applied to Sentinel-1 data and optimizes the inversion process by taking advantage of the high repeat frequency of the Sentinel observations. The second algorithm relies on a new method, based on the difference between backscattered Sentinel-1 radar signals observed on two consecutive days, expressed as a function of NDVI optical index. Both methods are applied to almost 1.5 years of satellite data (July 2015-November 2016), and are validated using field data acquired at a study site. This leads to an RMS error in volumetric moisture of approximately 0.087 m³/m³ and 0.059 m³/m³ for the first and second methods, respectively. No site calibrations are needed with these techniques, and they can be applied to any vegetation-covered area for which time series of SAR data have been recorded.
Soil moisture mapping at a high spatial resolution is very important for several applications in hydrology, agriculture and risk assessment. With the arrival of the free Sentinel data at high spatial ...and temporal resolutions, the development of soil moisture products that can better meet the needs of users is now possible. In this context, the main objective of the present paper is to develop an operational approach for soil moisture mapping in agricultural areas at a high spatial resolution over bare soils, as well as soils with vegetation cover. The developed approach is based on the synergic use of radar and optical data. A neural network technique was used to develop an operational method for soil moisture estimates. Three inversion SAR (Synthetic Aperture Radar) configurations were tested: (1) VV polarization; (2) VH polarization; and (3) both VV and VH polarization, all in addition to the NDVI information extracted from optical images. Neural networks were developed and validated using synthetic and real databases. The results showed that the use of a priori information on the soil moisture condition increases the precision of the soil moisture estimates. The results showed that VV alone provides better accuracy on the soil moisture estimates than VH alone. In addition, the use of both VV and VH provides similar results, compared to VV alone. In conclusion, the soil moisture could be estimated in agricultural areas with an accuracy of approximately 5 vol % (volumetric unit expressed in percent). Better results were obtained for soil with a moderate surface roughness (for root mean surface height between 1 and 3 cm). The developed approach could be applied for agricultural plots with an NDVI lower than 0.75.
In a perspective to develop an inversion approach for estimating surface soil moisture of crop fields from Sentinel-1/2 data (radar and optical sensors), the Water Cloud Model (WCM) was calibrated ...from C-band Synthetic Aperture Radar (SAR) data and Normalized Difference Vegetation Index (NDVI) values collected over crops fields and grasslands. The soil contribution that depends on soil moisture and surface roughness (in addition to SAR instrumental parameters) was simulated using the physical backscattering model IEM (Integral Equation Model). The vegetation descriptor used in the WCM is the NDVI because it can be directly calculated from optical images. A large dataset consisting of radar backscattered signal in Vertical transmit and Vertical receive (VV) and Vertical transmit and Horizontal receive (VH) polarizations with wide range of incidence angle, soil moisture, surface roughness, and NDVI-values was used. It was collected over two agricultural study sites. Results show that the soil contribution to the total radar backscattered signal is lower in VH than in VV because VH is more sensitive to vegetation cover. Thus, the use of VH alone or in addition to VV for retrieving the soil moisture is not advantageous in presence of well-developed vegetation cover.
This paper assesses the potential of Synthetic Aperture Radar (SAR) in the C and L bands to penetrate into the canopy cover of wheat, maize and grasslands. For wheat and grasslands, the sensitivity ...of the C and L bands to in situ surface soil moisture (SSM) was first studied according to three levels of the Normalized Difference Vegetation Index (NDVI < 0.4, 0.4 < NDVI < 0.7, and NDVI > 0.7). Next, the temporal evolution of the SAR signal in the C and L bands was analyzed according to SSM and the NDVI. For wheat and grasslands, the results showed that the L-band in HH polarization penetrates the canopy even when the canopy is well-developed (NDVI > 0.7), whereas the penetration of the C-band into the canopy is limited for an NDVI < 0.7. For an NDVI less than 0.7, the sensitivity of the radar signal to SSM is approximately 0.27 dB/vol.% for the L-band in HH polarization and approximately 0.12 dB/vol.% for the C-band (in both VV and VH polarizations). For highly developed wheat and grassland cover (NDVI > 0.7), the sensitivity of the L-band in HH polarization to SSM is approximately 0.19 dB/vol.%, whereas as the C-band is insensitive to SSM. For maize, only the temporal evolution of the C-band according to SSM and the NDVI was studied because the swath of SAR images in the L-band did not cover the maize plots. The results showed that the C-band in VV polarization is able to penetrate the maize canopy even when the canopy is well developed (NDVI > 0.7) due to high-order scattering along the soil-vegetation pathway that contains a soil contribution. According to results obtained in this paper, the L-band would penetrate a well-developed maize cover since the penetration depth of the L-band is greater than that of the C-band.
Radar data at C-band has shown great potential for the monitoring of soil and canopy hydric conditions of wheat crops. In this study, the C-band Sentinel-1 time series including the backscattering ...coefficients σ0 at VV and VH polarization, the polarization ratio (PR) and the interferometric coherence ρ are first analyzed with the support of experimental data gathered on three plots of irrigated winter wheat located in the Haouz plain in the center of Morocco covering five growing seasons. The results showed that ρ and PR are tightly related to the canopy development. ρ is also sensitive to soil preparation. By contrast, σ0 was found to be widely linked to changes in surface soil moisture (SSM) during the first growth stages when Leaf Area Index remains moderate (<1.5 m2/m2). In addition, drastic changes in the crop geometry associated to heading had a strong impact on the C-band σ0, in particular for VH polarization. The coupled water cloud and Oh models (WCM) were then calibrated and validated on the study sites. The comparison between the predicted and observed σ0 yielded a root mean square error (RMSE) values ranging from 1.50 dB to 2.02 dB for VV and between 1.74 dB to 2.52 dB for VH with significant differences occurring in the second part of the season after heading. Finally, new approaches based on the inversion of the WCM for SSM retrieval over wheat fields were proposed using Sentinel-1 radar data only. To this objective, the dry above-ground biomass (AGB) and the vegetation water content (VWC) were retrieved from the interferometric coherence and the PR. The relationships were then used as the vegetation descriptor in the WCM. The best retrieval results were obtained using the relationship between ρVV and the AGB (R and RMSE of 0.82, 0.05 m3/m3 respectively and no bias). The new retrieval approaches were then applied to a large database covering a rainfed field in Morocco and 18 plots of rainfed and irrigated wheat of the Kairouan plain (Tunisia) and compared to other classical techniques of SSM retrieval including simple linear relationships between SSM and σ0. The method based on the WCM and the ρVV-AGB relationships also provided with slightly better results than the others on the validation database (r = 0.75, RMSE = 0.06 m3/m3 and bias = 0.01 m3/m3 over the 18 plots of Tunisia) but the simple linear relationships performed also reasonably well (r = 0.62, RMSE = 0.07, bias = −0.01 in Tunisia for instance). This study opens perspectives for high resolution soil moisture mapping from Sentinel-1 data over south Mediterranean wheat crops and in fine, for irrigation scheduling and retrieval through the assimilation of these new products in an evapotranspiration model.
•Surface soil moisture retrieval on wheat crops from Sentinel-1 data only.•Use of interferometric coherence and polarization ratio as vegetation descriptors.•Comparison to classical approaches of SSM retrieval.•Evaluation on a large database covering rainfed and irrigated plots in North Africa.
The recently launched Sentinel-1 satellite with a Synthetic Aperture Radar (SAR) sensor onboard offers a powerful tool for irrigation monitoring under various weather conditions, with high spatial ...and temporal resolution. This research discusses the potential of different metrics calculated from the Sentinel-1 time series for mapping irrigated fields. A methodology for irrigation mapping using SAR data is proposed. The study is performed using VV (vertical–vertical) and VH (vertical–horizontal) polarizations over an agricultural site in Urgell, Catalunya (Spain). With field segmentation information from SIGPAC (the Geographic Information System for Agricultural Parcels), the backscatter intensities are averaged within each field. From the Sentinel-1 time series for each field, the statistics and metrics, including the mean value, the variance of the signal, the correlation length, and the fractal dimension, are analyzed. With the Support Vector Machine (SVM), the classification of irrigated crops, irrigated trees, and non-irrigated fields is performed with the metrics vector. The results derived from the SVM are validated with ground truthing from SIGPAC over the whole study area, with a good overall accuracy of 81.08%. Random Forest (RF) machine classification is also tested in this study, which gives an accuracy of around 82.2% when setting the tree depth at three. The methodology is based only on SAR data, which makes it applicable to all areas, even with frequent cloud cover, but this method may be less robust when irrigation is less dominated to soil moisture change.
Comprehensive knowledge about irrigation timing is crucial for water resource management. This paper presents a comparative analysis between C- and L-band Synthetic Aperture Radar (SAR) data for the ...detection of irrigation events. The analysis was performed using C-band time series data derived from the Sentinel-1 (S1) satellite and two L-band images from the PALSAR-2 (ALOS-2) sensor acquired over irrigated grassland plots in the Crau plain of southeast France. The S1 C-band time series was first analyzed as a function of rainfall and irrigation events. The backscattering coefficients in both the L and C bands were then compared to the time difference between the date of the acquired SAR image and the date of the last irrigation event occurring before the SAR acquisition (Δt). Sensitivity analysis was performed for 2 classes of the Normalized Difference Vegetation Index (NDVI ≤0.7 and NDVI >0.7). The main results showed that when the vegetation is moderately developed (NDVI ≤0.7), the C-band temporal variation remains sensitive to the soil moisture dynamics and the irrigation events could be detected. The C-VV signal decreases due to the drying out of the soil when the time difference between the S1 image and irrigation event increases. For well-developed vegetation cover (NDVI >0.7), the C-band sensitivity to irrigation events becomes dependent on the crop type. For well-developed Gramineae grass with longs stalks and seedheads, the C band shows no correlation with Δt due to the absence of the soil contribution in the backscattered signal, contrary to the legume grass type, where the C band shows a good correspondence between C-VV and Δt for NDVI > 0.7. In contrast, analysis of the L-band backscattering coefficient shows that the L band remains sensitive to the soil moisture regardless of the vegetation cover development and the vegetation characteristics, thus being more suitable for irrigation detection than the C band. The L-HH signal over Gramineae grass or legume grass types shows the same decreasing pattern with the increase in Δt, regardless of the NDVI-values, presenting a decrease in soil moisture with time and thus high sensitivity of the radar signal to soil parameters. Finally, the co-polarizations for both the C and L bands (L-HH and C-VV) tend to be more adequate for irrigation detection than the HV cross-polarization, as they show higher sensitivity to soil moisture values.
The main objective of this work was to retrieve surface soil moisture (SSM) by using scattering models and a support vector machine (SVM) technique driven by backscattering coefficients obtained from ...Sentinel-1 satellite images acquired over bare agricultural soil in the Tensfit basin of Morocco. Two backscattering models were selected in this study due to their wide use in inversion procedures: the theoretical integral equation model (IEM) and the semi-empirical model (Oh). To this end, the sensitivity of the SAR backscattering coefficients at V V ( σ v v ∘ ) and V H ( σ v h ∘ ) polarizations to in situ soil moisture data were analyzed first. As expected, the results showed that over bare soil the σ v v ∘ was well correlated with SSM compared to the σ v h ∘ , which showed more dispersion with correlation coefficients values (r) of about 0.84 and 0.61 for the V V and V H polarizations, respectively. Afterwards, these values of σ v v ∘ were compared to those simulated by the backscatter models. It was found that IEM driven by the measured length correlation L slightly underestimated SAR backscatter coefficients compared to the Oh model with a bias of about − 0.7 dB and − 1.2 dB and a root mean square (RMSE) of about 1.1 dB and 1.5 dB for Oh and IEM models, respectively. However, the use of an optimal value of L significantly improved the bias of IEM, which became near to zero, and the RMSE decreased to 0.9 dB. Then, a classical inversion approach of σ v v ∘ observations based on backscattering model is compared to a data driven retrieval technic (SVM). By comparing the retrieved soil moisture against ground truth measurements, it was found that results of SVM were very encouraging and were close to those obtained by IEM model. The bias and RMSE were about 0.28 vol.% and 2.77 vol.% and − 0.13 vol.% and 2.71 vol.% for SVM and IEM, respectively. However, by taking into account the difficultly of obtaining roughness parameter at large scale, it was concluded that SVM is still a useful tool to retrieve soil moisture, and therefore, can be fairly used to generate maps at such scales.
This paper presents a technique for the mapping of soil moisture and irrigation, at the scale of agricultural fields, based on the synergistic interpretation of multi-temporal optical and Synthetic ...Aperture Radar (SAR) data (Sentinel-2 and Sentinel-1). The Kairouan plain, a semi-arid region in central Tunisia (North Africa), was selected as a test area for this study. Firstly, an algorithm for the direct inversion of the Water Cloud Model (WCM) was developed for the spatialization of the soil water content between 2015 and 2017. The soil moisture retrieved from these observations was first validated using ground measurements, recorded over 20 reference fields of cereal crops. A second method, based on the use of neural networks, was also used to confirm the initial validation. The results reported here show that the soil moisture products retrieved from remotely sensed data are accurate, with a Root Mean Square Error (RMSE) of less than 5% between the two moisture products. In addition, the analysis of soil moisture and Normalized Difference Vegetation Index (NDVI) products over cultivated fields, as a function of time, led to the classification of irrigated and rainfed areas on the Kairouan plain, and to the production of irrigation maps at the scale of individual fields. This classification is based on a decision tree approach, using a combination of various statistical indices of soil moisture and NDVI time series. The resulting irrigation maps were validated using reference fields within the study site. The best results were obtained with classifications based on soil moisture indices only, with an accuracy of 77%.