Spatiotemporal Variation of Osmanthus fragrans Phenology in China in Response to Climate Change From 1973 to 1996

Climate change greatly affects spring and autumn plant phenology around the world consequently, and significantly impacts ecosystem function and the social economy. However, autumn plant phenology, especially autumn flowering phenology, has not been studied so far. In this study, we examined the spatiotemporal pattern of Osmanthus fragrans phenology, including both leaf phenology (the date of bud-bust, BBD; first leaf unfolding, FLD; and 50% of leaf unfolding, 50 LD) and flowering phenology (the date of first flowering, FFD; peak of flowering, PFD; and end of flowering, EFD). Stepwise multiple linear regressions were employed to analyze the relationships between phenophases and climatic factors in the long term phenological data collected by the Chinese Phenological Observation Network from 1973 to 1996. The results showed that spring leaf phenophases and autumn flowering phenophases were strongly affected by latitude. BBD, FLD, and 50LD of O. fragrans were delayed by 3.98, 3.93, and 4.40 days as per degree of latitude increased, while FFD, PFD and EFD in O. fragrans advanced 3.11, 3.26, and 2.99 days, respectively. During the entire study period, BBD was significantly delayed across the region, whereas no significant trends were observed either in FLD or 50LD. Notably, all flowering phenophases of O. fragrans were delayed. Both leaf and flowering phenophases negatively correlated with growing degree-days (GDD) and cold degree-days (CDD), respectively. BBD and FLD were negatively correlated with total annual precipitation. In addition to the effects of climate on autumn flowering phenology, we found that earlier spring leaf phenophases led to delayed autumn flowering phenophases. Our results suggest that future climate change and global warming might delay the phenological sequence of O. fragrans. Our findings also advanced the flowering mechanism study of autumn flowering plants, and facilitated the accurate prediction of future phenology and climate change.

Introduction

Phenology is the study of the timing of recurring life-cycle events in plants, which rely on various biotic and abiotic factors (Lieth, 1974) and are triggered by changes in environmental conditions (Cong et al., 2017; Liu et al., 2021). Plant phenology is one of the most reliable biological indicators of climate change (Fu et al., 2017), and changes in plant phenology have important impacts on ecosystem structure and function (Fu et al., 2020), including carbon, water and nutrient cycling, hydrology, demography and biological interactions (Estiarte and Peñuelas, 2015; Xie et al., 2018). These changes can also cause a feedback loop that further augmenting changes in the climate system (Thackeray et al., 2016; Zeng et al., 2017). Monitoring plant phenological processes is important to improve our understanding of the impacts of global warming on plant and ecosystem function (Richardson et al., 2013; Fu et al., 2017; Chen et al., 2019). Previous studies have mainly focused on spring phenology, but less on autumn phenology (Gallinat et al., 2015), which also controls nitrogen cycling, ecosystem functions, and the associated feedbacks to climate systems (Visser, 2016). Moreover, the delayed autumn phenology may have greater influences than advanced spring phenology (Liu et al., 2016a), for instance in the regulation of carbon balance in temperate plants (Richardson et al., 2009). Therefore, studying the change of phenological sequence and the corresponding climatic drivers can improve our understanding of the response of species-specific phenology to environmental factors, as well as the responses of plants and ecosystem to ongoing climate change.

The magnitude of plant phenological responses to climate change is diverse and varies significantly in terms of geographical location, time and species (Chen et al., 2017). According to the Hopkins’ Bioclimatic Law, increased latitude and altitude can lead to delayed spring phenology, advanced autumn phenology, and an overall shorter growing season (Wang et al., 2015; Zhu et al., 2018). Many previous studies have investigated the timing of phenological events and their connection to climate change by using long-term records and remote sensing data in temperate and cold regions. Earlier spring phenology and later autumn phenology have been confirmed by these studies (Jeong et al., 2011; Gerst et al., 2016). However, in contrast to the wide attention on the variation of spring phenology, autumn flowering phenology has not received as much attention, because of the difficulties in interpretation (Xie et al., 2018). As a result, there are only a few studies that have investigated the dynamics of autumn phenology and the related controlling factors (Estiarte and Peñuelas, 2015; Liu et al., 2016b), particularly at regional scales (Gallinat et al., 2015; Fu et al., 2018). Shifting flowering phenology is a key biological indicator of climate change. However, the response of autumn flowering species to climate change is unclear. More studies are needed on autumn phenology and its controlling factors to gain a more holistic understanding of the impact of global climate change on ecological processes (Peng et al., 2019).

Temperature, precipitation and photoperiod are considered the primary controls of phenology (Fu et al., 2020). Among which temperature is the most important indicators (Li et al., 2018; Cheng et al., 2021). Warming advances spring plant phenology in temperate and boreal zones (Richardson et al., 2010). In subtropical regions, increase in chilling can advance the leaf-out (Song et al., 2020a). Nonlinear relationships (Cook et al., 2012) have been reported between spring phenology and temperature (Fu et al., 2015). The response of plant spring phenology to precipitation was inconclusive (Hänel and Tielbörger, 2015). Higher preseason precipitation could also alter spring phenology (Dai et al., 2013; Piao et al., 2019; Zhou et al., 2019), because more precipitation may increase the heat demand of spring phenology (Fu et al., 2014). Photoperiod may also influence bud burst, but generally in a minor way (Chen et al., 2017; An et al., 2020). To date, the mechanisms by which photoperiod affects phenology remain unexplored experimentally (Liu et al., 2016a). Compared to spring phenology, the linkages between autumn phenology and climatic factors have not been better identified (Fu et al., 2017; Piao et al., 2019; An et al., 2020). A higher preseason temperature delays autumn phenology in temperate regions (Cong et al., 2017; Zhu et al., 2017; Li et al., 2018; Fu et al., 2020). Photoperiod (Way and Montgomery, 2015), drought and heavy precipitation (Xie et al., 2018), precipitation (Richardson et al., 2013; Estiarte and Peñuelas, 2015) have been considered to influence autumn phenology. The drivers influencing autumn phenology vary across regions and among plant species. Low temperature and short daylength are considered as the dominant cues for autumn phenology (Gunderson et al., 2012). At high latitude colder regions, leaf senescence may be more responsive to photoperiod rather than temperature (Way and Montgomery, 2015). In contrast, temperature is a key factor for autumn leaf color (Tanino et al., 2010). In addition, previous studies have indicated that spring phenology also affects autumn phenology (Fu et al., 2014, 2018; Liu et al., 2016b), where an earlier spring phenology leads to a later autumn phenology (Fu et al., 2018; Škrk et al., 2020). The mechanism between these climatic factors and phenology on large spatial and temporal scales is still unclear.

