Despite irrigation being one of the main sources of anthropogenic water consumption, detailed information about water amounts destined for this purpose are often lacking worldwide. In this study, a methodology which can be used to estimate irrigation amounts over a pilot area in Spain by exploiting remotely sensed soil moisture is proposed. Two high-resolution DISPATCH (DISaggregation based on Physical And Theoretical scale CHange) downscaled soil moisture products have been used: SMAP (Soil Moisture Active Passive) and SMOS (Soil Moisture and Ocean Salinity) at 1 km. The irrigation estimates have been obtained through the SM2RAIN algorithm, in which the evapotranspiration term has been improved to adequately reproduce the crop evapotranspiration over irrigated areas according to the FAO (Food and Agriculture Organization) model. The experiment exploiting the SMAP data at 1 km represents the main work analyzed in this study and covered the period from January 2016 to September 2017. The experiment with the SMOS data at 1 km, for which a longer time series is available, allowed the irrigation estimates to be extended back to 2011. For both of the experiments carried out, the proposed method performed well in reproducing the magnitudes of the irrigation amounts that actually occurred in four of the five pilot irrigation districts. The SMAP experiment, for which a more detailed analysis was performed, also provided satisfactory results in representing the spatial distribution and the timing of the irrigation events. In addition, the investigation into which term of the SM2RAIN algorithm plays the leading role in determining the amount of water entering into the soil highlights the importance of correct representation of the evapotranspiration process.

Due to their image-based nature, “contextual” approaches are very attractive to estimate evapotranspiration (ET) from remotely-sensed land surface temperature (LST) data. Their application is however limited to highly heterogeneous areas where the soil and vegetation temperature endmembers (Tends) can be observed at the thermal sensor resolution. This paper aims to develop a simple theoretical approach to estimate Tends independently from LST images. Soil Tends are simulated by a soil energy balance model forced by meteorological data. Vegetation Tends are obtained from soil Tends and air temperature. Model-derived soil Tends are first evaluated with in situ measurements made over an irrigated area in Morocco. The root mean square difference (RMSD) between modeled and ground-based soil Tends is estimated as 2.4 ºC. Model-derived soil Tends are next compared with the soil Tends retrieved from 90-m resolution ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data collected over two irrigated areas in Mexico and Spain. Such a comparison reveals a strong consistency between model-derived and high-resolution image-based soil Tends. A recent contextual ET model (SEB-1S) is then applied to 90-m resolution and to 1-km resolution (aggregated) ASTER data using the model-derived or image-based Tends as the input. The RMSD between 90-m resolution SEB-1S and in situ ET is estimated as 65 and 82 W·m^-2, and the RMSD between 1-km resolution SEB-1S and aggregated SEB-1S ET is estimated as 78 and 56 W·m^-2, for the image-based and model-derived Tends, respectively. In light of the above results, Tends should be estimated a priori when contextual models are applied to low resolution images. Moreover, the consistency over highly heterogeneous areas between model-derived and high-resolution image-based Tends provides a meaningful basis for developing mixed modeling observational approaches.

Root-zone soil moisture (RZSM) plays a key role for most water and energy budgets, as it is particularly relevant in controlling plant transpiration and hydraulic redistribution. RZSM data is needed for a variety of different applications, such as forecasting crop yields, improving flood predictions and monitoring agricultural drought, among others. Remote sensing provides surface soil moisture (SSM) retrievals, whose key advantage is the large spatial coverage on a systematic basis. This study tests a simple method to retrieve RZSM estimates from high-resolution SSM derived from SMAP (Soil Moisture Active Passive). A recursive exponential filter using a time constant $\tau$ is calibrated per land cover type, which uses as an intermediate step a long-term ISBA-DIF (Interaction Soil Biosphere Atmosphere—Diffusion scheme) dataset over an area located in Catalonia, NE of Spain. The $\tau$ values thus obtained are then used as an input to the same recursive exponential filter, to derive 1 km resolution RZSM estimates from 1 km SMAP SSM, which are obtained from the original data by downscaling to a 1 km resolution, through the DISPATCH (DISaggregation based on a Physical and Theoretical scale CHange) methodology. The results are then validated with scaled in situ observations at different depths, over two different areas, one representative of rainfed crops, and the other of irrigated crops. In general, the estimates agree well with the observations over the rainfed crops, especially at a 10 cm and 25 cm depth. Nash–Sutcliffe (NS) scores ranging between 0.33 and 0.58, and between 0.37 and 0.56 have been found, respectively. Correlation coefficients for these depths are high, between 0.76 and 0.91 (10 cm), and between 0.71 and 0.90 (25 cm). For the irrigated sites, results are poorer (partly due to the extremely high heterogeneity present), with NS scores ranging between −2.57 and 0.16, and correlations ranging between −0.56 and 0.48 at 25 cm. Given the strong correlations and NS scores found in the surface, the sensitivity of the filter to different $\tau$ values was investigated. For the rainfed site, it was found, as expected, with increasing $\tau$, increasing NS and correlations with the deeper layers, suggesting a better coupling. Nevertheless, a strong correlation with the surface (5 cm) or shallower depths (10 cm) observed over certain sites indicates a certain lack of skill of the filter to represent processes which occur at lower levels in the SM column. All in all, a calibration accounting for the vegetation was shown to be an adequate methodology in applying the recursive exponential filter to derive the RZSM estimates over large areas. Nevertheless, the relative shallow surface at which the estimates correlate in some cases seem to indicate that an effect of evapotranspiration in the profile is not well captured by the filter.

