Global soil moisture (SM) products are currently available from passive microwave sensors at typically 40 km spatial resolution. Although recent efforts have been made to produce 1 km resolution data from the disaggregation of coarse scale observations, the targeted resolution of available SM data is still far from the requirements of fine-scale hydrological and agricultural studies. To fill the gap, a new disaggregation scheme of Soil Moisture Active and Passive (SMAP) data is proposed at 100 m resolution by using the disaggregation based on physical and theoretical scale change (DISPATCH) algorithm. The main objectives of this paper is (i) to implement DISPATCH algorithm at 100 m resolution using SMAP SM and Landsat land surface temperature and vegetation index data and (ii) to investigate the usefulness of an intermediate spatial resolution (ISR) between the SMAP 36 km resolution and the targeted 100 m resolution. The sequential disaggregation approach from 36 km to ISR (ranging from 1 km to 30 km) and from ISR to 100 m resolution is evaluated over 22 irrigated field crops in central Morocco using in-situ SM measurements collected from January to May 2016. The lowest root mean square difference (RMSD) between the 100 m resolution disaggregated and in-situ SM is obtained when the ISR is around 10 km. Therefore, the two-step disaggregation is more efficient than the direct disaggregation from SMAP to 100 m resolution. Moreover, we propose a moving average window algorithm to increase the accuracy in the 100 m resolution SM as well as to reduce the low-resolution boxy artifacts on disaggregated images. The correlation coefficient between 100 m resolution disaggregated and in situ SM ranges between 0.5–0.9 for four out of the six extensive sampling dates. This methodology relies solely on remote sensing data and can be easily implemented to monitor SM at a high spatial resolution over irrigated regions.

Abstract. Thermal-based two-source energy balance modeling is essential to estimate the land evapotranspiration (ET) in a wide range of spatial and temporal scales. However, the use of thermal-derived land surface temperature (LST) is not sufficient to simultaneously constrain both soil and vegetation flux components. Therefore, assumptions (about either soil or vegetation fluxes) are commonly required. To avoid such assumptions, an energy balance model, TSEB-SM, was recently developed by Ait Hssaine et al. (2018b) in order to consider the microwave-derived near-surface soil moisture (SM), in addition to the thermal-derived LST and vegetation cover fraction (fc) normally used. While TSEB-SM has been successfully tested using in situ measurements, this paper represents its first evaluation in real life using 1 km resolution satellite data, comprised of MODIS (MODerate resolution Imaging Spectroradiometer) for LST and fc data and 1 km resolution SM data disaggregated from SMOS (Soil Moisture and Ocean Salinity) observations. The approach is applied during a 4-year period (2014–2018) over a rainfed wheat field in the Tensift basin, central Morocco. The field used was seeded for the 2014–2015 (S1), 2016–2017 (S2) and 2017–2018 (S3) agricultural seasons, while it remained unploughed (as bare soil) during the 2015–2016 (B1) agricultural season. The classical TSEB model, which is driven only by LST and fc data, significantly overestimates latent heat fluxes (LE) and underestimates sensible heat fluxes (H) for the four seasons. The overall mean bias values are 119, 94, 128 and 181 W m−2 for LE and −104, −71, −128 and −181 W m−2 for H, for S1, S2, S3 and B1, respectively. Meanwhile, when using TSEB-SM (SM and LST combined data), these errors are significantly reduced, resulting in mean bias values estimated as 39, 4, 7 and 62 W m−2 for LE and −10, 24, 7, and −59 W m−2 for H, for S1, S2, S3 and B1, respectively. Consequently, this finding confirms again the robustness of the TSEB-SM in estimating latent/sensible heat fluxes at a large scale by using readily available satellite data. In addition, the TSEB-SM approach has the original feature to allow for calibration of its main parameters (soil resistance and Priestley–Taylor coefficient) from satellite data uniquely, without relying either on in situ measurements or on a priori parameter values.

