Irrigation represents one of the most impactful human interventions in the terrestrial water cycle. Knowing the distribution and extent of irrigated areas as well as the amount of water used for irrigation plays a central role in modeling irrigation water requirements and quantifying the impact of irrigation on regional climate, river discharge, and groundwater depletion. Obtaining high-quality global information about irrigation is challenging, especially in terms of quantification of the water actually used for irrigation. Here, we review existing Earth observation datasets, models, and algorithms used for irrigation mapping and quantification from the field to the global scale. The current observation capacities are confronted with the results of a survey on user requirements on satellite-observed irrigation for agricultural water resources’ management. Based on this information, we identify current shortcomings of irrigation monitoring capabilities from space and phrase guidelines for potential future satellite missions and observation strategies.

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 $VV$ ( ${\sigma }_{vv}^{\circ }$ ) and $VH$ ( ${\sigma }_{vh}^{\circ }$ ) polarizations to in situ soil moisture data were analyzed first. As expected, the results showed that over bare soil the ${\sigma }_{vv}^{\circ }$ was well correlated with SSM compared to the ${\sigma }_{vh}^{\circ }$ , which showed more dispersion with correlation coefficients values (r) of about $0.84$ and $0.61$ for the $VV$ and $VH$ polarizations, respectively. Afterwards, these values of ${\sigma }_{vv}^{\circ }$ 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 ${\sigma }_{vv}^{\circ }$ 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.

Agricultural water use represents more than 70% of the world’s freshwater through irrigation water inputs that are poorly known at the field scale. Irrigation monitoring is thus an important issue for optimizing water use in particular with regards to the water scarcity that the semi-arid regions are already facing. In this context, the aim of this study is to develop and evaluate a new approach to predict seasonal to daily irrigation timing and amounts at the field scale. The method is based on surface soil moisture (SSM) data assimilated into a simple land surface (FAO-56) model through a particle filter technique based on an ensemble of irrigation scenarios. The approach is implemented in three steps. First, synthetic experiments are designed to assess the impact of the frequency of observation, the errors on SSM and the a priori constraints on the irrigation scenarios for different irrigation techniques (flooding and drip). In a second step, the method is evaluated using in situ SSM measurements with different revisit times (3, 6 and 12 days) to mimic the available SSM product derived from remote sensing observation. Finally, SSM estimates from Sentinel-1 are used. Data are collected on different wheat fields grown in Morocco, for both flood and drip irrigation techniques in addition to rainfed fields used for an indirect evaluation of the method performance. Using in situ data, accurate results are obtained. With an observation every 6 days to mimic the Sentinel-1 revisit time, the seasonal amounts are retrieved with R > 0.98, RMSE < 32 mm and bias < 2.5 mm. Likewise, a good agreement is observed at the daily scale for flood irrigation as more than 70% of the detected irrigation events have a time difference from actual irrigation events shorter than 4 days. Over the drip irrigated fields, the statistical metrics are R = 0.74, RMSE = 24.8 mm and bias = 2.3 mm for irrigation amounts cumulated over 15 days. When using SSM products derived from Sentinel-1 data, the statistical metrics on 15-day cumulated amounts slightly dropped to R = 0.64, RMSE = 28.7 mm and bias = 1.9 mm. The metrics on the seasonal amount retrievals are close to assimilating in situ observations with R = 0.99, RMSE = 33.5 mm and bias = −18.8 mm. Finally, among four rainfed seasons, only one false event was detected. This study opens perspectives for the regional retrieval of irrigation amounts and timing at the field scale and for mapping irrigated/non irrigated areas.

