Mechanisms of changes in spring phenology

Plant phenology, particularly in spring, is shifting with warming climate (Fitter and Fitter, 2002; Post et al., 2018; Primack and Gallinat, 2016). Numerous observational studies and controlled experiments, conducted over a range of climatic and phenological conditions (Huang et al., 2020; Parmesan, 2007; Prevéy et al., 2019; Root et al., 2003; Wolkovich et al., 2012), have produced varying results, with observed phenological changes differing by research approach (Wolkovich et al., 2012), climatic region (Parmesan, 2007; Post et al., 2018; Prevéy et al., 2017; Zhang et al., 2015), and functional group (Parmesan, 2007; Root et al., 2003; Willis et al., 2010; Wolkovich et al., 2013; Zettlemoyer et al., 2019). For example, observed changes are often smaller with experimental studies (Wolkovich et al., 2012), native species (Wolkovich et al., 2013), and warmer regions (Prevéy et al., 2017). Understanding these differences (also referred to as discrepancies or mismatches) is crucial for predicting changes in plant communities and ecosystems (Fitter and Fitter, 2002; Polgar et al., 2014; Primack and Gallinat, 2016) but difficult with sensitivity analysis that is not linked to drivers of spring phenology.

Plants in the northern hemisphere require an accumulation of cool winter temperatures (winter chilling) to break dormancy and an accumulation of warm spring temperatures (spring forcing, degree days) to initiate species event (e.g., leafing or flowering) (Piao et al., 2019; Way and Montgomery, 2015). This means that changes in spring phenology result primarily from changes in spring forcing prior to species events if chilling requirements are fully met and other phenological constraints (e.g., photoperiod effect and environmental stresses) remain constant. While the direction of phenological changes (advance or delay) is set by the sign of forcing changes (+ or -), the magnitude of changes depends on quantity of forcing changes, as well as level of growth temperatures at time of species event with a warmer climate, i.e., higher growth temperatures are associated with smaller observed changes and sensitivity and lower growth temperatures with larger observed changes and sensitivity (Chu et al., 2021, 2023; Prevéy et al., 2017; Shen et al., 2014). As species often differ in timing of spring events, both forcing change and growth temperature are species-specific, reflecting changes in spring temperatures that are relevant to phenological changes of individual species (Chu et al., 2021, 2023). The effects of photoperiod (day length) and environmental stresses (e.g., drought and spring freezing) are similar to those of insufficient winter chilling in constraining the advance of spring phenology (Chen et al., 2011; Huang et al., 2019; Körner and Basler, 2010; Ma et al., 2019; Way and Montgomery, 2015). The difference between observed phenological change and that expected from forcing change and growth temperature represents phenological difference due to effects of phenological constraints and can be thereafter named phenological lag. Separating effects of different constraints is possible if soil moisture/plant water status is monitored (Huang et al., 2019; Post et al., 2022), or species-specific chilling and photoperiod effects (Man et al., 2017a, 2021a; Fu et al., 2019) or cold hardiness (Man et al., 2017b, 2021b) are known from previous research. For example, phenological lag can represent the effect of insufficient chilling when photoperiod effect and environmental stresses are constant (Chen et al., 2011).

The changes in spring phenology are often quantified with sensitivity, a parameter that reflects rates of changes but not uneven warming and responses (Beaubien and Hamann, 2011; Keenan et al., 2020; Post et al., 2018; Rafferty et al., 2020) or magnitude and underlying mechanisms of changes in spring phenology (Chu et al., 2021, 2023). Here, we present a new method to partition observed phenological change into forcing change, growth temperature, and phenological lag, based on changes in spring phenology and temperatures with control (baseline) and warmer climates. We applied this analytical approach to meta-analyze phenological changes reported in the literature and examined the contributions of forcing change, growth temperature, and phenological lag to differential responses reported previously, i.e., research approach (observational or experimental study), species origin (native or exotic), climatic region (boreal or temperate), and growth form (tree, shrub, herb, or grass).

Partitioning observed changes