DISPATCH is a disaggregation algorithm of the low-resolution soil moisture (SM) estimates derived from passive microwave observations. It provides disaggregated SM data at typically 1 km resolution by using the soil evaporative efficiency (SEE) estimated from optical/thermal data collected around solar noon. DISPATCH is based on the relationship between the evapo-transpiration rate and the surface SM under non-energy-limited conditions and hence is well adapted for semi-arid regions with generally low cloud cover and sparse vegetation. The objective of this paper is to extend the spatio-temporal coverage of DISPATCH data by 1) including more densely vegetated areas and 2) assessing the usefulness of thermal data collected earlier in the morning. Especially, we evaluate the performance of the Temperature Vegetation Dryness Index (TVDI) instead of SEE in the DISPATCH algorithm over vegetated areas (called vegetation-extended DISPATCH) and we quantify the increase in coverage using Sentinel-3 (overpass at around 09:30 am) instead of MODIS (overpass at around 10:30 am and 1:30 pm for Terra and Aqua, respectively) data. In this study, DISPATCH is applied to 36 km resolution Soil Moisture Active and Passive SM data over three 50 km by 50 km areas in Spain and France to assess the effectiveness of the approach over temperate and semi-arid regions. The use of TVDI within DISPATCH increases the coverage of disaggregated images by 9 and 14% over the temperate and semi-arid sites, respectively. Moreover, including the vegetated pixels in the validation areas increases the overall correlation between satellite and in situ SM from 0.36 to 0.43 and from 0.41 to 0.79 for the temperate and semi-arid regions, respectively. The use of Sentinel-3 can increase the spatio-temporal coverage by up to 44% over the considered MODIS tile, while the overlapping disaggregated data sets derived from Sentinel-3 and MODIS land surface temperature data are strongly correlated (around 0.7). Additionally, the correlation between satellite and in situ SM is significantly better for DISPATCH (0.39–0.80) than for the Copernicus Sentinel-1-based (−0.03 to 0.69) and SMAP/S1 (0.37–0.74) product over the three studies (temperate and semi-arid) areas, with an increase in yearly valid retrievals for the vegetation-extended DISPATCH algorithm.

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 ‘SMO20’, 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 ‘SME16’; and the third approach disaggregates the soil moisture active and passive SM at 100 m resolution using Landsat optical/thermal data ‘SMO19’. 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/m2) and underestimates H (mean bias of −110 W/m2), 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 R2 of 0.72 and 0.92, and a mean RMSE of 31 and 36 W/m2 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.

Abstract. Thermal-based two-source energy balance modeling is very useful for estimating the land evapotranspiration (ET) at a wide range of spatial and temporal scales. However, the land surface temperature (LST) is not sufficient for constraining simultaneously both soil and vegetation flux components in such a way that assumptions (on either the soil or the vegetation fluxes) are commonly required. To avoid such assumptions, a new energy balance model (TSEB-SM) was recently developed in Ait Hssaine et al. (2018a) to integrate the microwave-derived near-surface soil moisture (SM), in addition to the thermal-derived LST and vegetation cover fraction (fc). Whereas, TSEB-SM has been recently tested using in-situ measurements, the objective of this paper is to evaluate the performance of TSEB-SM in real-life using 1 km resolution MODIS (Moderate resolution imaging spectroradiometer) LST and fc data and the 1 km resolution SM data disaggregated from SMOS (Soil Moisture and Ocean Salinity) observations by using DisPATCh. The approach is applied during a four-year period (2014–2018) over a rainfed wheat field in the Tensift basin, central Morocco, during a four-year period (2014–2018). The field was seeded for the 2014–2015 (S1), 2016–2017 (S2) and 2017–2018 (S3) agricultural season, while it was not ploughed (remained as bare soil) during the 2015–2016 (B1) agricultural season. The mean retrieved values of (arss, brss) calculated for the entire study period using satellite data are (7.32, 4.58). The daily calibrated αPT ranges between 0 and 1.38 for both S1 and S2. Its temporal variability is mainly attributed to the rainfall distribution along the agricultural season. For S3, the daily retrieved αPT remains at a mostly constant value (∼ 0.7) throughout the study period, because of the lack of clear sky disaggregated SM and LST observations during this season. Compared to eddy covariance measurements, TSEB driven only by LST and fc data significantly overestimates latent heat fluxes for the four seasons. The overall mean bias values are 119, 94, 128 and 181 W/m2 for S1, S2, S3 and B1 respectively. In contrast, these errors are much reduced when using TSEB-SM (SM and LST combined data) with the mean bias values estimated as 39, 4, 7 and 62 W/m2 for S1, S2, S3 and B1 respectively.

