Land Surface Phenology Research Articles

Revisiting the Past: Replicability of a Historic Long-Term Vegetation Dynamics Assessment in the Era of Big Data Analytics

Open and analysis-ready data, as well as methodological and technical advancements have resulted in an unprecedented capability for observing the Earth’s land surfaces. Over 10 years ago, Landsat time series analyses were inevitably limited to a few expensive images from carefully selected acquisition dates. Yet, such a static selection may have introduced uncertainties when spatial or inter-annual variability in seasonal vegetation growth were large. As seminal pre-open-data-era papers are still heavily cited, variations of their workflows are still widely used, too. Thus, here we quantitatively assessed the level of agreement between an approach using carefully selected images and a state-of-the-art analysis that uses all available images. We reproduced a representative case study from the year 2003 that for the first time used annual Landsat time series to assess long-term vegetation dynamics in a semi-arid Mediterranean ecosystem in Crete, Greece. We replicated this assessment using all available data paired with a time series method based on land surface phenology metrics. Results differed fundamentally because the volatile timing of statically selected images relative to the phenological cycle introduced systematic uncertainty. We further applied lessons learned to arrive at a more nuanced and information-enriched vegetation dynamics description by decomposing vegetation cover into woody and herbaceous components, followed by a syndrome-based classification of change and trend parameters. This allowed for a more reliable interpretation of vegetation changes and even permitted us to disentangle certain land-use change processes with opposite trajectories in the vegetation components that were not observable when solely analyzing total vegetation cover. The long-term budget of net cover change revealed that vegetation cover of both components has increased at large and that this process was mainly driven by gradual processes. We conclude that study designs based on static image selection strategies should be critically evaluated in the light of current data availability, analytical capabilities, and with regards to the ecosystem under investigation. We recommend using all available data and taking advantage of phenology-based approaches that remove the selection bias and hence reduce uncertainties in results.

Impact of High-Cadence Earth Observation in Maize Crop Phenology Classification

For farmers, policymakers, and government agencies, it is critical to accurately define agricultural crop phenology and its spatial-temporal variability. At the moment, two approaches are utilized to report crop phenology. On one hand, land surface phenology provides information about the overall trend, whereas weekly reports from USDA-NASS provide information about the development of particular crops at the regional level. High-cadence earth observations might help to improve the accuracy of these estimations and bring more precise crop phenology classifications closer to what farmers demand. The second component of the proposed solution requires the use of robust classifiers (e.g., random forest, RF) capable of successfully managing large data sets. To evaluate this solution, this study compared the output of a RF classifier model using weather, two different satellite sources (Planet Fusion; PF and Sentinel-2; S-2), and ground truth data to improve maize (Zea mays L.) crop phenology classification using two regions of Kansas (Southwest and Central) as a testbed during the 2017 growing season. Our findings suggests that high temporal resolution (PF) data can significantly improve crop classification metrics (f1-score = 0.94) relative to S-2 (f1-score = 0.86). Additionally, a decline in the f1-score between 0.74 and 0.60 was obtained when we assessed the ability of S-2 to extend the temporal forecast for crop phenology. This research highlights the critical nature of very high temporal resolution (daily) earth observation data for crop monitoring and decision making in agriculture.

Near Real-time Fine-resolution Land Surface Phenological Prediction Using Convolutional Neural Network and Data Fusion

Near real-time fine-resolution land surface phenology (LSP) prediction is essential for understanding surface attributes and ecosystem functions, and solving important ecological processes related to phenology at the landscape scale. In this paper, we applied the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM) to fuse image pairs of Landsat 8 and Moderate-resolution Imaging Spectroradiometer (MODIS) as train data, and then applied the first derivative method to retrieve phenophase transition dates from fused time series of satellite data as label data. The convolutional neural network (CNN) model was trained using fusion images as inputs and the label data as targets. The trained model was further used to predict LSP dates from individual Landsat images. As evaluated using the reference data, the predict land surface phenological dates and could match the reference well with the coefficient of determination of 0.77 and root mean squared errors of 3.535, and our study provides an alternative method to predict land surface phenological dates using individual Landsat images.

How Relevant is High-Cadence Earth Observation for Maize Crop Phenology Classification?