In this study, we estimate the spatiotemporal variations of spring leaf phenology and autumn-flowering phenology of Osmanthus fragrans in response to climate change in China from 1973 to 1996. Specifically, this study attempted to solve the following questions: (i) What are the spatiotemporal variations of the spring and autumn phenology of O. fragrans in China from 1973 to 1996; (ii) How do the spatiotemporal variations in the spring leaf and autumn flowering phenology of O. fragrans correlate with climate factors? and (iii) What is the relationship between spring leaf phenology and autumn flowering phenology? In view of about O. fragrans was an autumn flowering species, according to the results of Sherry et al. (2007), we hypothesize that O. fragrans will flower earlier at high latitude compared to low latitude because of the variation of temperature. In addition, considering the spatial variation of precipitation and photoperiod is disordered, we hypothesize that photoperiod did not regulate the variation of flowering time of O. fragrans.

Materials and Methods

Study Area and Plant Species

The study area extended from 24 to 34°N and from 108 to 118°E, and the altitude ranged from 23 to 720 m. This provided a broad geographical coverage of the cultivation distribution area of O. fragrans. Observations centered on five stations in three different climatic zones of China. According to China’s eco-regional classification, one site (Xi’an) was in the warm temperate zone. Another (Wuhu) was in the northern subtropical zone. One station (Changde) was located in the transition zone from the north subtropical zone to the mid-subtropical zone, and two sites (Guilin and Liuzhou) were in the mid-subtropical zone (Figure 1). The temperature regimes of these five sites are shown in Figure 2 and Supplementary Data 1.

FIGURE 1

Figure 1. Locations of observation sites, overlain on a map of China’s eco-geographical regions (issued by the Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences). The five study locations represent three different climatic regions for O. fragrans cultivation.

FIGURE 2

Figure 2. Mean monthly temperatures at the five observation sites.

Osmanthus fragrans belongs to the Oleaceae family. It is an evergreen tree or shrub that originated in southwest China. This species flowers in the autumn and is grown ornamentally for its attractive foliage and fragrant edible flowers. China is the distribution center of O. fragrans and is rich in germplasm resources for the species, which is a model species of the Oleaceae family. It is one of the top ten most famous flowers, having been cultivated in China for over 2,500 years. Now it is widely distributed, has been planted in open fields from Shandong (36° north latitude) to San’ya (18° north latitude), Hainan Province. O. fragrans is an excellent garden tree and an important spice plant. As ornamental plants, their flowering seasons have important aesthetic and economic benefits for some regions. For example, O. fragrans festivals are of major cultural importance in certain parts of China.

Phenological Data and Climate Data

The phenological data of O. fragrans collected at five locations in the study area from 1973 to 1996 was obtained from the Chinese Phenological Observation Network (CPON) (Supplementary Data 2). Details of the phenological observation method conform to standardized observation criteria and guidelines (Wan and Liu, 1979; Dai et al., 2014). Six phenophases were investigated, including three spring leaf phenophases and three autumn flowering phenophases. The BBD was defined as the date the appearance of the first green leaf tip on a few twigs of the observed tree (Murray et al., 1989; Basler and Körner, 2014). The FLD was defined as the date when the first batch of leaves is fully spread on a few twigs of the observed tree; whereas 50LD is defined as the date when leaves are fully spread on 50% of twigs of the observed tree (Chen et al., 2017). The FLD represented the date when the first leaf unfolded in at least one individual, while the 50LD represented the date when more than 50% of individuals had unfolded their leaves (Huang et al., 2019). The FFD was defined as the date when 10% flowers were open, the PFD was defined as the date when 90% flowers were in full bloom, and the EFD was the date on which the majority of small flowers (about 90%) had withered or dropped (Yang, 2013). To ensure the accuracy and validity of the statistical analysis, the phenological time series with more than 10 years’ worth of observations were chosen (Dai et al., 2014). The original phenological data was preprocessed according to the 30-day rule to remove outliers that possibly had been recorded incorrectly (Wang et al., 2018). Only five stations (Liuzhou, Guilin, Changde, Wuhu, and Xi’an) in the time period, from 1973 to 1996, were chosen for analysis. In these stations, O. fragrans was planted in batches and long observed. Generally, more than 5 trees of the same species were selected as the observation target to record the phenology, and a branch without shade was selected from the middle and upper four directions (East, South, West, and North) of each tree to observe its phenology. In total, 70, 80, 82, 77, 72, and 76 time series were analyzed for BBD, FLD, 50LD, FFD, PFD, and EFD, respectively. All onset data was converted to Julian day (DOY). The number of phenological samples (the total number of observations) observed at each site did differ (Table 1 and Supplementary Data 3). All phenological data were obtained from cultivated populations of O. fragrans.

TABLE 1

Table 1. Number of years of phenological period during 1973–1996 years at each location.

The data of daily climate factors, including daily mean air temperature, daily precipitation and daily sunshine hours were obtained from the China Meteorological Sharing Service Network. The meteorological stations are all located at or nearby the corresponding phenological stations. They were Liuzhou and Guilin stations in Guangxi Autonomous Region, Changde station in Hunan Province, and Wuhu station in Anhui Province, respectively. However, because we cannot directly download meteorological data of Xi’an in the China Meteorological Sharing Service Network, the data were collected from the Great Wild Goose Pagoda station in Xi’an Meteorological Bureau. These data are well checked and normalized by China Meteorological Administration before being issued.¹

Growing degree-days (GDD) is a heat accumulation index for analyzing the effect of temperature on spring phenophases (Fu et al., 2014; Liu et al., 2019). Cold degree-days (CDD) were applied as a corresponding metric to explain the variation in autumn phenology (Dragoni and Rahman, 2012; Xie et al., 2018).

According to McMaster and Wilhelm (1997), GDD was calculated as the thermal sum of the difference between daily mean temperature (T_m) and the base temperature (T_b) between DOY₁ and DOY₂ [formula (1)].

$$GDD=\sum _{DO{Y}_1}^{DO{Y}_2}\left(\genfrac{}{}{0.0pt}{}{\left(T_m-T_b\right)when{T}_m>T_b}{\mathrm{\hspace{0.25em}0}when{T}_m\le T_b}\right)(1)$$

Where, DOY₁ represents the starting date of the research phase. It was set as January 1st here. DOY₂ represents the ending date of the research phase. To analyze the correlation between spring phenology and GDD across different stations and years, the end dates of GDD were set to April 30th for BBD, FLD, and 50LD (Zhu et al., 2018).

Temperatures above 0°C (Hänninen, 1990) and a temperature of 5°C (Marchin et al., 2015) represented the minimum temperature threshold required for stimulating budburst. Experiments have shown that base temperature is generally accepted to be below 10°C, although it is species-specific (Harrington et al., 2010; Malyshev, 2020). Three base temperatures, 0, 5, and 10°C were tested in this study. Hereafter we only report results using the base temperature of 5°C. The temperature threshold of 0°C and 10°C were found similar results (Supplementary Table 1 and Supplementary Data 4).

Cold Degree-Days was calculated by summing the deviation between the base temperature (T_b) and daily mean temperature (T_m) between DOY₁ and DOY₂ [formula (2)].