This project explored the integrated use of satellite, ground observations and hydrological distributed models to support water resources assessment and monitoring in High Mountain Asia (HMA). Hydrological data products were generated taking advantage of the synergies of European and Chinese data assets and space-borne observation systems. Energy-budget-based glacier mass balance and hydrological models driven by satellite observations were developed. These models can be applied to describe glacier-melt contribution to river flow. Satellite hydrological data products were used for forcing, calibration, validation and data assimilation in distributed river basin models. A pilot study was carried out on the Red River basin. Multiple hydrological data products were generated using the data collected by Chinese satellites. A new Evapo-Transpiration (ET) dataset from 2000 to 2018 was generated, including plant transpiration, soil evaporation, rainfall interception loss, snow/ice sublimation and open water evaporation. Higher resolution data were used to characterize glaciers and their response to environmental forcing. These studies focused on the Parlung Zangbo Basin, where glacier facies were mapped with GaoFeng (GF), Sentinal-2/Multi-Spectral Imager (S2/MSI) and Landsat8/Operational Land Imager (L8/OLI) data. The geodetic mass balance was estimated between 2000 and 2017 with Zi-Yuan (ZY)-3 Stereo Images and the SRTM DEM. Surface velocity was studied with Landsat5/Thematic Mapper (L5/TM), L8/OLI and S2/MSI data over the period 2013–2019. An updated method was developed to improve the retrieval of glacier albedo by correcting glacier reflectance for anisotropy, and a new dataset on glacier albedo was generated for the period 2001–2020. A detailed glacier energy and mass balance model was developed with the support of field experiments at the Parlung No. 4 Glacier and the 24 K Glacier, both in the Tibetan Plateau. Besides meteorological measurements, the field experiments included glaciological and hydrological measurements. The energy balance model was formulated in terms of enthalpy for easier treatment of water phase transitions. The model was applied to assess the spatial variability in glacier melt. In the Parlung No. 4 Glacier, the accumulated glacier melt was between 1.5 and 2.5 m w.e. in the accumulation zone and between 4.5 and 6.0 m w.e. in the ablation zone, reaching 6.5 m w.e. at the terminus. The seasonality in the glacier mass balance was observed by combining intensive field campaigns with continuous automatic observations. The linkage of the glacier and snowpack mass balance with water resources in a river basin was analyzed in the Chiese (Italy) and Heihe (China) basins by developing and applying integrated hydrological models using satellite retrievals in multiple ways. The model FEST-WEB was calibrated using retrievals of Land Surface Temperature (LST) to map soil hydrological properties. A watershed model was developed by coupling ecohydrological and socioeconomic systems. Integrated modeling is supported by an updated and parallelized data assimilation system. The latter exploits retrievals of brightness temperature (Advanced Microwave Scanning Radiometer, AMSR), LST (Moderate Resolution Imaging Spectroradiometer, MODIS), precipitation (Tropical Rainfall Measuring Mission (TRMM) and FengYun (FY)-2D) and in-situ measurements. In the case study on the Red River Basin, a new algorithm has been applied to disaggregate the SMOS (Soil Moisture and Ocean Salinity) soil moisture retrievals by making use of the correlation between evaporative fraction and soil moisture.

The resolution of current satellite surface soil moisture (SM) estimates is very low, of tens of kilometers, which proves to be insufficient for various agricultural and hydrological applications. Amongst the existing downscaling approaches of remotely sensed SM, DISPATCH (DISaggregation based on a Physical And Theoretical scale CHange) improves the resolution of SMOS (Soil Moisture and Ocean Salinity) soil moisture data using soil evaporative efficiency (SEE) estimates at high resolution (HR) and a SEE(SM) model implemented at low resolution (LR). Defined as the ratio of actual to potential soil evaporation, SEE can be derived from the remotely sensed land surface temperature (LST) and normalized difference vegetation index (NDVI). The current version of DISPATCH uses a linear SEE(SM) model. This study aims at improving the SEE(SM) model and testing different calibration strategies, to ultimately have more robust and better downscaled SM products. A nonlinear SEE(SM) model is introduced and its influence on the derived HR SM products is studied over a range of conditions. Each model, linear and nonlinear, is calibrated from remote sensing data on a daily and a multi-date basis. The approaches were tested over two mixed dry and irrigated areas in Catalonia, Spain, and over one dry area in Morocco. When using the linear model, better statistical results were generally obtained using a daily calibration (current version of DISPATCH), most notably over one Spanish site. However, the best results were systematically obtained for an annually calibrated nonlinear model, in terms of all metrics considered: correlation coefficient, slope of the linear regression, bias, unbiased root mean square error. In particular, when using the annually calibrated nonlinear SEE (SM) model, the temporal slope of the linear regression between disaggregated and in situ soil moisture increased to 1.16 and 0.75 for one Spanish site and for the Moroccan site (as opposed to 0.44 and 0.58, respectively, when using the linear model with a daily calibration). The temporal correlation coefficient increased to 0.47 and 0.54 over the Spanish sites (as opposed to 0.18 and 0.27, respectively, when using the linear model with a daily calibration). Those contrasted results indicate compensation effects between the model type and the calibration strategy. Taking into account studies that report the strong nonlinear behavior of the SEE with respect to SM, the introduction of the nonlinear SEE(SM) model in DISPATCH, combined with a multi-date calibration, is proven to perform significantly better under various conditions, leading to more robust disaggregated SM products. The SEE modeling based on the nonlinear SM model, with a multi-date calibration, could be integrated into the CATDS—Centre Aval de Traitement des Données SMOS as a future product, as well as into existing evapotranspiration models, which are based on a combination of thermal and microwave data.