Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Temporal Trends and Spatial Variability of Vegetation Phenology over the Northern Hemisphere during 1982-2012

  • Siyuan Wang ,

    Contributed equally to this work with: Siyuan Wang, Bojuan Yang

    w_siyuan@126.com

    Affiliation Key Laboratory of Digital Earth, Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing, China

  • Bojuan Yang ,

    Contributed equally to this work with: Siyuan Wang, Bojuan Yang

    Affiliation Key Laboratory of Digital Earth, Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing, China

  • Qichun Yang ,

    ‡ These authors also contributed equally to this work.

    Affiliation Ecosystem Dynamics and Global Ecology Laboratory, School of Forestry and Wildlife Sciences, Auburn University, Auburn, Alabama, United States of America

  • Linlin Lu ,

    ‡ These authors also contributed equally to this work.

    Affiliation Key Laboratory of Digital Earth, Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing, China

  • Xiaoyue Wang ,

    ‡ These authors also contributed equally to this work.

    Affiliation Key Laboratory of Digital Earth, Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing, China

  • Yaoyao Peng

    ‡ These authors also contributed equally to this work.

    Affiliation Key Laboratory of Digital Earth, Institute of Remote Sensing and Digital Earth, Chinese Academy of Science, Beijing, China

Abstract

Satellite-derived vegetation phenology has been recognized as a key indicator for detecting changes in the terrestrial biosphere in response to global climate change. However, multi-decadal changes and spatial variation of vegetation phenology over the Northern Hemisphere and their relationship to climate change have not yet been fully investigated. In this article, we investigated the spatial variability and temporal trends of vegetation phenology over the Northern Hemisphere by calibrating and analyzing time series of the satellite-derived normalized difference vegetation index (NDVI) during 1982–2012, and then further examine how vegetation phenology responds to climate change within different ecological zones. We found that during the period from 1982 to 2012 most of the high latitude areas experienced an increase in growing period largely due to an earlier beginning of vegetation growing season (BGS), but there was no significant trend in the vegetation growing peaks. The spatial pattern of phenology within different eco-zones also experienced a large variation over the past three decades. Comparing the periods of 1982–1992, 1992–2002 with 2002–2012, the spatial pattern of change rate of phenology shift (RPS) shows a more significant trend in advancing of BGS, delaying of EGS (end of growing season) and prolonging of LGS (length of growing season) during 2002–2012, overall shows a trend of accelerating change. Temperature is a major determinant of phenological shifts, and the response of vegetation phenology to temperature varied across different eco-zones.

Introduction

Vegetation phenology is highly sensitive to climate change, so monitoring vegetation phenology is becoming an increasing important way to understand and quantify global environmental changes [1, 2]. Since the 1950s, the world’s mean surface temperature has increased by 0.6°C, and the warming trend has been faster in the high latitude areas, such as Eurasian continent [3]. As a result of the recent increase in surface temperature, earlier spring phenological events have been investigated in the northern hemisphere [1], [2], [4]. Shifts in vegetation phenology will affect ecosystem functions since they affect hydrological cycle, surface energy balance, and terrestrial carbon cycle [56]. Therefore, it is important to improve our understanding to quantify accurately the timing of phenological shifts for improving our understanding of terrestrial ecosystems response to global change.

In the past few decades, satellite remote sensing has been recognized as a valuable tool for acquiring and analyzing information about land surface development. In particular, time series of remote sensing data are an important data source for investigating vegetation phenological changes. The seasonal dynamics in a long-term vegetation index (VI) time series such as those from the AVHRR, SPOT, or MODIS sensors can be used for plant phenological detection. For example, time series of NDVI data derived from AVHRR or MODIS data have been used to estimate the inter-annual variations in vegetation phenology for the past decades [78]. Annual time series of Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) data have been used to identify key phenological events [9]. Recent, research efforts have been focused towards improving algorithms for calculating vegetation phenological shifts using other satellite-derived data (e.g. MODIS-enhanced vegetation index (EVI), leaf area index (LAI), FAPAR and Albedo) and modeling techniques.

Many previous studies investigated the timing of phenological events and their connection to climate change by using time-series of remote sensing data in the northern hemisphere [2],[1015]. Due to the increase in land surface temperature, strong vegetation phenology shifts have occurred in the northern hemisphere. For example, in North America, a delayed dormancy onset date (0.551 days yr-1) and an extended growing season length (0.683 days yr-1) were observed over the mid and high latitudes during 1982–2006 [16]. In China, the beginning of the growing season (BGS) has advanced in spring by 0.79 days yr-1 and the end of the growing season (EGS) delayed in autumn by 0.37 days yr-1 in temperate regions during 1982–1999 [17]. In Europe, an earlier BGS (0.54 days yr-1) and a prolonged length of vegetation growing season (LGS) (0.96 days yr-1) were found during 1982–2001 [18]. For the Northern Hemisphere as a whole, BGS advanced by 5.2 days during 1982–1999 but advance by 0.2 days during 2000–2008, and EGS was delayed by 4.3 days in the period of 1982–1999 and it was further delayed by another 2.3 days in the period of 2000–2008 [2]. However, although it took a great deal of effort to examine phenological shifts over the northern hemisphere, great uncertainties still exist. Just as our literature review revealed that many studies present different magnitudes of phenological shifts, including earlier or later phenology events and longer or shorter vegetation growing seasons [13], [1921]. Furthermore, some other studies only concentrated on earlier BGS or assessed later shifts of dormancy events [9], [2224].

In this paper, we intend to investigate the spatial variability and temporal trends of vegetation phenology shifts in the northern hemisphere during 1982–2012, and further examine the mechanisms of phenology response to climate change within different ecological zones. Our objectives in this study are to better understand the following key scientific questions: (1) What is the vegetation phenology pattern on the pixel level in the northern hemisphere? (2) How did the patterns of phenology shift within different ecological zones over the past three decades? And (3) what are the key factors affecting vegetation growing season change in recent years? Answers to these questions are essential to improve our understanding of phenological shifts caused by climate warming and to formulate measures for climate warming and environmental problems at the global level.

Materials and Methods

Data Source and Pre-Processing

NDVI dataset from GIMMS.

The NDVI3g data with a spatial resolution of 8 km×8 km and a temporal resolution of a 15-day interval were obtained from the NASA Global Inventory Modeling and Mapping Studies (GIMMS) group. Data were derived from the AVHRR instrument onboard the NOAA satellite series 7, 9, 11, 14, 16, and 17 for the time period of January 1982 to December 2011, and these data had been corrected for calibration, solar geometry, heavy aerosols, clouds and other effects not related to vegetation change [2527].

NDVI dataset from MODIS.

The Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra and Aqua satellites provide near-daily repeated coverage of the earth’s surface with 36 spectral bands and a swath width of approximately 2330 km. For this study we selected the MOD 13C1 for global vegetation information extraction. The MODIS products are distributed through NASA’s Earth Observing System Data and Information System (EOSDIS). The MODIS product (MOD 13C1), with a 0.05°spaital resolution and 16-d intervals from 2002~2012, was generated from atmospherically corrected by using bi-directional surface reflectance functions, and masked for water, ices, clouds, aerosols, and shadows [2829].

Holdridge life eco-zones data.

The Holdridge life eco-zone is one of the most widespread, quantitative schemes for classification of vegetation, land use, and environments.The concept of life zones was first published by Leslie Holdridge in 1947, and can be used to determine land use options and ecological types. The Holdridge life eco-zones data set was acquired from the Food and Agriculture Organization of the United Nations (FAO). The data set shows the Holdridge life eco-zones of the world, a combination of climate and vegetation types. The life eco-zones were devised using three indicators: mean annual biotemperature (based on the growing season length and temperature), mean annual precipitation (including rain and snow), and a potential evapotranspiration ratio [3031]. The data set has a spatial resolution of one-half degree latitude/longitude, and a total of 20 life eco-zones, ranging from the evergreen tropical rainforest zone to the boreal tundra woodland zone (Fig 1).

thumbnail
Fig 1.

(a) Spatial pattern of the Holdridge life eco-zones in the Northern Hemisphere. (b) Distribution of vegetation types in the Northern Hemisphere (adapted from UMD Global Land Cover Classification, 1 Kilometer, Department of Geography, University of Maryland).

https://doi.org/10.1371/journal.pone.0157134.g001

Global meteorological data set.

The global meteorological data set was acquired from the National Centers for Environmental Prediction/National Center for Atmospheric Sciences (NCEP/NCAR). The NCEP/NCAR reanalysis meteorological data have been widely used to study and evaluate vegetation phenology. The NCEP/NCAR Reanalysis project incorporates observations and numerical weather prediction (NWP) model output from 1948 to present. Results are available at 6 hour intervals (4-times daily), and two variables were acquired for this study: mean temperature and mean rainfall. Time series are composed of 10-day global images at 0.5 degree resolution.

Global vegetation cover data Set.

The global vegetation cover characteristics data were acquired from Department of Geography, the University of Maryland [3233]. The global vegetation classification data set, with 1 kilometer pixel resolution, was generated from 1-km AVHRR imagery data. Imagery from the AHHRR satellites acquired between 1981 and 2000 were analyzed to distinguish the following thirteen land cover types: evergreen needleleaf forest, evergreen broadleaf forest, deciduous needleleaf forest, deciduous broadleaf forest, mixed forest, woodland, wooded grassland, closed shurbland, open shrubland, grassland, crop land, bare ground, urban and built (Fig 1).

Preprocessing of NDVI time-series data.

GIMMS NDVI3g and MODIS NDVI datasets were processed to have uniform spatial projection, and to reduce the cloud and atmospheric effects before using them. The maximum value composition (MVC) was applied to the MODIS NDVI data and the AVHRR NDVI time series for cloud filtering. However, there was still some cloud contamination existing in some cases. To eliminate errors caused by clouds, snow and ice contamination, we applied a Savitzky-Golay filtering procedure to each annual NDVI cycle for smoothing and reconstrucing the NDVI time series, as described by Chen et al. [34] and Broich et al. [35]. For MODIS MOD13C1 product, a summary pixel reliability layer has been included in the MOD13C1, we first use the reliability layer information to locate the cloud or ice/snow pixels (where the pixel value equals 2 or 3 in the reliability layer) in the NDVI time-series. Then these cloud or ice/snow pixels can be replaced by a linear interpolation method using the adjacent NDVI data. The new NDVI values can be calculated as follows: (1) Where PR (Pixel Reliability) is the value of pixel reliability and NDVI(i, t) is the value of NDVI of the ith pixel in the tth time.

To further minimize the influence of soil and non-vegetated signals in NDVI time series, we chose targeted study regions with multiyear average NDVI greater than 0.1 in the Northern Hemisphere based on GIMMS (1982–2006) and MODIS (2002–2012) NDVI data.

Determination of Phenological Parameters

Several methods have been developed to detect the phenophases (greenup, maturity, senescence, and dormancy) of vegetation phenology by using NDVI data, such as NDVI thresholds [7, 36], largest NDVI increase [37], backward-looking moving averages [38], and fitting logistic functions [3839]. Comparative studies have shown that the Midpointpixel threshold method based on variations of an NDVIratio is the one of the most consistent methods for the estimation of the BGS, EGS, LGS and peak of growing season (PGS) over a variety of ecosystems and corroborates with ground-based phenology data [7]. Thus, for the determination of BGS, EGS based on NDVI, we used the Midpointpixel threshold method, where the value midpoint was set at fifty percent of the maximum NDVI [7]. Firstly we calculated the annual NDVI ratio time series curve during 1982–2012. Subsequently, the phenology date was determined for each pixel by using the midpoint pixel method, and BGS or EGS was determined according to the slope of the fitted curve. Fig 2 shows the general scheme of the method for determining the phenological parameters. The NDVI ratio was calculated by using the following formula: (2) Where NDVIratio is the output ratio, ranging from 0–1, NDVI is the daily NDVI, NDVImax is the annual maximum NDVI, and NDVImin is the annual minimum NDVI. Thus the phenological parameters can be extracted with a Midpointpixel thresholds method for each pixel.

thumbnail
Fig 2. Schematic figures illustrating method identifying the beginning of growing season (BGS), end of growing season (EGS), length of growing season (LGS) and peak of growing season (PGS).

https://doi.org/10.1371/journal.pone.0157134.g002

Phenology Trend Analysis

The spatial distribution of the phenology shift trends in the Northern Hemisphere from 1982–2012 based on the GIMMS and MODIS NDVI datasets was characterized by the change rate of phenology shift (RPS), which was calculated by the following relationship: (3) Where i is the year for 1 to n; n is the total number of years; and Dayi is the seasonal phenology (BGS, EGS, LGS and PGS) of year i. A positive value of RPS indicates an increase of phenology in a time series. The timings of these parameters were calculated to the day because we took the day of the satellite observation into account.

Correlation coefficients and significance tests were used to examine the ralationship betweeen climate and phenology shift over the period of analysis, 1982–2012. We analysed the effects of seasonal temperature and precipitation on the trends of phenology shift in different Holdridge life ecozones by Pearson correlation analysis. In addition, a t-test determined the significance of the correlation, with p values <0.1, 0.05 and 0.01 are considered to be low, middle and high significance levels, respectively.

Results and Analyses

