A robust and unified land surface phenology algorithm for diverse biomes and growth cycles in China by using harmonized Landsat and Sentinel-2 imagery

Land surface phenology (LSP) is beneficial to understand ecosystem response to climate change, vegetation and crop type discrimination, and ecological modeling. However, the existing efforts based on coarse resolution data (≥500 m) cannot perform well in regions with higher spatial heterogeneity and multi-cropping system, such as China. Given that the majority of 10 m/30 m-based phenological research has focused on North America and Europe, developing spatiotemporally explicit LSP data in China is imperative. More importantly, the existing 30 m LSP products are mainly suitable for vegetation types with a single vegetation cycle, but cannot work well for biomes with complex seasonality (e.g., multiple growth cycles). Here we first harmonized three vegetation indices, i.e., the normalized difference vegetation index (NDVI), two-band enhanced vegetation index (EVI2), and land surface water index (LSWI) from Landsat-7/8 and Sentinel-2 imagery on the Google Earth Engine (GEE) platform. We then developed a new 30 m LSP algorithm that unified different phenological cycle-seeking processes per vegetation type and improved the existing algorithm. Furthermore, we used the algorithm to estimate the LSP product (LSP30CHN) for 2016–2019 across China, suitable for all vegetation types. The validation results showed a reasonably high accuracy (R2 > 0.6, RMSE < 15 days, mostly) of the LSP30CHN data against multi-sources in-situ observational (e.g., PhenoCam) and satellite-retrieved vegetation phenology data. Moreover, LSP30CHN data showed a consistent pattern but finer spatial details with the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) phenology product (MCD12Q2) at the homogenous area. We also found that phenological differences between LSP30CHN and MCD12Q2 increased with surface fragmentation, suggesting the potential of LSP30CHN to delineate phenological information on more fragmented landscapes. In contrast, the 500 m LSP data cannot provide such details in the regions with mixed cropping structures (e.g., corn, rice, and soybean) and multiple cropping index (e.g., single- and double-cropping systems). This study offers high accuracy of the LSP map for China, valuable for finer phenology-based services such as field-level crop management and agricultural phenology monitoring. It opens up new insights about exploring large-scale refined agricultural management and ecological assessment for other regions with complicated, fragmented landscapes and vegetation seasonality.