Crop phenology can be defined as the study of biological processes such as emergence, flowering, and senescence that are associated with and affected by environmental growing conditions. The ability to reliably detect crop phenology and its spatial-temporal variability is critical for farmers, policymakers, and government agencies, since it has implications for the entire food chain. Currently, two methods are the most used to report crop phenology. Land surface phenology provides insight into the overall trend, whereas USDA-NASS weekly reports provide insight into the development of specific crops at the regional level. High-cadence earth observations may be able to improve the accuracy of these estimations and bring more precise crop phenology classifications closer to what farmers need. The use of robust classifiers (e.g., random forest, RF) to manage large data sets is required to successfully achieve this goal. This study compared the output of an RF classifier model using weather, two different satellite sources (Planet Fusion; PF and Sentinel-2; S-2), and ground truth data to improve maize (Zea mays L.) crop phenology classification during the 2017 growing season in Kansas. Our findings indicate that high-cadence (PF) data can enhance crop classification metrics (f1-score = 0.94) as compared to S-2 (f1-score = 0.86). This study emphasizes the significance of very high temporal resolution (daily) earth observation data for agricultural crop monitoring and decision-making tools.

Monitoring Cropland Phenology on Google Earth Engine Using Gaussian Process Regression.

Monitoring cropland phenology from optical satellite data remains a challenging task due to the influence of clouds and atmospheric artifacts. Therefore, measures need to be taken to overcome these challenges and gain better knowledge of crop dynamics. The arrival of cloud computing platforms such as Google Earth Engine (GEE) has enabled us to propose a Sentinel-2 (S2) phenology end-to-end processing chain. To achieve this, the following pipeline was implemented: (1) the building of hybrid Gaussian Process Regression (GPR) retrieval models of crop traits optimized with active learning, (2) implementation of these models on GEE (3) generation of spatiotemporally continuous maps and time series of these crop traits with the use of gap-filling through GPR fitting, and finally, (4) calculation of land surface phenology (LSP) metrics such as the start of season (SOS) or end of season (EOS). Overall, from good to high performance was achieved, in particular for the estimation of canopy-level traits such as leaf area index (LAI) and canopy chlorophyll content, with normalized root mean square errors (NRMSE) of 9% and 10%, respectively. By means of the GPR gap-filling time series of S2, entire tiles were reconstructed, and resulting maps were demonstrated over an agricultural area in Castile and Leon, Spain, where crop calendar data were available to assess the validity of LSP metrics derived from crop traits. In addition, phenology derived from the normalized difference vegetation index (NDVI) was used as reference. NDVI not only proved to be a robust indicator for the calculation of LSP metrics, but also served to demonstrate the good phenology quality of the quantitative trait products. Thanks to the GEE framework, the proposed workflow can be realized anywhere in the world and for any time window, thus representing a shift in the satellite data processing paradigm. We anticipate that the produced LSP metrics can provide meaningful insights into crop seasonal patterns in a changing environment that demands adaptive agricultural production.

Earth observation method for an operational assessment of crop phenology metric “start of the season”

Vegetation phenology is one of the key indicators for assessing inter annual changes of ecosystems while crop phenology has crucial and practical importance for precise timing of farming practices of today’s smart agriculture. The aim of this study is to develop and implement a procedure for filtering and processing daily satellite signals from Moderate Resolution Imaging Spectroradiometer (MODIS) satellites for the operational derivation of vegetation indices and land surface phenology (LSP) metrics. Second, this work investigates and selects precisely adjusted parameters for the calculation of a specific LSP metric – the SOS (start of the season). The SOS obtained from the satellite data was compared with available in situ ground measurements at 80 sites with records from 2000-2012 across the Czech Republic. The coefficients of determination, R2, for the land surface phenology based on Earth observation and traditional ground phenology indicated differences between the two methods. However, the average R2 at all sites for agricultural land was 0.31 for winter wheat onset of leaf sheath elongation; the R2 exceeded 0.5 at 21 sites and was approximately 0.7 at some sites. The developed procedure proved suitable for operational monitoring of actual crop conditions and other influence factors that are important for drought monitoring.

Comparing Time-Lapse PhenoCams with Satellite Observations across the Boreal Forest of Quebec, Canada

