Assessment of Agricultural Water Requirements for Semi-Arid Areas: A Case Study of the Boufakrane River Watershed (Morocco)

: This work was undertaken to develop a low-cost but reliable assessment method for agricultural water requirements in semi-arid locations based on remote sensing data/techniques. In semi-arid locations, water resources are often limited, and long-term water consumption may exceed the natural replenishment rates of groundwater reservoirs. Sustainable land management in these locations must include tools that facilitate assessment of the impact of potential future land use changes. Agricultural practices in the Boufakrane River watershed (Morocco) were used as a case study application. Land use practices were mapped at the thematic resolution of individual crops, using a total of 13 images generated from the Sentinel-2 satellites. Using a supervised classiﬁcation scheme, crop types were identiﬁed as cereals, other crops followed by cereals, vegetables, olive trees, and fruit trees. Two classiﬁers were used, namely Support vector machine (SVM) and Random forest (RF). A validation of the classiﬁed parcels showed a high overall accuracy of 89.76% for SVM and 84.03% for RF. Results showed that cereal is the most represented species, covering 8870.43 ha and representing 52.42% of the total area, followed by olive trees with 4323.18 ha and a coverage rate of 25%. Vegetables and other crops followed by cereals cover 1530.06 ha and 1661.45 ha, respectively, representing 9.4% and 9.8% of the total area. In the last rank, fruit trees occupy only 3.67% of the total area, with 621.06 ha. The Food and Agriculture Organization (FAO) free software was used to overlay satellite data images with those of climate for agricultural water resources management in the region. This process facilitated estimations of irrigation water requirements for all crop types, taking into account total potential evapotranspiration, effective rainfall, and irrigation water requirements. Results showed that olive trees, fruit trees, and other crops followed by cereals are the most water demanding, with irrigation requirements exceeding 500 mm. The irrigation requirements of cereals and vegetables are lower than those of other classes, with amounts of 300 mm and 150 mm, respectively. irrigation requirements, other crops followed by cereals, olive trees, and fruit trees were the three types of crops that required a very large quantity of water, exceeding 500 mm. Vegetables required about 450 mm; the demand of cereals did not exceed 200 mm. The potential evapotranspiration was strongly correlated with the water demand; the crop types that demanded a lot of water were those that recorded high values of evapotranspiration. For the ER, the class that recorded the lowest values was vegetables, the other classes recorded values higher than 300