$$CDD=\sum _{DO{Y}_1}^{DO{Y}_2}\left(\genfrac{}{}{0.0pt}{}{\left(T_b-T_m\right)when{T}_b>T_m}{\mathrm{\hspace{0.25em}0}when{T}_b\le T_m}\right)(2)$$

The length of preseason was chosen to be 2 months before mean phenological events following previous studies (Cong et al., 2013; Sakkir et al., 2015). Temperatures during the preseason is the most important for flowering phenology for most species in Europe and Asia (Fu et al., 2015; Zhang et al., 2015). In the previous researches, DOY₁ was set as 1st August, DOY₂ was set to 31st October (Xie et al., 2018; Zhu et al., 2018). The earliest flowering phenology of O. fragrans occurs in Xi’an was in late September, and the average flowering phenology was in the early October across all the stations and years. In this manuscript the DOY₁ was set as 1st July, DOY₂ was set to 31st October to calculate the CDD.

The threshold temperature is accepted to be between 20 and 26°C (Delpierre et al., 2009). The base temperature for CDD is 25°C in this paper. We also tested the base temperature threshold of 20 and 30°C, and very similar results were observed (Supplementary Table 2 and Supplementary Data 5).

The cumulative precipitation (PPT) and cumulative sunshine duration (SSD) during the same periods of GDD or CDD were calculated for each station and each year. Considering photoperiod may also regulate phenology, we took the relationships between daily sunshine hours (i.e., day length) and phenological variables as the photoperiod effect.

Statistical Analyses

A stepwise multiple linear regression model was used to estimate the effect of geographical factors and years on the spatiotemporal variation of phenology in the five study areas from 1973 to 1996. Geographical factors (longitude, latitude, and altitude) and year were used as independent variables, and phenological metrics were dependent variable. The variance inflation factor (VIF) was used to quantify the multi-collinearity among variables in the models. The model was accepted when the VIF of individual predictors was less than three, which indicated a lack of multi-issues with collinearity (Zuur et al., 2010). Squared semi-partial correlation coefficients were used to test the relative contribution of each independent variable in a given model, which were determined as the reduction in R² with removing a given predictor from the set of independent variables (Watson et al., 2011). Values were reported as fractions of the original R² (Zhu et al., 2018).

Stepwise multiple linear regression analysis was used to analyze the climate controls on the phenophases across all stations and years. Simple linear regression analyses were conducted between climatic factors and phenophases to quantify the effects of climatic factors on the phenophases, One-way ANOVA was used to test the significance of climate factors including GDD, CDD, TTP, and SSD in five stations. The correlation between spring and autumn phenophases was assessed using simple linear correlation. P-values less than 0.05 were considered significant. All analyses were performed using SPASS 22.0 (SPASS Inc., Chicago, IL, United States).

Results

Spatiotemporal Variation of Spring Leaf Phenophases in Osmanthus fragrans

Based on the results from the multiple stepwise regression analysis, the phenophases showed significant correlation with geographical factors (Table 2 and Supplementary Data 6). Latitude has the largest squared semi-partial correlation coefficients of all the factors, suggesting that it was the most profound factor that influencing phenology, moreover, it also was the factor that contributes most to the spatial variation in all phenological metrics. Due to the strong dependence of phenophases on latitude, BBD, FLD and 50LD showed significant positive correlations with latitude in O. fragrans (P < 0.001) (Figures 3A–C and Supplementary Data 7).

TABLE 2

Table 2. Stepwise multiple linear regression coefficients of each independent variable related to phenological metrics.

FIGURE 3

Figure 3. Relationship between O. fragrans phenophase and latitude. (A) BBD, (B) FLD, (C) 50LD, (D) FFD, (E) PFD, and (F) EFD. BBD, date of bud-burst; FLD, date of first leaf unfolding; 50LD, date of 50% of leaf unfolding; FFD, first flowering day; PFD, peak flowering day; EFD, end of flowering day. DOY indicates the Julian day of year.

The impact of geographic factors on phenophases varied (Table 2). BBD, FLD and 50LD were delayed in O. fragrans by 3.98, 3.93 and 4.40 days on average per degree of latitude, respectively. The geographical factors explained more than 60% of the spatial variation in the spring phenology. No significant relationship was found between leaf phenophases and longitude or altitude, except for the correlation between 50LD and altitude (p < 0.01) (Table 2), where 50LD was delayed 0.018 days per m increasing in altitude.

Across all the stations, the multiple stepwise regression analysis shows that BBD was delayed by an average of about 0.69 day per year from 1973 to 1996 (P < 0.01) (Table 2), whereas FLD and 50LD didn’t show significant temporal trends (Table 2). The data from the stations at Changde and Xi’an had significant positive trends in BBD, but none of other stations (Figure 4A). In addition, the only delayed trends were observed at Changde station in FLD and 50LD (Figures 4B,C and Supplementary Data 8).

FIGURE 4

Figure 4. Interannual variations of (A) BBD, (B) FLD, (C) 50LD, (D) FFD, (E) PFD, and (F) EFD for O. fragrans at the five observation sites and the mean observation of five sites from 1973 to 1996. BBD, date of bud-burst; FLD, date of first leaf unfolding; 50LD, date of 50% of leaf unfolding; FFD, first flowering day; PFD, peak flowering day; EFD, end of flowering day. DOY indicates the Julian day of year.

Spatiotemporal Variation of Autumn Flowering Phenophases in Osmanthus fragrans

In contrast to spring leaf phenophases, autumn FFD, PFD, and EFD in O. fragrans had significant and negative correlations with latitude across the whole area from 1973 to 1996 (Figures 3D–F and Supplementary Data 7). FFD, PFD, and EFD in O. fragrans advanced 3.11, 3.26, and 2.99 days, respectively, by every degree northward. In each of these measures, the geographical factors accounted for 61, 61, and 55% of the spatial variation, respectively. There was no significant variation among autumn flowering phenophases regarding longitude and altitude (Table 2 and Supplementary Data 6).

The flowering phenophases of O. fragrans showed a delayed trend in the study area (Table 2). FFD was delayed by 0.67 day (P < 0.001) (Table 2) per year, PFD was delayed by 0.52 day (P < 0.01) per year (Table 2), and EFD was delayed by 0.43 day (P < 0.05) (Table 2) per year. However, no significant temporal trends were observed in flowering phenophases in each station (Figures 4D–F and Supplementary Data 8).

Effects of Climatic Factors on Osmanthus fragrans Phenology

Spring phenophases showed negative correlation with GDD and PPT, and positive correlation with SSD by simple linear regression analyses. However, autumn flowering phenophases were only negatively correlated with CCD, but not with PPT and SSD (Supplementary Table 3 and Supplementary Data 9). SSD was no significantly different during CCD among the five stations by one-way ANOVA (P > 0.05). Stepwise multiple linear regression analysis demonstrated that GDD and CDD were the major factors influencing O. fragrans phenophases (Table 3). No significant correlations were found between other phenophases and PPT. Only BBD (P < 0.01) and FLD (P < 0.01) were significantly negatively correlated with PPT. The correlation was extremely weak between phenophases and SSD (Table 3 and Supplementary Data 10). BBD, FLD and 50LD advanced 4.9, 5.2, and 5.2 days, respectively, with every GDD increase of 100°C-days (Figures 5A–C). FFD, PFD and 50LD advanced 4.0, 4.4, and 3.8 days, respectively, with every CDD increase of 100°C-days (Figures 5D–F and Supplementary Data 10).