No abstract

The use of soil moisture (SM) measurements from satellites has grown in recent years, fostering the development of new products at high resolution. This opens the possibility of using them for certain applications that were normally carried out using in situ data. We investigated this hypothesis through two main analyses using two high-resolution satellite-based soil moisture (SBSM) products that combined microwave with thermal and optical data: (1) The Disaggregation based on Physical And Theoretical scale Change (DISPATCH) and, (2) The Soil Moisture Ocean Salinity-Barcelona Expert Center (SMOS-BEC Level 4). We used these products to analyse the SM differences among pixels with contrasting vegetation. This was done through the comparison of the SM measurements from satellites and the measurements simulated with a simple antecedent precipitation index (API) model, which did not account for the surface characteristics. Subsequently, the deviation of the SM from satellite with respect to the API model (bias) was analysed and compared for contrasting land use categories. We hypothesised that the differences in the biases of the varied categories could provide information regarding the water retention capacity associated with each type of vegetation. From the satellite measurements, we determined how the SM depended on the tree cover, i.e., the denser the tree cover, the higher the SM. However, in winter periods with light rain events, the tree canopy could dampen the moistening of the soil through interception and conducted higher SM in the open areas. This evolution of the SM differences that depended on the characteristics of each season was observed both from satellite and from in situ measurements taken beneath a tree and in grass on the savanna landscape. The agreement between both types of measurements highlighted the potential of the SBSM products to investigate the SM of each type of vegetation. We found that the results were clearer for DISPATCH, whose data was not smoothed spatially as it was in SMOS-BEC. We also tested whether the relationships between SM and evapotranspiration could be investigated using satellite data. The answer to this question was also positive but required removing the unrealistic high-frequency SM oscillations from the satellite data using a low pass filter. This improved the performance scores of the products and the agreement with the results from the in situ data. These results demonstrated the possibility of using SM data from satellites to substitute ground measurements for the study of land–atmosphere interactions, which encourages efforts to improve the quality and resolution of these measurements.

Soil moisture (SM) data are required at high spatio-temporal resolution—typically the crop field scale every 3–6 days—for agricultural and hydrological purposes. To provide such high-resolution SM data, many remote sensing methods have been developed from passive microwave, active microwave and thermal data. Despite the pros and cons of each technique in terms of spatio-temporal resolution and their sensitivity to perturbing factors such as vegetation cover, soil roughness and meteorological conditions, there is currently no synergistic approach that takes advantage of all relevant (passive, active microwave and thermal) remote sensing data. In this context, the objective of the paper is to develop a new algorithm that combines SMAP L-band passive microwave, MODIS/Landsat optical/thermal and Sentinel-1 C-band radar data to provide SM data at the field scale at the observation frequency of Sentinel-1. In practice, it is a three-step procedure in which: (1) the 36 km resolution SMAP SM data are disaggregated at 100 m resolution using MODIS/Landsat optical/thermal data on clear sky days, (2) the 100 m resolution disaggregated SM data set is used to calibrate a radar-based SM retrieval model and (3) the so-calibrated radar model is run at field scale on each Sentinel-1 overpass. The calibration approach also uses a vegetation descriptor as ancillary data that is derived either from optical (Sentinel-2) or radar (Sentinel-1) data. Two radar models (an empirical linear regression model and a non-linear semi-empirical formulation derived from the water cloud model) are tested using three vegetation descriptors (NDVI, polarization ratio (PR) and radar coherence (CO)) separately. Both models are applied over three experimental irrigated and rainfed wheat crop sites in central Morocco. The field-scale temporal correlation between predicted and in situ SM is in the range of 0.66–0.81 depending on the retrieval configuration. Based on this data set, the linear radar model using PR as a vegetation descriptor offers a relatively good compromise between precision and robustness all throughout the agricultural season with only three parameters to set. The proposed synergistical approach combining multi-resolution/multi-sensor SM-relevant data offers the advantage of not requiring in situ measurements for calibration.