Intercomparison of satellite-derived vegetation phenology is scarce in remote locations because of the limited coverage area and low temporal resolution of field observations. By their reliable near-ground observations and high-frequency data collection, PhenoCams can be a robust tool for intercomparison of land surface phenology derived from satellites. This study aims to investigate the transition dates of black spruce (Picea mariana (Mill.) B.S.P.) phenology by comparing fortnightly the MODIS normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI) extracted using the Google Earth Engine (GEE) platform with the daily PhenoCam-based green chromatic coordinate (GCC) index. Data were collected from 2016 to 2019 by PhenoCams installed in six mature stands along a latitudinal gradient of the boreal forests of Quebec, Canada. All time series were fitted by double-logistic functions, and the estimated parameters were compared between NDVI, EVI, and GCC. The onset of GCC occurred in the second week of May, whereas the ending of GCC occurred in the last week of September. We demonstrated that GCC was more correlated with EVI (R2 from 0.66 to 0.85) than NDVI (R2 from 0.52 to 0.68). In addition, the onset and ending of phenology were shown to differ by 3.5 and 5.4 days between EVI and GCC, respectively. Larger differences were detected between NDVI and GCC, 17.05 and 26.89 days for the onset and ending, respectively. EVI showed better estimations of the phenological dates than NDVI. This better performance is explained by the higher spectral sensitivity of EVI for multiple canopy leaf layers due to the presence of an additional blue band and an optimized soil factor value. Our study demonstrates that the phenological observations derived from PhenoCam are comparable with the EVI index. We conclude that EVI is more suitable than NDVI to assess phenology in evergreen species of the northern boreal region, where PhenoCam data are not available. The EVI index could be used as a reliable proxy of GCC for monitoring evergreen species phenology in areas with reduced access, or where repeated data collection from remote areas are logistically difficult due to the extreme weather.

Comprehensive Assessment of Performances of Long Time-Series LAI, FVC and GPP Products over Mountainous Areas: A Case Study in the Three-River Source Region, China

Vegetation biophysical products offer unique opportunities to examine long-term vegetation dynamics and land surface phenology (LSP). It is important to understand the time-series performances of various global biophysical products for global change research. However, few endeavors have been dedicated to assessing the performances of long-term change characteristics or LSP extraction derived from different satellite products, especially in mountainous areas with highly fragmented and rugged surfaces. In this paper, we assessed the time-series characteristics and LSP detections of Global LAnd Surface Satellite (GLASS) leaf area index (LAI), fractional vegetation cover (FVC), and gross primary production (GPP) products across the Three-River Source Region (TRSR). The performances of products’ temporal agreements and their statistical relationship as a function of topographic indices and heterogeneous pixels, respectively, were investigated through intercomparison among three products during the period 2000 to 2018. The results show that the phenological differences between FVC and two other products are beyond 10 days over more than 35% of the pixels in TRSR. The long-term trend of FVC diverges significantly from GPP and LAI for 13.96% of the total pixels, and the percentages of mismatched pixels between FVC and two other products are 33.24% in the correlation comparison. Moreover, good agreements are observed between GPP and LAI, both in terms of LSP and interannual variations. Finally, the LSP and long-term dynamics of the three products exhibit poor performances on heterogeneous surfaces and complex topographic areas, which reflects the potential impacts of environmental factors and algorithmic imperfections on the quality and performances of different products. Our study highlights the spatiotemporal disparities in detections of surface vegetation activity in mountainous areas by using different biophysical products. Future global change studies may require multiple high-quality satellite products with long-term stability as data support.

Greenup dates change across a temperate forest-grassland ecotone in northeastern China driven by spring temperature and tree cover