Spatial Patterns of Vegetation Phenological Metrics

Fig 3 shows the spatial distributions of the phenological metrics BGS, EGS, LGS and PGS over the northern hemisphere. The BGS dates ranged widely from approximately 15 days (15 January) at low latitudes to approximately 225 days (15 July) at high latitudes. The mean BGS date was progressively delayed with increasing latitude. Concurrently, the progressive patterns that vegetation turn green gradually along latitudinal gradient are mainly attributed to the climatological temperature [2], [40]. The BGS in Eurasia was much earlier than in North America at the same latitude, and the latest dates of BGS occured in northern Canada and northern Siberia owing to low temperatures. The EGS dates ranged from 15 May to approximately 30 November, the LGS value ranged from 64 days to 240 days, and the PGS value ranged from 64 days to approximately 240 days. There were opposite spatial patterns between BGS and EGS over the northern hemipshere. The dates of vegetation senescence occured in reverse order of the green-up onset. Compared with other regions in the hemisphere, east Asia, Europe and southern North America had earlier BGS and later EGS, and thus had longer growing seasons, with some regions lasting up to 10 months. In contrast, in northern Canada and northern Siberia the vegetation growing seasons lasted only 1 month.

thumbnail
Fig 3. Spatial distributions of the mean (a) BGS, (b) EGS, (c) LGS and (d) PGS in the northern hemisphere during 2002–2012.

The results were calculated by 11-year MODIS NDVI data.

https://doi.org/10.1371/journal.pone.0157134.g003

Spatial and Temporal Variations of Phenological Metrics over the Northern Hemisphere

The spatial pattern of the changes of BGS, EGS, LGS and PGS for periods of 1982–1992, 1992–2002 and 2002–2012 are mapped in Fig 4. There were clear temporal shifts of phenology metrics (BGS, EGS, LGS and PGS) from 1982 to 2012. During the entire period of 1982–2012, many high-latitude regions showed a trend towards earlier BGS and delayer EGS, resulting in a longer LGS, but trends in PGS show a more fragmented pattern than BGS, EGS and LGS. During 1982–1992, more than 80% of earlier BGS (over 0.5 days yr-1) mainly occurred in northern America and Eurasia, with exceptions being parts of Russia, which showed a delayed BGS (Fig 4a). A slight pattern change occurred in EGS, with more than 65% of the study region experiencing an advancing trend in EGS (Fig 4b). A notable change occurred in LGS, which became about 80% longer (over 0.5 days yr-1) in northern America and Eurasia (Fig 4c), which was similar to BGS trend. Results indicated that a longer LGS trend is mainly caused by a significant advance in spring phenology during 1982–1992, which is consistend with previous studies [17], [41]. During 1992–2002, BGS advanced throughout most of the hemisphere, with exceptions being parts of central northern America and southern Asia, while EGS was delayed throughout almost the entire hemisphere, resulting in a longer LGS (Fig 4e). During 2002–2012, a marked advancing trend (over 0.8 days yr-1) in BGS occurred form 55°N to 75°N (Fig 4i), correspondingly, the LGS can also be found to be siginificantly prolonged in this region (Fig 4m). From 1982 to 2012, unlike BGS, EGS and LGS, PGS showed a more fragmented pattern, but overall there was not a significant trend in PGS. From the point of spatial patterns, PGS was more relvant with BGS than EGS or LGS (Fig 4d, 4h and 4m).

thumbnail
Fig 4. Spatial distribution of the changes in dates of BGS, EGS, LGS and PGS over the past three decades in the Northern Hemisphere.

Positive values (red colors) represent later onset (BGS), later finish (EGS), longer duration (LGS) and later peak of the growing season. (a), (e), (i): changes in dates of vegetation green-up for 1982–1992, 1992–2002 and 2002–2012, respectively (days yr-1). (b), (f), (k): changes in dates of vegetation senescence for 1982–1992, 1992–2002 and 2002–2012, respectively (days yr-1). (c), (g), (l): changes in dates of vegetation growing season length for 1982–1992, 1992–2002 and 2002–2012, respectively (days yr-1). (d), (h), (m): changes in dates of vegetation growing peak for 1982–1992, 1992–2002 and 2002–2012, respectively (days yr-1).

https://doi.org/10.1371/journal.pone.0157134.g004

Fig 5 shows the overall trends of phenological metrics changes during 1982–2012 and interannual variability of vegetation phenology along the latitudinal gradient. As seen in Fig 5a and Table 1, changes in BGS, EGS, LGS and PGS showed high interannual variations. Long-term variations in BGS, EGS and LGS showed ‘earlier green-up onset’, ‘delayed dormancy onset’, and ‘prolonged growing season’ features, which are consistent with previous studies [2],[10],[17]. Only using NOAA NDVI data, linear regression between phenological timing and time indicates that the rate of change is -0.19 days per year (R2 = 0.33, P < 0.001) for BGS, 0.21 days per year (R2 = 0.41, P < 0.001) for EGS, 0.32 days per year (R2 = 0.44, P < 0.001) for LGS and 0.01 days per year (R2 = 0.00, P < 0.001) for PGS for the period of 1982–2002 over the Northern Hemisphere. For using MODIS NDVI data, the rate of change is -0.50 days per year (R2 = 0.25, P < 0.001) for BGS, 0.35 days per year (R2 = 0.16, P < 0.001) for EGS, 0.59 days per year (R2 = 0.38, P < 0.001) for LGS and 0.05 days per year (R2 = 0.01, P < 0.001) for PGS for the period of 2002–2012. As seen in Fig 5b, the vegetation phenology was strongly dependent on latitude, which indicates temperature is the primary controlling factor on vegetation phenology. However, it should be noted that phenological metrics have large variations along the longitudinal gradient (see the interannual variability in Fig 5b), which reflects the vegetation phenology response to global warming for different latitude. Linear regression between phenological timing and latitude indicates that the rate of change is 1.24 days per degree of latitude (R2 = 0.96, P < 0.001) for BGS, -0.51 days (R2 = 0.57, P < 0.001) for EGS, -2.02 days (R2 = 0.97, P < 0.001) for LGS and 0.54 days (R2 = 0.76, P < 0.001) for PGS in Northern Hemisphere between 30° N and 75° N.

thumbnail
Fig 5. Overall trends of phenological transitional dates (DOY, day of year) along the longitudinal gradient in Northern Hemisphere during 1982–2012.

(a) overall trends in dates of vegetation phenology for 1982–2012. (b) variability of vegetation phenology along the latitudinal gradient from 1982 to 2002.The dark solid line indicates the mean and the shaded zone indicates the range of interannual variability.

https://doi.org/10.1371/journal.pone.0157134.g005

