Introduction

Plant phenology has proven to be a very sensitive indicator for climate change impacts1,2. The observed change of plant phenology across the globe became notable in responses to recent climate change3,4,5,6. The change in plant spring phenological events and life cycles of animals could significantly impact the structure and function of ecosystem, e.g. interspecific relation, carbon uptake and nutrient cycling7,8. Furthermore, the change in phenology can feed back on atmosphere through physical process (e.g. surface albedo and energy partitioning) or biogeochemical circulation (e.g. carbon cycle)9,10. A better understanding of phenology—climate interactions are crucial for climate change research.

Many studies in middle and high latitudes demonstrate that the temperature is the main driving force and interannual modulator of phenological change, while other factors (e.g. photoperiod) only play a secondary role as limiting factors1. Since the temperature sensitivity of plant phenological stages (phenophases change in days per degree Celsius) determines the magnitude of phenological shifts in response to future climate warming, more attention has been paid to it, both in observational records and warming experiment studies11.

Observed evidence shows that the temperature sensitivity of plant spring phenophases exhibited interspecific and spatial pattern5. Early-season plant species often exhibit stronger temperature sensitivity or faster rate of phenological change at a fixed site or across the same region12,13,14,15,16,17. With respect to the spatial pattern, the response of spring phenophases to temperature is stronger in warmer than in colder countries across Europe4. Likewise, in China and the United States, the plant phenophases at lower latitudes (warmer sites) show stronger temperature sensitivity18,19,20,21. An interesting question here is why phenological responsiveness to temperature is greater in early-season plant species or warmer regions. According to recent studies, the interspecific difference in temperature sensitivity may be attributed to the difference in life strategy of plant species22,23 or drivers (abiotic versus biotic) of selection on plant phenology24. The spatial difference of temperature sensitivity may be caused by the genetic difference due to the adaptation to local climate25,26. However, the detailed mechanisms about the interspecific and spatial pattern of phenological response to temperature remain unclear.

In the present study, we developed a model for estimating temperature sensitivity of spring plant phenophases through parameterizing temperature sensitivity from an existing phenological model. Second, using the derived model, we calculated the temperature sensitivity of first leaf date (FLD) for 20 woody plant species at 43 sites in China. Third, we discussed the interspecific variation and spatial patterns of temperature sensitivity across eastern China. Our overall objectives are to derive a generalized model for calculating temperature sensitivity of spring phenophases and use it to interpret the observed interspecific variation and spatial pattern in phenological response to temperature.

Results

Using a phenological model based on degree day accumulation, the sensitivity in space (SpaceSens) model for calculating temperature sensitivity of spring plant phenophases was derived (equation (4), see Methods section). According to the SpaceSens model, the temperature sensitivity of the spring plant phenophase equals the length of period (L, which is used for performing linear regression between phenophases and temperature) divided by the difference between temperature at the onset date of phenophase (Tp) and base temperature threshold (Tb).

The temperature sensitivity of FLD for 20 plants species in Eastern China was simulated by using the SpaceSens model. The results show that the temperature sensitivity of FLD varies among different species and locations (Figure 1a). For example, on average, Salix babylonica FLD is most responsive with 3.68 days earlier for every 1°C increase in temperature of the previous month, while the response of Firmiana platanifolia FLD is minimum, with only 2.10 days earlier per 1°C increase (Figure 1b). Overall, the temperature sensitivity of FLD estimated from the SpaceSens model is −2.98 ± 0.61 days °C−1 (mean ± standard deviation, n = 668) for all site/species.

Figure 1
figure 1

Temperature sensitivity (Tsen) for first leaf date of 20 plant species plotted according to their onset date.

(a) Tsen for all site/species; (b) Tsen for all species averaged from their distribution sites (y = 0.02x − 4.7, R2 = 0.15, P = 0.08). Two plant species showing strongest and weakest temperature sensitivity are marked with different colors; (c) Tsen for all species at site 18 (y = 0.03x − 5.46, R2 = 0.26, P = 0.02); (d) Tp: y = 0.21x − 7.16; R2 = 0.98, P = 4.4*10−16; Tb: y = 0.13x − 9.63, R2 = 0.50, P = 5.0*10−4.