Over semi-arid agricultural areas, the surface energy balance and its components are largely dependent on the soil water availability. In such conditions, the land surface temperature (LST) retrieved from the thermal bands has been commonly used to represent the high spatial variability of the surface evaporative fraction and associated fluxes. In contrast, however, the soil moisture (SM) retrieved from microwave data has rarely been used thus far due to the unavailability of high-resolution (field scale) SM products until recent times. Soil evaporation is controlled by the surface SM. Moreover, the surface SM dynamics is temporally related to root zone SM, which provides information about the water status of plants. The aim of this work was to assess the gain in terms of flux estimates when integrating microwave-derived SM data in a thermal-based energy balance model at the field scale. In this study, SM products were derived from three different methodologies: the first approach inverts SM, labeled hereafter as ‘SM_O20’, from the backscattering coefficient and the interferometric coherence derived from Sentinel-1 products in the water cloud model (WCM); the second approach inverts SM from Sentinel-1 and Sentinel-2 data based on machine learning algorithms trained on a synthetic dataset simulated by the WCM noted ‘SM_E16’; and the third approach disaggregates the soil moisture active and passive SM at 100 m resolution using Landsat optical/thermal data ‘SM_O19’. These SM products, combined with the Landsat based vegetation index and LST, are integrated simultaneously within an energy balance model (TSEB-SM) to predict the latent (LE) and sensible (H) heat fluxes over two irrigated and rainfed wheat crop sites located in the Haouz Plain in the center of Morocco. H and LE were measured over each site using an eddy covariance system and their values were used to evaluate the potential of TSEB-SM against the classical two source energy balance (TSEB) model solely based on optical/thermal data. Globally, TSEB systematically overestimates LE (mean bias of 100 W/m²) and underestimates H (mean bias of −110 W/m²), while TSEB-SM significantly reduces those biases, regardless of the SM product used as input. This is linked to the parameterization of the Priestley Taylor coefficient, which is set to α_PT = 1.26 by default in TSEB and adjusted across the season in TSEB-SM. The best performance of TSEB-SM was obtained over the irrigated field using the three retrieved SM products with a mean R² of 0.72 and 0.92, and a mean RMSE of 31 and 36 W/m² for LE and H, respectively. This opens up perspectives for applying the TSEB-SM model over extended irrigated agricultural areas to better predict the crop water needs at the field scale.

A better understanding of the hydrological functioning of irrigated crops using remote sensing observations is of prime importance in semi-arid areas where water resources are limited. Radar observations, available at high resolution and with a high revisit time since the launch of Sentinel-1 in 2014, have shown great potential for the monitoring of the water content of the upper soil and of the canopy. In this paper, a complete set of data for radar signal analysis is shared with the scientific community for the first time to our knowledge. The data set is composed of Sentinel-1 products and in situ measurements of soil and vegetation variables collected during three agricultural seasons over drip-irrigated winter wheat in the Haouz plain in Morocco. The in situ data gather soil measurements (time series of half-hourly surface soil moisture, surface roughness and agricultural practices) and vegetation measurements collected every week/2 weeks including aboveground fresh and dry biomasses, vegetation water content based on destructive measurements, the cover fraction, the leaf area index, and plant height. Radar data are the backscattering coefficient and the interferometric coherence derived from Sentinel-1 GRDH (Ground Range Detected High Resolution) and SLC (Single Look Complex) products, respectively. The normalized difference vegetation index derived from Sentinel-2 data based on Level-2A (surface reflectance and cloud mask) atmospheric-effects-corrected products is also provided. This database, which is the first of its kind made available open access, is described here comprehensively in order to help the scientific community to evaluate and to develop new or existing remote sensing algorithms for monitoring wheat canopy under semi-arid conditions. The data set is particularly relevant for the development of radar applications including surface soil moisture and vegetation variable retrieval using either physically based or empirical approaches such as machine and deep learning algorithms. The database is archived in the DataSuds repository and is freely accessible via the following DOI: https://doi.org/10.23708/8D6WQC (Ouaadi et al., 2020a).

The difficulty of calculating the daily water budget of irrigated fields is often due to the uncertainty surrounding irrigation amounts and timing. The automated detection of irrigation events has the potential to greatly simplify this process, and the combination of high-resolution SAR (Sentinel-1) and optical satellite observations (Sentinel-2) makes the detection of irrigation events feasible through the use of a surface soil moisture (SSM) product. The motivation behind this study is to utilize a large irrigation dataset (collected during the ESA Irrigation + project over five sites in three countries over three years) to analyze the performance of an established algorithm and to test potential improvements. The study’s main findings are (1) the scores decrease with SSM observation frequency; (2) scores decrease as irrigation frequency increases, which was supported by better scores in France (more sprinkler irrigation) than in Germany (more localized irrigation); (3) replacing the original SSM model with the force-restore model resulted in an improvement of about 6% in the F-score and narrowed the error on cumulative seasonal irrigation; (4) the Sentinel-1 configuration (incidence angle, trajectory) did not show a significant impact on the retrieval of irrigation, which supposes that the SSM is not affected by these changes. Other aspects did not allow a definitive conclusion on the irrigation retrieval algorithm: (1) the lower scores obtained with small NDVI compared to large NDVI were counter-intuitive but may have been due to the larger number of irrigation events during high vegetation periods; (2) merging different runs and interpolating all SSM data for one run produced comparable F-scores, but the estimated cumulative sum of irrigation was around −20% lower compared to the reference dataset in the best cases.