thumbnail
Table 1. Shifts of phenology metrics based on NOAA AVHRR data (1982–2002) and MODIS data (2002–2012) (days per decade).

https://doi.org/10.1371/journal.pone.0157134.t001

The variation in mean BGS, EGS and PGS for specific vegetation classes were estimated by using AVHRR data and MODIS data. Fig 6 provides a temporal depiction of the variation in BGS, EGS, LGS and PGS for specific vegetation classes during 1982–2012. It was found that deciduous broadleaf forest and woodland exhibited a higher rate of change for both BGS (−0.24 and −0.25 days yr−1, respectively) and EGS (+0.48 and +0.33 days yr−1, respectively) than the other classes. Shrubland showed the least variation in BGS (-0.08 days yr−1) and EGS (0.05 days yr−1) compared to the other classes. Interestingly, including deciduous broadleaf forest, woodland and deciduous needleleaf forest, the delay in EGS is consistently larger than the advancement in BGS,which leads to the inference that autumnal response of vegetation is greater and faster than that in spring, which is consistent with the findings of previous studies [42]. All these specific vegetation classes revealed a prolonging growth cycle with earlier BGS and delayed EGS.

thumbnail
Fig 6. Interannual variations of (a) BGS, (b) EGS, (c) LGS and (d) PGS for specific vegetation classes over Northern Hemisphere based on AVHRR data (back lines, 1982–2006) and MODIS data (red lines, 2002–2012).

https://doi.org/10.1371/journal.pone.0157134.g006

Variation in Phenology within Different Holdridge Life Eco-Zones

In order to examine the phenological trends and their characteristic differences, the index RPS was calculated for the Holdridge life eco-zones in the study area. The index RPS was divided arbitrarily into four groups based on the value of the change rate of phenology shift. The first group represented the slight mode of phenological shift, in which the value of RPS was from 0.2 to 0.5 or from -0.5 to -0.2. The second group represented the moderate mode where the RPS was from 0.5 to 0.8 or from -0.8 to -0.5. The third group belonged to the fast mode with the RPS from 0.8 to 1.0 or from -1.0 to -0.8, and the group four was for the acute mode with the vulue of RPS above 1.0 or below -1.0.

Fig 7 indicated that during 1982–1992, for BGS, most of the eco-zones, including boreal coniferous forest (Ba), boreal mountain system (BM), boreal tundra woodland (Bb), temperate continental forest (TeDc) and temperate steppe (TeBSK), were in the first RPS group (Fig 7a, a light blue area shows a slight advance trend in spring phenology). For EGS, no significant change occurred in many eco-zones (Fig 7b, yellow areas). Due to most of the eco-zones experienced an advancing trend in BGS, resulting in a longer growing season length. So for LGS, boreal coniferous forest belonged to the group one, boreal mountain system and temperate steppe belonged to the group two (Fig 7c, light red areas), all eco-zones experienced an extend trend in LGS. During 1992–2002, the same spatial pattern of RPS as the period of 1982–1992, but we should note that a remarkable extend occurred in the vegetation growing length, due to the ‘earlier green-up onset’, ‘delayed dormancy onset’ features in eco-zones. During 2002–2012, compared with the period of 1982–1992 and 1992–2002, the spatial pattern of RPS showed a more significant trend in advancing of BGS, delaying of EGS and prolonging of LGS during 2002–2012, overall shows a trend of accelerating change. Such as boreal coniferous forest eco-zone, belonged to the group three, group two and group four for BGS, EGS and LGS, respectively (Fig 7g, 7h and 7i). It should be noted that in some eco-zones, such as polar, temperate desert and tropical desert, due to poor vegetation cover and snow (ice) and soil contamination, the spatial pattern of RPS shows a meaningless trend for the whole period of 1982–2012.

thumbnail
Fig 7. Spatial distribution of the change rate of phenology shift (RPS) within different Holdridge life eco-zones during the growing season from 1982 to 2012 in the Northern Hemisphere.

Positive values (red colors) represent later onset (BGS), later finish (EGS), longer duration (LGS) of phenology in a time series. (a), (d), (g): RPS of BGS for 1982–1992, 1992–2002 and 2002–2012, respectively. (b), (e), (h): RPS of EGS for 1982–1992, 1992–2002 and 2002–2012, respectively. (c), (f), (i): RPS of LGS for 1982–1992, 1992–2002 and 2002–2012, respectively.

https://doi.org/10.1371/journal.pone.0157134.g007

Analysis of the Relationship “Phenology-Climate” within Different Eco-Zones

In order to investigate the relationship between vegetation phenology and climate factors, correlation coefficients were calculated over Holdridge life eco-zones from 1982 to 2012 based on NOAA-derived phenology, MODIS-derived phenology and NCEP/NCAR reanalysis meteorological data. As we previously mentioned, in the Northern Hemisphere, phenology is primarily determined by temperature in the months period preceding the event (e.g., bud-burst, flowering, and leaf out), and higher pre-season temperatures may advance the BGS [2], [11], [17], [43], [44]. In this study, we determining the correlation analyses between the BGS and spring temperature (March, April and May), EGS and autumn temperature (August, September and October), BGS and precipitation, EGS and precipitation in the main ecozones of phenological changes. Fig 8 presents the statistical distributions of the correlation coefficients derived from Pearson correlation analyses between the interannual variations of BGS, EGS and of the climatic data (spring temperature, autumn temperature and mean annual precipitation). The presentation of the pixels were restricted to those eco-zones (Ba, Bb, BM, TeBSK, TeDc) only where the RPS exhibited a significant trend in the period 1982–2012. About 75% of all pixels showed a negative correlation with significance at p<0.05 level between the interannual variations of BGS and of the spring temperature, with Pearson’s correlation coefficient ranging from -0.97 to -0.20, and histogram peak ranging from 0.4 to 0.6, which indicates that there is a strong negative correlation between BGS and spring temperature. About 45% of all pixels exhibited a positive correlation with significance at p<0.05 level between EGS and autumn temperature, with Pearson’s correlation coefficient ranging from 0.2 to 0.95, and histogram peak ranging from 0.1 to 0.4, which indicates that there is a strong positive correlation between EGS and autumn temperature.This is consistent with previous studies that the earlier BGS and delayer EGS are the function of mean or accumulated temperatures over the Northern Hemisphere [2], [11], [45], [46]. Like the distribution pattern of the correlation coefficients between BGS, EGS and temperature, about 65% of all pixels showed a negative correlation between BGS and mean annual precipitation, with Pearson’s correlation coefficient ranging from 0.2 to 0.90, and histogram peak ranging from 0.2 to 0.4. However, overall EGS and mean annual precipitation did not exhibit a significant correlation. It should be noted that comparing the correlation coefficient of vegetation phenology and temperature with that of phenology and precipitation, we can see the correlation between temperature and phenology is much higher than that between precipitation and phenology, displaying that temperature is a major determinant of changes in vegetation phenology. This is consistent with previous studies reporting the cycles of vegetation in temperate zones are primarily controlled by temperature [2], [17], [44], [4748].