Changes in climate and land cover are potential drivers of land surface phenology changes. Here, we investigate whether subpixel percent tree cover (PTC) change is an important driver of trends in satellite derived vegetation spring greenup date (GUD) across the Hulunbuir temperate forest-grassland ecotone in northeastern China. GUD was estimated using the MODIS-derived enhanced vegetation index time series during 2001–2020 with a spatial resolution of 500 m. To understand the influential mechanisms of PTC on GUD, we examined relationships between the spatial variations in GUD and PTC at multiple spatio-temporal extents. Forested pixels with greater PTC were found to have generally earlier GUDs for all forest types. The GUD of forests was also generally earlier than that of grassland. On the other hand, we observed approximately 23.7% and 1.2% significantly earlier and later trends in GUD across the region, respectively. Meanwhile, widespread increases in preseason land surface temperature (LST) and PTC were detected. Both increases in LST and PTC contributed to the earlier GUD in the forested region. Specifically, we found negative correlations (Spearman correlation coefficient -0.17 to -0.55) between the change slopes of GUD and PTC in every forest and grassland type. The results highlight the important impacts of subpixel PTC on GUD variations, and improve the understanding of ecosystem changes under the effects of climate and human activities (e.g., afforestation) over the Hulunbuir temperate forest-grassland ecotone.

Global Vegetation Photosynthetic Phenology Products Based on MODIS Vegetation Greenness and Temperature: Modeling and Evaluation

Land surface phenology (LSP) products that are derived from different data sources have different definitions and biophysical meanings. Discrepancies among these products and their linkages with carbon fluxes across plant functional types and climatic regions remain somewhat unclear. In this study, to differentiate LSP related to gross primary production (GPP) from LSP related to remote sensing data, we defined the former as vegetation photosynthetic phenology (VPP), including the starting and ending days of GPP (SOG and EOG, respectively). Specifically, we estimated VPP based on a combination of observed VPP from 145 flux-measured GPP sites together with the vegetation index and temperature data from MODIS products using multiple linear regression models. We then compared VPP estimates with MODIS LSP on a global scale. Our results show that the VPP provided better estimates of SOG and EOG than MODIS LSP, with a root mean square error (RMSE) for SOG of 12.7 days and a RMSE for EOG of 10.5 days. The RMSE was approximately three weeks for both SOG and EOG estimates of the non-forest type. Discrepancies between VPP and LSP estimates varied across plant functional types (PFTs) and climatic regions. A high correlation was observed between VPP and LSP estimates for deciduous forest. For most PFTs, using VPP estimates rather than LSP improved the estimation of GPP. This study presents a useful method for modeling global VPP, investigates in detail the discrepancies between VPP and LSP, and provides a more effective global vegetation phenology product for carbon cycle modeling than the existing ones.

Widespread Mismatch Between Phenology and Climate in Human‐Dominated Landscapes

AbstractPlants track changing climate partly by shifting their phenology, the timing of recurring biological events. It is unknown whether these observed phenological shifts are sufficient to keep pace with rapid climate changes. Phenological mismatch, or the desynchronization between the timing of critical phenological events, has long been hypothesized but rarely quantified on a large scale. It is even less clear how human activities have contributed to this emergent phenological mismatch. In this study, we used remote sensing observations to systematically evaluate how plant phenological shifts have kept pace with warming trends at the continental scale. In particular, we developed a metric of spatial mismatch that connects empirical spatiotemporal data to ecological theory using the “velocity of change” approach. In northern mid‐to high‐latitude regions (between 30–70°N) over the last three decades (1981–2014), we found evidence of a widespread mismatch between land surface phenology and climate where isolines of phenology lag behind or move in the opposite direction to the isolines of climate. These mismatches were more pronounced in human‐dominated landscapes, suggesting a relationship between human activities and the desynchronization of phenology dynamics with climate variations. Results were corroborated with independent ground observations that indicate the mismatch of spring phenology increases with human population density for several plant species. This study reveals the possibility that not even some of the foremost responses in vegetation activity match the pace of recent warming. This systematic analysis of climate‐phenology mismatch has important implications for the sustainable management of vegetation in human‐dominated landscapes under climate change.

Mapping grassland mowing events across Germany based on combined Sentinel-2 and Landsat 8 time series

