In this study, the potential of phenological indicators derived from Sentinel-2A (S2) time series were evaluated to explore the key variables that allow identifying both cropland and crop types. Based on the derived S2 phenological metrics and fitted vegetation indices (VI), ten feature sets were developed and assessed to discriminate different crop types via Random Forest (RF) classifier. The comparison between VI data-based classifications has shown that NDVI and EVI2 phenological sets could delineate and identify crop types more accurately compared to RENDVI data. Overall, the combined use of fitted VI and phenological features rather than being used separately achieved the best performances. Further, the result of using optimum features was the most accurate among the ten feature sets, with an overall accuracy of 88% and Kappa of 0.84. This study constitutes a substantial improvement in crop-types identification, which gives a valuable tool to monitor agricultural areas.

Changing land use patterns is of great importance in environmental studies and critical for land use management decision making over farming systems in arid and semi-arid regions. Unfortunately, ground data scarcity or inadequacy in many regions can cause large uncertainties in the characterization of phenological changes in arid and semi-arid regions, which can hamper tailored decision making towards best agricultural management practices. Alternatively, state-of-the-art methods for phenological metrics’ extraction and long time-series analysis techniques of multispectral remote sensing imagery provide a viable solution. In this context, this study aims to characterize the changes over farming systems through trend analysis. To this end, four farming systems (fallow, rainfed, irrigated annual, and irrigated perennial) in arid areas of Morocco were studied based on four phenological metrics (PhM) (i.e., great integral, start, end, and length of the season). These were derived from large Moderate resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) time-series using both a machine learning algorithm and a pixel-based change analysis method. Results showed that during the last twenty-year period (i.e., 2000–2019), a significant dynamism of the plant cover was linked to the behavior of farmers who tend to cultivate intensively and to invest in high-income crops. More specifically, a relevant variability in fallow and rainfed areas, closely linked to the weather conditions, was found. In addition, significant lag trends of the start (−6 days) and end (+3 days) were found, which indicate that the length of the season was related to the spatiotemporal variability of rainfall. This study has also highlighted the potential of multitemporal moderate spatial resolution data to accurately monitor agriculture and better manage land resources. In the meantime, for operationally implementing the use of such work in the field, we believe that it is essential consider the perceptions, opinions, and mutual benefits of farmers and stakeholders to improve strategies and synergies whilst ensuring food, welfare, and sustainability.