TABLE 3

Table 3. Partial correlation coefficients between the phenological metrics and climatic factors across all stations and years.

FIGURE 5

Figure 5. The relationship between the dates of phenophases and growing degree-days/cold degree-days (GDD/CDD) across all stations and years for (A) BBD, (B) FLD, (C) 50LD, (D) FFD, (E) PFD, and (F) EFD. BBD, date of bud-burst; FLD, date of first leaf unfolding; 50LD, date of 50% of leaf unfolding; FFD, first flowering day; PFD, peak flowering day; EFD, end of flowering day. DOY indicates the Julian day of year.

Relationship Between Spring Phenology and Autumn Phenology

A significantly negative correlation between leaf and flowering phenophases dominated the entire study area from 1973 to 1996 (R² = 0.22 to 0.35, P < 0.001, Figure 6 and Supplementary Data 11). Flowering phenophases were delayed by 0.42–0.48 day by every single day of advanced BBD. FFD, PFD, and EFD were delayed 0.52, 0.51, and 0.51 day, respectively, by every day of advanced FLD. Autumn flowering phenophases were delayed 0.52–0.55 day for every day of advanced 50LD.

FIGURE 6

Figure 6. The relationship of (A) FFD with BBD, (B) PFD with BBD, (C) EFD with BBD, (D) FFD with FLD, (E) PFD with FLD, (F) EFD with FLD, (G) FFD with 50FLD, (H) PFD with 50FLD, (I) EFD with 50FLD. BBD, date of bud-burst; FLD, date of first leaf unfolding; 50LD, date of 50% of leaf unfolding; FFD, first flowering day; PFD, peak flowering day; EFD, end of flowering day. DOY indicates the Julian day of year. Numbers in brackets indicate the squared correlation coefficients.

Discussion

The Effect of Geographical Factors on the Spatial Distributions of Phenology

In this study, there was a clear pattern of the geographical distribution of O. fragrans phenology in China. All spring and autumn phenological metrics of O. fragrans showed strong dependence on geographical factors (Table 2). We found that spring phenophases were delayed with increased latitude, from south to north, while autumn flowering phenophases were advanced at higher latitudes. In other words, the period between leaf phenology and flowering phenology shortened from south to north (Figure 3). In line with the previous study, we focused on the geographical distribution of leaf phenophases in China (Dai et al., 2014). The spring phenophases were delayed about 4 days per degree of latitude across the studied areas (P < 0.001), which is consistent with hypothesis that spring and summer phenologies are delayed 4 days per degree in latitude (Hopkins, 1918; Richardson et al., 2019). These latitudinal patterns are consistent with the previous studies that had examined widespread spring species along latitudinal gradients (Dai et al., 2014; Wang et al., 2015; Gerst et al., 2016; Song et al., 2021). Similar latitudinal patterns were also detected in autumn phenology (Gill et al., 2015; Zhu et al., 2018). Many spring plant phenological events require a certain amount of chilling accumulation for dormancy release (Cook et al., 2012; Fu et al., 2014; Du et al., 2019). Chilling requirements vary among populations of the same species, depending on latitude. Low-latitude populations may require less chilling than that of high-latitude populations (Harrington et al., 2010; Liang, 2016). This latitudinal pattern might be attributed to the long-term adaptation to local climate (Vitasse et al., 2018), and it is likely connected to geographic variation in population genetics (Liang, 2016).

Changes in Phenology of Osmanthus fragrans

This study demonstrated that BBD was significantly delayed at an average rate of 0.69 day per year, while FLD and 50LD displayed marginal trends (Table 2). These findings are inconsistent with most of previous studies (Ge et al., 2015; Wang et al., 2018), which have demonstrated that spring plant phenophases have occurred earlier over the past several decades with global warming. These discrepancies are likely attributed to climate zone. The advance has largely been observed in temperate, boreal, alpine, or in subalpine climates (Shen et al., 2014; Wang et al., 2015). Notably, some studies indeed shown delayed spring bud burst due to climate warming in some warm temperate to subtropical regions (Yu et al., 2010; Zhang et al., 2021), consistent with our finding of spring phenophases. The delay might be due to the warm living habits and the distributional ranges of O. fragrans. The insufficient accumulation of chilling temperature during dormancy counteracted the force of temperature accumulation induced by climate warming, which in total delayed leaf phenophases (Xie et al., 2015; Chen et al., 2017). Meanwhile we found that autumn flowering phenophases were significantly delayed in O. fragrans (Table 2), which is consistent with the patterns in other autumn phenological studies (Ge et al., 2015; Shen et al., 2020). Previous reports on fall phenophases, such as leaf senescence and the end of season (EOS), have revealed substantially delayed trends associated with increasing temperatures (Gallinat et al., 2015; Cong et al., 2017; Fu et al., 2017). Control and observation experiments have also shown that the flowering of fall flowering species is often delayed (Sherry et al., 2007; Pearson, 2019). The delays in reproductive phases generally could be attributed to the fact that warming delays the fulfillment of chilling requirements and thus leads to later reproduction (Yu et al., 2010; Peng et al., 2019).

In addition, different temporal trends in O. fragrans phenology have been found at each site, suggesting that plant phenology may also be influenced by other complex factors. In addition to temperature, other environmental factors such as precipitation, photoperiod and nutrient availability, as well as biological factors, may have certain impacts on spring and autumn phenology (Piao et al., 2019). Photoperiod and temperature may co-regulate plant phenology (Fu et al., 2020), and the effect of precipitation on plant phenology can be explained by its indirect impacts on the thermal requirement (GDD or CDD) for both spring and autumn phenological events (Fu et al., 2014; Hänninen, 2016). Further, the interplay between soil water content and nutrients also has been linked to plant phenology (Estiarte and Peñuelas, 2015). Our results reinforce the need for phenological analysis of the interaction of multiple environmental factors. In order to fully understand the physiological and molecular genetic mechanisms underlying the phenological patterns and processes, controlled physiological experiments are needed.

Climatic Factors Driving the Spatiotemporal Variation in Phenology