Spatially explicit knowledge on grassland extent and management is critical to understand and monitor the impact of grassland use intensity on ecosystem services and biodiversity. While regional studies allow detailed insights into land use and ecosystem service interactions, information on a national scale can aid biodiversity assessments. However, for most European countries this information is not yet widely available. We used an analysis-ready-data cube that contains dense time series of co-registered Sentinel-2 and Landsat 8 data, covering the extent of Germany. We propose an algorithm that detects mowing events in the time series based on residuals from an assumed undisturbed phenology, as an indicator of grassland use intensity. A self-adaptive ruleset enabled to account for regional variations in land surface phenology and non-stationary time series on a pixel-basis. We mapped mowing events for the years from 2017 to 2020 for permanent grassland areas in Germany. The results were validated on a pixel level in four of the main natural regions in Germany based on reported mowing events for a total of 92 (2018) and 78 (2019) grassland parcels. Results for 2020 were evaluated with combined time series of Landsat, Sentinel-2 and PlanetScope data. The mean absolute percentage error between detected and reported mowing events was on average 40% (2018), 36% (2019) and 35% (2020). Mowing events were on average detected 11 days (2018), 7 days (2019) and 6 days (2020) after the reported mowing. Performance measures varied between the different regions of Germany, and lower accuracies were found in areas that are revisited less frequently by Sentinel-2. Thus, we assessed the influence of data availability and found that the detection of mowing events was less influenced by data availability when at least 16 cloud-free observations were available in the grassland season. Still, the distribution of available observations throughout the season appeared to be critical. On a national scale our results revealed overall higher shares of less intensively mown grasslands and smaller shares of highly intensively managed grasslands. Hotspots of the latter were identified in the alpine foreland in Southern Germany as well as in the lowlands in the Northwest of Germany. While these patterns were stable throughout the years, the results revealed a tendency to lower management intensity in the extremely dry year 2018. Our results emphasize the ability of the approach to map the intensity of grassland management throughout large areas despite variations in data availability and environmental conditions.

Quantifying latitudinal variation in land surface phenology of Spartina alterniflora saltmarshes across coastal wetlands in China by Landsat 7/8 and Sentinel-2 images

Spartina alterniflora is an aggressive invasive plant spreading along the coastal China, spanning a latitudinal range of 20°N - 39°N, and its invasion resulted in dramatic decline in both native plant diversity and ecosystem functioning. Phenology of S. alterniflora saltmarshes is a critical feature to elucidate the invasion dynamics over geographical regions but has not been well understood yet. In this study, we examined the variation of S. alterniflora saltmarsh phenology across coastal China during 2018–2020 by using time series Landsat 7/8 and Sentinel-2 images. Combined Landsat 7/8 and Sentinel-2 images provided more good-quality observations in a year, which could facilitate phenological retrieval. We applied and assessed three widely used phenology retrieval methods (i.e., NDVI-based pixel-specific statistical threshold, NDVI-based double logistic mathematical equation, and LSWI-based biological threshold) for retrieving the start and end of season (SOS and EOS) as well as the length of season (LOS) of S. alterniflora saltmarshes. The SOS and EOS dates derived from three phenology retrieval methods showed similar patterns in latitudinal phenology variation: SOS became later and LOS became shorter as latitude increased, and the latitudinal trend of EOS was not as large as that of SOS. This study shows the potential of Landsat 7/8 and Sentinel-2 to quantify land surface phenology of S. alterniflora saltmarshes, which not only enhances our understanding of the spatial-temporal dynamics of coastal saltmarshes in China but also improves the management of this plant invader that threatens native saltmarshes in the world.

Mapping burn severity in the western Italian Alps through phenologically coherent reflectance composites derived from Sentinel-2 imagery