The early leaf unfolding time is associated with greater responsiveness to temperature variation across all sites/species, but this higher sensitivity for early-season species usually breaks down well before late-April (Figure 1a). If the FLD and its temperature sensitivity for each species is averaged over all sites, early-season species is also connected with stronger temperature sensitivity (Figure 1b, R = 0.38, P < 0.1). When considering all species at a single site, the temperature sensitivity of FLD correlated positively with their onset date in most of 43 sites (e.g. at site 18, R = 0.51, P < 0.05, Figure 1c). According to the SpaceSens model, if the L is constant, the temperature at the onset date of phenophase (Tp) and base temperature threshold (Tb) are the two factors determining the temperature sensitivity. Both Tp and Tb increase linearly with the onset date of FLD (taking Site 18 as an example, Figure 1d). However, the slope of Tb is less than Tp, which leads to the stronger temperature sensitivity for earlier phenophases.

The temperature sensitivity of FLD for a specific plant species highly depends upon latitude, although the relationships are not linear (Figure S1). Generally, the temperature sensitivity of FLD is greater at lower latitude than at higher latitude. For most species, the temperature sensitivity remains stable between 40°N and 50°N (taking Sophora japonica as an example, Figure 2). From 40°N to 28°N, the temperature sensitivity gradually becomes stronger. According to equation (4), the spatial variation of temperature sensitivity is determined by the spatial pattern of Tp, since Tb and L remain constant among locations for a specific species. The change of Tp along a latitudinal gradient is consistent with the latitudinal pattern of temperature sensitivity.

Figure 2
figure 2

The temperature at the onset date of phenophase (Tp) and temperature sensitivity (Tsen) for first leaf date of Sophora japonica plotted according to latitude.

The overall dependence of Tp and Tsen on latitude is high: for Tp: R2 = 0.93, y = −59.98 + 4.922*x − 0.1059x2 + 7.452*10−4x3, n = 41, P < 0.001; for Tsen: R2 = 0.94, y = −38.37 + 2.450*x − 0.059x2 + 4.204*10−4x3, n = 41, P < 0.001.

Discussion

The magnitude of phenological response to climate change obtained by using the SpaceSens model (−2.98 ± 0.61 days °C−1) closely matches the results of previous studies in Europe and China4,18,27,28,29. All these studies showed 2.6 to 6 days advancement in spring phenological events per °C increase5. Therefore, the SpaceSens model gives a reliable estimation for temperature sensitivity of spring phenology.

In early spring, the plant species with earlier FLD exhibit the strongest reactions to temperature change based on the SpaceSens model. This result coincides with the results of previous studies5,11. Across late spring, sensitivities are relatively stable (Figure 1a). This pattern was also documented in a study based on a phenological dataset in Europe30. Some studies give several possible explanations for these findings from an ecological perspective. First, the phenological adaptation related to the life strategy of plant species could cause the interspecific difference in phenological response22. The short-lived, early successional species are opportunistic and usually adopt a more risky life strategy31. Therefore, the phenological responsiveness of opportunistic species to temperature change is stronger than long-lived, late successional species32. Second, the earlier phenophases are more adapted to extreme low temperature and thus follow more closely climate warming, while the later phenophases are more cautious to advance their spring phenology because they are more easily to be injured by late frost events33,34. Third, in mid-latitude temperate systems, biotic factors (e.g. competition between species for soil resources, light or pollinators) rather than abiotic factors (e.g. temperature) might dominate phenology during mid-season periods, when the climate is more stable than the early season24. This study gives an explanation from the perspective of biometeorology: Tp-Tb correlated positively with the onset date of FLD (because the slope of Tb is less than Tp, see Figure 1d), which leads to the stronger temperature sensitivity for earlier phenophases.

Since the early-season species are more responsive to temperature change, the early-flowering species may flower even earlier than the late-flowering species if the climate warming continues in the future. More dispersed flowering times across species would have important implications. For example, the pollinators, especially which rely on continuously pollen recourse, may face a challenge to survive17. For plant species, their competition for pollinators and patterns of hybridization may be altered13. However, these consequences of interspecific difference in phenological change need continued and comprehensive monitoring in the future.

The spatial patterns of temperature sensitivity in the present study supports the results that the response of spring phenology to temperature in warmer locations or at lower latitudes is stronger than in colder location or at higher latitudes4,18,19. The local adaptation of plant phenology to regional climate may be one reason19. For example, the within-species variations of phenological response to temperature in Japan were less in the plant populations with lower genetic diversity26. More studies comparing the phenological responses between cloned plants and natural plants will be useful for quantifying the degree of genotypic difference affecting phenological responses. Our study indicated, even though there is no genotypic difference within species, the lower Tp at lower latitude (warmer locations) also can explain the spatial variation of phenological responsiveness to temperature (Figure 2). However, another study found that the observational sensitivities of phenology were not associated with latitude11. We noticed that the phenological observation data in reference 11 was mainly from the area between 40°N and 60°N, where the Tp is stable (Figure 2). If their study area is expanded to more southern area, the temperature sensitivity of spring phenophases would increase.

