The quantity of land covered by various crops in a specific time span, referred to as a cropping pattern, dictates the level of agricultural production. However, retrieval of this information at a landscape scale can be challenging, especially when high spatial resolution imagery is not available. This study hypothesized that utilizing the unique advantages of multi-date and medium spatial resolution freely available Sentinel-2 (S2) reflectance bands (S2 bands), their vegetation indices (VIs) and vegetation phenology (VP) derivatives, and Sentinel-1 (S1) backscatter data would improve cropping pattern mapping in heterogeneous landscapes using robust machine learning algorithms, i.e., the guided regularized random forest (GRRF) for variable selection and the random forest (RF) for classification. This study’s objective was to map cropping patterns within three sub-counties in Murang’a County, a typical African smallholder heterogeneous farming area, in Kenya. Specifically, the performance of eight classification scenarios for mapping cropping patterns was compared, namely: (i) only S2 bands; (ii) S2 bands and VIs; (iii) S2 bands and VP; (iv) S2 bands and S1; (v) S2 bands, VIs, and S1; (vi) S2 bands, VP, and S1; (vii) S2 bands, VIs, and VP; and (viii) S2 bands, VIs, VP, and S1. Reference data of the dominant cropping patterns and non-croplands were collected. The GRRF algorithm was used to select the optimum variables in each scenario, and the RF was used to perform the classification for each scenario. The highest overall accuracy was 94.33% with Kappa of 0.93, attained using the GRRF-selected variables of scenario (v) S2, VIs, and S1. Furthermore, McNemar’s test of significance did not show significant differences (p ≤ 0.05) among the tested scenarios. This study demonstrated the strength of GRRF in selecting the most important variables and the synergetic advantage of S2 and S1 derivatives to accurately map cropping patterns in small-scale farming-dominated landscapes. Consequently, the cropping pattern mapping approach can be used in other sites of relatively similar agro-ecological conditions. Additionally, these results can be used to understand the sustainability of food systems and to model the abundance and spread of crop insect pests, diseases, and pollinators.

The proportion of area under various crops at a given point in time, known as a cropping pattern, plays an essential role in determining the level of agricultural production. In this study, cropping patterns of three sub-counties in Murang’a County, a typical African smallholder farming area in Kenya, were mapped. Specifically, we compared the performance of eight classification scenarios for mapping cropping patterns; namely using (i) only Sentinel-2 reflectance bands (S2), (ii) S2 and S2 derived vegetation indices (VIs); (iii) S2 and S2 vegetation phenology (VP); (iv) S2 and Sentinel-1 radar backscatter data (S1); (v) S2, VIs, and S1; (vi) S2, VP, and S1; (vii) S2, VIs and VP, and (viii) S2, VIs, VP and S1. Reference data of the dominant cropping patterns and non-croplands were collected. The guided regularized random forest (GRRF) algorithm was used to select the optimum variables and to perform the respective classification for each scenario. The most accurate result of the overall accuracy of 93.16% was attained from the scenario (viii) S2, VIs, VP, and S1. The McNemar’s test of significance did not show significant differences (p≤0.05) among the tested scenarios. Our study demonstrated the strength of GRRF and the synergetic advantage of S2 and S1 derivatives to map cropping patterns in a heterogeneous landscape where high resolution imagery are inaccessible. Our cropping pattern mapping approach can be used in other sites of relatively similar agro-ecological conditions.