Previous studies have shown that the spatiotemporal variation of plant phenology was strongly driven by climate, and that temperature was the primary factor influencing plant phenology (Jeong et al., 2011; Xu and Chen, 2013). This study observed significant negative partial correlations between spring/autumn phenophases and GDD/CDD in O. fragrans across all stations and years. More thermal accumulation tended to advance the occurrence of spring phenophases, while less heat deficit resulted in delayed autumn phenophases (Figure 5 and Table 3). Similarly, negative relationships have been discovered between temperature and phenological events, such as flowering and leaf-out (Willis et al., 2017). The relationships between plant autumn phenology and various climatic factors are more complicated (Ge et al., 2015; An et al., 2020). Temperature, expressed in terms of CDD, was supposed to be the most important controlling factor during the end of season (EOS) in deciduous forests (Dragoni and Rahman, 2012). Air temperature and leaf senescence have been linked at the bio-molecular scale (Allona et al., 2008). Previous studies have reported that cumulative cold temperatures below a certain threshold trigger autumn phenology (Delpierre et al., 2009; Dragoni and Rahman, 2012; Gill et al., 2015). Warmer temperatures in autumn lower CDD and postpone autumn phenology (Ge et al., 2015; Liu et al., 2016b; Shen et al., 2020). For autumn flowering species, flowering dates were delayed by warm preseason temperature (Sherry et al., 2007; Pearson, 2019), most of which were consistent with our studies.

Our results showed that BBD and FLD advanced with increased precipitation, while other phenological metrics remain unchanged (Table 3), which are consistent with those reported by previous studies (Hänel and Tielbörger, 2015; Zhou et al., 2019). Water availability played a dominant role in regulating plant growth (Liu et al., 2009). In other words, plants grow earlier with increased precipitation, accordingly, increased soil water content. In addition, the temporal distribution of precipitation, rather than the total amount, may be more related to plant phenology (Dragoni and Rahman, 2012).

In our study, the spring phenology of O. fragrans was positively correlated with SSD by simple linear regression (Supplementary Table 3), but weak relationship between spring phenology and SSD by stepwise multiple linear regression analysis (Table 3). Our results are consistent with these previous studies (Polgar et al., 2014; Du et al., 2019; Song et al., 2020b). For all species, the responses of phenology to photoperiod are species-specific (Basler and Körner, 2014). Most woody plant species are not sensitive to photoperiod, even for the photoperiod sensitive species, many studies suggest their budburst dates were not sensitive to photoperiod, when the chilling demand was met (Zohner et al., 2016).

The effect of one climatic factor could be exacerbated or mitigated by another (Mitchell et al., 2015), so accurate prediction of climate change on plant phenology requires a holistic consideration of the joint and interactive effects of all factors. However, few experimental studies have quantified the individual and combined effects of different environmental cues on plant phenology (Pérez-Ramos et al., 2020).

Conclusion

In our study, we demonstrated that the spring phenophases of O. fragrans were delayed and the autumn phenophases were advanced as latitude increased. Further analysis found that temperature rather than precipitation was the most important driver of the spatiotemporal variation in the phenologies of O. fragrans. In addition to the effects of climate on autumn flowering phenology, earlier spring leaf phenophases could cause delayed autumn flowering phenophases. Overall, these results suggest that temperature affects the autumn phenology of O. fragrans through both direct and indirect factors. Therefore, the future climate warming may further delay the autumn phenological sequence of O. fragrans, thus, it is necessary to conduct deeper studies of the autumn phenology of O. fragrans.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

XW and FS conceived the ideas. XW, YL, XL, and MZ analyzed the data. XW, XL, and SH drew Figure. XW drafted the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Henan Province Major Research Fund of Public Welfare (No. 201300110900) and the National Natural Science Foundation of China (No. 31600375).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank China phenological observation network for providing the phenological observation data of Osmanthus fragrans. The authors are grateful to Elizabeth Tokarz at Yale University and Ding Liu at Harvard University for their assistances with English language and grammatical editing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.716071/full#supplementary-material

Footnote

1. ^ http://data.cma.cn/

References

Allona, I., Ramos, A., Ibáñez, C., Contreras, A., Casado, R., and Aragoncillo, C. (2008). Molecular control of winter dormancy in establishment in trees. Span. J. Agric. Res 6, 201–210.

Google Scholar

An, S., Chen, X. Q., Zhang, X. Y., and Lang, W. G. (2020). Precipitation and minimum temperature are primary climatic controls of alpine grassland autumn phenology on the Qinghai-Tibet Plateau. Remote Sens. 12:431. doi: 10.3390/rs12030431

CrossRef Full Text | Google Scholar

Basler, D., and Körner, C. (2014). Photoperiod and temperature responses of bud swelling and bud burst in four temperate forest tree species. Tree Physiol. 34, 377–388. doi: 10.1093/treephys/tpu021

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Huang, J. G., Ma, Q. Q., Hänninen, H., Tremblay, F., and Bergeron, Y. (2019). Long-term changes in the impacts of global warming on leaf phenology of four temperate tree species. Glob. Chang. Biol. 25, 997–1004. doi: 10.1111/gcb.14496

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X. Q., Wang, L. X., and Inouye, D. (2017). Delayed response of spring phenology to global warming in subtropics and tropics. Agric. For. Meteorol. 234-235, 222–235. doi: 10.1016/j.agrformet.2017.01.002

CrossRef Full Text | Google Scholar

Cheng, W., Li, Z., and Yan, L. (2021). Uniforming spring phenology under non-uniform climate warming across latitude in China. Sci. Total Environ. 762:143177. doi: 10.1016/j.scitotenv.2020.143177

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, N., Shen, M. G., and Piao, S. L. (2017). Spatial variations in responses of vegetation autumn phenology to climate change on the Tibetan Plateau. J. Plant Ecol. 10, 744–752. doi: 10.1093/jpe/rtw084

CrossRef Full Text | Google Scholar

Cong, N., Wang, T., Nan, H., Ma, Y., Wang, X., Myneni, R. B., et al. (2013). Changes in satellite-derived spring vegetation green-up date and its linkage to climate in China from 1982 to 2010: a multimethod analysis. Glob. Change Biol. 19, 881–891. doi: 10.1111/gcb.12077

PubMed Abstract | CrossRef Full Text | Google Scholar

Cook, B. I., Wolkovich, E. M., and Parmesan, C. (2012). Divergent responses to spring and winter warming drive community level flowering trends. Proc. Natl. Acad. Sci. U.S.A. 109, 9000–9005. doi: 10.1073/pnas.1118364109

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, J. H., Wang, H. J., and Ge, Q. (2014). The spatial pattern of leaf phenology and its response to climate change in China. Int. J. Biometeorol. 58, 521–528. doi: 10.1007/s00484-013-0679-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, J. H., Wang, H. J., and Ge, Q. S. (2013). Multiple phenological responses to climate change among 42 plant species in Xi’an, China. Int. J. Biometeorol. 57, 749–758. doi: 10.1007/s00484-012-0602-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Delpierre, N., Dufrêne, E., Soudani, K., Ulrich, E., Cecchini, S., Boé, J., et al. (2009). Modelling ainterannual and spatial variability of leaf senescence for three deciduous tree species in France. Agric. For. Meteorol. 149, 938–948. doi: 10.1016/j.agrformet.2008.11.014

CrossRef Full Text | Google Scholar

Dragoni, D., and Rahman, A. F. (2012). Trends in fall phenology across the deciduous forests of the Eastern USA. Agric. For. Meteorol. 157, 96–105. doi: 10.1016/j.agrformet.2012.01.019

