Characterizing spatiotemporal patterns of crop phenology across North America during 2000–2016 using satellite imagery and agricultural survey data.

• Detected the major crop phenological stages using an improved EVI-curve-based approach over North America during 2000–2016. • Offered the first attempt to characterize the spatial and temporal patterns in crop phenology at a continental scale. • Crop planting and harvesting dates are highly correlated with the min/maximum temperatures. Crop phenology represents an integrative indicator of climate change and plays a vital role in terrestrial carbon dynamics and sustainable agricultural development. However, spatiotemporal variations of crop phenology remain unclear at large scales. This knowledge gap has hindered our ability to realistically quantify the biogeochemical dynamics in agroecosystems, predict future climate, and make informed decisions for climate change mitigation and adaptation. In this study, we improved an EVI-curve-based approach and used it to detect spatiotemporal patterns in cropping intensity and five major phenological stages over North America during 2000–2016 using vegetation index in combination with agricultural survey data and other ancillary maps. Our predicted crop phenological stages showed strong linear relationships with the survey-based datasets, with R2, RMSEs, and MAEs in the ranges of 0.35 –0.99, three to ten days, and two to eight days, respectively. During the study period, the planting dates were advanced by 0.60 days/year (p < 0.01), and harvesting dates were delayed by 0.78 days/year (p < 0.01) over North America. A minimum temperature increase by 1 °C caused a 4.26-day planting advance (r = −0.50, p < 0. 01) or a 0.66-day harvest delay (r = 0.10, p < 0.01). While, a higher maximum temperature resulted in a planting advance by 4.48 days/°C (r = −0.62, p < 0.01) or a harvest advance by 2.22 days/°C (r = −0.40, p < 0.01). Our analysis illustrated evident spatiotemporal variations in crop phenology in response to climate change and management practices. The derived crop phenological datasets and cropping intensity maps can be used in regional climate assessments and in developing adaptation strategies. [ABSTRACT FROM AUTHOR]