The spatial pattern of temperature sensitivity is a considerable factor impacting temporal trends of plant phenophases across space. If the increase in temperature is constant across space, the phenological trends would be weaker in higher latitude, since the temperature sensitivity is weaker there. Actually, higher latitude has warmed more than the lower latitude over the past 100 years35. Thus the effect of weaker temperature sensitivity and stronger warming trend in high latitude could partially offset each other, which make the latitude not an important factor influencing the phenological trends36.

Furthermore, the observed evidence demonstrates that temperature sensitivity changes with time for specific plant phenophases37,38. For example, the response to a 1°C increase in spring (March–May) temperature during a 30-year shifting time window was −2.5 to −5 days in Switzerland and −2.5 to −15 days in the UK from 1753 to 195839. In Beijing, China, the temperature sensitivity of first flowering date for 20 species during the 1990–2007 are shown to be greater than during the previous period (1963–1989), with a difference of −0.87 days per 1°C on average15. In this study, we use the average climate from 1961 to 1990 as the normal scenario to calculate temperature sensitivity. If using climate in other periods, the estimated temperature sensitivity is expected to be different. So the temporal shift of temperature sensitivity also can be explained by employing our SpaceSens model.

A prevalent decline in the temperature sensitivity of spring phenology with increasing local spring temperature variance was found at the species level from ground observations40. Through our analysis, local spring temperature variance (the multiyear averaged standard deviation of daily detrended air temperature from March to May) correlated significantly and positively with latitude (R = 0.49, P < 0.01, Figure S2). Therefore, higher latitude was associated with greater spring temperature variance and lower temperature sensitivity. This explains the result that temperature sensitivity declines with increasing local spring temperature variance40.

Methods

Model development

We adapt an existing phenological model to derive a generalized description of temperature sensitivity of spring plant phenophases. The degree day model (also called thermal time model) based on the linearity between development rate of plant and temperature has been well used in phenological studies41,42,43. The equation of the degree day model is:

where xt is the daily mean temperature at the day t. Tb is the base temperature threshold for temperature accumulation. Fi is the critical value of degree days for onset of the plant phenophase at site i. ty is the onset date of the plant phenophase. to is the start date of the degree days calculation. The degree day model implicitly assumes that environmental conditions (e.g. chilling requirement or day length) required to break dormancy have been met before t041.

The degree day model of a specific plant can only be applied at a single site, since the critical value of degree days sum for onset of the plant phenophase is different among different locations44. The previous experiments showed that degree day required to budburst decreases exponentially with increased amount of previous chilling43,45,46,47. Therefore, Fi can be described as an exponential function of coldness of the winter48:

where d, f, e are constant parameters. Ci is the mean winter temperature (December to February) at site i during a reference period (1961–1990). Equation (1) and (2) constitute the spatio-temporal phenological model.

In phenological studies, the temperature sensitivity of phenophases is usually measured as the slope coefficient of a linear regression model with phenophase as the dependent variable and temperature during a specific period before the mean phenophase as the independent variable. According to the principle of the spatio-temporal model (Figure 3), if the temperature during above-mentioned specific period (with a length of L) increased by ΔT (ΔT → 0), the increased degree days due to the temperature rise and the decreased degree days due to phenological advance should be equal (i.e. the areas of the two blue regions shown in Figure 3 should be equal). The equation could be expressed as follows:

where ty is the onset date of the plant phenophase in the normal scenario. TP is the temperature on date ty. ty1 is the onset date of plant phenophases in the warmer scenario. L is length of period, which is used for performing linear regression between phenophases and temperature (note that L is usually less than tyts, where ts is the day when the temperature exceed Tb).

Figure 3
figure 3

Scheme of sensitivity in space (SpaceSens) model.

Bold gray line represent normal temperature scenario in reference period and thin red line represent warmer scenario (a specific period with a length of L before phenophases increased by ΔT, ΔT → 0). t0: the staring date of the degree day calculation; Tb: base temperature threshold; ty: onset date of the plant phenophase in the normal scenario; ts: the day when temperature exceed Tb; Tp: temperature at the day ty; ty1: the onset date of plant phenophase in the warmer scenario.