Deriving burn severity from multispectral satellite data is a widely adopted approach to infer the degree of environmental change caused by fire. Burn severity maps obtained by thresholding bi-temporal indices based on pre- and post-fire Normalized Burn Ratio (NBR) can vary substantially depending on temporal constraints such as matched acquisition and optimal seasonal timing. Satisfying temporal requirements is crucial to effectively disentangle fire and non-fire induced spectral changes and can be particularly challenging when only a few cloud-free images are available. Our study focuses on 10 wildfires that occurred in mountainous areas of the Piedmont Region (Italy) during autumn 2017 following a severe and prolonged drought period. Our objectives were to: (i) generate reflectance composites using Sentinel-2 imagery that were optimised for seasonal timing by embedding spatial patterns of long-term land surface phenology (LSP); (ii) produce and validate burn severity maps based on the modelled relationship between bi-temporal indices and field data; (iii) compare burn severity maps obtained using either a pair of cloud-free Sentinel-2 images, i.e. paired images, or reflectance composites. We proposed a pixel-based compositing algorithm coupling the weighted geometric median and thematic spatial information, e.g. long-term LSP metrics derived from the MODIS Collection 6 Land Cover Dynamics Product, to rank all the clear observations available in the growing season. Composite Burn Index data and bi-temporal indices exhibited a strong nonlinear relationship (R2 > 0.85) using paired images or reflectance composites. Burn severity maps attained overall classification accuracy ranging from 76.9% to 83.7% (Kappa between 0.61 and 0.72) and the Relative differenced NBR (RdNBR) achieved the best results compared to other bi-temporal indices (differenced NBR and Relativized Burn Ratio). Improvements in overall classification accuracy offered by the calibration of bi-temporal indices with the dNBR offset were limited to burn severity maps derived from paired images. Reflectance composites provided the highest overall classification accuracy and differences with paired images were significant using uncalibrated bi-temporal indices (4.4% to 5.2%) while they decreased (2.8% to 3.2%) when we calibrated bi-temporal indices derived from paired images. The extent of the high severity category increased by ~19% in burn severity maps derived from reflectance composites (uncalibrated RdNBR) compared to those from paired images (calibrated RdNBR). The reduced contrast between healthy and burnt conditions associated with suboptimal seasonal timing caused an underestimation of burnt areas. By embedding spatial patterns of long-term LSP metrics, our approach provided consistent reflectance composites targeted at a specific phenological stage and minimising non-fire induced inter-annual changes. Being independent from the multispectral dataset employed, the proposed pixel-based compositing approach offers new opportunities for operational change detection applications in geographic areas characterised by persistent cloud cover.

Multi-year monitoring land surface phenology in relation to climatic variables using MODIS-NDVI time-series in Mediterranean forest, Northeast Tunisia

The Mediterranean region is one of the most vulnerable regions to climate change. The majority of climate models forecast a rise in temperatures and less rainfall, which have been observed in recent decades. These changes will affect several vegetation properties, especially phenological dynamics and traits, by increasing drought intensity and recurrence. In this climate change context, the present study aims to assess the evolution of vegetation state and its relation with the climate dynamics in the Mediterranean forest region of northeast Tunisia using Land Surface Phenology (LSP) metrics and the vegetation index (NDVI) analysis from 2000 to 2017. To conduct this work, we used precipitation and temperature data from the two closest weather stations and 16-day NDVI composite images from the MODIS satellite source, with 250-m spatial resolution. Three phenological metrics— start of season (SOS), end of season (EOS), and length of season (LOS) — were obtained and compared for different vegetation types. The LSP variation in response to climatic metrics was also analyzed. The results showed that the LSP in our study area changed significantly during the 2000–2017 period, which includes an average 7.8 days delay in the SOS, an average advance in the EOS by 5 days, and LOS shortened by an average 12.8 days. Autumn (Pr_9) and spring (Pr_3 and P3_4) precipitations, as well as maximum temperature (Tx9+10), represent the best climate parameters to explain the changes in LSP. Both the NDVI and SPEI showed a significant high correlation (p < 0.001) on longer time scales. LSP and NDVI proved useful tools for monitoring the vegetation state according to climate for better planning territory purposes.

Exploring the Applicability and Scaling Effects of Satellite-Observed Spring and Autumn Phenology in Complex Terrain Regions Using Four Different Spatial Resolution Products