CrossRef Full Text | Google Scholar

Du, Y. J., Pan, Y. Q., and Ma, K. P. (2019). Moderate chilling requirement controls budburst for subtropical species in China. Agric. For. Meteorol. 278:107693. doi: 10.1016/j.agrformet.2019.107693

CrossRef Full Text | Google Scholar

Estiarte, M., and Peñuelas, J. (2015). Alteration of the phenology of leaf senescence and fall in winter deciduous species by climate change: effects on nutrient proficiency. Glob. Chang. Biol. 21, 1005–1017. doi: 10.1111/gcb.12804

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y. H., Piao, S., Delpierre, N., Hao, F., Hänninen, H., Liu, Y., et al. (2017). Larger temperature response of autumn leaf senescence than spring leaf-out phenology. Glob. Chang. Biol. 24, 2159–2168. doi: 10.1111/gcb.14021

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y. H., Piao, S., Zhao, H., Jeong, S. J., Wang, X., Vitasse, Y., et al. (2014). Unexpected role of winter precipitation in determining heat requirement for spring vegetation green-up at northern middle and high latitudes. Glob. Chang. Biol. 20, 3743–3755. doi: 10.1111/gcb.12610

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y. H., Zhao, H., Piao, S., Peaucelle, M., Peng, S., Zhou, G., et al. (2015). Declining global warming effects on the phenology of spring leaf unfolding. Nature 526, 104–107. doi: 10.1038/nature1540

CrossRef Full Text | Google Scholar

Fu, Y. S., Li, X. X., Zhou, X. C., Geng, X. J., Guo, Y. H., and Zhang, Y. R. (2020). Progress in plant phenology modeling under global climate change. Sci. China Earth Sci. 63, 1237–1247. doi: 10.1007/s11430-019-9622-2

CrossRef Full Text | Google Scholar

Fu, Y. Y., He, H., Zhao, J. J., Larsen, D., Zhang, H. Y., Sunde, M. G., et al. (2018). Climate and spring phenology effects on autumn phenology in the Greater Khingan Mountains, northeastern China. Rem. Sens. 10:449. doi: 10.3390/rs10030449

CrossRef Full Text | Google Scholar

Gallinat, A. S., Primack, R. B., and Wagner, D. L. (2015). Autumn, the neglected season in climate change research. Trends Ecol. Evol. 30, 169–176. doi: 10.1016/j.tree.2015.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, Q. S., Wang, H. J., Rutishauser, T., and Dai, J. H. (2015). Phenological response to climate change in China: a meta-analysis. Glob. Chang. Biol. 21, 265–274. doi: 10.1111/gcb.12648

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerst, K. L., Kellermann, J. L., Enquist, C. A. F., Rosemartin, A. H., and Denny, E. G. (2016). Estimating the onset of spring from a complex phenology database: trade-offs across geographic scales. Int. J. Biometeorol. 60, 391–400. doi: 10.1007/s00484-015-1036-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill, A. L., Gallinat, A. S., Sanders-Demott, R., Rigden, A. J., Short Gianotti, D. J., Mantooth, J. A., et al. (2015). Changes in autumn senescence in northern hemisphere deciduous trees: a meta-analysis of autumn phenology studies. Ann. Bot. 116, 875–888. doi: 10.1093/aob/mcv055

PubMed Abstract | CrossRef Full Text | Google Scholar

Gunderson, C. A., Edwards, N. T., Walker, A. V., O’hara, K. H., Campion, C. M., and Hanson, P. J. (2012). Forest phenology and a warmer climate–growing season extension in relation to climatic provenance. Glob. Chang. Biol. 18, 2008–2025. doi: 10.1111/j.1365-2486.2011.02632.x

CrossRef Full Text | Google Scholar

Hänel, S., and Tielbörger, K. (2015). Phenotypic response of plants to simulated climate change in a long-term rain-manipulation experiment: a multi-species study. Oecologia 177, 1015–1024. doi: 10.1007/s00442-015-3231-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hänninen, H. (1990). Modelling bud dormancy release in trees from cool and temperate regions. Acta For. Fennica 213, 1–47. doi: 10.14214/aff.7660

CrossRef Full Text | Google Scholar

Hänninen, H. (2016). Boreal and Temperate Trees in a Changing Climate. Netherlands: Springer, doi: 10.1007/978-94-017-7549-6

CrossRef Full Text | Google Scholar

Harrington, C. A., Gould, P. J., and Clair, J. B. (2010). Modeling the effects of winter environment on dormancy release of Douglas-fir. For. Ecol. Manag. 259, 798–808. doi: 10.1016/j.foreco.2009.06.018

CrossRef Full Text | Google Scholar

Hopkins, A. D. (1918). Periodical events and natural law as guides to agricultural research and practice. Mon. Weather Rev. 9:42.

Google Scholar

Huang, W. J., Ge, Q. S., Wang, H. J., and Dai, J. H. (2019). Effects of multiple climate change factors on the spring phenology of herbaceous plants in Inner Mongolia, China: Evidence from ground observation and controlled experiments. Int. J. Climatol. 39, 5140–5153. doi: 10.1002/joc.6131

CrossRef Full Text | Google Scholar

Jeong, S. J., Ho, C. H., Gim, H. J., and Brown, M. E. (2011). Phenology shifts at start vs. end of growing season in temperate vegetation over the Northern Hemisphere for the period 1982-2008. Glob. Chang. Biol. 17, 2385–2399. doi: 10.1111/j.1365-2486.2011.02397.x

CrossRef Full Text | Google Scholar

Li, P., Peng, C. H., Wang, M., Luo, Y. P., Li, M. X., Zhang, K. R., et al. (2018). Dynamics of vegetation autumn phenology and its response to multiple environmental factors from 1982 to 2012 on Qinghai-Tibetan Plateau in China. Sci. Total Environ. 855, 637–638. doi: 10.1016/j.scitotenv.2018.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, L. (2016). Beyond the bioclimatic law: geographic adaptation patterns of temperate plant phenology. Prog. Phys. Geog. 40, 811–834. doi: 10.1177/0309133316656558

CrossRef Full Text | Google Scholar

Lieth, H. (1974). “Purposes of a phenology book,” in Phenology and Seasonality Modeling, ed. H. Lieth (New York, NY: Springer), 3–19. doi: 10.1007/s40656-021-00428-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Fu, Y. H., Zeng, Z. Z., Huang, M. T., Li, X., and Piao, S. L. (2016a). Temperature, precipitation, and insolation effects on autumn vegetation phenology in temperate China. Glob. Chang. Biol. 22, 644–655. doi: 10.1111/gcb.13081

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Fu, Y. H., Zhu, Z. C., Liu, Y. W., Liu, Z., Huang, M. T., et al. (2016b). Delayed autumn phenology in the Northern Hemisphere is related to change in both climate and spring phenology. Glob. Chang. Biol. 22, 3702–3711. doi: 10.1111/gcb.13311

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W. X., Zhang, Z., and Wan, S. Q. (2009). Predominant role of water in regulating soil and microbial respiration and their responses to climate change in a semiarid grassland. Glob. Change Biol 15, 184–195. doi: 10.1111/j.1365-2486.2008.01728.x

