Remote sensing of plant phenology as an indicator of climate change and for mapping land cover has received significant scientific interest in the past two decades. The advancing of spring events, the lengthening of the growing season, the shifting of tree lines, the decreasing sensitivity to warming and the uniformity of spring across elevations are a few of the important indicators of trends in phenology. The Sentinel-2 satellite sensors launched in June 2015 (A) and March 2017 (B), with their high temporal frequency and spatial resolution for improved land mapping missions, have contributed significantly to knowledge on vegetation over the last three years. However, despite the additional red-edge and short wave infra-red (SWIR) bands available on the Sentinel-2 multispectral instruments, with improved vegetation species detection capabilities, there has been very little research on their efficacy to track vegetation cover and its phenology. For example, out of approximately every four papers that analyse normalised difference vegetation index (NDVI) or enhanced vegetation index (EVI) derived from Sentinel-2 imagery, only one mentions either SWIR or the red-edge bands. Despite the short duration that the Sentinel-2 platforms have been operational, they have proved their potential in a wide range of phenological studies of crops, forests, natural grasslands, and other vegetated areas, and in particular through fusion of the data with those from other sensors, e.g., Sentinel-1, Landsat and MODIS. This review paper discusses the current state of vegetation phenology studies based on the first five years of Sentinel-2, their advantages, limitations, and the scope for future developments.

Several methods exist for extracting plant phenological information from time series of satellite data. However, there have been only a few successful attempts to temporarily match satellite observations (Land Surface Phenology or LSP) with ground based phenological observations (Ground Phenology or GP). The classical pixel to point matching problem along with the temporal and spatial resolution of remote sensing data are some of the many issues encountered. In this study, MODIS-sensor’s Normalised Differenced Vegetation Index (NDVI) time series data were smoothed using two filtering techniques for comparison. Several start of season (SOS) methods established in the literature, namely thresholds of amplitude, derivatives and delayed moving average, were tested for determination of LSP-SOS for broadleaf forests at a site in southwestern Germany using 2001–2013 time series of NDVI data. The different LSP-SOS estimates when compared with species-rich GP dataset revealed that different LSP-SOS extraction methods agree better with specific phases of GP, and the choice of data processing or smoothing strongly affects the LSP-SOS extracted. LSP methods mirroring late SOS dates, i.e., 75% amplitude and 1st derivative, indicated a better match in means and trends, and high, significant correlations of up to 0.7 with leaf unfolding and greening of late understory and broadleaf tree species. GP-SOS of early understory leaf unfolding partly were significantly correlated with earlier detecting LSP-SOS, i.e., 20% amplitude and 3rd derivative. Early understory SOS were, however, more difficult to detect from NDVI due to the lack of a high resolution land cover information.

Phenology serves as a major indicator of ongoing climate change. Long-term phenological observations are critically important for tracking and communicating these changes. The phenological observation network across Germany is operated by the National Meteorological Service with a major contribution from volunteering activities. However, the number of observers has strongly decreased for the last decades, possibly resulting in increasing uncertainties when extracting reliable phenological information from map interpolation. We studied uncertainties in interpolated maps from decreasing phenological records, by comparing long-term trends based on grid-based interpolated and station-wise observed time series, as well as their correlations with temperature. Interpolated maps in spring were characterized by the largest spatial variabilities across Bavaria, Germany, with respective lowest interpolated uncertainties. Long-term phenological trends for both interpolations and observations exhibited mean advances of −0.2 to −0.3 days year⁻¹ for spring and summer, while late autumn and winter showed a delay of around 0.1 days year⁻¹. Throughout the year, temperature sensitivities were consistently stronger for interpolated time series than observations. Such a better representation of regional phenology by interpolation was equally supported by satellite-derived phenological indices. Nevertheless, simulation of observer numbers indicated that a decline to less than 40% leads to a strong decrease in interpolation accuracy. To better understand the risk of declining phenological observations and to motivate volunteer observers, a Shiny app is proposed to visualize spatial and temporal phenological patterns across Bavaria and their links to climate change–induced temperature changes.