An accurate assessment of evapotranspiration (ET) is crucially needed at the basin scale for studying the hydrological processes and water balance especially from upstream to downstream. In the mountains, this term is poorly understood because of various challenges, including the vegetation complexity, plant diversity, lack of available data and because the in situ direct measurement of ET is difficult in complex terrain. The main objective of this work was to investigate the potential of a Two-Source-Energy-Balance model (TSEB) driven by the Landsat and MODIS data for estimating ET over a complex mountain region. The complexity is associated with the type of the vegetation canopy as well as the changes in topography. For validating purposes, a large-aperture scintillometer (LAS) was set up over a heterogeneous transect of about 1.4 km to measure sensible (H) and latent heat (LE) fluxes. Additionally, two towers of eddy covariance (EC) systems were installed along the LAS transect. First, the model was tested at the local scale against the EC measurements using multi-scale remote sensing (MODIS and Landsat) inputs at the satellite overpasses. The obtained averaged values of the root mean square error (RMSE) and correlation coefficient (R) were about 72.4 Wm⁻² and 0.79 and 82.0 Wm⁻² and 0.52 for Landsat and MODIS data, respectively. Secondly, the potential of the TSEB model for evaluating the latent heat fluxes at large scale was investigated by aggregating the derived parameters from both satellites based on the LAS footprint. As for the local scale, the comparison of the latent heat fluxes simulated by TSEB driven by Landsat data performed well against those measured by the LAS (R = 0.69, RMSE = 68.0 Wm⁻²), while slightly more scattering was observed when MODIS products were used (R = 0.38, RMSE = 99.8 Wm⁻²). Based on the obtained results, it can be concluded that (1) the TSEB model can be fairly used to estimate the evapotranspiration over the mountain regions; and (2) medium- to high-resolution inputs are a better option than coarse-resolution products for describing this kind of complex terrain.

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.

In this study, a simple model, based on a light-use-efficiency model, was developed in order to estimate growth and yield of the irrigated winter wheat under semi-arid conditions. The originality of the proposed method consists in (1) the modifying of the expression of the conversion coefficient (ε_conv) by integrating an appropriate stress threshold (k_sconv) for triggering irrigation, (2) the substitution of the product of the two maximum coefficients of interception (ε_imax) and conversion (ε_{conv_max}) by a single parameter ε_max, (3) the modeling of ε_max as a function of the Cumulative Growing Degree Days (CGDD) since sowing date, and (4) the dynamic expression of the harvest index (HI) as a function of the CGDD and the final harvest index (HI₀) depending on the maximum value of the Normalized Difference Vegetation Index (NDVI). The calibration and validation of the proposed model were performed based on the observations of wheat dry matter (DM) and grain yield (GY) which were collected on the R3 irrigated district of the Haouz plain (center of Morocco), during three agricultural seasons. Further, the outputs of the simple model were also evaluated against the AquaCrop model estimates. The model calibration allowed the parameterization of ε_max in four periods according to the wheat phenological stages. By contrast, a linear evolution was sufficient to represent the relationship between HI and CGDD. For the model validation, the obtained results showed a good agreement between the estimated and observed values with a Root Mean Square Error (RMSE) of about 1.07 and 0.57 t/ha for DM and GY, respectively. These correspond to a relative RMSE of about 19% for DM and 20% for GY. Likewise, although of its simplicity, the accuracy of the proposed model seems to be comparable to that of the AquaCrop model. For GY, R², and RMSE values were respectively 0.71 and 0.44 t/ha for the developed approach and 0.88 and 0.37 t/ha for AquaCrop. Thus, the proposed simple light-use-efficiency model can be considered as a useful tool to correctly reproduce DM and GY values.