CrossRef Full Text | Google Scholar

Liu, Y., Li, G., Wu, X., Niklas, K. J., and Sun, S. (2021). Linkage between species traits and plant phenology in an alpine meadow. Oecologia 195, 409–419. doi: 10.1007/s00442-020-04846-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. J., Qin, Y., and Ge, Q. (2019). Spatiotemporal differentiation of changes in maize phenology in China from 1981 to 2010. J. Geogr. Sci. 29, 351–362. doi: 10.1007/s11442-019-1602-5

CrossRef Full Text | Google Scholar

Malyshev, A. V. (2020). Warming events advance or delay spring phenology by affecting bud dormancy depth in trees. Front. Plant Sci. 11:856. doi: 10.3389/fpls.2020.00856

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchin, R. M., Salk, C. F., Hoffmann, W. A., and Dunn, R. R. (2015). Temperature alone does not explain phenological variation of diverse temperate plants under experimental warming. Glob. Chang. Biol. 21, 3138–3151. doi: 10.1111/gcb.12919

PubMed Abstract | CrossRef Full Text | Google Scholar

McMaster, G. S., and Wilhelm, W. W. (1997). Growing degree-days: one equation, two interpretations. Agric. For. Meteorol. 87, 291–300. doi: 10.1016/S0168-1923(97)00027-0

CrossRef Full Text | Google Scholar

Mitchell, P., Wardlaw, T., and Pinkard, L. (2015). “Combined stresses in forests,” in Combined Stresses in Plants: Physiological, Molecular, and Biochemical Aspect, ed. R. Mahalingam (Switzerland: Springer International Publishing), doi: 10.1007/978-3-319-07899-1

CrossRef Full Text | Google Scholar

Murray, M. B., Cannell, M. G. R., and Smith, R. I. (1989). Date of budburst of 15 tree species in Britain following climatic warming. J. Appl. Ecol 26, 693–700. doi: 10.2307/2404093

CrossRef Full Text | Google Scholar

Pearson, K. D. (2019). Spring- and fall-flowering species show diverging phenological responses to climate in the Southeast USA. Int. J. Biometeorol. 63, 481–492. doi: 10.1007/s00484-019-01679-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, J., Wu, C., Zhang, X., Wang, X., and Gonsamo, A. (2019). Satellite detection of cumulative and lagged effects of drought on autumn leaf senescence over the Northern Hemisphere. Glob. Chang. Biol. 25, 2174–2188. doi: 10.1111/gcb.14627

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Ramos, I. M., Cambrollé, J., Hidalgo-Galvez, M. D., Matías, L., Montero-Ramírez, A., Santolaya, S., et al. (2020). Phenological responses to climate change in communities of plants species with contrasting functional strategies. Environ. Exp. Bot. 170:103852. doi: 10.1016/j.envexpbot.2019.103852

CrossRef Full Text | Google Scholar

Piao, S., Liu, Q., Chen, A., Janssens, I. A., Fu, Y., Dai, J., et al. (2019). Plant phenology and global climate change: current progresses and challenges. Glob. Chang. Biol. 25, 1922–1940. doi: 10.1111/gcb.14619

PubMed Abstract | CrossRef Full Text | Google Scholar

Polgar, C., Gallinat, A., and Primack, R. B. (2014). Drivers of leaf-out phenology and their implications for species invasions: insights from Thoreau’s Concord. New Phytol. 202, 106–115. doi: 10.1111/nph.12647

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, A. D., Black, T. A., Ciais, P., Delbart, N., Friedl, M. A., Gobron, N., et al. (2010). Influence of spring and autumn phenological transitions on forest ecosystem productivity. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 365, 3227–3246. doi: 10.1098/rstb.2010.0102

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, A. D., Hollinger, D. Y., Dail, D. B., Lee, J. T., Munger, J. W., and O’keefe, J. (2009). Influence of spring phenology on seasonal and annual carbon balance in two contrasting New England forests. Tree Physiol. 29, 321–331. doi: 10.1093/treephys/tpn040

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, A. D., Hufkens, K., Li, X., and Ault, T. R. (2019). Testing hopkins’ bioclimatic law with phenocam data. Appl. Plant Sci. 7:e1228. doi: 10.1002/aps3.1228

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, A. D., Keenan, T. F., Migliavacca, M., Ryu, Y. L., Sonnentag, O., and Toomey, M. (2013). Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. Agric. For. Meteorol. 169, 156–173. doi: 10.1016/j.agrformet.2012.09.012

CrossRef Full Text | Google Scholar

Sakkir, S., Shah, J. N., Cheruth, A. J., and Kabshawi, M. (2015). Phenology of desert plants from an arid gravel plain in eastern United Arab Emirates. J. Arid Land 7, 54–62. doi: 10.1007/s40333-014-0036-2

CrossRef Full Text | Google Scholar

Shen, M., Jiang, N., Peng, D., Rao, Y., and Tang, Y. (2020). Can changes in autumn phenology facilitate earlier green-up date of northern vegetation? Agr. For. Meteorol. 291:108077. doi: 10.1016/j.agrformet.2020.108077

CrossRef Full Text | Google Scholar

Shen, M., Tang, Y., Chen, J., Yang, X., Wang, C., Cui, X., et al. (2014). Earlier-season vegetation has greater temperature sensitivity of spring phenology in northern hemisphere. PLoS One 9:e88178. doi: 10.1371/journal.pone.0088178

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherry, R. A., Zhou, X., Gu, S., Arnone, J. A. III, Schimel, D. S., Verburg, P. S., et al. (2007). Divergence of reproductive phenology under climate warming. Proc. Natl. Acad. Sci. U.S.A. 104, 198–202. doi: 10.1073/pnas.0605642104

PubMed Abstract | CrossRef Full Text | Google Scholar

Škrk, N., Balzano, A., Repinek, Z., and Ufar, K. (2020). Inter-tree variability of autumn leaf phenology of european beech (fagus sylvatica) on a site in ljubljana, slovenia. Les 69, 5–20. doi: 10.26614/les-wood.2020.v69n02a09

CrossRef Full Text | Google Scholar

Song, Z. Q., Fu, Y. S., Du, Y. J., and Huang, Z. L. (2021). Global warming increases latitudinal divergence in flowering dates of a perennial herb in humid regions across eastern Asia. Agric. For. Meteorol. 296:108209. doi: 10.1016/j.agrformet.2020.108209

CrossRef Full Text | Google Scholar

Song, Z. Q., Fu, Y. S., Du, Y. J., Li, L., Ouyang, X. J., Ye, W. H., et al. (2020a). Flowering phenology of a widespread perennial herb shows contrasting responses to global warming between humid and non-humid regions. Funct. Ecol. 34, 1870–1881. doi: 10.1111/1365-2435.13634