The information on land surface phenology (LSP) was extracted from remote sensing data in many studies. However, few studies have evaluated the impacts of satellite products with different spatial resolutions on LSP extraction over regions with a heterogeneous topography. To bridge this knowledge gap, this study took the Loess Plateau as an example region and employed four types of satellite data with different spatial resolutions (250, 500, and 1000 m MODIS NDVI during the period 2001–2020 and ~10 km GIMMS3g during the period 1982–2015) to investigate the LSP changes that took place. We used the correlation coefficient (r) and root mean square error (RMSE) to evaluate the performances of various satellite products and further analyzed the applicability of the four satellite products. Our results showed that the MODIS-based start of the growing season (SOS) and end of the growing season (EOS) were highly correlated with the ground-observed data with r values of 0.82 and 0.79, respectively (p &lt; 0.01), while the GIMMS3g-based phenology signal performed badly (r &lt; 0.50 and p &gt; 0.05). Spatially, the LSP that was derived from the MODIS products produced more reasonable spatial distributions. The inter-annual averaged MODIS SOS and EOS presented overall advanced and delayed trends during the period 2001–2020, respectively. More than two-thirds of the SOS advances and EOS delays occurred in grasslands, which determined the overall phenological changes across the entire Loess Plateau. However, both inter-annual trends of SOS and EOS derived from the GIMMS3g data were opposite to those seen in the MODIS results. There were no significant differences among the three MODIS datasets (250, 500, and 1000 m) with regard to a bias lower than 2 days, RMSE lower than 1 day, and correlation coefficient greater than 0.95 (p &lt; 0.01). Furthermore, it was found that the phenology that was derived from the data with a 1000 m spatial resolution in the heterogeneous topography regions was feasible. Yet, in forest ecosystems and areas with an accumulated temperature ≥10 °C, the differences in phenological phase between the MODIS products could be amplified.

Specific Drivers and Responses to Land Surface Phenology of Different Vegetation Types in the Qinling Mountains, Central China

Land surface phenology (LSP), as a precise bio-indicator that responds to climate change, has received much attention in fields concerned with climate change and ecology. Yet, the dynamics of LSP changes in the Qinling Mountains (QMs)—A transition zone between warm-temperate and north subtropical climates with complex vegetation structure—under significant climatic environmental evolution are unclear. Here, we analyzed the spatiotemporal dynamics of LSP for different vegetation types in the QMs from 2001 to 2019 and quantified the degree of influence of meteorological factors (temperature, precipitation, and shortwave radiation), and soil (temperature and moisture), and biological factors (maximum of NDVI and middle date during the growing season) on LSP changes using random forest models. The results show that there is an advanced trend (0.15 days/year) for the start of the growing season (SOS), a delayed trend (0.24 days/year) for the end of the growing season (EOS), and an overall extended trend (0.39 days/year) for the length of the growing season (LOS) in the QMs over the past two decades. Advanced SOS and delayed EOS were the dominant patterns leading to a lengthened vegetation growing season, followed by a joint delay of SOS and EOS, and the latter was particularly common in shrub and evergreen broadleaved forests. The growth season length increased significantly in western QMs. Furthermore, we confirmed that meteorological factors are the main factors affecting the interannual variations in SOS and EOS, especially the meteorological factor of preseason mean shortwave radiation (SWP). The grass and crop are most influenced by SWP. The soil condition has, overall, a minor influence the regional LSP. This study highlighted the specificity of different vegetation growth in the QMs under warming, which should be considered in the accurate prediction of vegetation growth in the future.

Fusing Geostationary Satellite Observations with Harmonized Landsat-8 and Sentinel-2 Time Series for Monitoring Field-Scale Land Surface Phenology

Accurate and timely land surface phenology (LSP) provides essential information for investigating the responses of terrestrial ecosystems to climate changes and quantifying carbon and surface energy cycles on the Earth. LSP has been widely investigated using daily Visible Infrared Imaging Radiometer Suite (VIIRS) or Moderate Resolution Imaging Spectroradiometer (MODIS) observations, but the resultant phenometrics are frequently influenced by surface heterogeneity and persistent cloud contamination in the time series observations. Recently, LSP has been derived from Landsat-8 and Sentinel-2 time series providing detailed spatial pattern, but the results are of high uncertainties because of poor temporal resolution. With the availability of data from Advanced Baseline Imager (ABI) onboard a new generation of geostationary satellites that observe the earth every 10–15 min, daily cloud-free time series could be obtained with high opportunities. Therefore, this study investigates the generation of synthetic high spatiotemporal resolution time series by fusing the harmonized Landsat-8 and Sentinel-2 (HLS) time series with the temporal shape of ABI data for monitoring field-scale (30 m) LSP. The algorithm is verified by detecting the timings of greenup and senescence onsets around north Wisconsin/Michigan states, United States, where cloud cover is frequent during spring rainy season. The LSP detections from HLS-ABI are compared with those from HLS or ABI alone and are further evaluated using PhenoCam observations. The result indicates that (1) ABI could provide ~3 times more high-quality observations than HLS around spring greenup onset; (2) the greenup and senescence onsets derived from ABI and HLS-ABI are spatially consistent and statistically comparable with a median difference less than 1 and 10-days, respectively; (3) greenup and senescence onsets derived from HLS data show sharp boundaries around the orbit-overlapped areas and shifts of ~13 days delay and ~15 days ahead, respectively, relative to HLS-ABI detections; and (4) HLS-ABI greenup and senescence onsets align closely to PhenoCam observations with an absolute average difference of less than 2 days and 5 days, respectively, which are much better than phenology detections from ABI or HLS alone. The result suggests that the proposed approach could be implemented the monitor of 30 m LSP over regions with persistent cloud cover.