Species-specific forcing change (FC) can be calculated from the difference between degree days (> 0 °C, see Man and Lu, 2010) of baseline (or control) and warmer temperatures at the time of species event (leafing or flowering) with baseline temperatures (OCi) (see Figure 1).

A diagram showing the relationships among forcing change (FC, difference in spring forcing (degree days) between baseline and warmer temperatures at species event with baseline temperatures OCi), expected response (NE, difference between climatic conditions of baseline and warmer temperatures in reaching species forcing needs with baseline temperatures, i.e., spring forcing at OCi), growth temperature (FC/NE, average temperature at species event with warmer temperatures), and phenological lag (NC, difference between expected NE and observed NO responses with warmer temperatures).

The phenological response expected under the null hypothesis that climate warming does not induce changes in phenological constraints can be determined from forcing change and species phenology with baseline temperatures (Figure 1):

where NE is the difference in the number of days (n) between baseline and warmer temperatures in reaching species forcing needs with baseline temperatures (spring forcing at OCi). Under this hypothesis, spring phenology is expected to shift from OCi with baseline temperatures to EWi with warmer temperatures. NE is typically positive, meaning earlier leafing and flowering, but can be negative if temperatures decrease over time in observational studies.

Growth temperature (TN) is the average temperature or rate of forcing accumulation within the window of expected phenological response (between OCi and EWi) and can be determined from forcing change and expected response.

The difference between expected (NE) and observed (NO) responses is phenological lag (NC) reflecting phenological difference due to effects of phenological constraints.

NC > 0 indicates increasing phenological constraints that may result from increasing chilling deficiency, photoperiod restriction, moisture stress, or risks of spring frosts; NC = 0 indicates no changes; and NC < 0 indicates decreasing constraints with warming climate.

Meta-analysis of global data

Phenological Data

We followed the guidelines of PRISMA (preferred reporting items for systematic reviews and meta-analyses) (O’Dea et al., 2021; Page et al., 2020) to build up the dataset for responses of plant spring phenology to experimental and climatic warming. We used Web of Science and Google Scholar to search experimental and observational studies on warming-induced changes in spring phenology. The following criteria were used to select papers: (a) accessible peer-reviewed articles published in scientific journals for boreal (MAT < 6 °C) and temperate (MAT ≥ 6 °C) regions in the northern hemisphere, (b) studies reporting spring leafing or flowering before Julian day 213 (July 31), (c) studies containing different climatic conditions of baseline (lower temperatures in early period of observations or unwarmed control treatment) and warmer (higher temperatures in later period of observations or warmer treatment) climates, (d) access to local temperature data (only one study was excluded due to lack of nearby temperature data within 50 km), and (e) reported phenological changes in unique locations, species, and periods of observations/warmer treatments. When publications did not contain sufficient phenological information, we accessed data directly from phenological databases (e.g., China’s National Earth System Science Data Center). When results were reported in graphical format, we extracted the required data using the distance measuring tool in Adobe Acrobat Reader DC. In experimental studies with multiple warming treatments, we selected the treatment with the smallest temperature increases to represent a likely scenario of climate change and to maintain data independence. We also excluded studies where the timing of phenological events, total temperature change, or observational time periods were not available.

In total, 66 studies were identified from 87 locations, containing 1,527 phenological responses for 980 species (Tables S1–S3). For studies reporting multiple phenological stages, we only used data for the first occurrence of the event. To account for variation in phenological responses, each study location was characterized by research approach (observation or experiment), climatic region (boreal or temperate), biogeographic origin of species (exotic or native to study area), and growth form (tree, shrub, herb, or grass).

Climate Data

Baseline and warmer temperature data for calculating forcing change, expected response, growth temperature, phenological lag, spring warming (average temperature change from Julian day 1 to 182), and long-term MAT and MAP were obtained from 9 different sources (Table S4). When temperature data were available from multiple sources, we selected the data set with the fewest missing values. For each phenological response, the number of weather stations within 50 km of data collection ranged from 1 for studies at a single location to 14 for observational studies across large geographic areas (Table S2).

Data analysis

