High quality and long-term precipitation data are required to study the variability and trends of rainfall and the impact of climate change. In developing countries like Morocco, the quality of climate data collected from various weather stations faces numerous obstacles. This paper presents methods for collecting, correcting, reconstructing, and homogenizing precipitation series of Morocco’s Fez-Meknes region from 1961 to 2019. Data collected from national specialized agencies based on 83 rain gauge stations was processed through an algorithm specially designed for the homogenization of climatic data (Climatol). We applied the Mann-Kendall test and Sen’s slope estimator to raw and homogenized data to calculate rainfall trend magnitudes and significance. The homogenization process allows for the detection of a larger number of stations with statistically significant negative trends with 95% and 90% confidence levels, particularly in the mountain ranges, that threatens the main sources of water in the largest watershed in the country. The regionalization of our rain gauge stations is highlighted and compared to previous studies. The monthly and annual means of raw and homogenized data show minor differences over the three main climate zones of the region.

This paper explores the temporal and spatial patterns of annual, seasonal, and monthly rainfall series during the period of 1961–2018 at 15 stations in the agriculturally important Sebou river basin, northern Morocco. Trends were investigated using the classical non-parametric Mann–Kendall test and the Theil–Sen approach at 90%, 95% and 99% confidence levels. A general decreasing trend was found at the annual scale, significant at the 95% confidence level at 8 stations out of 15 (53%). A particularly large decreasing trend between −30 mm and −50 mm per decade was found in the north and eastern parts of the basin. Autumn rainfall tended to increase, but this was not statistically significant. During the winter months, rainfall tended to decrease sharply (−27 mm and −40 mm per decade) in the northern slopes of the Rif mountains, while in spring, the mountainous area of the basin recorded decreases ranging between −12 mm and −16 mm per decade. During winter and spring, negative trends were significant at ten stations (66%). Summer rainfall tends toward a decrease, but the absolute change is small. These results help to understand the rainfall variability in the Sebou river basin and allow for improved mitigation strategies and water resource plans based on a prospective view of the impact of climate change on the river basin.

Exploring the relationship between cereal yield and the remotely sensed normalized difference vegetation index (NDVI) is of great importance to decision-makers and agricultural stakeholders. In this study, an approach based on the Pearson correlation coefficient and linear regression is carried out to reveal the relationship between cereal yield and Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI data in the Fez-Meknes region of Morocco. The results obtained show strong correlations reaching 0.70 to 0.89 between the NDVI and grain yield. The linear regression model explains 58 to 79% of the variability in yield in regional provinces marked by the importance of cereal cultivation, and 51 to 53% in the mountainous provinces with less agricultural land devoted to major cereals. The regression slopes indicate that a 0.1 increase in the NDVI results in an expected increase in grain yield of 4.9 to 8.7 quintals (q) per ha, with an average of 6.8 q/ha throughout the Fez-Meknes region. The RMSE ranges from 2.12 to 4.96 q/ha. These results are promising in terms of early yield forecasting based on MODIS-NDVI data, and consequently, in terms of grain import planning, especially since the national grain production does not cover the demand. Such remote sensing data are therefore essential for administrations that are in charge of food security decisions.

Morocco’s meteorological observation network is quite old, but the spatial coverage is insufficient to conduct studies over large areas, especially in mountainous regions, such as the Fez-Meknes region, where spatio-temporal variability in precipitation depends on altitude and exposure. The lack of station data is the main reason that led us to look for alternative solutions. TerraClimate (TC) reanalysis data were used to remedy this situation. However, reanalysis data are usually affected by a bias in the raw values. Bias correction methods generally involve a procedure in which a “transfer function” between the simulated and corrected variable is derived from the cumulative distribution functions (CDFs) of these variables. We explore the possibilities of using TC precipitation data for the Fez-Meknes administrative region (Morocco). This examination is of great interest for the region whose mountain peaks constitute the most important reservoir of water in the country, where TC data can overcome the difficulty of estimating precipitation in mountainous regions where the spatio-temporal variability is very high. Thus, we carried out the validation of TC data on stations belonging to plain and mountain topographic units and having different bioclimatic and topographic characteristics. Overall, the results demonstrate that the TC data capture the altitudinal gradient of precipitation and the average rainfall pattern, with a maximum in November and a minimum in July, which is a characteristic of the Mediterranean climate. However, we identified quasi-systematic biases, negative in mountainous regions and positive in lowland stations. In addition, summer precipitation is overestimated in mountain regions. It is considered that this bias comes from the imperfect representation of the physical processes of rainfall formation by the models. To reduce this bias, we applied the quantile mapping (QM) method. After correction using five QM variants, a significant improvement was observed for all stations and most months, except for May. Validation statistics for the five bias correction variants do not indicate the superiority of any particular method in terms of robustness. Indeed, results indicate that most QM methods lead to a significant improvement in TC data after monthly bias corrections.

The aim of this paper was to present a precipitation trend analysis using gridded data at annual, seasonal and monthly time scales over the Fez-Meknes region (northern Morocco) for the period 1961–2019. Our results showed a general decreasing trend at an annual scale, especially over the mountain and the wetter parts of the region, which was statistically significant in 72% of the grid points, ranging down to −30 mm per decade. A general upward trend during autumn, but still non-significant in 95% of the grid points, was detected, while during winter, significant negative trends were observed in the southwest (−10 to −20 mm per decade) and northeast areas (more than −20 mm per decade) of the region. Spring rainfall significantly decreased in 86% of the grid points, with values of this trend ranging between 0 and −5 mm per decade in the upper Moulouya and −5 to −10 mm per decade over the rest of the region (except the northwest). At a monthly time scale, significant negative trends were recorded during December, February, March and April, primarily over the northeast Middle Atlas and the northwest tip of the region, while a significant upward trend was observed during the month of August, especially in the Middle Atlas. These results could help decision makers understand rainfall variability within the region and work out proper plans while taking into account the effects of climate change.