Annual Cropping Intensity Dynamics in China from 2001 to 2023

Spatial and temporal information about cropping patterns of single and multiple crops is important for monitoring crop production and land-use intensity. We used time-series MODIS NDVI 8-day composite data to develop annual cropping pattern products at a 250 m spatial resolution for China, covering the period from 2001 to 2023. To address the potential impacts of varying parameters in both data pre-processing and the peak detection algorithm on the accuracy of cropping pattern mapping, we employed a grid-search method to fine-tune these parameters. This process focused on optimizing the Savitzky–Golay smoothing window size and the peak width parameters using a calibration dataset. The results highlighted that an optimal combination of a five to seven MODIS composite window size in Savitzky–Golay smoothing and a peak width of four MODIS composites achieved good overall mapping accuracy. Pixel-wise accuracy assessments were conducted for the selected mapping years of 2001, 2011, and 2021. Overall accuracies were between 89.7% and 92.0%, with F1 scores ranging from 0.921 to 0.943. Nationally, this study observed a fluctuating trend in multiple cropping percentages, with a notable increase after 2013, suggesting shifts toward more intensive agricultural practices in recent years. At a finer spatial scale, the combination of Mann–Kendall and Sen’s slope analyses revealed that approximately 12.9% of 3 km analytical windows exhibited significant changes in cropping intensity. We observed spatial clusters of increasing and decreasing crop intensity trends across provinces such as Hebei, Shandong, Shaanxi, and Gansu. This study underscores the importance of data smoothing and peak detection methods in analyzing high temporal resolution remote sensing data. The generation of annual single/multiple cropping pattern maps at a 250 m spatial resolution enhances our comprehension of agricultural dynamics through time and across different regions.