This study aims to evaluate a remote sensing-based approach to allow estimation of the temporal and spatial distribution of crop evapotranspiration (ET) and irrigation water requirements over irrigated areas in semi-arid regions. The method is based on the daily step FAO-56 Soil Water Balance model combined with a time series of basal crop coefficients and the fractional vegetation cover derived from high-resolution satellite Normalized Difference Vegetation Index (NDVI) imagery. The model was first calibrated and validated at plot scale using ET measured by eddy-covariance systems over wheat fields and olive orchards representing the main crops grown in the study area of the Haouz plain (central Morocco). The results showed that the model provided good estimates of ET for wheat and olive trees with a root mean square error (RMSE) of about 0.56 and 0.54 mm/day respectively. The model was then used to compare remotely sensed estimates of irrigation requirements (RS-IWR) and irrigation water supplied (WS) at plot scale over an irrigation district in the Haouz plain through three growing seasons. The comparison indicated a large spatio-temporal variability in irrigation water demands and supplies; the median values of WS and RS-IWR were 130 (175), 117 (175) and 118 (112) mm respectively in the 2002–2003, 2005–2006 and 2008–2009 seasons. This could be attributed to inadequate irrigation supply and/or to farmers’ socio-economic considerations and management practices. The findings demonstrate the potential for irrigation managers to use remote sensing-based models to monitor irrigation water usage for efficient and sustainable use of water resources.

In mountainous regions, the scarcity of air temperature (Ta) measurements is a major limitation for hydrological and crop monitoring. An alternative to in situ measurements could be to downscale the reanalysis Ta data provided at high-temporal resolution. However, the relatively coarse spatial resolution of these products (i.e., 9 km for ERA5-Land) is unlikely to be directly representative of actual local Ta patterns. To address this issue, this study presents a new spatial downscaling strategy of hourly ERA5-Land Ta data with a three-step procedure. First, the 9 km resolution ERA5 Ta is corrected at its original resolution by using a reference Ta derived from the elevation of the 9 km resolution grid and an in situ estimate over the area of the hourly Environmental Lapse Rate (ELR). Such a correction of 9 km resolution ERA5 Ta is trained using several machine learning techniques, including Multiple Linear Regression (MLR), Support Vector Regression (SVR), and Extreme Gradient Boosting (Xgboost), as well as ancillary ERA5 data (daily mean, standard deviation, hourly ELR, and grid elevation). Next, the trained correction algorithms are run to correct 9 km resolution ERA5 Ta, and the corrected ERA5 Ta data are used to derive an updated ELR over the area (without using in situ Ta measurements). Third, the updated hourly ELR is used to disaggregate 9 km resolution corrected ERA5 Ta data at the 30-meter resolution of SRTM’s Digital Elevation Model (DEM). The effectiveness of this method is assessed across the northern part of the High Atlas Mountains in central Morocco through (1) k-fold cross-validation against five years (2016 to 2020) of in situ hourly temperature readings and (2) comparison with classical downscaling methods based on a constant ELR. Our results indicate a significant enhancement in the spatial distribution of hourly local Ta. By comparing our model, which included Xgboost, SVR, and MLR, with the constant ELR-based downscaling approach, we were able to decrease the regional root mean square error from approximately 3 ${}^{\circ }$C to 1.61 ${}^{\circ }$C, 1.75 ${}^{\circ }$C, and 1.8 ${}^{\circ }$C, reduce the mean bias error from −0.5 ${}^{\circ }$C to null, and increase the coefficient of determination from 0.88 to 0.97, 0.96, and 0.96 for Xgboost, SVR, and MLR, respectively.