thumbnail
Fig 8. Statistical distribution of Correlation coefficients of the interannual variations of BGS, EGS and of the climatic data (spring temperature, autumn temperature and mean annual precipitation) in the Northern Hemisphere during 1982–2012.

https://doi.org/10.1371/journal.pone.0157134.g008

Even though temperature is considered as the major determinant for vegetation phenological change, little is known on vegetation phenological responses to temperature changes [44], [4950]. Moreover, what is the phenological change per unit temperature. To further examine the temperature-sensitivity of the BGS and EGS, we analyzed the statistical distribution pattern of BGS and EGS changes in relation to mean annnual temperature within the different Holdridge life eco-zones (Fig 9). In those eco-zones whose RPS showed a significant trend during 1982–2012, the vegetation phenology was strongly dependent on temperature, in general, which indicates a marked advancing trend in BGS and a delaying trend in EGS along increasing temperature gradient. Interestingly, in all eco-zones, the BGS and EGS exhibited the largest variation in the temperature range of 5–12°C, which shows the vegetation phenological response to global warming is more sensitive in this temperature zone. Furthermore, in Ba, TeDC and TeBSK eco-zones, the advancement in BGS was larger than the delay in EGS, which leads to the inference that vernal vegetation response to temperature is greater and faster than that in autumn in these eco-zones. In the BM eco-zone, by contrast, the delay in EGS is consistently larger than the advancement in BGS, which also leads to the inference that autumnal vegetation response to temperature is greater and faster than that in spring in this zone. These inferences require detailed ground verification.

thumbnail
Fig 9. Distribution of BGS, EGS changes in relation to the mean annual temperature (upper), and the frequency of pixel in relation to mean annual temperature (bottom) within the different Holdridge life eco-zones during 1982–2012.

We used the 30-year mean annual temperature to represent the climatic temperature condition. Only pixels that BGS and EGS related to temperature for statistically significant at 90% confidence level are included. (Holdridge life eco-zones including boreal coniferous forest (Ba), boreal tundra woodland (Bb), boreal mountain system (BM), temperate continental forest (TeDc), and temperate steppe (TeBSK)).

https://doi.org/10.1371/journal.pone.0157134.g009

We next examined the vegetation phenological changes associated with temperature change within the different Holdridge life eco-zones (Fig 10). We found in all cases, higher mean temperature change associated with earlier mean BGS, such as during 2002–2012, the two largest advancement in BGS were found within the Ba (−0.84 days yr−1) and Bb (−0.55 days yr−1) zones, whose standard deviations of temperature change were 0.58 and 0.56, respectively. During 1992–2002, the largest advancement in BGS was found within the TeDc zone (−0.46 days yr−1), whose standard deviation of temperature change was 0.53.During 1982–1992, the largest advancement in BGS was found within the TeDc zone (−0.71 days yr−1), whose standard deviation of temperature change was 0.52. This is consistent with the results of previous studies showing that higher temperature-sensitivity associated with earlier onset of spring phenological events [44], [5152]. Furthermore, in all eco-zones, overall vegetation phenology revealed a longer phenological cycle with earlier BGS and delayed EGS, which represents an extension of the growth cycle. Specifically, when comparing BGS with EGS, although the advanced BGS and delayed EGS have been reported to be the key factors that cause prolonged vegetation growing seasons, but in our study, the earlier BGS was found to be a more important factor that causes the LGS changes during 1982–2012, which differs from the results of Jeong et al.[2].

thumbnail
Fig 10. Sensitivity of vegetation phenological changes to temperature change within the different Holdridge life eco-zones during 1982–2012.

The 30-year mean annual temperature was used to represent the climatic temperature condition. Only pixels that BGS and EGS related to temperature for statistically significant at 90% confidence level are included.

https://doi.org/10.1371/journal.pone.0157134.g010

Discussion

Uncertainties in Remote Sensing Estimates on Phenological Metrics

Satellite-based techniques have been widely used for phenological monitoring and eco-environmental evaluation. However, there were many uncertainties in the results due to the temporal, spatial, and ecological complexity of vegetation biochemical processes. In this study, the uncertainties in the phenological results are twofold: uncertainties related to the different datasets and uncertainties related to the methods for determing phenological metrics.

  1. Uncertainties related to the different datasets. In this paper, the uncertainties related to datasets include five ways: satellite drift, sensor differences, calibration uncertainties, atmopsheric effects and sub-pixel clouds [10], [2527]. For example, sensor differences, such as NOAA AVHRR sensor or MODIS sensor, with different sensor characteristics and spatial resolution (e.g., AVHRR and MODIS have different bandwidths and central wavelengths, different number of bands and radiometric resolution; and MODIS sensor has more better onboard geometric and radiometric calibrations than AVHRR sensor), can affect the quantification of phenology shifts. Some previous studies also show that the phenology paremeters (e.g., BGS, EGS) from the different datasets were significantly different over the nothern high lattitudes [20], [2324], [53]. Furthermore, some studies have found that the newer generation sensors such as MODIS and SPOT performed better than AVHRR for investigating the phenology in alpine vegetation [20], [24], [54]. We found similar results with AVHRR and MODIS sensor, for instance, the BGS and EGS values appear big differences in both satellite sensors (see Figs 6 and 7), but the trends of BGS and EGS show generally similar spatio-temporal patterns and good agreement between the two sensors. To reduce the uncertainty of the results, in this paper the whole study period was divided into three periods (1982–1992, 1992–2002 and 2002–2012, respectively), in each period we only use the same dataset for trend analysis, thus these disagreements between AVHRR and MODIS sensors did not influence the results reported in this paper.
  2. Uncertainties related to the methods for determing phenological metrics. On a hemispherical scale, satellite-based studies have been demonstrated to be useful for phenological monitoring. However, there are many uncertainties in the results when determing phenological metrics due to the temporal, spatial, and ecological complexity of phenological processes. Over the past few years, several models have been developed for detecting the phenological metrics by using satellite data, including NDVI thresholds, largest NDVI increase, backward-looking moving averages, and fitting logistic functions, but each model has its uncertainty for estimating the phenological metrics. Overall, comparative studies have shown that the Midpointpixel threshold model based on variations of an NDVIratio is the one of the most consistent methods for determing phenological metrics over a variety of ecosystems and corroborates with field-based observational data [7], [55]. Thus, in order to determine the global phenological metrics based on AVHRR and MODIS data, we used the Midpointpixel threshold method in this study. Furthermore, in response to concerns over the validity of satellite derived phenology, our research adopted a strict approach to the selection of pixels for analysis through statistical tests of significance at different levels.