W calculated forcing change, expected response, growth temperature, and phenological lag for each of the 1527 responses compiled. We used linear mixed-effects models to explore the differences in observed responses and phenological lags between research approaches (observations vs. experiments), species origins (native vs. exotic), and climatic regions (boreal vs. temperate), and among growth forms (trees, shrubs, herbs, and grasses). Separate analysis was conducted for each of the four fixed effects and two species events (leafing and flowering). Location and species were treated as random effects.

To assess the influences of climatic, phenological, and biological variables, we used a stepwise regression with automated combined forward and backward selection by Akaike information criterion (AIC) (Burnham and Anderson, 2002) to select the best combination of variables for predicting observed responses (NO).

The variables included were altitude (AL), latitude (LA), MAT, MAP, spring phenology (DT), growth temperature (TN), forcing change (FC), spring warming (SW), and phenological lag (NC).

Observed changes in spring phenology

Across 66 studies with 980 species and 1527 responses, observed responses (NO) ranged from −23.8 to +26.8 days for leafing and −36.3 to +49.9 days for flowering, with an average of +6.0 days for both events. These changes differed from expected responses (NE) which ranged from −3 to +27 days for leafing and −3 to +26 days for flowering and averaged +7.8 and +9.5 days, respectively. Observed phenological responses were thus smaller than the expected by 1.8 and 3.5 days, on average, for leafing and flowering, respectively.

In experimental studies, observed flowering response was non-significantly smaller than those from observational studies (1.7 days, p=0.232, Figure 2a), while phenological lag was significantly longer for both leafing (3.0 days, p=0.028) and flowering (2.7 days, p=0.027) (Figure 2b). Observed flowering response in exotic species was 2 days greater than in native species (p=0.016), likely due to differences in phenological lag (1.9 days, p=0.019). Observed leafing and flowering responses in boreal region were 4.3 and 1.9 days smaller (p=0.004, p=0.088), respectively, than in temperate region, likely due to smaller forcing changes (equivalent to 2.8 days in leafing and 1.5 days in flowering based on forcing changes and growth temperatures, Table 1 note a). Observed leafing response in grasses was 3.8 days smaller than in trees and shrubs (p=0.033), likely due to differences in forcing changes (6.0 days, Table 1 note b).

Observed responses and phenological lags (least square means ± standard errors) in leafing and flowering (left and right set of bars in each pair, respectively) by research approach (observation and experiment), species origin (native and exotic), climatic region (boreal and temperate), and growth form (tree, shrub, herb, and grass) extracted from reported plant phenological changes in spring.

Phenological lags are calculated from the differences between observed responses and those expected from species-specific changes in spring temperatures. Different letters indicate means differ significantly (p<0.05).

Mean values of climatic, phenological, and biological variables by species event (leafing and flowering), research approach (observation and experiment), species origin (native and exotic), climatic region (boreal and temperate), and growth form (tree, shrub, herb, and grass).

Growth temperature=temperatures at species event with warmer climate, forcing change =warming-induced changes in spring forcing experienced by individual species, spring warming=average spring temperature changes, and phenological lag= phenological difference due to effects of phenological constraints.

Other than phenological lag, forcing change and growth temperature were the most influential variables affecting observed responses (Table 2). The remaining variables combined (altitude, latitude, MAT, and spring warming) explained <2.5% variations in observed leafing and flowering responses.

Final stepwise regression coefficients and variable influence on observed phenological responses by Akaike information criterion (AIC).

Differential phenological responses

Observations vs. experiments