Introduction
According to World Bank data, the global rural population declined from 66.39% of the total population in 1960 to 44.72% in 2016. Meanwhile, agricultural added value, as a percentage of GDP (gross domestic product), decreased from 7.59 to 3.43% during the period 1994-2017. Despite these changes, the continued growth of the global human population has resulted in phases of deforestation together with innovations, increasing the efficiency of agricultural systems, and subsequent crises of unprecedented demographic, economic and urban expansion [1,2].
In semi-arid regions like Morocco, agricultural practices are facing a series of challenges not limited to climate change, which is reflected in an increasingly warmer and drier climate [3], coupled with the increasingly random spatiotemporal variability of rainfall, and associated droughts and floods [4][5][6][7][8]. These issues are exacerbated by increasingly complex land use and land cover practices that in turn adversely affect socio-economic development [9]. Agriculture, as an important economic sector, is therefore deeply impacted, given crop production dependence on the annual rainfall distribution [4][5][6][7][8]. It is clear that humans have created and are now witnessing a great agricultural ecosystem disturbance [8].
Recently, geospatial technologies have been used extensively for spatiotemporal monitoring of environmental phenomena, including land use/land cover changes [6,[9][10][11][12][13][14][15], understanding the ecosystem functions [16,17], identifying agricultural systems and crop mapping [3,8,18,19], estimating fractional crop cover and crop residue [20], estimating the impacts of urbanization on agricultural dynamics [3], identifying the karst cavities in agricultural areas [21][22][23], and water balance assessments at regional and local scales [6,24]. Many investigations have shown a great deal of potential in terms of different machine learning approaches in imagery classification, such as vector machine support [3,[25][26][27][28][29] and random forest [30][31][32]. These innovative data interrogation and modeling approaches are critically important for estimating agricultural crop water use. This is important in order to implement effective strategies for advanced water resource management for agriculture in response to contemporary water budget challenges. Several works have been published supporting these needs in recent years [33][34][35][36][37][38][39][40][41], particularly in semi-arid zones [42]. Ofentse Moseki et al. [42] used the CROPWAT model to determine the irrigation needs of the Jatropha crop in Botswana. They used the CROPWAT model to estimate baseline evapotranspiration (ETo), evapotranspiration (ETc), irrigation water requirements (IWR) and yield response to irrigation scheduling in Botswana. The results showed that the annual ETo from 2014 to 2016 at the station was 1456 mm. The lowest monthly ETo (50.10 mm) was observed in June and the highest (182.59 mm) in January.
This model is widely used, especially in understanding the changes in crop water requirements [38], which are defined as the depth of water required to meet the evapotranspiration water loss (ETc) of a disease-free crop growing in large fields; this parameter is important for promoting sustainable development. In particular, the model is used in the determination of crop water and of the effects of irrigation programming on the crop [40,41]. This model allows calculation of the water requirements of the different crops using soil data, climatic data, and data on the crops themselves. Therefore, to determine the crop's water requirement, several parameters were calculated.
Calculation of Potential Evaporation of Crop ETc: Before calculating the ETc, specific studies on the water requirements of crops in the area should be examined; the meteorological and research stations and the environment should also be visited. The calculation of this parameter is done by the following two main steps: Irrigation requirements: Part of the crop water requirements is met by rainfall (Pe), groundwater (Ge) and stored soil water (Wb); or and is determined on a monthly basis.
The specific objectives of this work were to (i) use high spatial resolution Sentinel-2 images to map crop types in the Boufakrane watershed; (ii) evaluate machine learning (ML) methods such as support vector machine (SVM) classifier and random forest (RF) in crop species' mapping; (iii) use the CROPWAT 8.0 model to estimate water demand for agriculture in the study area through the calculation of potential evaporation and effective rainfall. Finally, irrigation water requirements were estimated.

Study Area
The Boufakrane River watershed is located in the headwaters region of the Great Sebou Basin between longitudes 5 • 25 46.13 and 5 • 37 49.71 W and between latitudes 33 • 28 54.40 and 33 • 58 32.93 N ( Figure 1). Locally, the area is part of the Fez-Meknes region, which is one of the most important and productive areas for agriculture in the region, given its relatively high water availability and good quality soils. The regional Useful Agricultural Area (UAA) is estimated to be approximately 1,340,826 hectares, representing 15% of the national total area. The UAA is dominated by cereals (816,000 ha) and olive trees (350,000 ha), and 14% of the area is irrigated (184,162 ha). Climatically, the region is characterized by a semi-arid climate, with a mean annual rainfall of 500 mm, a mean annual reference evapotranspiration of 907 mm, and a dry season extending from June to October ( Figure 2).

Data
The regions of interest used for land use classification were developed based on a series of field missions throughout the study area along with high-resolution Google Earth visualizations. A total of 88,546 pixels with a resolution of 10 m 2 was used for classification; 65% were used as reference data and 35% for the validation. The sampling was done in a random way; it was chosen to be representative and well distributed in the area. Figure 3 shows the field data used for this work. Weather data were collected for the study period of September 2018 to August 2019. These data included monthly mean rainfall,

Data
The regions of interest used for land use classification were developed based on a series of field missions throughout the study area along with high-resolution Google Earth visualizations. A total of 88,546 pixels with a resolution of 10 m 2 was used for classification; 65% were used as reference data and 35% for the validation. The sampling was done in a random way; it was chosen to be representative and well distributed in the area. Figure 3 shows the field data used for this work. Weather data were collected for the study period of September 2018 to August 2019. These data included monthly mean rainfall, monthly average minimum temperature, monthly average maximum temperature, humidity (%), wind speed (m/s), and sunshine intensity (hours). Figure 4 shows climate data collected by the climate station at the Faculty of Science of Meknes, Moulay Ismail University, coordinates: Latitude: 33 • 52 11.12 N, Longitude: 5 • 32 35.11 W, Z = 554 m.
For the satellite data, a total of 13 satellite images covering a whole crop year were used to carry out this work. All these images were obtained from the Sentinel-2 sensor of the European Space Agency (https://sentinel.esa.int/web/sentinel/sentinel-data-access (accessed on 15 August 2021)). This mission was launched in June 2015, with a revisit time (i.e., image interval) of 10 days and image spatial resolution of 10 m to 60 m in thirteen spectral bands from visible to mid-infrared. Images were downloaded from https://scihub.copernicus.eu/dhus/#/home (accessed on 15 August 2021) for the period August 2018-August 2019 (Table 1). the European Space Agency (https://sentinel.esa.int/web/sentinel/sentinel-data-access (accessed on 15 August 2021)). This mission was launched in June 2015, with a revisit time (i.e., image interval) of 10 days and image spatial resolution of 10 m to 60 m in thirteen spectral bands from visible to mid-infrared. Images were downloaded from https://scihub.copernicus.eu/dhus/#/home (accessed on 15 August 2021) for the period August 2018-August 2019 (Table 1).

Methodolgy
Normalized difference vegetation index (NDVI) was calculated to construct NDVI time-series images. In parallel, several surveys and field missions were conducted in the region to collect reference (validation) points for the different agricultural crops. These data were combined with the high spatial resolution Google Earth images. Spectral profiles were then constructed and used as input data for a machine learning approach to map different crop species in the region. As a final step, and to associate satellite data with observed water management practices, CROPWAT 8.0 software was used to estimate crop water requirements ( Figure 5). The simulation was carried out on three parameters: potential evapotranspiration for the different crop types (ETc), effective rainfall (ER), and irrigation water requirements (IWR). These three parameters are dependent on each other. The ETc is the amount of water that should be transpired in a given time by the crop, while the ER is defined as the rainfall fraction that responds to the crops' water needs [43]. During the rainy months, rainfall covers the water requirements of the crops, while during the dry months, rainfall must be supplemented by irrigation water to cover water requirements.

Support Vector Machine (SVM)
SVMs belong to a family of algorithms that use supervised learning and are specialized in solving mathematical discrimination and regression problems. They were developed in 1998 by Vladimir Vapnik [44]. Support vector machines (SVMs) represent a group of theoretically superior machine learning algorithms. The development of this method was initially triggered by the exploration and formalization of machine learning capacity control and over-fitting problems [44] and represents an efficient technique, with reduced data and processing demands. The method avoids the problems of over-adjustment and does not require any hypothesis on the type of data. Although non-parametric, the method is capable of developing efficient decision limits and can therefore minimize classification errors. This is done by searching for the optimal separation between classes [45]. Their work was quickly adopted because of their ability to work with large data, their theoretical guarantees and the good results achieved in practice. Requiring a small number of parameters, SVMs are appreciated for their simplicity of use.

Random Forest (RF) Classifier
Random forest (RF) was developed by [46]. It is a supervised non-parametric method applicable for both classification and prediction [47,48]. Model subroutines are composed of a combination of decision trees used independently to assign the most frequent class to the input data, and the majority vote of the trees determines the class prediction. The part of the data not used in tree training is used for performance evaluation.
For the current investigation, after extracting the crop type characteristic based on the data collected in the field, twelve decision trees were constructed and were the basis for RF classifier ( Figure 6). These decision trees make it possible to predict the different classes.

Support Vector Machine (SVM)
SVMs belong to a family of algorithms that use supervised learning and are specialized in solving mathematical discrimination and regression problems. They were developed in 1998 by Vladimir Vapnik [44]. Support vector machines (SVMs) represent a group of theoretically superior machine learning algorithms. The development of this method was initially triggered by the exploration and formalization of machine learning capacity control and over-fitting problems [44] and represents an efficient technique, with reduced data and processing demands. The method avoids the problems of over-adjustment and does not require any hypothesis on the type of data. Although non-parametric, the method is capable of developing efficient decision limits and can therefore minimize classification errors. This is done by searching for the optimal separation between classes [45]. Their work was quickly adopted because of their ability to work with large data, their theoretical guarantees and the good results achieved in practice. Requiring a small number of parameters, SVMs are appreciated for their simplicity of use.

Random Forest (RF) Classifier
Random forest (RF) was developed by [46]. It is a supervised non-parametric method applicable for both classification and prediction [47,48]. Model subroutines are composed of a combination of decision trees used independently to assign the most frequent class to the input data, and the majority vote of the trees determines the class prediction. The part of the data not used in tree training is used for performance evaluation.
For the current investigation, after extracting the crop type characteristic based on the data collected in the field, twelve decision trees were constructed and were the basis for RF classifier ( Figure 6). These decision trees make it possible to predict the different classes. Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 18

NDVI Time-Series Spectral Profile Curves
Field data and data collected from the Regional Directorate of Agriculture showed the presence of five main cropping systems in the region, including cereals, other crops followed by cereals, vegetables (onion, potatoes, tomatoes), olive trees and fruit trees. More than two hundred profiles were developed for the different crop types. These profiles were associated with the field data and the visualization of high spatial resolution Google Earth images in order to collect input data for classification (Figures 7 and 8).

NDVI Time-Series Spectral Profile Curves
Field data and data collected from the Regional Directorate of Agriculture showed the presence of five main cropping systems in the region, including cereals, other crops followed by cereals, vegetables (onion, potatoes, tomatoes), olive trees and fruit trees. More than two hundred profiles were developed for the different crop types. These profiles were associated with the field data and the visualization of high spatial resolution Google Earth images in order to collect input data for classification (Figures 7 and 8).
From Figure 8, it is possible to discriminate the spectral characteristics of the different crops in relation to NDVI values during the year. This index, proposed for the first time by Rouse et al., 1973 [49], is widely used and provides information on the quantity and vigor of vegetation, taking into account the near infrared (NIR) and visible red bands of the electromagnetic spectrum [49,50] calculated by the following equation: where ρ N IR : the reflectance in the near − infrared reflectance, ρ RED : the reflectance in the red band.
For example, for an olive pixel, the NDVI value did not change significantly throughout the year, with an increase around February. While for a pixel of cereals, the NDVI values did increase with the crop growth cycle. Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 18 From Figure 8, it is possible to discriminate the spectral characteristics of the different crops in relation to NDVI values during the year. This index, proposed for the first time by Rouse et al. 1973 [49], is widely used and provides information on the quantity and vigor of vegetation, taking into account the near infrared (NIR) and visible red bands of the electromagnetic spectrum [49,50] calculated by the following equation: where : the reflectance in the near − infrared reflectance, : the reflectance in the red band.
For example, for an olive pixel, the NDVI value did not change significantly throughout the year, with an increase around February. While for a pixel of cereals, the NDVI values did increase with the crop growth cycle.

Overall Accuracy
The classification approaches used in this work were selected based on the confusion matrix and the Kappa index [51][52][53], whose the overall accuracy is the proportion of the area mapped correctly. It provides the user of the map with the probability that a randomly selected location on the map is correctly classified. The Kappa coefficient measures the agreement between the resulting classes of the classifier and the true values [52,54], with values ranging from 0 to 1, where 0 represents no agreement and 1 represents perfect agreement. The table 2 shows the confusion matrix calculated from the reference data and the map classes; the reference data is presented in the row, and the map classes in the columns. The results showed that the two approaches showed high classification accuracies. The overall accuracy for the SVM classifier exceeded 89.76%, and a significant agreement by the Kappa index of 0.79 was obtained. The overall accuracy for RF was 84.03%, with a Kappa index of 0.68. For the validation, 57,554 pixels were used. In most cases, it

Overall Accuracy
The classification approaches used in this work were selected based on the confusion matrix and the Kappa index [51][52][53], whose the overall accuracy is the proportion of the area mapped correctly. It provides the user of the map with the probability that a randomly selected location on the map is correctly classified. The Kappa coefficient measures the agreement between the resulting classes of the classifier and the true values [52,54], with values ranging from 0 to 1, where 0 represents no agreement and 1 represents perfect agreement. The Table 2 shows the confusion matrix calculated from the reference data and the map classes; the reference data is presented in the row, and the map classes in the columns. The results showed that the two approaches showed high classification accuracies. The overall accuracy for the SVM classifier exceeded 89.76%, and a significant agreement by the Kappa index of 0.79 was obtained. The overall accuracy for RF was 84.03%, with a Kappa index of 0.68. For the validation, 57,554 pixels were used. In most cases, it showed a high accuracy of this classification for most crop species. Few confusions between classes were recorded using the RF approach (e.g., crops followed by cereals, cereals, vegetables, and olive trees). In addition to the overall accuracy and the Kappa index, other types of accuracy and errors were calculated for both classifiers, including the producer's accuracy, the user's accuracy, the commission error, and the omission error (Table 3) [53].
(i) Producer's accuracy is defined as the probability that a value in a reference dataset was correctly classified. Producer's accuracy is the complement to the probability of omission error. (ii) User's accuracy represents the probability that a resulting value in a certain class is really that class. User's accuracy is the complement to the probability of commission error. (iii) Commission errors represent the fraction of the resulting values in a class that does not belong to that class. (iv) Omission errors represent the fraction of values that belongs to one class but was predicted in a different class. Table 3. Commission, omission, producer's and user's accuracy for the SVM and RF classifiers. In terms of accuracy, the producer's accuracy and user's accuracy confirm the results found for the overall accuracy. The producer's accuracy showed very high values, with more than 80% for all classes in the both classifiers SVM and RF, except the class other crops followed by cereals, which presented a value of 41.86% for the RF classifier. The user's accuracy showed very high values, with more than 79% for all classes in the SVM classifier, except for other crops followed by cereals, which presented a value 17.07%. For the RF classifier, this accuracy shows high values for the three classes of cereals, vegetables, and olive trees: 98.09%, 80.62%, and 66.05%, respectively. The other two classes, fruit trees and other crops followed by cereals, represent low precision, with values of 46.71% and 7.74%, respectively.

Classes Commission (%) Omission (%) Producer's Accuracy (%) User's Accuracy (%)
In terms of commission error, the results showed that the other crops followed by cereals represents the highest values for the two classifiers SVM and RF, with values of 82.93% and 92.26%, respectively, followed by vegetables and fruit trees. The other two classes represented a low value for this error. For the omission error, cereals and other crops followed by cereals represented the highest value for the SVM, whereas for RF, cereals and other crops followed by cereals represent a high value of 58.14%, followed by fruit trees with 30.54% and olive trees with 29.95%.

Crop Mapping
Crop mapping was performed based on field data combined with a detailed study of the chlorophyll response (NDVI) for each of crop type. The crops determined in the region were cereals, other crops followed by cereals, vegetables (onions, potatoes, tomatoes), olive trees, and fruit trees. In order to calculate the areas of each class, pixel size was used. Therefore, after obtaining the classes, the area of each class was obtained by multiplying the number of pixels and the pixel size

CROPWAT for Water Crop Requirements
This section includes calculation of the crop water requirements using the FAO free software CROPWAT 8.0, based on climate, soil, and crop data. Thus, three main variables were estimated by the units of water depth (mm): ETc, ER, and IWR. For preliminary planning, monthly data are frequently used, and the total of the data of the different crops over the area constitutes the basis for determining the supply.
The climatic data were used to calculate the reference evapotranspiration ET0, and by determining the timing of planting or sowing, the rate of crop development, the duration of crop development stages, and the growing season kc for a given crop were chosen. Then, the ETc for each crop type was calculated for each 10 day period. Figure 10 shows the different crop types' evapotranspiration from October 2018 to August 2019. The curves show that this parameter increases in the driest months (July and August) for all crop types, while it decreases in the rainiest months (December to February). For the class of other crops followed by cereals, for example, the potential evapotranspiration reached up

CROPWAT for Water Crop Requirements
This section includes calculation of the crop water requirements using the FAO free software CROPWAT 8.0, based on climate, soil, and crop data. Thus, three main variables were estimated by the units of water depth (mm): ETc, ER, and IWR. For preliminary planning, monthly data are frequently used, and the total of the data of the different crops over the area constitutes the basis for determining the supply.
The climatic data were used to calculate the reference evapotranspiration ET 0 , and by determining the timing of planting or sowing, the rate of crop development, the duration of crop development stages, and the growing season kc for a given crop were chosen. Then, the ETc for each crop type was calculated for each 10 day period. Figure 10 shows the different crop types' evapotranspiration from October 2018 to August 2019. The curves show that this parameter increases in the driest months (July and August) for all crop types, while it decreases in the rainiest months (December to February). For the class of other crops followed by cereals, for example, the potential evapotranspiration reached up to 70 mm. Not all precipitation is effective, and in the most cases some of this precipitation can be lost through surface runoff, deep percolation, or evaporation. Only part of the highintensity rain can penetrate and be stored in the root zone. These rains can be 100% effective when the vegetation cover is complete, while they can be only 60% effective with a low percentage of vegetation cover. The relationship between the average monthly ER and the average monthly rainfall is shown for different values of the average monthly ETc [55]. Figure 11 shows the evolution of ER for the different types of crops. The evolution curves of this parameter show that the crops' needs were met in the rainy months. However, in the dry months, these crops suffered from water stress. Irrigation water requirements are calculated using the field water balance, based on ETc and ER. They allow for optimal production in a given growing environment. Figure 12 shows the evolution of water requirements for agriculture for the different types of crops. In the rainy months, rainfall covers the water requirements of the crops; Not all precipitation is effective, and in the most cases some of this precipitation can be lost through surface runoff, deep percolation, or evaporation. Only part of the high-intensity rain can penetrate and be stored in the root zone. These rains can be 100% effective when the vegetation cover is complete, while they can be only 60% effective with a low percentage of vegetation cover. The relationship between the average monthly ER and the average monthly rainfall is shown for different values of the average monthly ETc [55]. Figure 11 shows the evolution of ER for the different types of crops. The evolution curves of this parameter show that the crops' needs were met in the rainy months. However, in the dry months, these crops suffered from water stress. Not all precipitation is effective, and in the most cases some of this precipitation can be lost through surface runoff, deep percolation, or evaporation. Only part of the highintensity rain can penetrate and be stored in the root zone. These rains can be 100% effective when the vegetation cover is complete, while they can be only 60% effective with a low percentage of vegetation cover. The relationship between the average monthly ER and the average monthly rainfall is shown for different values of the average monthly ETc [55]. Figure 11 shows the evolution of ER for the different types of crops. The evolution curves of this parameter show that the crops' needs were met in the rainy months. However, in the dry months, these crops suffered from water stress. Irrigation water requirements are calculated using the field water balance, based on ETc and ER. They allow for optimal production in a given growing environment. Figure 12 shows the evolution of water requirements for agriculture for the different types of crops. In the rainy months, rainfall covers the water requirements of the crops; Figure 11. Effective rain.
Irrigation water requirements are calculated using the field water balance, based on ETc and ER. They allow for optimal production in a given growing environment. Figure 12 shows the evolution of water requirements for agriculture for the different types of crops. In the rainy months, rainfall covers the water requirements of the crops; during this time, the water requirements for irrigation are expected to be very low. This parameter is strongly related to the climatic conditions and is directly influenced by variations in conditions. It is inversely correlated with the daily rainfall. Thus, it is high in the dry months and low in the rainy months.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 18 during this time, the water requirements for irrigation are expected to be very low. This parameter is strongly related to the climatic conditions and is directly influenced by variations in conditions. It is inversely correlated with the daily rainfall. Thus, it is high in the dry months and low in the rainy months.  Figure 13 shows the total for the year of the three estimated parameters (ETc, ER, and IWR) for the different types of crops. For water irrigation requirements, other crops followed by cereals, olive trees, and fruit trees were the three types of crops that required a very large quantity of water, exceeding 500 mm. Vegetables required about 450 mm; the demand of cereals did not exceed 200 mm. The potential evapotranspiration was strongly correlated with the water demand; the crop types that demanded a lot of water were those that recorded high values of evapotranspiration. For the ER, the class that recorded the lowest values was vegetables, while the other classes recorded values higher than 300 mm.

Discussion
Using the classification approach described earlier, and based on spectral analysis results, a map of the different agricultural species was produced for the study area ( Figure  9). Results generally showed that the classification using both classifiers was satisfactory,  Figure 13 shows the total for the year of the three estimated parameters (ETc, ER, and IWR) for the different types of crops. For water irrigation requirements, other crops followed by cereals, olive trees, and fruit trees were the three types of crops that required a very large quantity of water, exceeding 500 mm. Vegetables required about 450 mm; the demand of cereals did not exceed 200 mm. The potential evapotranspiration was strongly correlated with the water demand; the crop types that demanded a lot of water were those that recorded high values of evapotranspiration. For the ER, the class that recorded the lowest values was vegetables, while the other classes recorded values higher than 300 mm. during this time, the water requirements for irrigation are expected to be very low. This parameter is strongly related to the climatic conditions and is directly influenced by variations in conditions. It is inversely correlated with the daily rainfall. Thus, it is high in the dry months and low in the rainy months. Figure 12. Irrigation requirements. Figure 13 shows the total for the year of the three estimated parameters (ETc, ER, and IWR) for the different types of crops. For water irrigation requirements, other crops followed by cereals, olive trees, and fruit trees were the three types of crops that required a very large quantity of water, exceeding 500 mm. Vegetables required about 450 mm; the demand of cereals did not exceed 200 mm. The potential evapotranspiration was strongly correlated with the water demand; the crop types that demanded a lot of water were those that recorded high values of evapotranspiration. For the ER, the class that recorded the lowest values was vegetables, while the other classes recorded values higher than 300 mm.

Discussion
Using the classification approach described earlier, and based on spectral analysis results, a map of the different agricultural species was produced for the study area ( Figure  9). Results generally showed that the classification using both classifiers was satisfactory,

Discussion
Using the classification approach described earlier, and based on spectral analysis results, a map of the different agricultural species was produced for the study area ( Figure 9). Results generally showed that the classification using both classifiers was satisfactory, with the exception of some confusion between a few classes, which is likely due to the spectral similarity of the crops.
Previous investigations using geospatial techniques have shown a strong efficiency in land use/land cover change monitoring [6,11], crop mapping [3,18], and identification of agricultural systems [8]. Ouzemou et al. [3] carried out an insightful study in the plains of Tadla, Morocco. The objectives included mapping different agricultural species using high-resolution satellite data and machine learning approaches and comparing the different used approaches. Their study showed an overall accuracy of 89.26%, 85.27%, and 57.17%, respectively, for random forest, support vector machine, and spectral angle mapper, with a Kappa index of 0.85, 0.80, and 0.4, respectively. Comparing our result with this one, our study also showed a high overall accuracy of 89.76% for SVM and 84.03% for RF, and a Kappa index of 0.79 and 0.68, respectively.
The CROPWAT model was used to determine water demand for the different agricultural species in the region. Three parameters were calculated, namely, the crops' potential evaporation, the effective rain, and the irrigation requirements. The results obtained by [40] showed that the irrigation requirements varied according to the location, whereas the required water quantity per palm varied between 115 and 200 liters per day. Comparing our result with this one, our study showed that olive trees, fruit trees, and other crops followed by cereals are the most water demanding, with needs exceeding 500 mm. The water demands of cereals and vegetables are lower than that of other classes, with amounts of 300 mm and 150 mm, respectively.
As explained, water consumption is increasing in the Saiss plain. This is mainly due to excessive exploitation. According to the 1939-2002 groundwater data record, there is a constant deficit of approximately 100 Mm 3 /year, with an estimated inflow of 242 Mm 3 /year and an estimated outflow of 342 Mm 3 /year. While the output includes abstraction with 260 Mm 3 /year, and rivers and springs with 82 Mm 3 /year, 22% of the water balance is dedicated to human drinking water supplies and 78% to private irrigation [56].

Conclusions
The assessment and estimation of water demand for agriculture is crucial to improve water resource management in a given region. The final objective was to determine the water demand for agriculture in the Boufakrane River watershed through several steps. First, a map of the different crop types was produced using the SVM and RF machine learning algorithms, based on field data combined with the high spatial resolution Google Earth images. Five crop types were mapped, including cereals, other crops followed by cereals, vegetables, olive trees, and fruit trees. Then, the evaluation of the classification map was made based on the Kappa index and the overall accuracy. Finally, the satellite data were combined with climate, soil, and crop data before being used as inputs for CROPWAT software to estimate the water requirements for agriculture.
The mapping results showed a strong potential of high-resolution satellite data in agricultural species mapping. The evaluation of the two classifiers used (RF and SVM) showed a Kappa index higher than 0.67 and an overall accuracy exceeding 83%. The irrigation requirements showed that the other crops followed by cereals, olive trees, and fruit trees were the three types of crops that required a very large quantity of water, exceeding 500 mm. Vegetables required an amount of about 450 mm; the demand of cereals did not exceed 200 mm.
The method developed in this work facilitates estimations of the agriculture water demand in the study area, thereby promoting sustainable water resource management. Through this study, we recommend a combination of these methods with existing real data for the implementation of a system for the quantification of water resources for crops throughout Morocco, which would allow validation using global crop yield data. It was also an opportunity to see the link between water demand and known groundwater reserves and existing data on actual evapotranspiration in this area.