Temperature sensitivity of phenophases (Tsen) can be regarded as the number of days changed in phenophases (Δt) divided by the change in temperature in degrees Celsius (ΔT). According to equation (3), the equation for Δt/ΔT is derived:

where Tsen is the temperature sensitivity of the plant phenophase; the meanings of three parameters (L, Tb and Tp) are same with equation (3). Hereafter, we call equation (4) sensitivity in space (SpaceSens) model.

Data source and statistical analysis

A total number of 5004 FLD observations for 20 woody plants from 43 sites with altitude less than 350 meters a.s.l. in eastern China were selected for calibration of the phenological models (Figure 4). The phenological data (1963–2009) was from the China Phenological Obsevation Network (CPON, http://www.cpon.ac.cn/). Phenology of each species was observed at 11 to 23 sites. The phenological time series have various lengths at different sites ranging from 1 to 37 years. According to the obsevation criteria of CPON, FLD was defined as the date when a individual formed a full first leaf. Usually, for each plant species, only one individual was observed at each site. If more than one individual is observed, the reported date is the average among all individuals of the same species. The sites span a wide latitudinal range (between 28°N and 49°N) and have a temperate or subtropical monsoon climate49. The daily mean temperature data from 1961 to 2009 at these sites were downloaded from the China Meteorological Data Sharing Service System (http://cdc.cma.gov.cn/cdc_en/home.dd).

Figure 4
figure 4

Locations of 43 phenological observation sites for studying the spatial pattern of temperature sensitivity.

The software ArcGis 9.3 was used to create the map.

Because the dataset from CPON included more than 600 plant species, we selected a number of typical woody species following these specific criteria: (1) the species must be a widespread plant which has been observed at more than 10 sites, because we want to study the spatial pattern of phenological responses; (2) the data of a species have as many observations as possible in order to calibrate the phenological models more accurately. As a result, 20 widespread plant species which have most abundant data were selected to investigate the interspecific and spatial difference of temperature sensitivity (Table 1). The current distribution range of each plant species was taken from a comprehensive atlas of woody plants in China50. Model parameters (t0, Tb, d, e, f) of spatio-temporal phenological model for simulating first leaf dates (FLD) were calibrated by Ge et al.48. External validation with independent observations in Ge et al.48 indicates that the model simulates FLD with a mean root-mean-square error (RMSE) of 6 days (see Table S1). The annual FLD of these 20 woody plants at sites within their present distribution range were simulated applying the model parameters in Table S1. The mean FLD of these species across their distribution sites all occurred in spring (from 20 March to 21 April, Table 1).

Table 1 Twenty plant species, their distribution ranges and first leaf date (FLD) averaged from sites within respective distribution range in eastern China. N: the latitude at the northern range limit; S: the latitude at the southern range limit. SD: standard deviation of FLD across sites within respective distribution range

For caculating temperature sensitivities, daily mean temperature time series were averaged over the reference period 1961–1990 for each site. Such daily mean temperature series can be regarded as temperature in the normal scenario (see Figure 3). The methods of estimating three variables in SpaceSens model were the following:

  1. 1

    Tp is a variable, which depends on species, sites and climatological mean. Firstly, we calculated the onset date of FLD (ty) in the normal scenario for each species/site by using the spatio-temporal phenological model described in Ge et al48 (see above). Then TP was measured as the temperature at ty (Figure 3).

  2. 2

    Tb is a species-specific variable, which was given through calibration of model parameters (Table S1).

  3. 3

    Because our goal was to compare the pattern of temperature sensitivity derived from our SpaceSens model with that derived from regression analysis, the value of L should follow the convention in phenological studies based on regression analysis. For most studies, the length of specific period (L) considered for regressing phenology against temperature is one (or two) month or a number of days before the mean phenophases and remained the same for all species and at all sites4,13,14,51. Thus we set L to 30 days as commonly used. In addition, since L did not vary with species and sites, the value of L would not affect the spatial and interspecific patterns of temperature sensitivity.

Consequently, the temperature sensitivities of FLD for each species/site could be estimated by using equation (4). At last, the relationship between the temperature sensitivity and FLD of 20 plant species at a fixed site and the relationship between the temperature sensitivity of FLD and latitude for a specific plant species were interpreted through regression analysis.