This study aims to investigate the effects of an optimized sowing calendar for wheat over a surface irrigation scheme in the semi-arid region of Haouz (Morocco) on irrigation water requirements, crop growth and development and on yield. For that, a scenario-based simulation approach based on the covariance matrix adaptation–evolution strategy (CMA-ES) was proposed to optimize both the spatiotemporal distribution of sowing dates and the irrigation schedules, and then evaluate wheat crop using the 2011–2012 growing season dataset. Six sowing scenarios were simulated and compared to identify the most optimal spatiotemporal sowing calendar. The obtained results showed that with reference to the existing sowing patterns, early sowing of wheat leads to higher yields compared to late sowing (from 7.40 to 5.32 t/ha). Compared with actual conditions in the study area, the spatial heterogeneity is highly reduced, which increased equity between farmers. The results also showed that the proportion of plots irrigated in time can be increased (from 40% to 82%) compared to both the actual irrigation schedules and to previous results of irrigation optimization, which did not take into consideration sowing dates optimization. Furthermore, considerable reduction of more than 40% of applied irrigation water can be achieved by optimizing sowing dates. Thus, the proposed approach in this study is relevant for irrigation managers and farmers since it provides an insight on the consequences of their agricultural practices regarding the wheat sowing calendar and irrigation scheduling and can be implemented to recommend the best practices to adopt.

In semi-arid regions of Morocco, where the majority of water withdrawals are devoted to irrigation, optimizing irrigation practices in agriculture is a national priority in the face of recurring droughts and growing pressure on groundwater resources. However, the hydrological impacts of different drip-irrigation systems in the soil–plant–atmosphere continuum remain insufficiently understood. We monitored the stable isotope composition (δ²H, δ¹⁸O) across the two agricultural plots in Marrakech (Morocco) with surface drip and subsurface drip irrigation treatments for a complete hydrologic year (June 2022 to June 2023). Weekly to daily samples of rainfall, irrigation water, groundwater, and soil at various depths (5–50 cm) were sampled, and water from branch xylem was extracted using the cryogenic vacuum distillation method. We found that the subsurface irrigation treatment, which delivered water directly to the root zone, maintained narrow isotopic ranges in water of soils beyond 30 cm, as well as in branch xylem and leaf water. By contrast, surface irrigation treatment plots showed pronounced evaporative isotopic enrichment: summer topsoil water δ¹⁸O peaked at −1.1‰ (vs. −8.7‰ in subsurface irrigation treatment), and leaf water reached +13‰ (vs. +8‰ in subsurface). Despite this larger isotopic heterogeneity in surface irrigation site, branch xylem water δ¹⁸O remained within −6 to 2.5‰ across all soil depth, similar to subsurface irrigation treatment, which ranged between −5 and 0‰. This suggests that olive roots accessed soil water uniformly from the upper 50 cm under both irrigation treatments. Seasonal xylem isotopic enrichment in spring and midsummer mirrored shifts towards shallow, evaporatively altered soil water under surface irrigation, but not under the subsurface. The results suggest that subsurface drip irrigation can significantly improve drought resilience and water-use efficiency in the expanding olive sector of the Maghreb, while continuous isotope monitoring serves as a practical approach to enhance sustainable and adaptive water management in water-limited regions.

Evapotranspiration (ET) is challenging to measure directly, motivating the use of remote sensing products as alternatives. We evaluated five high-resolution (≤1 km) global ET products (SSEBop, MOD16, ETMonitor, PMLv2, and FAO’s WaPOR) against five eddy covariance (EC) measurements in Morocco’s semi-arid Tensift Al-Haouz region, with observations spanning from 2006 to 2019. These five products were selected because they offer the finest spatial resolution (around 1 km or less) among freely downloadable global ET datasets, making them well-suited for comparison with local EC flux tower data. The study area was chosen for its reliable ground-truth EC stations, extensive knowledge of local irrigation practices, and a semi-arid climate that provides a rigorous testbed for ET model evaluation in water-limited conditions. Precipitation observations were included to assess each product’s sensitivity to soil moisture and precipitation-driven ET variations, particularly to identify which models respond to rainfall and irrigation inputs (i.e., differences between rainfed and irrigated fields). Results indicate that PMLv2 achieved the best agreement with EC (R² up to 0.65, RMSE as low as 0.4 mm/day, and PBIAS under 10% at most sites), followed by WaPOR and SSEBop which captured seasonal ET patterns (R² ~0.3–0.5) with moderate bias (~20–30%). In contrast, ETMonitor and MOD16 underperformed, showing larger errors (RMSE ~1–2.5 mm/day) and substantial underestimation biases (e.g., MOD16 PBIAS ~50–80% in irrigated sites). These findings underscore the impact of algorithmic differences and highlight PMLv2, SSEBop, and WaPOR as more reliable options for estimating ET in semi-arid agricultural regions lacking in situ measurements.