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