By partitioning observed changes based on drivers of spring phenology, we clarified some of the uncertainty regarding the underlying mechanisms of differences in plant phenological response to climate warming. By synthesizing global data, we found that observed flowering responses in experimental studies were non-significantly smaller than that in observational studies, but this was not the case for leafing responses; these findings differ from those based on temperature sensitivity (Wolkovich et al., 2012). The greater warming in experimental studies did not produce greater observed responses, due to longer phenological lags, which may have resulted in the reduced sensitivity (Wolkovich et al., 2012). First, experimental studies are often conducted at higher altitude sites with lower MAT, later spring phenology, and therefore longer accumulation of winter chilling. The longer phenological lags are therefore unlikely the result of insufficient winter chilling. Second, artificial warming has typically provided higher forcing changes and growth temperatures than natural climate warming (see Table 1). Greater forcing changes produce larger expected responses (observed response + phenological lag), whereas higher growth temperatures reduce expected responses (Chu et al., 2021; Prevéy et al., 2017; Wolkovich et al., 2021). The high temperatures in the artificial warming structures (Marion et al., 1997) can also reduce humidity and soil moisture (Ettinger et al., 2019; Huang et al., 2019), slowing spring development (Ganjurjav et al., 2020; Huang et al., 2019; Moore et al., 2015), particularly in early spring when temperatures in the surrounding environment are low, thus restricting soil water movement. While the different forcing changes adequately accounts for the differences in expected leafing responses (2.5 out of 3.3 days, Table 1 note c), the large differences in forcing changes but similar expected responses in flowering (8.63 and 9.07 days, Figure 2) suggests a strong effect of growth temperature (Table 1 note d) (Chu et al., 2023). Thus, warming-induced moisture stress and high growth temperatures may interact to alter observed phenological responses and temperature sensitivity in experimental studies (Wolkovich et al., 2012). Future experiments should consider humidity, soil moisture, and plant water status to better understand disparities between observational and experimental studies and relate plant phenological changes from experimental warming to those from natural climate change (Wolkovich et al., 2012).

Native vs. exotic species

Exotic species are often reported with more noticeable observed responses to warming, a trend that is generally based on flowering (Calinger et al., 2013; Willis et al., 2010; Wolkovich et al., 2013). The ecological implications of this response are well described (Polgar et al., 2014; Zettlemoyer et al., 2019), but the underlying mechanisms are not well understood. Lower growth temperatures faced by early-start exotics are likely partially the cause (Chu et al., 2021, 2023), as shown by a reverse trend for late-start exotics (Zohner and Renner, 2014). In this synthesis, phenological lag did not increase from early leafing to late flowering for exotic species, a trend that is consistent with dormancy release commonly reported on leaf and flower buds (Campoy et al., 2013; Gariglio et al., 2006; Hussain et al., 2015; Wall et al., 2008), contrary to that for native species (Figure 2b). The smaller observed response in flowering for native species is associated with longer phenological lags and accumulations of winter chilling (Table 1), indicatingg more stressful environments later in the season or higher stress sensitivity with reproductive events. Comparatively, plants starting growth early in the growing season may be less restricted by soil moisture that often recharges from winter precipitation and is depleted with increased moisture consumption over time (Sherry et al., 2007; Wolkovich et al., 2013; Zettlemoyer et al. 2019), particularly in dry climates (Man and Greenway, 2013). Moisture stress can progressively develop over the season (Piao et al. 2019; Wolkovich et al., 2013; Zettlemoyer et al., 2019) and differentially affect exotic and native species with differing spring phenology (Stuble et al., 2021; Willis et al., 2010; Wolkovich et al., 2013) The smaller observed responses in flowering of native species suggests reproductive disadvantage (IPCC, 2014; Zettlemoyer et al., 2019). The consistent leafing response, however, does not support suggestions that late-start native species would have relatively shorter active growing seasons and, therefore, a competitive growth disadvantages relative to early-start exotic species with climate warming (Fitter and Fitter, 2002; Morin et al., 2009; Primack and Gallinat, 2016). Total thermal benefits from climate warming do not differ among species with differing spring phenology (Chu et al., 2021). Given the projected increases in temperature and decreases in precipitation for many parts of the world (IPCC, 2014), future studies should assess the differences in drought sensitivity between vegetative and reproductive events and among functional groups and climatic regions.

Boreal vs. temperate region