Influencing Factors of Phenological Shifts

Many factors affect phenological shift patterns, such as topographic conditions, meteorological conditions, human disturbance and structural geology, among other factors. As far the global vegetation phenology, as reported by the IPCC’s AR4 [56], many natural environmental factors strongly influence phenological shifts, such as temperature, precipitation/soil moisture, photoperiod, winter chilling. As we previously mentioned, in the Northern Hemisphere, phenology is primarily determined by temperature, specifically, alpine plant species are considered to be the most sensitive vegetation ecosystem in response to rising temperature. It should be noted that different plant species have different temperature requiements for development, also different plant species have different responses to temperature change. For instance, Jin et al. [57] reported that on the Tibetan Plateau, 1°C increase in soil temperature will advance the BGS by 4.6 to 9.9 days, and postpone the EGS by 7.3 to 10.5 days among the different plant species. Meanwhile, our study also indicates that there were different plant phenological responses to temperature within different eco-zones. In particular, our study further revealed that higher temperature warming related to earlier onset of BGS.

Even though temperature clearly played the most important role in vegetation phenology at the high latitudes, precipitation/soil moisture was also an important factor for vegetation phenology. Especially in arid and semi-arid regions, precipitation/moisture was a more important factor in vegetation phenology. For example, Prieto et al. [58] found that flowering time was strongly influenced by precipitation patterns in Mediterranean shurbs. Lesica et al. [59], Broich et al. [60] and Krishnaswamy et al. [15] also reported that precipitation and moisture may play a more important role in phenological cycles in arid and semi-arid systems. In the high-latitude regions snowmelt is also an important factor affecting alpine vegetation phenology, especially in some cold regions, the time of snowmelt is often considered to be the observable indicators of phenological events in tundra plants [61]. Some recent studies have found that the spatial pattern of vegetation growing season mainly corresponded to the snowmelt gradient and plant species in the alpine ecosystem [6263]. Jeganathan et al.[42] also found some link between the trend in the onset of snowmelt and inter-annual variation in BGS and EGS. Interestingly, in some snow-rich regions, such as the Kola Peninsula in Russia, and outer Troms County in northern Norway, mountain birch and herbaceous plants may have green leaves before snowmelt in spring when there is enough light through the snow [6465]. Overall, in some high-latitude alpine areas, timing of snowmelt and temperature are the two most important factors controlling phenology.

In addition to temperature and precipitation, the photoperiod and winter chilling also have an influence on vegetation life cycle at the high latitudes. Although it is commonly assumed that winter chilling plays a important role in regulation of phenological phases, in fact, the interactions between chilling and photoperiod might play a more important role in controlling phenological phases at the high latitudes. For example, some recent experiments suggest that phtotoperiod might play an important role in controlling budburst and scenescence [66]. Some previous studies found that warmer temperatures during the winter will slow the dormancy breaking process by delaying the fulfillment of the chilling requirement. However, warm temperatures during the spring will advance the plant life growth by accelerating the accumulation of heat. [45], [52], [67]. Especially for humid extratropical areas, phenology is more controlled by winter chilling, photoperiod and temperature [6869].

Conclusions

During the period from 1982 to 2012, most of the high-latitude regions experienced significant trends in advancing of BGS, delaying of EGS and prolonging of LGS, but no significant trends in PGS. Comparing the period of 2002–2012 with the periods of 1982–1992 and 1992–2002, the spatial pattern of RPS shows a more significant trend in advancing of BGS, delaying of EGS and prolonging of LGS during 2002–2012, overall shows a trend of accelerating change. The correlation coefficients between climate and phenological events for the Holdridge life eco-zones from 1982 to 2012 showed that temperature was a major determinant of changes in vegetation phenology, with a significant advanced trend in BGS and a delayed trend in EGS along an increased temperature gradient. However, the phenological response of vegetation to temperature varied across different eco-zones. This study has shown that Midpointpixel threshold method is useful, and the trends in this study are generally in agreement with trends in previous studies. However, there are considerable uncertainties associated with the data and methods. To reduce uncertainies, it is clearly needed to design a biogeochemical model that integrates satallite-derived parameters, plant types, natural environmental factors and ground-based observations in future research.

Acknowledgments

The authors wish to thank the NASA and Global Inventory Modeling and Mapping Studies (GIMMS) group for making the 25 years of AVHRR NDVI3g data freely available to researchers, and the authors are very grateful to Prof. Hanqin Tian and the anonymous reviewers for their constructive and critical comments on this manuscript.

Author Contributions

Conceived and designed the experiments: SW BY. Performed the experiments: SW BY QY XW. Analyzed the data: SW BY QY XW. Contributed reagents/materials/analysis tools: SW BY QY XW. Wrote the paper: SW BY QY LL YP.