A global increase in tree cover extends the growing season length as observed from satellite records

Plant phenology provides information on the seasonal dynamics of plants, and changes herein are important for understanding the impact of climate change and human management on the biosphere. Land surface phenology is the study of plant phenology across large spatial scales estimated by satellite observations. However, satellite observations (pixels) are often composed of a mixture of vegetation types, like woody vegetation and herbaceous vegetation, having different phenological characteristics. Therefore, any changes in tree cover presumably impact land surface phenology, as trees usually have a different seasonal cycle compared to herbaceous vegetation. On the other hand, changes in land surface phenology are often interpreted as a result of climate change-induced impacts on the photosynthetic activity of vegetation. Therefore, it is important to better understand the role of changes in vegetation cover (here, the proportion between tree and short vegetation cover) in satellite-derived land surface phenology analysis. We studied the impact of changes in tree cover on satellite observed land surface phenology at a global scale over the past three decades. We found an extension of the growing season length in 36.6% of the areas where tree cover increased, whereas only 20.1% of the areas where tree cover decreased showed an increase in growing season length. Furthermore, the ratio between tree cover and short vegetation cover was found to affect changes in the length of the growing season, with the denser tree cover showing a more pronounced extension of the growing season length (especially in boreal forests). These results highlight the importance of changes in tree cover when analyzing the impact of climate change on vegetation phenology. Our study thereby addresses a critical knowledge gap for an improved understanding of changes in land surface phenology during recent decades in the context of climate and human-induced global land cover change.

GCI30: a global dataset of 30 m cropping intensity using multisource remote sensing imagery

Abstract. The global distribution of cropping intensity (CI) is essential to our understanding of agricultural land use management on Earth. Optical remote sensing has revolutionized our ability to map CI over large areas in a repeated and cost-efficient manner. Previous studies have mainly focused on investigating the spatiotemporal patterns of CI ranging from regions to the entire globe with the use of coarse-resolution data, which are inadequate for characterizing farming practices within heterogeneous landscapes. To fill this knowledge gap, in this study, we utilized multiple satellite data to develop a global, spatially continuous CI map dataset at 30 m resolution (GCI30). Accuracy assessments indicated that GCI30 exhibited high agreement with visually interpreted validation samples and in situ observations from the PhenoCam network. We carried out both statistical and spatial comparisons of GCI30 with six existing global CI estimates. Based on GCI30, we estimated that the global average annual CI during 2016–2018 was 1.05, which is close to the mean (1.09) and median (1.07) CI values of the existing six global CI estimates, although the spatial resolution and temporal coverage vary significantly among products. A spatial comparison with two satellite-based land surface phenology products further suggested that GCI30 was not only capable of capturing the overall pattern of global CI but also provided many spatial details. GCI30 indicated that single cropping was the primary agricultural system on Earth, accounting for 81.57 % (12.28×106 km2) of the world's cropland extent. Multiple-cropping systems, on the other hand, were commonly observed in South America and Asia. We found large variations across countries and agroecological zones, reflecting the joint control of natural and anthropogenic drivers on regulating cropping practices. As the first global-coverage, fine-resolution CI product, GCI30 is expected to fill the data gap for promoting sustainable agriculture by depicting worldwide diversity of agricultural land use intensity. The GCI30 dataset is available on Harvard Dataverse: https://doi.org/10.7910/DVN/86M4PO (Zhang et al., 2020).