Spatial variations along altitudinal and latitudinal gradients are often reported but not well understood (Ge et al., 2015; Morin et al., 2009; Parmesan, 2007; Zhang et al., 2015). Our analysis indicates smaller observed responses in boreal region (MAT<6 °C) due to less forcing changes (Tables 1), but not to longer phenological lag or chilling deficiency (Figure 2). The similar spring warming but less forcing changes suggests that the temperature increases in boreal region occur most frequently in the winter (Beaubien and Hamann, 2011; Shen et al., 2015) when temperatures are below freezing and don’t contribute to spring forcing (Man and Lu, 2010). The uneven warming has been reported in both natural and experimental settings (Beaubien and Hamann, 2011; Prevéy et al., 2017; Shen et al., 2015; Slaney et al., 2007; Yang et al., 2020) and supports the claim that average temperature changes do not represent forcing changes (Chu et al., 2021; Keenan et al., 2020; Wolkovich et al., 2021). When forcing changes are held constant, however, observed responses increase with altitude and latitude, as shown by the positive regression coefficients in both leafing and flowering models (Table 2). The possible mechanisms behind this may be lower growth temperature and less moisture stress at higher altitude and latitude (Rafferty et al., 2020), as shown by the negative collinearity of altitude and latitude with growth temperature and phenological lag in both leafing and flowering (data not shown). Therefore, boreal region, depending on forcing changes, can have smaller (Ge et al., 2015; Menzel et al., 2006; Rafferty et al., 2020; Shen et al., 2015) or greater (Ge et al., 2015; Parmesan, 2007; Post et al., 2018; Prevéy et al., 2017) phenological responses, spatial trends that may be different from those suggested by phenological sensitivity (Liu et al., 2019; Post et al., 2018). Consequently, climate change may not result in general phenological convergence across altitudes and latitudes (Rafferty et al., 2020; Shen et al., 2015; Tao et al., 2021), contrary to suggestions by others for small scale studies (Prevéy et al., 2017; Vitasse et al., 2018; Ziello et al., 2009).

Trees, shrubs, herbs, and grasses

Smaller observed responses in herb and grass leafing support findings from some early syntheses (König et al, 2018; Parmesan, 2007), but not those that show a greater response from non-woody plants (Ge et al., 2015; Post and Stenseth, 1999; Root et al., 2003) or no differences (Stuble et al., 2021). Herbs and grasses often leaf out earlier than woody shrubs and trees (Badeck et al., 2004; Heberling et al., 2019) and should have lower growth temperatures and therefore greater phenological responses to climate warming (Chu et al., 2021; Prevéy et al., 2017; Wolkovich et al., 2021). However, different growth forms are often studied separately on sites in different climatic regions, making comparisons challenging. Furthermore, most studies on herb and grass leafing are conducted in boreal region at high altitudes, low MAT, late spring phenology, low forcing changes, and low spring warming (Table 1); the extremely small forcing changes resulted in small observed responses despite cooler growth temperatures (Figure 2a). The observed responses with herbs and grasses would not be smaller if converted to temperature sensitivity (see Table 1 for large differences in spring warming for leafing). In contrast to leafing, herb and grass flowering studies have greater forcing changes and warmer growth temperatures (Table 1), which lowers both expected and observed responses (Table 1 note d) (Chu et al., 2021; Prevéy et al., 2017; Wolkovich et al., 2021) and potentially induces moisture stress and increases phenological lags (Ganjurjav et al., 2020; Huang et al., 2019; Moore et al., 2015; Post et al., 2022), particularly for grasses (Sherry et al., 2007; Zettlemoyer et al., 2019) that have shallower root systems (Kulmatiski and Beard, 2013; Schenk and Jackson, 2002).

Mechanistic assessment of changes in spring phenology