CrossRef Full Text | Google Scholar

Song, Z. Q., Song, X. Q., Pan, Y. Q., Dai, K., Shou, J. J., Chen, Q., et al. (2020b). Effects of winter chilling and photoperiod on leaf-out and flowering in a subtropical evergreen broadleaved forest in China. For. Ecol. Manag. 458:117766. doi: 10.1016/j.foreco.2019.117766

CrossRef Full Text | Google Scholar

Tanino, K. K., Kalcsits, L., Silim, S., Kendall, E., and Gray, G. R. (2010). Temperature-driven plasticity in growth cessation and dormancy development in deciduous woody plants: a working hypothesis suggesting how molecular and cellular function is affected by temperature during dormancy induction. Plant Mol. Biol. 73, 49–65. doi: 10.1007/s11103-010-9610-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Thackeray, S. J., Henrys, P. A., Hemming, D., Bell, J. R., Botham, M. S., Burthe, S., et al. (2016). Phenological sensitivity to climate across taxa and trophic levels. Nature 535, 241–245. doi: 10.1038/nature18608

PubMed Abstract | CrossRef Full Text | Google Scholar

Visser, M. E. (2016). Phenology: Interactions of climate change and species. Nature 535, 236–237. doi: 10.1038/nature18905

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitasse, Y., Signarbieux, C., and Fu, Y. H. (2018). Global warming leads to more uniform spring phenology across elevations. Proc. Natl. Acad. Sci. U.S.A. 115, 1004–1008. doi: 10.1073/pnas.1717342115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, M., and Liu, X. (1979). Chinese Phenology Observation Method [M]. Beijing: Science Press.

Google Scholar

Wang, H., Ge, Q., Dai, J., and Tao, Z. (2015). Geographical pattern in first bloom variability and its relation to temperature sensitivity in the USA and China. Int. J. Biometeorol. 59, 961–969. doi: 10.1007/s00484-014-0909-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H. J., Dai, J. H., Rutishauser, T., Gonsamo, A., Wu, C. Y., and Ge, Q. (2018). Trends and variability in temperature sensitivity of lilac flowering phenology. J. Geophys. Res. Biogeogr. 123, 807–817. doi: 10.1002/2017JG004181

CrossRef Full Text | Google Scholar

Watson, J. R., Hays, C. G., Raimondi, P. T., Mitarai, S., Dong, C., McWilliams, J. C., et al. (2011). Currents connecting communities: nearshore community similarity and ocean circulation. Ecology 92, 1193–1200. doi: 10.2307/23034990

CrossRef Full Text | Google Scholar

Way, D. A., and Montgomery, R. A. (2015). Photoperiod constraints on tree phenology, performance and migration in a warming world. Plant Cell Environ. 38, 1725–1736. doi: 10.1111/pce.12431

PubMed Abstract | CrossRef Full Text | Google Scholar

Willis, C. G., Law, E., Williams, A. C., Franzone, B. F., Bernardos, R., Bruno, L., et al. (2017). CrowdCurio: an online crowd-sourcing platform to facilitate climate change studies using herbarium specimens. New Phytol. 215, 479–488. doi: 10.1111/nph.14535

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y. Y., Wang, X. J., and Silander, J. A. (2015). Deciduous forest responses to temperature, precipitation, and drought imply complex climate change impacts. Proc. Natl. Acad. Sci. U.S.A. 112, 13585–13590. doi: 10.1073/pnas.1509991112

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y. Y., Wang, X. J., Wilson, A. M., and Silander, J. A. (2018). Predicting autumn phenology: how deciduous tree species respond to weather stressors. Agric. For. Meteorol. 250-251, 127–137. doi: 10.1016/j.agrformet.2017.12.259

CrossRef Full Text | Google Scholar

Xu, L., and Chen, X. Q. (2013). Regional unified model-based leaf unfolding prediction from1960 to 2009 across northern China. Glob. Chang. Biol. 19, 1275–1284. doi: 10.1111/gcb.12095

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, K. M. (2013). Chinese Osmanthus [M]. China: China Forestry Press.

Google Scholar

Yu, H., Luedeling, E., and Xu, J. (2010). Winter and spring warming result in delayed spring phenology on the Tibetan Plateau. Proc. Natl. Acad. Sci. U.S.A. 107, 22151–22156. doi: 10.1073/pnas.1012490107

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Z. Z., Piao, S. L., Li, L. Z. X., Zhou, L. M., Ciais, P., Wang, T., et al. (2017). Climate mitigation from vegetation biophysical feedbacks during the past three decades. Nat. Clim. Chang. 7, 432–436. doi: 10.1038/nclimate3299

CrossRef Full Text | Google Scholar

Zhang, F. C., Wang, J. B., Zheng, J. H., Lin, H. K., and Hänninen, J. (2021). Chilling accumulation and photoperiod regulate rest break and bud burst in five subtropical tree species. For. Ecol. Manag. 485:118813. doi: 10.1016/j.foreco.2020.118813

CrossRef Full Text | Google Scholar

Zhang, H., Yuan, W., Liu, S., Dong, W., and Fu, Y. (2015). Sensitivity of flowering phenology to changing temperature in China. J. Geophys. Res. 120, 1658–1665. doi: 10.1002/2015JG003112

CrossRef Full Text | Google Scholar

Zhou, Z. X., Li, Y., Song, J., Ru, J. Y., Lei, L. J., Zhong, M. X., et al. (2019). Growth controls over flowering phenology response to climate change in three temperate steppes along a precipitation gradient. Agric. For. Meteorol. 274, 51–60. doi: 10.1016/j.agrformet.2019.04.011

CrossRef Full Text | Google Scholar

Zhu, W. Q., Jiang, N., Chen, G. S., Zhang, D. H., Zheng, Z. T., and Fan, D. Q. (2017). Divergent shifts and responses of plant autumn phenology to climate change on the Qinghai-Tibetan Plateau. Agric. For. Meteorol. 239, 166–175. doi: 10.1016/j.agrformet.2017.03.013

CrossRef Full Text | Google Scholar

Zhu, W. Q., Zheng, Z. T., Jiang, N., and Zhang, D. H. (2018). A comparative analysis of the spatio-temporal variation in the phenologies of two herbaceous species and associated climatic driving factors on the Tibetan Plateau. Agric. For. Meteorol. 248, 177–184. doi: 10.1016/j.agrformet.2017.09.021

CrossRef Full Text | Google Scholar

Zohner, C. M., Benito, B. M., Svenning, J., and Renner, S. S. (2016). Day length unlikely to constrain climate-driven shifts in leaf-out times of northern woody plants. Nat. Clim. Chang. 6:1120. doi: 10.1038/nclimate3138

CrossRef Full Text | Google Scholar

Zuur, A. F., Ieno, E. N., and Elphick, C. S. (2010). A protocol for data exploration to avoid common statistical problems. Methods Ecol. Evol. 1, 3–14. doi: 10.1111/j.2041-210X.2009.00001.x

CrossRef Full Text | Google Scholar