References

  1. 1. Piao SL, Friedlingstein P, Ciais P, Viovy N, Demarty J. Growing season extension and its effects on terrestrial carbon flux over the last two decades. Global Biogeochemical Cycles 2007; 21, GB3018,
  2. 2. Jeong SJ, Ho CH, GIM HJ, Brown M. Phenology shifts at start vs end of growing season in temperate vegetation over the Northern Hemisphere for the period 1982–2008. Global Change Biology 2011; 17: 2385–2399.
  3. 3. IPCC (Ed.) Climate change 2007: Impacts, adaptation and vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK Cambridge University Press 2007.
  4. 4. Geerken Roland A.. An algorithm to classify and monitor seasonal variations in vegetation phenologies and their inter-annual change. ISPRS Journal of Photogrammetry and Remote Sensing 2009; 64: 422–431.
  5. 5. Richardson AD, Black TA, Ciais P, Delbart N, Friedl MA, Gobron N, et al. Influence of spring and autumn phenological transitions on forest ecosystem productivity. Phil. Trans. R. Soc. B 2010; 365: 3227–3246. pmid:20819815
  6. 6. Steven R, Garrity G, Bohrer KD, Maurer KL, Mueller CS, Vogel , et al. A comparison of multiple phenology data sources for estimating seasonal transitions in deciduous forest carbon exchange. Agricultural and Forest Meteorology 2011; 151: 1741–1752.
  7. 7. White MA, Thornton PE, Running SW. A continental phenology model for monitoring vegetation responses to interannual climatic variability. Global Biogeochemical Cycles 1997; 11(2), 217–234.
  8. 8. Ivits E, Cherlet M, Sommer S, Mehl W. Addressing the complexity in non-linear evolution of vegetation phonological change with time-series of remote sensing images. Ecological Indicators 2013; 26: 49–60.
  9. 9. Verstrate MM, Gobron N, Aussedat O. An automatic procedure to identify key vegetation phenology events using the JRC-FARAR products. Advances in Space Research 2008; 41: 1773–1783.
  10. 10. Zhou L, Tucker CJ, Kaufmann RK, Slayback D, Shabanov NV, Myneni RB. Variations in northern vegetation activity inferred from satellite data of vegetation index during 1981–1999. Journal of Geophysical Research 2001; 106, 20069–20083.
  11. 11. De Beurs KM, Henebry GM. Land surface phenology and temperature variation in the International Geosphere-Biosphere Program high-latitude transects. Global Change Biology 2005; 11: 779–790.
  12. 12. De Beurs KM, Henebry GM. Northern annular mode effects on the land surface phenologies of northern Eurasia. Journal of Climate 2008; 21: 4257–4279.
  13. 13. Stockli R, Rutishauser T, Baker I, Liniger MA, Denning AS. A global reanalysis of vegetation phenology. J.Geophys.Res. 2011; 116, G03020,
  14. 14. Shrestha UB, Gautam S, Bawa KS. Widespread Climate Change in the Himalayas and Associated Changes in Local Ecosystems. PLoS ONE 2012; 7(5): e36741. pmid:22615804
  15. 15. Krishnaswamy J, John R, Joseph S. Consistent response of vegetation dynamics to recent climate change in tropical mountain regions. Global Change Biology 2014; 20(1): 203–215. pmid:23966269
  16. 16. Zhu WQ, Tian HQ, Xu XF, Pan YZ, Chen GS, Lin WP. Extension of the growing season due to delayed autumn over mid and high latitudes in North America during 1982–2006. Global Ecology and Biogeography 2011; 21(2): 260–271.
  17. 17. Piao SL, Fang JY, Zhou LM, Ciais P, Zhu B. Variations in satellite-derived phenology in China’s temperate vegetation. Global Change Biology 2006; 12: 672–685.
  18. 18. Stockli R, Vidale PL. European plant phenology and climate as seen in a 20-year AVHRR land-surface parameter dataset. INT.J. Remote Sensing 2004; 25(17): 3303–3330.
  19. 19. Cleland EE, Chuine I, Menzel A, Mooney HA. Shifting plant phenology in response to global change. Trends in Ecology and Evolution 2007; 22(7): 357–365. pmid:17478009
  20. 20. Zeng HQ, Jia GS, Epstein H. Recent changes in phenology over the northern high latitudes detected from multi-satellite data. Environmental Research Letters 2011; 6(4): 1–11.
  21. 21. Horion S, Cornet Y, Erpicum M, Tychon B. Studying interactions between climate variability and vegetation dynamic using a phenology based approach. Int. J. Appl. Earth Observ. Geoinf. 2012;
  22. 22. Goward SN, Tucker CJ, Dye DG. North American vegetation patterns observed with the NOAA advance very high resolution radiometer. Vegetation 1985; 64(1): 3–14.
  23. 23. Delbart N, Toan LT, Kergoat L, and Fedotova V. Remote sensing of spring phenology in boreal regions: A free of snow-effect method using NOAA-AVHRR and SPOT-VGT data (1982–2004). Remote Sensing of Environment 2006; 101: 52–62.
  24. 24. Zhang GL, Zhang YJ, Dong JW, Xiao XM. Green-up dates in the Tibetan Plateau have continuously advanced from 1982–2011. PNAS 2013; 110(11): 4309–4314. pmid:23440201
  25. 25. Tucker CJ, Pinzon JE, Brown ME. Global Inventory Modeling and Mapping Studies, NA94apr15b.n11-VIg, 2.0, Global Land Cover Facility, University of Maryland, College Park, Maryland, 2004.
  26. 26. Tucker CJ, Pinzon JE, Brown ME, Slayback D, Pak EW, Mahoney R, Vermote E, Saleous N. An Extended AVHRR 8-km NDVI Data Set Compatible with MODIS and SPOT Vegetation NDVI Data. International Journal of Remote Sensing 2005; 26(20): 4485–4498.
  27. 27. Pinzon J, Brown ME, Tucker CJ. Satellite time series correction of orbital drift artifacts using empirical mode decomposition. In: Huang N. (Editor), Hilbert-Huang Transform: Introduction and Applications 2005; pp. 167–186.
  28. 28. Huete A, Justice C, Leeuwen W. MODIS Vegetation Index (MOD13). Algorithm Theoretical Basis Document 1999; http://modis.gsfc.nasa.gov/data/atbd/atbd_mod13.pdf.
  29. 29. Huete A, Didan K, Miura T, Rodriguez EP, Gao X, Ferreira LG. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensing of Environment 2002; 83: 195–213.
  30. 30. Holdridge LR. Life Zone Ecology. San Jose, Costa Rica, Tropical Science Center. 1967.
  31. 31. FAO. Global Ecological Zones Mapping. Workshop Report, Cambridge, 28–30 July 1999. Rome, 2000. FAO-Forest Resources Assessment Programme, Working Paper 2000; 26.
  32. 32. Hansen M, DeFries R, Townshend JRG., Sohlberg R. UMD Global Land Cover Classification, 1 Kilometer, 1.0, Department of Geography, University of Maryland, College Park, Maryland 1998; 1981–1994.
  33. 33. Hansen M, DeFries R, Townshend JRG, Sohlberg R. Global land cover classification at 1km resolution using a decision tree classifier. International Journal of Remote Sensing 2000; 21: 1331–1365.
  34. 34. Chen J, Jonsson P, Tamura M, Gu ZH, Matsushita B, Eklundh L. A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter. Remote Sensing of Environment 2004; 91: 332–344.
  35. 35. Broich M, Huete A, Paget M, Ma XL, Tulbure M, Coupe N R, Evans B, Beringer J, Devadas R, Davies K, Held A. A spatially explicit land surface phenology data product for science, monitoring and natural resources management applications. Environmental Modelling & Software 2015; 64:191–204.
  36. 36. Lloyd D. A phenological classification of terrestrial vegetation cover using shortwave vegetation index imagery. International Journal of Remote Sensing 1990; 11: 2269–2279.
  37. 37. Kaduk J, Heimann M. A prognostic phenology model for global terrestrial carbon cycle models. Climate Research 1996; 6: 1–19.
  38. 38. Zhang XY, Friedl MA, Schaaf CB, Strahler AH, hodges JCF, Gao F, Reed BC, Huete A., Monitoring vegetation phenology using MODIS. Remote Sensing of Environment 2003; 84(3): 471–475.
  39. 39. Zhang XY, Friedl MA, Schaaf CB. Global vegetation phenology from Moderate Resolution Imaging Spectroradiometer (MODIS): evaluation of global patterns and comparison with in situ measurements. Journal of Geophysical Research 2006; 111, G04017,
  40. 40. Schwartz MD, Ahas R, Aasa A. Onset of spring starting earlier across the northern hemisphere. Global Change Biology 2006; 12: 343–351.
  41. 41. Walther GR, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, et al. Ecological responses to recent climate change. Nature 2002; 416: 389–395. pmid:11919621
  42. 42. Jeganathan C, Dash J, Atkinson PM. Remotely sensed trends in the phenology of northern high latitude terrestrial vegetation controlling for land cover change and vegetation type. Remote sensing of Environment 2014; 143: 154–170.
  43. 43. Hunter AF, Lechowicz MJ. Predicting the timing of budburst in temperate trees. Journal of Applied Ecology 1992; 29: 597–604.
  44. 44. Shen M, Tang Y, Chen J, Yang X, Wang C, Cui X, Yang Y, Han L, Li L, Du J, Zhang G, Cong N. Earlier-Season Vegetation Has Greater Temperature Sensitivity of Spring Phenology in Northern Hemisphere. PLoS ONE 2014; 9(2), e88178. pmid:24505418
  45. 45. Zhang XY, Tarpley D, Sullivan JT. Diverse responses of vegetation phenology to a warming climate. Geophysical Research Letters 2007; 34: 1–5.
  46. 46. Ivits E, Cherlet M, Sommer S, Mehl W, Vogt J, Micale F. Combing satellite derived phenology with climate data for climate change impact assessment. Global and Planetary Change 2012; 88: 85–97.
  47. 47. Barber VA, Juday GP, Finney BP. Reduced growth of Alaskan white spruce in the twentieth century from temperature-induced drought stress. Nature. 2000; 405: 668–673. pmid:10864320
  48. 48. Gong DY, and Ho C.H. Detection of large-scale climate signals in spring vegetation index (normalized difference vegetation index) over the Northern Hemisphere. Journal of Geophysical Research 2003; 108, 4066,
  49. 49. Ibanez I, Primack RB, Miller-Rushing AJ, Ellwood E, Higuchi H, Lee SD, et al. Forecasting phenology under global warming. Philosophical Transactions of the Royal Society B-Biological Sciences 2010; 365:3247–3260.
  50. 50. Vitasse Y, Francois C, Delpierre N, Dufrene E, Kremer A, Chuine I, et al. Assessing the effects of climate change on the phenology of European temperate trees. Agricultural and Forest Meteorology 2011;
  51. 51. Fu YH, Compioli M, Deckmyn G, Janssens IA. The impact of winter and spring temperatures on temperate tree budburst dates: results from an experimental climate manipulation. PLoS ONE 2012; 7, e47324. pmid:23071786
  52. 52. Cook BI, Wolkovich EM, Parmesan C. Divergent responses to spring and winter warming drive community level flowering trends. PANS 2012; 109: 9000–9005.
  53. 53. Fisher JI, Mustard JF. Cross-scalar satellite phenology from ground, Landsat, and MODIS data. Remote Sensing of Environment 2007; 109: 261–273.
  54. 54. Fontana F, Rixen C, Jonas T, Aberegg G, Wunderle S. Alpine grassland phenology as seen in AVHRR, VEGETATION, and MODIS NDVI time series—-a comparison with in situ measurements. Sensors 2008; 4: 2833–2853.
  55. 55. De Berus KM, and Henebry GM. Spatio-temporal statistical methods for modelling land surface phenology, Phenological Research: Methods for Environmental and Climate Change Analysis ed Hudson I L and Keatley M R (Dordrecht: Springer) 2010; pp 177–208.
  56. 56. Rosenzweig C, Casassa G, Karoly DJ, Imeson A, Liu C, Menzel A, et al. Assessment of observed changes and responses in natural and managed systems. In: Parry M.L., Canziani O.F., Palutikof J.P., van der Linden P.J. 2007.
  57. 57. Jin ZN, Zhuang QL, He JS, Luo TX, Shi Y. Phenology shift from 1989 to 2008 on the Tibetan Plateau: an analysis with a process-based soil physical model and remote sensing data. Climate change 2013; 119: 435–449.
  58. 58. Prieto P, Peñuelas J, Ogaya R, Estiarte M. Precipitation-dependent flowering of Globularia alypum and Erica multiflora in Mediterranean shrubland under experimental drought and warming and its interannual variability. Annals of Botany 2008; 102:275–285. pmid:18565983
  59. 59. Lesica P, Kitelson PM. Precipitation and temperature are associated with advanced flowering phenology in a semi-arid grassland. Journal of Arid Environments 2010; 74:1013–1017.
  60. 60. Broich M, Huete A. Tulbure MG Ma X Paget M Restrepo-Coupe N, et al. Land surface phenological response to decadal climate variability across Australia using satellite remote sensing. Boigeosciences 2014; 11: 5181–5198.
  61. 61. Odland A. Estimation of the growing season length in alpine areas: effects of snow and temperatures. In: Scmidt JG (ed) Alpine environment: geology, ecology and conservation. Nova Science Publication, New York 2011.
  62. 62. Buus-Hinkler J, Hansen BU, Tamstorf M, and Pedersen SB. Snow-vegetation relations in a high Arctic ecosystem: inter-annual variability inferred from new monitoring and model concepts. Remote Sensing of Environment 2006; 105:237–247.
  63. 63. Reiko I, Hiroyuki O. A cost-effective monitoring method using digital time-lapse cameras for detecting temporal and spatial variations of snowmelt and vegetation phenology in alpine ecosystems. Eclological Informatics 2013; 16: 25–34.
  64. 64. Shutova E, Wielgolaski FE, Karlsen SR, Makarova O, Berlina N, Filimonova T, et al. Growing seasons of Nordic mountain birch in northernmost Europe as indicated by long-term field studies and analyses of satellite images. Int J Biometeorol. 2006; 51: 155–166. pmid:16832653
  65. 65. Wielgolaske FE, Inouye DW. Phenology: An Integrative Environmental Science, Springer Netherlands 2013; pp.225–247.
  66. 66. Basler D, and Korner C. Photoperiod sensitivity of bud burst in 14 temperate forest tree species. Agricultural and Forest Meteorology 2012; 165(15): 73–81.
  67. 67. Yu HY, Luedeling E, Xu JC. Winter and spring warming result in delayed spring phenology on the Tibetan Plateau. PNAS 2010; 107(51): 22151–22156. pmid:21115833
  68. 68. Korner C, Basler D. Phenology under global warming. Science 2010; 327: 1461–1462. pmid:20299580
  69. 69. Vitasse Y, Basler D. What role for photoperiod in the bud burst phenology of European beech. Enr J Forest Res. 2012;