Averaged across all observational and experimental studies, observed responses and phenological lags are both positive, suggesting that climate warming advances spring phenology but not at rates expected from changes in spring temperatures. The positive lags are likely due to more stressful environments with warmer and drier climate (Huang et al., 2019; Sherry et al., 2007; Zettlemoyer et al., 2019). We quantified phenological lags from changes in spring phenology and temperatures, an approach that does not require species-specific chilling or forcing needs that are often unavailable (Fitter and Fitter, 2002; Morin et al., 2009) or variable in methods of estimation across species and studies (Man and Lu, 2010; Ettinger et al., 2020; Zhang et al., 2018). The relatively early stage of climate warming and association of longer phenological lags with longer accumulation of winter chilling (Table 1) suggest that climate warming is not likely to induce a general chilling shortage (Chu et al., 2023; Ettinger et al., 2020; Tao et al., 2021; Yang et al., 2020), although incidental effects on individual species or in particular conditions can occur (Fu et al., 2015). Similarly, the variable and minor effects of photoperiod (Basler and Körner, 2012; Chuine et al., 2010; Ettinger et al., 2020; Zohner et al., 2016) (insignificant spring phenology; see Table 2) often reported in studies with extreme warming scenarios (Caffarra and Donnelly, 2011; Fu et al., 2019; Laube et al., 2014) or compounded with that of growth temperatures (Chu et al., 2021), minor effects by nutrient availability (Piao et al., 2019), or incidental spring freezing (Man et al., 2009; 2021b) are not likely to cause the systemic differences in phenological lag between observational and experimental studies or between native and exotic species. Phenological lag indicates the overall lag effect and the need to investigate individual constraints.

Our analysis indicates that both quantity (forcing change) and rate (growth temperature) of spring forcing accumulation strongly influence the magnitude of change in spring phenology, consistent with the findings by Chu et al. (2021) that smaller phenological response with late season species are largely due to their higher growth temperatures. In phenological research, the role of growth temperature is not adequately recognized (Chu et al., 2021, 2023) and can be confounded with progressive increases of phenological lag with greater chilling requirements or insufficient chilling (Asse et al., 2018; Morin et al., 2009; Primack and Gallinat, 2016), photoperiod constraints (Körner and Basler, 2010; Shen et al., 2014; Way and Montgomery, 2015), or moisture stress (Piao et al., 2019; Wolkovich et al., 2013; Zettlemoyer et al., 2019) suggested for late season species. Compared to average spring temperature change that explained <0.5% of the variation in observed responses (Table 2), forcing change quantifies climate warming effect relevant to phenological change of individual species and reduces uneven warming effects among species with different spring phenology (Chu et al., 2021) and areas in different climatic regions (Post et al., 2018; Yang et al., 2020). Both forcing change and growth temperature can be readily extracted from phenological and temperature data, as demonstrated in this synthesis, but not examined in assessment of plant phenology changes in spring.

Conclusions

Differential phenological responses reported previously are mainly due to different forcing changes, growth temperatures, and environmental stresses from warmer and drier climates. The overall effects of insufficient winter chilling and photoperiod restriction are likely limited but may change with continuous progress of climate warming and advance of spring phenology (IPCC, 2014; Ettinger et al., 2020; Tao et al., 2021). Our method allows for mechanistic partitioning of observed changes by drivers of spring phenology and helps to understand changes in spring phenology and differences in plant phenological responses across different functional groups or climatic regions.

Acknowledgements

This work was supported by Ontario Ministry of Natural Resources and Forestry (Canada), Guangxi Normal University (China), Jilin Provincial Academy of Forestry Sciences (China), Shanghai Botanical Garden (China), Research Institute of Subtropical Forestry (China), and Lakehead University (Canada). Lisa Buse and Gillian Muir of Ontario Ministry of Natural Resources and Forestry and two anonymous reviewers provided constructive suggestions for improving earlier drafts of the manuscript. Phenological and temperature data provided by China’s National Earth System Science Data Center are greatly appreciated.

Additional information

Author Contributions

R.M. and J.T. conceived the general idea of the method; Y.J., X.C., X.Y., R.M., and J.T. compiled data; and Y.J., S.J.M., X.C., X.Y., R.M., J.T., and Q.L.D. contributed to the interpretations of the results and to the writing of the paper.

Data availability

Temperature data and R codes for calculating forcing changes, expected responses, growth temperature, spring warming, and phenological lags are deposited at Dryad https://doi.org/10.5061/dryad.dncjsxm9x.

Additional files

Supplementary data. Table S1. A list of references for all studies included. Table S2. A summary of study information. Table S3. MetaData containing all information used in the analysis. Table S4. Sources of weather data