USGS ScienceSearch

USGS · 70227642

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

Abstract

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 R 2 , 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.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 8.407168163601076° to 75.05035357407698° latitude; -174.7265625° to -52.734375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Yanjun Yang, Wei Ren, Bo Tao, Lei Ji, Li Li, Alex C. Ruran, J. B. Fisher, Jiangui Liu, Michael Sama, Zhe Li, Qingjiu Tian. 2020. Characterizing spatiotemporal patterns of crop phenology across North America during 2000–2016 using satellite imagery and agricultural survey data. https://doi.org/10.1016/j.isprsjprs.2020.10.005

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

National Land Cover Database 2019: A new strategy for creating clean leaf-on and leaf-off Landsat composite images

National Land Cover Database (NLCD) 2019 is a new epoch of national land cover products for the conterminous United States. Image quality is fundamental to the quality of any land cover product. Image preprocessing has often taken a considerable proportion of overall time and effort for this kind of national project. An approach to prepare image inputs for NLCD 2019 production was developed to ensure efficiency and quality of operational production. Here, we introduce a new and comprehensive strategy to produce clear Landsat composite images for NLCD 2019 production. First, we developed a new median-value compositing method. Second, we designed parameter settings for selecting images and pixels to generate 4 composite images (leaf-on, leaf-off, primary reference, and complementary reference) for a target year based on the US Landsat Analysis Ready Data surface reflectance dataset. Third, we developed a method, referred to as Detection and Filling with Simulated Image, to detect and replace clouds and cloud shadow pixels to produce the final clean leaf-on and leaf-off image composites. This image compositing and processing strategy was implemented for the entire conterminous United States to produce images for NLCD 2019. Our image results and NLCD 2019 change detection and land cover products, which were released in July 2021, showed this new strategy to be effective and efficient.

ISPRS Journal of Photogrammetry and Remote Sensing

A novel regression method for harmonic analysis of time series

Harmonic analysis of time series is an important technique in remote sensing to reveal seasonal land surface dynamics. However, frequency selection in the harmonic analysis is often difficult because high-frequency components are useful for delineating seasonal dynamics but sensitive to noise and gaps in time series. On the other hand, it is challenging to obtain temporally continuous satellite data with high quality because of atmospheric contamination. We developed a novel regression method named Harmonic Adaptive Penalty Operator (HAPO) for harmonic analysis of unevenly distributed time series. We introduced a new penalty function to minimize unexpected fluctuations in the model, which can substantially reduce the overfitting issue of regression in time series with temporal gaps. Specifically, the new penalty function minimizes the length of the model curve and the value range difference between the model and the time series observations. We compared HAPO with three widely used regression methods (OLS: Ordinary Least Squares; LASSO: Least Absolute Shrinkage and Selection Operator; and Ridge) in different scenarios using Landsat time series data across the United States. First, we evaluated methods using the Landsat surface reflectance time series within a single year. HAPO showed low and consistent monthly Root Mean Square Deviation (RMSD) values, in which most of the time RMSD of predicted reflectance were less than 0.04. More importantly, HAPO showed consistent and less bias given varying density and irregularity of time series. Second, we evaluated methods using multi-year time series. HAPO was a better predictor of relatively short time series (< 4 years) with steady low RMSD values. When a longer time series ( 4 years) was used, all four methods showed similar RMSD values, but HAPO outperformed the other methods if there were temporal gaps. Therefore, for places with large seasonal observation gaps or for time series that are relatively short (less than 4 years), HAPO can provide more consistent and accurate results in harmonic analysis of time series.

ISPRS Journal of Photogrammetry and Remote Sensing

Sentinel-2 and WorldView-3 atmospheric correction and signal normalization based on ground-truth spectroradiometric measurements

Remote sensing satellite Earth Observing Systems (EOS) provide a variety of products for monitoring Earth surface processes at varying spatial and spectral resolutions. Combining information from high and medium spatial resolution images is valuable for monitoring ground cover and vegetation status in cropland, grassland, forests, and other natural settings. However, coupling information from different EOS requires compensating for atmospheric and view angle effects before integrating comparable surface reflectance (SR) values. The objectives of this study were i) to assess how different atmospheric constituents affect the atmospheric correction results in Sentinel-2 and WorldView-3 imagery, ii) to establish a relationship with field spectra measurements, and iii) to develop an empirical approach to ensure that SR values extracted from different EOS can be normalized for use in monitoring vegetation and land cover status. We compared surface reflectance values derived from Sentinel-2 images corrected with Sen2Cor, MODTRAN or FLAASH atmospheric correction approaches for the visible-to-near infrared regions. Additionally, this information was compared to SR values extracted from WorldView-3 imagery acquired from the same dates and location (Central Spain) and corrected with MODTRAN and FLAASH approaches. Assessment of the atmospheric correction was conducted by comparing satellite image SR with ground-truth spectra acquired with a FieldSpec hand-held spectroradiometer. The results emphasized the importance of using common atmospheric parameters collected from ancillary data sources (i.e. MODIS Atmosphere & Land products) to ensure a reliable SR comparison. When compared to field-collected spectral data, SR from corrected Sentinel-2 push-broom imagery showed a reliable match (<4% difference in the visible bands and <0.52% difference in the near infrared bands). However, SR imagery from the pointable WorldView-3 instrument showed significant deviation, likely resulting from the effects of steep off-nadir acquisition angles (24.6° to 39.1°) combined with surface anisotropy. The magnitude and sign of the deviation in SR differed depending on the vegetation type, wavelength and sun-surface-sensor geometry. Therefore, it was necessary to account for angular effects to ensure reliable comparisons of imagery from the different EOS. In this study, an empirical angular correction approach was developed based on calibrating each WorldView-3 band against the ground-truth spectra. This correction allowed for the accurate signal normalization of WorldView-3 and Sentinel-2 imagery SR in the visible-to-near infrared regions.

ISPRS Journal of Photogrammetry and Remote Sensing