1. Introduction
Olive trees are a drought-tolerant [
1] species of high economic importance, traditionally grown in the Mediterranean Region, where 95 percent of global olive cultivation is located [
2]. Over the last three decades, super-high density (SHD) planting systems have been taking over the olive sector worldwide, since their smaller varietal sizes and tree-to-tree distances allow for the complete mechanization of pruning and harvest processes, lowering production costs in the long term and increasing yields [
3]. However, olive trees also suffer under increasing droughts and higher temperatures [
4], and the recommended climate change adaptation would consist of increasing spacing between trees [
5]. Besides contrarily reducing tree-to-tree distance, an increasing number of plantations are also adopting (often backup) irrigation measures that are strongly encouraged by governmental subsidies on drip irrigation [
6]. Hereby, the limited water resources of many drought-prone areas dedicated to olive cultivation risk being more quickly exploited and in extreme cases depleted, increasing the environmental risk of desertification but also affecting the local population dependent on such water resources [
7,
8,
9].
Over the last ten years, Moroccan olive production has more than doubled, placing the country among the five main producers worldwide [
2]. This development has been strongly encouraged by a national strategy called the Green Morocco Plan (GMP) [
10], in which important subsidies were allocated to efficient irrigation systems [
11] and in the conversion to arboriculture [
12]. Over 30 percent of olive cultivation occurs in the Fès and Meknès Region (FMR) [
13]. Here, the Saïss plain, a very fertile area of economic importance, is one of the regions with the most intensive use of groundwater in Morocco [
8]. It is a source that will be depleted within the coming 25 years [
14], as a 100 Mm³ deficit is registered yearly [
15]. A spatio-temporal assessment of SHD olive plantations can thus help evaluate the impact of olive orchard intensification on local water resources. For this, using remote sensing (RS) data and methods can be a cost and time-efficient approach for large-scale mapping.
Among RS studies on tree crop mapping, the use of different spatial resolutions of satellite imagery has had varying levels of success depending on the tree spacing and tree crown size that were mapped [
16,
17]. Sparsely planted orchards are mostly mapped using object-based image analysis (OBIA) on very high spatial resolution (VHR) imagery [
18,
19]. However, VHR data have high acquisition costs [
20] and high computational power requirements [
21]. For homogeneous and intensive orchards, some studies achieved promising results using 10 m spatial resolution [
17,
22].
The adoption of multi-temporal approaches in land use and land cover (LULC) classification studies allows for capturing crop-specific phenology patterns and helps increase classification accuracies. However, the trade-offs between spatial and temporal resolution concerning commercial high spatial resolution data remain. The high acquisition costs of such imagery often limit the feasibility of their application, especially in less-developed countries [
23].
Besides spatial and temporal resolution, spectral resolution also plays a crucial role in discriminating crop types. Spectral properties of leaves can provide information about their physical and chemical composition. There are many factors that influence a leaf’s optical properties, such as the structure of the epidermis, waxes, and the cutin [
24]. Evergreen leaves, for example, have lower reflectance in the infrared wavelengths than deciduous leaves [
25]. Given these characteristics, different vegetation indices (VIs), which reflect the relationship between combinations of spectral bands, have been developed through the years and allow for determining the greenness of vegetation and their water content, and thus, they have different capacities for tree crop discrimination based on their leaves and planting characteristics.
The most widely used VI in crop classification studies is the Normalized Difference Vegetation Index (NDVI) [
26,
27]. However, it possesses limitations such as soil background brightness [
28] and saturation [
29,
30,
31]. The Modified Soil Adjustment Vegetation Index (MSAVI-2) [
32] allows one to simultaneously increase the vegetation signal while decreasing the soil-induced variations, thus minimizing the influences of soil on vegetation spectra. Specifically, studies on tree crop classification, such as [
33], concluded that NDVI performs worse on the discrimination of tree types than of annual crops. Furthermore, due to the high similarity amongst temporal reflectance profiles of fruit trees, recent work shifted from VI- to full-band-based multi-temporal approaches that consider the entire spectral resolution available in the image data set [
34].
In the case of olive trees, their evergreen phenology presents a discriminatory advantage against deciduous trees and other land use classes with high seasonal variability. The lack of seasonal changes in olive trees’ phenology makes the use of multi-temporal data redundant. As [
16] concluded, multi-temporal datasets did not improve accuracy in olive detection itself, but rather, they improved the classification accuracy of surrounding classes, reducing misclassification errors. We thus argue that, for olive orchard detection, the use of two images per year allows us to eliminate surrounding ground classes with changing phenology. This would further ease the use of higher spatial resolution with lower temporal availability or associated with high acquisition costs.
In LULC studies, besides both higher data availability and accessibility, increasing computational capacity has allowed for the development and adoption of machine learning algorithms (MLAs). MLAs can be subcategorized into supervised methods, which consist of training a classifying algorithm using reference data [
35,
36], and unsupervised methods. While supervised approaches are preferred due to their higher robustness and precision [
37], when no labelled data is available, as in the present study, unsupervised methods are required.
Unsupervised methods recognize unclassified data without any prior knowledge of LULC types, with an interpreter assigning a class to each cluster of pixels via visual interpretation [
38]. They group pixels based on the similarity of their values into clusters or spectral classes [
37]. There is the possibility of applying the clustering algorithm in a hierarchical form, either divisive or agglomerative, starting with a large cluster, which is iteratively subdivided, or many clusters that are iteratively merged [
35,
39].
k-means [
40,
41] is one of the simplest and most well-known clustering algorithms and is often used in RS studies [
42,
43].
In this regard, the overall objective of this study was to develop a two-step classification approach using k-means and high-resolution satellite data to enable a spatio-temporal assessment of olive orchard intensification in the Saïss Plain (Morocco). Specifically, we aimed:
- (i)
To map the extent of SHD olive plantations in the Saïss plain;
- (ii)
To assess the development of this planting pattern from the implementation of the GMP in 2010 until 2020.
Given the limited availability of HR imagery and the lack of reference data over the past years, we developed an unsupervised approach for single-class detection, which uses the advantage of olive trees’ evergreen phenology and the semi-arid climatic conditions in the study area. We hereby aimed to prove that super-intensive olive plantations are a suitable crop for mapping with satellite imagery of low temporal resolution and unsupervised classification techniques. This can be valuable when reference data is scarce and image availability is low or acquisition costs are high.
4. Discussion
4.1. Methodological Achievements and Shortcomings
In our work, super-intensive and intensive olive plantations were successfully mapped in a highly heterogeneous study area in Morocco. For this, we implemented an unsupervised classification approach based on the principles of hierarchical clustering, which allowed us to detect SHD and intensive olive plantations on 5 m resolution imagery. Hereby, two frequent cost-related problems in remote sensing were addressed: (i) while in most land use analyses the cost- and time-intensive gathering of reference data is necessary, e.g., to train supervised classification algorithms, in our approach no labelled data were required; (ii) the design of the clustering and masking sequence and the phenological characteristics of the target crop allowed us to use a very low temporal resolution, requiring only two images per year of analysis. This can be especially advantageous to reducing acquisition costs when using commercial imagery of higher spatial resolution.
Compared to supervised approaches, which can map intra-class variability, our method could only detect dense and vigorous olive plantations, including mainly super-intensive, but also a large share of intensive, and some traditional olive plantations. Younger plantations with smaller tree crown sizes could not be detected or led to unmapped areas within the same plantation. However, this is a problem that other studies using supervised classification for tree crop mapping with high-resolution data of 3 m pixel size, such as [
16], also faced during their analyses. On the other hand, high commission errors were scored on intensive olive plantations due to their similar reflectance values, something that supervised classification may overcome by training the classifying algorithm on different planting intensities as separate land use classes.
Previous work on tree crop mapping, specifically on olive trees, found that using multi-temporal information would not improve olive detection accuracy itself, but rather improve the accuracy of other crop classes, and hereby overall accuracy [
16]. Additionally, the difficulty of mapping perennial crops merely based on their phenological profiles was already discussed in [
17] and [
34], who found that using all reflectance bands helped discriminate different tree crop types. Our study confirmed these findings, since using the low temporal resolution of only two images per year and making use of the NIR band allowed us to reach a remarkably high accuracy in detecting olive plantations. However, limitations of NIR wavelengths were found in the discrimination between some annual crops and olive plantations, which confirms previous findings from the crop separability analysis conducted in
Section 2.4 (
Figure 9).
Based on previous studies on semi-arid environments, we tested different indices (NDVI vs. MSAVI-2) to see whether they helped improve mapping on sites with higher soil background reflectance. While MSAVI-2 successfully reduced soil background reflectance, to map super-intensive olive plantations, which are less affected by this phenomenon, we found that NDVI performed better, and using MSAVI-2 led to higher commission errors among traditional olive orchards.
Our sample-wise verification of 2010 results using VHR Google Earth Imagery revealed that a higher number of smaller plots of annual crops were misclassified in 2010 than in 2020, and that main commission errors in 2020 results were deciduous tree plantations. This may confirm the expected shift from cereal-based agriculture in 2010 towards arboriculture in 2020, incentivized by the GMP. Hereby, k-means cluster centers may also have shifted, creating such a difference among commission errors, both qualitatively and quantitatively. While applying the GDAL sieving tool removed most of the smaller misclassified spots in both years and improved accuracy, further research may investigate whether other unsupervised clustering methods are more appropriate for our proposed approach.
As opposed to most studies on olive mapping, which were mostly performed on homogeneous environments with a low diversity of surrounding ground classes [
19], our study succeeded in extracting olive plantations, mostly of super-intensive and to a great extent also of intensive planting patterns, on a very heterogeneous ground. However, it should be noted that no other evergreen tree crops are known to be cultivated in the study area. Thus, for replicability, this approach needs to be adapted if other evergreen tree crops are also present.
4.2. Implications for Agricultural Management and Policy Monitoring in the Study Region
This study addressed a subject that is currently affecting the olive sector not only in Morocco but also globally and will be of growing concern in times of climate change. Although irrigated olive cultivation can be an important carbon sink and can fight erosion while improving the profitability of the sector, the limited water resources of many drought-prone areas dedicated to olive cultivation may be also exploited and depleted more quickly, presenting not only an environmental risk but also affecting the local population who depend on such water resources.
Mapping super-intensive olive plantations can help locate important actors in the olive oil sector. However, although our study successfully mapped the extent and evolution of an important actor of the groundwater economy in the Saïss plain, we found that super-intensive olive plantations only represent a share of different land uses in the study area and that other actors also need to be addressed to improve water resource management.
Finally, our approach was helpful to analyze policy-induced land use change in the study area. It revealed that the average plantation size increased between 2010 and 2020, reflecting large agribusinesses taking over the sector to the detriment of smaller family farms and medium-size investors. Another remarkable finding was the detection of several larger SHD olive farms that were given up, despite the encouragement of the GMP to adopt intensive forms of arboriculture, especially olive trees. The reasons behind this can be attributed to the short lifecycle of SHD olive plantations. Other contributing factors, such as access to water or better revenues from other crops, should also be considered and could be further investigated.
5. Conclusions
Detecting super-intensive olive plantations may help to quantify the potential risks to local water resource management to address stakeholders and implement climate-smart agriculture technologies when needed.
While the use of remote sensing data and methods for mapping olive and tree crops has long been challenged by high spatial resolution requirements, we found that mapping super-intensive olive plantations comes with advantages compared to other tree crops and planting patterns. Using 5 m resolution allowed us to map all super-intensive olive plantations in the study area, except for areas with younger trees with smaller tree crown sizes.
Furthermore, the olive’s evergreen nature allowed us to use a reduced number of images per year. Especially when the aim was to extract super-intensive olive plantations within a single-class classification approach, as was the case in our study, we concluded that two satellite images were enough to suppress all other surrounding land use classes with higher variability in their phenology. This was made possible by using an unsupervised approach, which also allowed us to overcome the challenge of gathering labelled data, usually necessary for supervised classification methods.
The approach we developed, based on the principles of hierarchical clustering, consisted of applying a two-cluster k-means algorithm alternatingly on one image from the dry season and one from the winter season, extracting the vegetated cluster and creating a mask that was applied on the subsequent image to cluster. Using the NIR band was crucial to refine results and remove commission errors such as deciduous orchards with undercropping practices during the winter season, confirming findings from previous studies on tree crop mapping. In addition, we also concluded that NDVI was better to map SHD olive plantations than MSAVI-2, since they are less affected by soil background reflectance than other orchards of lower tree density.
Finally, this study leaves three open questions that may be addressed by future research: (i) while the performance of the approach we developed was quite promising, sample-wise verification of 2010 results revealed a larger share of commission errors among annual crops, suggesting that k-means may not be the best clustering algorithm for the proposed method, which further research may investigate; (ii) our study revealed that water-intensive SHD olive plantations, among other forms of land use, are a major player in the groundwater economy of the Saïss plain; however, further research is required to investigate which other land use classes are putting the Saïss aquifer at risk of depletion; (iii) lastly, our results showed that, despite the considerable commission errors among the small-sized plots in 2010, there was also a considerable number of larger plantations that appeared to have been given up and converted to other land uses. Future research may investigate the reasons for some areas converting from arboriculture to annual crops, despite the political framework encouraging the opposite.