Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Mixed Coniferous Broad-Leaved Forests as Road Shelter Forests: Increased Urban Traffic Noise Reduction Effects and Economic Benefits
Next Article in Special Issue
Chilling and Forcing Requirements of Wintersweet (Chimonanthus praecox L.) Flowering in China
Previous Article in Journal
Herbaceous Plant Diversity and Soil Physicochemical Properties under Different Artificial Forests in the Bulianta Core Mine, Inner Mongolia, China
Previous Article in Special Issue
Influence of Early-Season Drought on the Peak of Growing Season in China Varies by Drought Timing and Biomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shifts in Plant Phenology and Its Responses to Climate Warming in Three Temperate Cities of China during 1963–2020

1
Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
2
Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
3
Akesu National Station of Observation and Research for Oasis Agro-Ecosystem, Akesu 843017, China
4
University of Chinese Academy of Sciences, Beijing 100049, China
5
China–Pakistan Joint Research Center on Earth Sciences, Chinese Academy of Sciences—Higher Education Commission of Pakistan, Islamabad 45320, Pakistan
*
Authors to whom correspondence should be addressed.
Forests 2024, 15(10), 1712; https://doi.org/10.3390/f15101712
Submission received: 28 August 2024 / Revised: 20 September 2024 / Accepted: 26 September 2024 / Published: 27 September 2024
(This article belongs to the Special Issue Woody Plant Phenology in a Changing Climate)

Abstract

:
The advance of spring phenology and the delay of autumn phenology caused by global warming have been documented by many studies. However, most research has focused on natural areas, with limited studies conducted on phenological observations in urban environments. Here, we selected the first flowering date (FFD), first leaf date (FLD), and leaf coloring date (LCD) at three sites (Beijing, Harbin, and Mudanjiang) from the China Phenological Observation Network. We analyzed the phenological changes of 84 species between 1963–1991 and 1992–2020 to examine their response to urban warming. We then quantified the correlations and regressions between phenological events and preseason temperature. The results show the following: (1) Among the three sites, the mean FFD and FLD were earliest in Beijing, while the mean LCD occurred earliest in Harbin and latest in Beijing. (2) FFD and FLD exhibited a significant trend towards earlier occurrences at all three sites, while LCD showed a significant delay trend except for the Mudanjiang site. Specifically, at the Beijing, Harbin, and Mudanjiang sites, the mean FFD advanced by 8.32 days, 6.11 days, and 2.60 days in the latter period (p < 0.05), whereas the mean FLD advanced by 11.30 days, 7.21 days, and 5.02 days (p < 0.05), respectively. (3) In Beijing, Harbin, and Mudanjiang, both FFD and FLD were significantly negatively correlated with preseason temperature. However, no consistent relationship was observed between LCD and preseason temperature. These results enhance our understanding of the response of plant phenology to urban warming.

1. Introduction

Vegetation phenology is a sensitive indicator of climate change [1,2]. Phenological changes can alter ecosystem structure and function by leading to phenological mismatches of interacting species [3], helping the success of invasive plant species [4], and increasing the effective duration of plant assimilation [5]. Meanwhile, phenological changes also provide feedback to the climate system by influencing the water, energy, and carbon cycle [2,6]. Understanding of the phenological changes under historical climate variations could provide insights into how ecosystems might respond to future climate change. Various environmental factors have been identified as influencing vegetation phenology, including temperature [7], photoperiod [8], moisture [9], and CO2 [10]. Among all these factors, temperature reveals a primary control. Before initiating spring phenological events, plants first experience a period of endodormancy that needs to be broken by cold temperature and then ecodormancy that is accelerated by warm temperature [11]. Elevated autumn temperature was also found to slow down the senescence rate because plants receive more time for photosynthesis and nutrient reallocation under a suitable environment [12,13].
Phenological events are critical timings in plant developmental stages and are sensitive to climate change [6]. In the context of global warming, advanced spring phenological events (such as the leaf-out and flowering date) and delayed autumn phenological events (such as leaf coloring and leaf fall) have been generally reported by previous studies using ground-based phenology observation data across continents [7,14,15,16]. For spring phenological events, a comprehensive study based on more than 125,000 data series from phenological networks in 21 European countries during 1971–2000 demonstrated that the average advance of spring (summer) phenology was 2.5 d/decade [7]. In East Asia, the start of the growing season for Ginkgo biloba L. was found to occur 4 days earlier during 1953–2000 according to the phenological observation data [14]. For autumn phenological events, the leaf senescence dates were delayed by 2.6 d/decade in 1982–2011 on the basis of a meta-analysis in China [15]. Similarly, in the USA, leaf coloring was also found to be delayed by an average of 3.6 d/decade [17]. However, these observations and results were mainly derived from natural ecosystems. Considering that urban warming tends to be stronger than natural vegetation, phenological changes in cities may be more significant than those in natural ecosystems. For example, the length of the growing season in urban core areas was found to be increased by 9.9 ± 0.77 days compared to rural counterparts according to a satellite-based study in 343 Chinese cities [18]. Although there have been studies using remote sensing data to investigate the trends of phenological metrics in major cities of China [19,20], in situ observational studies in cities are still lacking.
In our study, we selected three phenological observation sites located within the urban area of Beijing, Harbin, and Mudanjiang as the study area to investigate the temporal changes of plant phenology and its response to urban warming. These sites record a long and continuous time series of first leaf dates (FLD), first flowering dates (FFD), and leaf coloring dates (LCD) from 1963 to 2020. We first divided the 1963–2020 period into two 29-year periods (1963–1991 and 1992–2020), and calculated the average FLD, FFD, and LCD for each period, as well as the differences between the two 29-year periods. Then, we conducted the correlation and regression analysis between phenophases (FLD, FFD, and LCD) and preseason temperature. In summary, our aims were (1) to identify the distinctions between the phenophases (FLD, FFD, and LCD) in different time periods (1963–1991 and 1992–2020); (2) to quantify the response of phenophases (FLD, FFD, and LCD) to preseason temperature.

2. Data and Methods

2.1. Study Site

We selected Beijing (40°2′ N, 116°33′ E), Harbin (45°70′ N, 126°63′ E), and Mudanjiang (44°43′ N, 129°67′ E) as the study area (Figure 1a) according to the principle of longer observation and better data continuity. The three stations have a temperate monsoon climate with a hot, rainy summer and a cold, dry winter. The coldest month is January, and the hottest month is July. The average temperature from 1963 to 2020 is 12.5 °C, 4.5 °C, and 4.4 °C for Beijing, Harbin, and Mudanjiang, respectively (Figure 1b). The average annual precipitation is 546 mm, 517 mm, and 533 mm for Beijing, Harbin, and Mudanjiang, respectively. The main vegetation type is deciduous broad-leaved forest.

2.2. Data Sources

All phenological data for this study are from the China Phenological Observation Network (CPON), which was established in 1963. We investigated three phenophases, including the first flowering date (FFD), first leaf date (FLD), and leaf coloring date (LCD). Observations of each phenophase conform to uniform observation criteria and guidelines [21]. FFD was defined as the date when a fixed individual has at least 3 flowers fully opens. FLD was defined as the date when a fixed individual with at least 3 leaves fully opens. LCD was defined as the date when an individual shows yellow leaves on 95%–100% of its crowns. FFD and FLD are considered representative indicators of spring phenology, while LCD represents autumn phenology [22,23]. Additionally, FLD and LCD can be interpreted as the start and end dates of the photosynthesis period, respectively [23,24]. The phenological data span from 1963 to 2020, and the effective observation years ranged between 36 and 52 years among sites because there are a few missing data in certain years due to discontinuous observations (Table A1, Table A2 and Table A3). A total of 84 species (belonging to 58 genera), including 67 species of trees, 27 species of shrubs, and 1 vine, were analyzed in three cities.
The climate data, including daily mean temperatures and precipitation, for the period from 1963 to 2020 were sourced from the China Meteorological Data Sharing Service System (https://data.cma.cn/, accessed on 27 September 2024).

2.3. Methods

2.3.1. Statistical Method of Phenological Changes

In order to analyze the impact of climate change on phenology at three sites over the past 60 years, the study period was divided into two periods (1963–1991 and 1992–2020) with the same duration. A comparison of temperature data from the Beijing, Harbin, and Mudanjiang sites between the two periods revealed that the mean monthly average temperatures at all three locations during 1992–2020 were higher than those during 1963–1991 by 1.53 °C, 1.43 °C, and 1.23 °C, respectively. The number of months with significant differences (p < 0.05) was eleven in Beijing and ten for both Harbin and Mudanjiang. Harbin and Mudanjiang experienced higher precipitation during 1992–2020 compared to that in the earlier period by an increase of 2.64 mm and 11.85 mm, respectively. However, Beijing station witnessed a decrease in precipitation by about 5.04 mm during the later period as opposed to its earlier period. The number of months with significant differences (p < 0.05) in monthly precipitation was three for Beijing station, five for Harbin station, and six for Mudanjiang (Figure 2). The significant difference in the monthly mean temperatures between the two periods provides a basis for analyzing the response of phenology to global warming.
The phenophase is represented by the day of the year (DOY). We only analyzed the species whose cumulative observation years are greater than or equal to 10 years. The average FFD, FLD, and LCD of the two periods and the whole period (1963–2020) were counted, the deviation between the two periods (later period minus earlier period) was analyzed, and the frequency distribution diagram and box plot were drawn. An unequal variance t-test was performed to test whether the deviations differed significantly from zero.

2.3.2. Analysis of the Relationship between Phenology and Preseason Temperature

Plant phenology is influenced by fluctuations in mean temperature during the period preceding phenological events [25]. Numerous studies have demonstrated a significant correlation between plant phenology and the mean temperature of the month in which the plant’s multi-year mean phenological onset occurs, as well as the two preceding months [23,26]. However, the optimum period (OP) for temperature influence cannot be limited to a single day. It is generally accepted that plant phenology is affected by temperatures from at least two weeks prior. Therefore, in this study, the OP was calculated using temperatures from 15 to 120 days before the phenological event. The OP was utilized in this study to examine the correlation between phenology and temperature, as defined by [14,23]
OP = [BP,EP],
where EP is the end date of OP (in DOY), defined as the average date of the phenophase from 1963 to 2020. BP is the beginning date of the period (in DOY). Then, the Pearson’s correlation coefficient (R) between the phenophase (FFD, FLD, or LCD) and average temperature during [BP, EP] was calculated by moving each BP from EP-15 to EP-120 days by a step length of 1 day. The period [BP, EP] with the highest absolute value of correlation coefficient was designated as the OP. The BP of LCD was restricted to occur after the summer solstice [13]. The regression slope for each phenophase (FFD, FLD, LCD) was calculated in relation to the corresponding average temperature within the OP, allowing for the quantification of the temperature sensitivity of each phenophase. The Adjusted R2 and RMSE for each phenophase (FFD, FLD, LCD) was calculated in relation to the corresponding average temperature within the OP using the Sigmoidal model (Equation (2), Figure A1) and Linear model (Equation (3), Figure A2). The significance of the correlation coefficients, regression slopes, Adjusted R2, and RMSE values were assessed using a two-sample t-test.
f ( x ) = c 1 + e ( a + b x ) + d
where c and d represent the horizontal asymptotes, while b is the growth rate parameter. The inflection point of the curve is located at x = −a/b.
y = i x + j
where i is the slope, and j is the y-intercept.

3. Results

3.1. Phenological Variation across Sites

There are significant differences in the mean FFD, FLD, and LCD of each site from 1963 to 2020 (Figure 3). The FFD range across three sites spanned 75 to 198 days, and more than 50% of the species (52.75%) were distributed in DOY 105–150. The FFD of the Beijing site was 24.77 days and 26.58 days earlier than Harbin and Mudanjiang. The FLD range was 75–197 days, and over 80% of species (80.90%) were distributed in DOY 89–125. The mean FLD was earliest in Beijing, followed by Mudanjiang, and the latest in Harbin. The mean LCD of Harbin was slightly earlier than that of Mudanjiang, and the mean LCD of Beijing was 31.14 days later than that of Harbin.

3.2. Phenological Change between 1963–1991 and 1992–2020

Based on the phenological differences observed across all species between the two periods (1963–1991 and 1992–2020), both FFD and FLD advanced at all three sites, while LCD was delayed at all sites except for Mudanjiang (Figure 4). Specifically, the median FFD at Beijing, Harbin, and Mudanjiang in 1992–2020 was earlier by 8.10 days, 2.99 days, and 2.18 days compared to 1963–1991, respectively. The median FLD in 1992–2020 was earlier by 11.15 days, 7.54 days, and 1.15 days, respectively, at these sites compared to the period 1963–1991. The median LCD at Beijing and Harbin in 1992–2020 was later by 5.10 and 7.05 days, respectively, compared to 1963–1991, while it was earlier by 0.17 days at Mudanjiang.
Changes in the frequency distribution of the three phenological phases from 1992–2020 and 1963–1991 reveal an overall advancement in FFD (Figure 5) and FLD (Figure 6), while a delay was observed in LCD at Beijing, Harbin, and Mudanjiang (Figure 7). Specifically, Beijing (85.11%) recorded the highest proportion of species with significantly advanced FFD, followed by Mudanjiang (57.69%), and the lowest proportion was found in Harbin (22.22%). The magnitude of advancement followed the order of Beijing (8.32 days), Harbin (6.11 days), and Mudanjiang (2.60 days). Similarly, for FLD, Beijing (2.60 days) had the highest proportion of significantly advanced species, followed by Harbin (75.00%), and Mudanjiang (8.70%) had the lowest. The magnitude of advancement also followed the same order of Beijing (11.30 days), Harbin (7.21 days), and Mudanjiang (5.02 days). In contrast, Harbin (81.82%) had the highest proportion of species with significantly delayed LCD, followed by Beijing (46.34%), while Mudanjiang showed no clear pattern of change. The magnitude of delay in LCD was greatest in Harbin (11.66 days), followed by Beijing (8.43 days) and Mudanjiang (0.92 days). These findings further underscore the notable advancing trend in spring and summer phenophases and the delaying trend in autumn phenophases.

3.3. Relationship between Phenophases and Temperature

Figure 8 shows the OP corresponding to vegetation phenology at different sites, revealing that FLD and FFD are generally significantly correlated with OP temperature (p < 0.05). In contrast, the correlation between LCD and OP temperature is not significant for most species (Figure 6). The OP duration of FFD ranged from 17 days to 120 days across the three sites, with the mean OP duration at Harbin (50.81 days), Beijing (51.55 days), and Mudanjiang (59.10 days) in ascending order. For FLD, the OP duration varied from 16 days to 120 days across all sites, with a mean OP of 42.43 days, 42.98 days, and 51.20 days at Mudanjiang, Beijing, and Harbin, respectively. For LCD, the OP duration ranged from 16 days to 111 days, with a mean OP duration of 38.33 days, 46.00 days, and 46.80 days at Harbin, Mudanjiang, and Beijing, respectively.
A comparison of the Adjusted R2 values between the Sigmoidal and Linear model revealed that, with the exception of the significant difference (p < 0.05) in Adjusted R2 for FFD and OP temperatures at the Mudanjiang site, the differences in Adjusted R2 between the two models were not statistically significant in other cases (Figure A3). Additionally, no significant differences were observed between the RMSEs derived from the Linear and Sigmoidal models. Notably, the RMSEs for all phenological periods (FFD, FLD, and LCD) in relation to OP temperatures were consistently smaller for the Linear model across all three sites (Figure A4). Consequently, the results of the Linear model were employed in subsequent analyses.
An analysis of the correlation between vegetation phenology and OP temperature revealed that at Beijing, Harbin, and Mudanjiang, FFD and FLD exhibited negative correlations with OP temperature, whereas the correlation between LCD and OP temperature varied significantly among sites (Figure 7). Specifically, the mean correlation coefficients between FFD and OP temperature were −0.76 ± 0.23 (mean ± SD), −0.69 ± 0.22, and −0.53 ± 0.27, respectively, with 89.36%, 88.89%, and 38.46% of species showing significance at the 5% level. For FLD, the mean correlation coefficients were −0.79 ± 0.13, −0.71 ± 0.12, and −0.81 ± 0.10, with 97.83%, 100.00%, and 100.00% of species significant at a 0.05 level. In contrast, the mean correlation coefficients between LCD and OP temperature were 0.21 ± 0.18, 0.49 ± 0.22, and −0.43 ± 0.19, respectively, with only 12.20%, 54.55%, and 15.38% of species displaying a significant correlation coefficient (p < 0.05).
The temperature sensitivity (regression slope) of vegetation phenology was significantly different among sites and species (Figure 9). Significant advancements predominantly characterized the responses of FFD and FLD to OP temperatures. At the Beijing, Harbin, and Mudanjiang sites, the temperature sensitivity of FFD ranged from −0.52 to −6.99 days °C−1, −1.94 to −7.17 days °C−1, and −2.33 to −18.02 days °C−1, respectively. For FLD, the temperature sensitivity ranged from −1.44 to −6.99 days °C−1, −2.42 to −4.47 days °C−1, and −2.44 to −3.95 days °C−1, respectively. In contrast, the response of LCD to OP temperatures did not exhibit a consistent pattern among sites. At Beijing, Harbin, and Mudanjiang, the proportions of species with positive slopes (indicating advancement with increasing OP temperatures) were 71.43%, 91.67%, and 28.57%, respectively. The magnitudes of these responses ranged from 0.0037 to 4.16, 0.73 to 6.36, and 0.57 to 1.35 days °C−1, respectively. For the remaining species, LCD advanced by 0.35–2.31 days, 2.79–2.79 days, and 0.85–7.50 days, respectively, with a 1 °C increase in OP temperatures. Overall, the mean temperature sensitivity across all species at the three sites was −5.90, −3.28, and 0.35 days °C−1 for FFD, FLD, and LCD, respectively.

4. Discussion

In terms of spatial pattern, our results support the bioclimatic law proposed by Hopkins, which states that spring phenology would be delayed but autumn phenology would be advanced by latitude north [27]. Concretely, the median FLD and FFD in Beijing were earlier than in Mudanjiang, followed by Harbin in both observation periods. This latitudinal gradient of spring phenology might be attributed to the response of plants to varying environmental conditions. However, the gradient of FLD and FFD across latitudes between Beijing and Harbin (4.36 days for FFD and 4.53 days for FLD per 1° N) was higher than Hopkins’ hypothesis (4 days per 1° N). Compared to the other cities of China (2.46 to 3.48 days per 1° N), North America (2.62 days per 1° N), or the whole Northern Hemisphere (0.7 to greater than 1.0 days per 1° N), the spring phenology has been more dramatically affected by latitude gradient in this study [28,29,30]. In addition, the LCD in Beijing was significantly earlier than in the other two cities, but the disparity between the other two cities was not conspicuous. These inconsistent quantifications in geographic gradient might first be caused by the difference across observed species themselves [31]. Another reason is that spring phenology could be dramatically advanced and autumn phenology could be delayed by high urbanization intensity [19].
Consistent with previous studies, widespread advancements in the FLD and FFD and inapparent delay in LCD were found in this study [15,32,33,34]. However, the magnitude of the shifts in spring and autumn phenology varied among stations. For example, our results showed that the advancing trends of spring phenology in Beijing across the two periods (deviation in phenophases/29 years) were the most significant (0.28 days year−1 for FFD and 0.38 days year−1 for FLD), which were close to a previous study in China (0.22 to 0.57 days year−1) but slightly higher than those observed in Europe (0.24 days year−1) [15,34]. Nevertheless, the greatest changes in LCD were observed in Harbin (0.40 days year−1 delay), which was much higher than both the previous study of China (0.19–0.25 delay per year) and Europe (0.036 delay per year) [15,34]. This is because the phenological changes in urban areas would be more dramatic than in natural environments, with 0.50 days year−1 much earlier in spring and 0.78 days year−1 much later in autumn [35]. It is worth noting that the phenological shifts in Mudanjiang were more slight than in other cities (0.08 days year−1 for FLD, 0.04 days year−1 for FFD, and 0.01 days year−1 for LCD). This may be attributed to the weaker urban heat island effect in Mudanjiang compared to the other two cities.
Based on the linear regression between the preseason temperature and phenophases, we calculated the temperature sensitivity of spring and autumn phenology. Both temperature sensitivities of FLD and FFD in all three sites were negative overall, which is consistent with the results in Europe and North America [15,36]. In a previous study, spring temperature sensitivity in most of the regions of the mid- and high-latitude Northern Hemisphere showed positive clines across latitudes [37]. However, in this study, there is no evidence to suggest the temperature sensitivity of spring phenology changes along the latitudinal gradient. It indicates that spring temperature sensitivity for a single phenology site might be affected by multiple factors such as biodiversity, radiation, precipitation, soil nitrogen, etc. [38]. In our study area, the Beijing site is the Summer Palace, and the Harbin site is the Heilongjiang Forest Botanical Garden. Precipitation and drought events were not considered, as both sites are artificially irrigated. Furthermore, Figure 2 indicates that there are few significant differences in mean precipitation between 1963–1991 and 1992–2020. Consequently, it is not feasible to analyze the effect of precipitation on phenology in our study. Additionally, some studies suggest that the spring phenology of most vegetation is not sensitive to photoperiod [39,40], and the start of growing season for vegetation in the Northern Hemisphere is not significantly affected by radiation intensity [41]. Therefore, this study did not incorporate the effects of photoperiod or radiation into the analysis of phenological patterns.
In addition, different from the previous study, i.e., plants that have earlier phenophases are more sensitive to temperature than later plants, our results did not show any significant relationship between temperature sensitivities of spring phenology and the timing of phenophases [42]. Unlike spring phenology, the direction of LCD response to preseason temperature was not consistent across sites, where the sensitivity of LCD in Beijing and Harbin was positive, but in Mudanjiang, it was negative. Due to the different urbanization processes, autumn phenology in North China cities exhibited a greater delay (3.1 days to less than 1 day) than in Northeast China [18]. The other reason could be the inconsistent standard in LCD recording caused by subjective observations, which introduced potential errors in the linear regression model [43]. Further, autumn phenological changes are affected by complex factors, and the leaf coloring response may not be linear [44].

5. Conclusions

By utilizing FFD, FLD, and LCD data from the China Phenological Observation Network and daily temperature data, we analyzed the phenological changes in three cities (Beijing, Harbin, and Mudanjiang) between 1963–1991 and 1992–2020 to examine their response to urban warming. The results show the following: (1) FFD and FLD occurred earliest in Beijing, while LCD occurred earliest at Mudanjiang and latest in Beijing. (2) In Beijing, Harbin, and Mudanjiang, the median FFD (FLD) advanced by 8.10 (11.15), 2.99 (7.54), and 2.18 (1.15) days, respectively, during the period 1992–2020 compared to 1963–1991. Meanwhile, the median LCD was delayed by 5.10 days in Beijing and 7.05 days in Harbin but advanced by 0.17 days in Mudanjiang. (3) In Beijing, Harbin, and Mudanjiang, both FFD and FLD were significantly negatively correlated with preseason temperature. However, no consistent relationship was observed between LCD and preseason temperature. Our study suggests that vegetation phenology in cities will experience more pronounced changes than in natural environments due to the urban heat island effect.

Author Contributions

Conceptualization, M.Z.; software, W.L.; validation, C.G. and W.B.; writing—original draft preparation, L.C.; writing—review and editing, S.L.; visualization, Y.H.; supervision, J.D. and X.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant no. 2023YFF1303804), National Natural Science Foundation of China (grant no. 42271062), and the Natural Science Fund for Distinguished Young Scholars of Xinjiang Uygur Autonomous Region (grant no. 2022D01E02).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Summary of phenological data from Beijing stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
Table A1. Summary of phenological data from Beijing stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
SpeciesFamilyLFNFFD1NFFD2MFFDNFLD1NFLD2MFLDNLCD1NLCD2MLCD
1Fraxinus chinensisFraxinusT18244/1924264/14222510/19
2Castanea mollissimaCastaneaT26255/3026264/19252511/2
3Syringa reticulata subsp. amurensisSyringaT22195/1824194/2211811/5
4Prunus persicaDuplexPrunusT18214/817204/11121611/3
5Platycladus orientalisPlatycladusT22233/28------
6Ailanthus altissimaAilanthusT24205/2425224/18172110/22
7Robinia pseudoacaciaRobiniaT25265/225264/14232310/30
8Pyrus betulifoliaPyrusT17214/1417194/8162010/28
9Juniperus chinensisJuniperusT23214/1------
10Salix matsudanaSalixT26194/224193/29201811/12
11Albizia julibrissinAlbiziaT26186/926185/1221411/1
12Styphnolobium japonicumStyphnolobiumT23207/1423204/25181311/14
13Rosa xanthinaRosaS14204/2314214/2101811/10
14Cotinus coggygria var. cinereusCotinusT17194/2116184/12171910/25
15Rhodotypos scandensRhodotyposT17184/1617174/7161111/13
16Populus × canadensisPopulusT23263/2625244/7232110/31
17Diospyros lotusDiospyrosT22255/1723254/17162010/28
18Forsythia suspensaForsythiaS18213/2818194/8151911/4
19Styphnolobium japonicumPendulaStyphnolobiumT16197/1118214/13---
20Koelreuteria paniculataKoelreuteriaT16195/3119214/11172010/27
21Populus tomentosaPopulusT22263/1721264/9141811/6
22Paeonia × suffruticosaPaeoniaS26254/2025263/31192310/29
23Paulownia fortuneiPaulowniaT16164/2015124/23---
24Catalpa bungeiCatalpaT18184/2719194/12---
25Prunus × yedoensisPrunusT26204/1426194/11251910/31
26Morus albaMorusS24234/2224234/20212110/30
27Prunus davidianaPrunusT27263/2227264/1181911/6
28Punica granatumPunicaT16205/2717214/23131910/30
29Diospyros kakiDiospyrosT23235/1423234/15212210/29
30Metasequoia glyptostroboidesMetasequoiaT---21204/11192011/3
31Viburnum farreriViburnumT16113/3014103/29---
32Salix matsudanaPendulaSalixT24204/626213/30242011/13
33Firmiana simplexFirmianaT23236/2024254/27212410/29
34Malus × micromalusMalusT26214/1326213/31201711/11
35Populus simoniiPopulusT13154/314154/6111411/4
36Prunus armeniacaPrunusT26264/224254/13222311/6
37Ginkgo bilobaGinkgoT---16204/10162010/30
38Jasminum nudiflorumJasminumT10193/17------
39Pinus tabuliformisPinusT25214/2625224/26---
40Ulmus pumilaUlmusT25183/1626184/8191811/11
41Prunus trilobaPrunusT21264/620254/8202410/30
42Yulania denudataYulaniaT26213/3026214/12251911/1
43Sorbaria sorbifoliaSorbariaS15186/816203/30141411/8
44Syringa oblataSyringaS24264/1024263/31202111/3
45Cercis chinensisCercisS24264/1423264/14222510/28
46Wisteria sinensisWisteriaV24264/2124264/21222311/8
47Lagerstroemia indicaLagerstroemiaS22257/524254/26212410/28
48Yulania liliifloraYulaniaS26204/425204/11241810/29
49Prunus × subhirtellaPrunusT16164/1015154/11151510/30
Table A2. Summary of phenological data from Harbin stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
Table A2. Summary of phenological data from Harbin stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
SpeciesFamilyLFNFFD1NFFD2MFFDNFLD1NFLD2MFLDNLCD1NLCD2MLCD
1Syringa reticulata subsp. amurensisSyringaT12186/712184/2312109/10
2Forsythia mandschuricaForsythiaS11184/1911185/211129/28
3Philadelphus schrenkiiPhiladelphusS14185/2814184/29121010/5
4Pinus koraiensisPinusT---10176/1---
5Juglans mandshuricaJuglansT17175/1817175/10---
6Viburnum opulus subsp. calvescensViburnumS13185/2513184/28131710/8
7Lonicera maackiiLoniceraS16185/2416184/30---
8Viburnum burejaeticumViburnumS13185/1713184/26131310/9
9Morus albaMorusS11185/2012185/1812159/27
10Prunus sibiricaPrunusT15185/115185/10---
11Caragana arborescensCaraganaS12185/1912185/8---
12Fraxinus mandshuricaFraxinusT11185/314185/1314159/29
13Xanthoceras sorbifoliumXanthocerasT11185/1911185/17---
14Flueggea suffruticosaFlueggeaS10176/1910185/18101610/2
15Ulmus pumilaUlmusT17184/1617185/3161210/8
16Prunus trilobaPrunusT16185/216185/7141410/2
17Pinus sylvestris var. mongolicaPinusT---10165/28---
18Sorbaria sorbifoliaSorbariaS15167/917184/18---
19Catalpa ovataCatalpaT12166/2413185/21---
20Syringa oblataSyringaS18185/818185/2141110/13
Table A3. Summary of phenological data from Mudanjiang stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
Table A3. Summary of phenological data from Mudanjiang stations investigated in this study. NFFD1, NFLD1, NLCD1, number of observation years during the period 1963–1991; NFFD2, NFLD2, NLCD2, number of observation years during the period 1992–2020; MFLD, mean timing of first leaf dates; MLCD, mean timing of leaf coloring date; MFFD, mean timing of first flowering date (all from 1963 to 2020); LF, life form; T, tree species; S, shrub species; V, vine species. The slash denotes that the data are less than 10 years old and have been excluded.
SpeciesFamilyLFNFFD1NFFD2MFFDNFLD1NFLD2MFLDNLCD1NLCD2MLCD
1Acer tataricum subsp. ginnalaAcerT15285/2414275/214299/24
2Prunus padusPrunusT15285/714284/2211289/28
3Salix babylonicaSalixT13294/3013274/23112610/21
4Picea koraiensisPiceaT---11285/8---
5Lespedeza bicolorLespedezaT14287/6---10239/20
6Styphnolobium japonicumStyphnolobiumT10107/711175/1710119/24
7Phellodendron amurensePhellodendronT11296/7------
8Prunus salicinaPrunusT16275/612255/4111810/3
9Spiraea salicifoliaSpiraeaS13226/1812184/21---
10Rubus komaroviiRubusS11266/411274/30---
11Malus baccataMalusT15285/1214264/24112210/3
12Crataegus pinnatifida var. majorCrataegusT16285/2715284/2911269/28
13Rhamnus davuricaRhamnusS11275/2610255/1---
14Caragana arborescensCaraganaS11285/16------
15Acer saccharumAcerT12294/2611284/27---
16Euonymus alatusEuonymusS14295/1913284/2911199/29
17Populus pseudosimoniiPopulusT10244/2710255/2---
18Populus simoniiPopulusT14204/3014265/5---
19Rhododendron dauricumRhododendronS10294/29------
20Larix gmeliniiLarixT11295/1111294/26102910/10
21Prunus pseudocerasusPrunusT15205/413185/1111610/12
22Ulmus pumilaUlmusT12254/2211215/8---
23Prunus trilobaPrunusT13125/312114/30111210/5
24Quercus dentataQuercusT14235/1314235/312199/30
25Sorbaria sorbifoliaSorbariaS12187/1612194/13---
26Corylus heterophyllaCorylusS---10285/1---
27Catalpa ovataCatalpaT11276/3011275/16---
28Tilia amurensisTiliaT11286/28------

Appendix B

Figure A1. Sigmoidal model analysis of FFD and OP temperature relationship in Mudanjiang.
Figure A1. Sigmoidal model analysis of FFD and OP temperature relationship in Mudanjiang.
Forests 15 01712 g0a1
Figure A2. Linear model analysis of FFD and OP temperature relationship in Mudanjiang.
Figure A2. Linear model analysis of FFD and OP temperature relationship in Mudanjiang.
Forests 15 01712 g0a2

Appendix C

Figure A3. Adjusted R2 comparison of Sigmoidal and Linear models for vegetation phenology (FFD, FLD, LCD) in relation to OP temperature in Beijing, Harbin, and Mudanjiang.
Figure A3. Adjusted R2 comparison of Sigmoidal and Linear models for vegetation phenology (FFD, FLD, LCD) in relation to OP temperature in Beijing, Harbin, and Mudanjiang.
Forests 15 01712 g0a3
Figure A4. RMSE comparison of Sigmoidal and Linear models for vegetation phenology (FFD, FLD, LCD) in relation to OP temperature in Beijing, Harbin, and Mudanjiang.
Figure A4. RMSE comparison of Sigmoidal and Linear models for vegetation phenology (FFD, FLD, LCD) in relation to OP temperature in Beijing, Harbin, and Mudanjiang.
Forests 15 01712 g0a4

References

  1. Cleland, E.E.; Chuine, I.; Menzel, A.; Mooney, H.A.; Schwartz, M.D. Shifting plant phenology in response to global change. Trends Ecol. Evol. 2007, 22, 357–365. [Google Scholar] [CrossRef] [PubMed]
  2. Richardson, A.D.; Keenan, T.F.; Migliavacca, M.; Ryu, Y.; Sonnentag, O.; Toomey, M. Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. Agric. For. Meteorol. 2013, 169, 156–173. [Google Scholar] [CrossRef]
  3. Visser, M.E.; Gienapp, P. Evolutionary and demographic consequences of phenological mismatches. Nat. Ecol. Evol. 2019, 3, 879–885. [Google Scholar] [CrossRef] [PubMed]
  4. Wolkovich, E.M.; Cleland, E.E. The phenology of plant invasions: A community ecology perspective. Front. Ecol. Environ. 2011, 9, 287–294. [Google Scholar] [CrossRef]
  5. Piao, S.; Friedlingstein, P.; Ciais, P.; Viovy, N.; Demarty, J. Growing season extension and its impact on terrestrial carbon cycle in the Northern Hemisphere over the past 2 decades. Glob. Biogeochem. Cycles 2007, 21, GB3018. [Google Scholar] [CrossRef]
  6. Piao, S.; Liu, Q.; Chen, A.; Janssens, I.A.; Fu, Y.; Dai, J.; Liu, L.; Lian, X.; Shen, M.; Zhu, X. Plant phenology and global climate change: Current progresses and challenges. Glob. Chang. Biol. 2019, 25, 1922–1940. [Google Scholar] [CrossRef]
  7. Menzel, A.; Sparks, T.H.; Estrella, N.; Koch, E.; Aasa, A.; Ahas, R.; Alm-Kübler, K.; Bissolli, P.; Braslavská, O.g.; Briede, A.; et al. European phenological response to climate change matches the warming pattern. Glob. Chang. Biol. 2006, 12, 1969–1976. [Google Scholar] [CrossRef]
  8. Laube, J.; Sparks, T.; Estrella, N.; Höfler, J.; Ankerst, D.; Menzel, A. Chilling outweighs photoperiod in preventing precocious spring development. Glob. Chang. Biol. 2014, 20, 170–182. [Google Scholar] [CrossRef]
  9. Peng, J.; Wu, C.; Zhang, X.; Wang, X. Satellite detection of cumulative and lagged effects of drought on autumn leaf senescence over the Northern Hemisphere. Glob. Chang. Biol. 2019, 25, 2174–2188. [Google Scholar] [CrossRef]
  10. Reyes-Fox, M.; Steltzer, H.; Trlica, M.J.; McMaster, G.S.; Andales, A.A.; LeCain, D.R.; Morgan, J.A. Elevated CO2 further lengthens growing season under warming conditions. Nature 2014, 510, 259–262. [Google Scholar] [CrossRef]
  11. Lang, G.; Early, J.; Martin, G.; Darnell, R. Endo-, Para-, and Ecodormancy: Physiological Terminology and Classification for Dormancy Research. HortScience 1987, 22, 371–377. [Google Scholar] [CrossRef]
  12. Guo, Y.; Ren, G.; Zhang, K.; Li, Z.; Miao, Y.; Guo, H. Leaf senescence: Progression, regulation, and application. Mol. Hortic. 2021, 1, 5. [Google Scholar] [CrossRef] [PubMed]
  13. Zohner, C.; Mirzagholi, L.; Renner, S.; Mo, L.; Rebindaine, D.; Bucher, R.; Palouš, D.; Vitasse, Y.; Fu, Y.; Stocker, B.; et al. Effect of climate warming on the timing of autumn leaf senescence reverses after the summer solstice. Science 2023, 381, eadf5098. [Google Scholar] [CrossRef] [PubMed]
  14. Matsumoto, K.; Ohta, T.; Irasawa, M.; Nakamura, T. Climate change and extension of the Ginkgo biloba L. Growing season in Japan. Glob. Chang. Biol. 2003, 9, 1634–1642. [Google Scholar] [CrossRef]
  15. Ge, Q.; Wang, H.; Rutishauser, T.; Dai, J. Phenological response to climate change in China: A meta-analysis. Glob. Chang. Biol. 2015, 21, 265–274. [Google Scholar] [CrossRef]
  16. Pearse, W.D.; Stemkovski, M.; Lee, B.R.; Primack, R.B.; Lee, S.D. Consistent, linear phenological shifts across a century of observations in South Korea. New Phytol. 2023, 239, 824–829. [Google Scholar] [CrossRef]
  17. Jeong, S.J.; Medvigy, D. Macroscale prediction of autumn leaf coloration throughout the continental United States. Glob. Ecol. Biogeogr. 2014, 23, 1245–1254. [Google Scholar] [CrossRef]
  18. Jia, W.; Zhao, S.; Zhang, X.; Liu, S.; Henebry, G.M.; Liu, L. Urbanization imprint on land surface phenology: The urban-rural gradient analysis for Chinese cities. Glob. Chang. Biol. 2021, 27, 2895–2904. [Google Scholar] [CrossRef]
  19. Yang, L.; Zhao, S.; Liu, S. Urban environments provide new perspectives for forecasting vegetation phenology responses under climate warming. Glob. Chang. Biol. 2023, 29, 4383–4396. [Google Scholar] [CrossRef]
  20. Zhou, D.; Zhao, S.; Zhang, L.; Liu, S. Remotely sensed assessment of urbanization effects on vegetation phenology in China’s 32 major cities. Remote Sens. Environ. 2016, 176, 272–281. [Google Scholar] [CrossRef]
  21. Wan, M.; Liu, X. China’s National Phenological Observational Criterion; Science: Beijing, China, 1979. [Google Scholar]
  22. Fu, Y.H.; Geng, X.; Chen, S.; Wu, H.; Hao, F.; Zhang, X.; Wu, Z.; Zhang, J.; Tang, J.; Vitasse, Y.; et al. Global warming is increasing the discrepancy between green (actual) and thermal (potential) seasons of temperate trees. Glob. Chang. Biol. 2023, 29, 1377–1389. [Google Scholar] [CrossRef] [PubMed]
  23. Dai, J.; Wang, H.; Ge, Q. Multiple phenological responses to climate change among 42 plant species in Xi’an, China. Int. J. Biometeorol. 2013, 57, 749–758. [Google Scholar] [CrossRef] [PubMed]
  24. Menzel, A. Plant phenological anomalies in Germany and their relation to air temperature and NAO. Clim. Chang. 2003, 57, 243–263. [Google Scholar] [CrossRef]
  25. Xu, Y.J.; Zhong, S.Y.; Dai, J.H.; Tao, Z.X.; Wang, H.J. Changes in flowering phenology of plants and their model simulation in Mudanjiang, China. Geogr. Res. 2017, 36, 779–789. [Google Scholar]
  26. Askeyev, O.V.; Sparks, T.H.; Askeyev, I.V.; Tishin, D.V.; Tryjanowski, P. East versus West: Contrasts in phenological patterns? Glob. Ecol. Biogeogr. 2010, 19, 783–793. [Google Scholar] [CrossRef]
  27. Hopkins, A.D. Bioclimatics: A Science of Life and Climate Relations (No. 280); US Department of Agriculture: Washington, DC, USA, 1938. [Google Scholar]
  28. Zhu, W.; Tian, H.; Xu, X.; Pan, Y.; Chen, G.; Lin, W. Extension of the growing season due to delayed autumn over mid and high latitudes in North America during 1982–2006. Glob. Ecol. Biogeogr. 2012, 21, 260–271. [Google Scholar] [CrossRef]
  29. Dai, J.; Wang, H.; Ge, Q. The spatial pattern of leaf phenology and its response to climate change in China. Int. J. Biometeorol. 2014, 58, 521–528. [Google Scholar] [CrossRef]
  30. Shen, M.; Cong, N.; Cao, R. Temperature sensitivity as an explanation of the latitudinal pattern of green-up date trend in Northern Hemisphere vegetation during 1982–2008. Int. J. Climatol. 2015, 35, 3707–3712. [Google Scholar] [CrossRef]
  31. Zeng, Z.A.; Wolkovich, E.M. Weak evidence of provenance effects in spring phenology across Europe and North America. New Phytol. 2024, 242, 1957–1964. [Google Scholar] [CrossRef]
  32. Fu, Y.H.; Piao, S.; Op de Beeck, M.; Cong, N.; Zhao, H.; Zhang, Y.; Menzel, A.; Janssens, I.A. Recent spring phenology shifts in western Central Europe based on multiscale observations. Glob. Ecol. Biogeogr. 2014, 23, 1255–1263. [Google Scholar] [CrossRef]
  33. Keenan, T.F.; Gray, J.; Friedl, M.A.; Toomey, M.; Bohrer, G.; Hollinger, D.Y.; Munger, J.W.; O’Keefe, J.; Schmid, H.P.; Wing, I.S.; et al. Net carbon uptake has increased through warming-induced changes in temperate forest phenology. Nat. Clim. Chang. 2014, 4, 598–604. [Google Scholar] [CrossRef]
  34. Menzel, A.; Yuan, Y.; Matiu, M.; Sparks, T.; Scheifinger, H.; Gehrig, R.; Estrella, N. Climate change fingerprints in recent European plant phenology. Glob. Chang. Biol. 2020, 26, 2599–2612. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, L.; Zhao, S. A stronger advance of urban spring vegetation phenology narrows vegetation productivity difference between urban settings and natural environments. Sci. Total Environ. 2023, 868, 161649. [Google Scholar] [CrossRef] [PubMed]
  36. Seyednasrollah, B.; Young, A.; Li, X.; Milliman, T.; Ault, T.; Frolking, S.; Friedl, M.; Richardson, A. Sensitivity of Deciduous Forest Phenology to Environmental Drivers: Implications for Climate Change Impacts Across North America. Geophys. Res. Lett. 2020, 47, e2019GL086788. [Google Scholar] [CrossRef]
  37. Gao, M.; Wang, X.; Meng, F.; Liu, Q.; Li, X.; Zhang, Y.; Piao, S. Three-dimensional change in temperature sensitivity of northern vegetation phenology. Glob. Chang. Biol. 2020, 26, 5189–5201. [Google Scholar] [CrossRef] [PubMed]
  38. Shen, P.; Wang, X.; Zohner, C.M.; Peñuelas, J.; Zhou, Y.; Tang, Z.; Xia, J.; Zheng, H.; Fu, Y.; Liang, J.; et al. Biodiversity buffers the response of spring leaf unfolding to climate warming. Nat. Clim. Chang. 2024, 14, 863–868. [Google Scholar] [CrossRef]
  39. Zohner, C.M.; Benito, B.M.; Svenning, J.-C.; Renner, S.S. Day length unlikely to constrain climate-driven shifts in leaf-out times of northern woody plants. Nat. Clim. Chang. 2016, 6, 1120–1123. [Google Scholar] [CrossRef]
  40. Wang, H.; Wang, H.; Ge, Q.; Dai, J. The Interactive Effects of Chilling, Photoperiod, and Forcing Temperature on Flowering Phenology of Temperate Woody Plants. Front. Plant Sci. 2020, 11, 443. [Google Scholar] [CrossRef]
  41. Descals, A.; Verger, A.; Yin, G.; Filella, I.; Fu, Y.H.; Piao, S.; Janssens, I.A.; Penuelas, J. Radiation-constrained boundaries cause nonuniform responses of the carbon uptake phenology to climatic warming in the Northern Hemisphere. Glob. Chang. Biol. 2023, 29, 719–730. [Google Scholar] [CrossRef]
  42. Kopp, C.W.; Neto-Bradley, B.M.; Lipsen, L.P.J.; Sandhar, J.; Smith, S. Herbarium records indicate variation in bloom-time sensitivity to temperature across a geographically diverse region. Int. J. Biometeorol. 2020, 64, 873–880. [Google Scholar] [CrossRef]
  43. Gallinat, A.S.; Primack, R.B.; Wagner, D.L. Autumn, the neglected season in climate change research. Trends Ecol. Evol. 2015, 30, 169–176. [Google Scholar] [CrossRef] [PubMed]
  44. Gao, C.; Wang, H.; Ge, Q. Interpretable machine learning algorithms to predict leaf senescence date of deciduous trees. Agric. For. Meteorol. 2023, 340, 109623. [Google Scholar] [CrossRef]
Figure 1. Summary of study area. (a) Map of observation sites. (b) Monthly mean temperature and total precipitation of the 3 cities (1963–2020). The curves represent the monthly mean temperature, and the bars represent the monthly total precipitation.
Figure 1. Summary of study area. (a) Map of observation sites. (b) Monthly mean temperature and total precipitation of the 3 cities (1963–2020). The curves represent the monthly mean temperature, and the bars represent the monthly total precipitation.
Forests 15 01712 g001
Figure 2. Monthly mean temperature and precipitation of the 3 stations in 1963–1991 (blue bar and curves) and 1992–2020 (red bar and curves). (a) Beijing; (b) Harbin; (c) Mudanjiang. The curves stand for monthly mean temperature, and the bars stand for monthly mean precipitation. * indicates that mean temperature or precipitation was significantly different (p < 0.05) between 1963–1991 and 1992–2020.
Figure 2. Monthly mean temperature and precipitation of the 3 stations in 1963–1991 (blue bar and curves) and 1992–2020 (red bar and curves). (a) Beijing; (b) Harbin; (c) Mudanjiang. The curves stand for monthly mean temperature, and the bars stand for monthly mean precipitation. * indicates that mean temperature or precipitation was significantly different (p < 0.05) between 1963–1991 and 1992–2020.
Forests 15 01712 g002
Figure 3. Frequency distribution of the mean FLD, LCD, and FFD of the 3 stations in 1963–2020. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Figure 3. Frequency distribution of the mean FLD, LCD, and FFD of the 3 stations in 1963–2020. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Forests 15 01712 g003
Figure 4. Boxplot of mean FLD, LCD, and FFD of the 3 stations in 1963–1991 and 1992–2020. (a) Beijing; (b) Harbin; (c) Mudanjiang.
Figure 4. Boxplot of mean FLD, LCD, and FFD of the 3 stations in 1963–1991 and 1992–2020. (a) Beijing; (b) Harbin; (c) Mudanjiang.
Forests 15 01712 g004
Figure 5. (af) Frequency distributions of the mean FFD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in FFD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Figure 5. (af) Frequency distributions of the mean FFD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in FFD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Forests 15 01712 g005
Figure 6. (af) Frequency distributions of the mean FLD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in FLD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Figure 6. (af) Frequency distributions of the mean FLD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in FLD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Forests 15 01712 g006
Figure 7. (af) Frequency distributions of the mean LCD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in LCD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Figure 7. (af) Frequency distributions of the mean LCD of the 3 stations in 1963–1991 (ac) and 1992–2020 (df). (gi) Frequency distributions of the deviations in LCD in the period 1992–2020 from the mean over the period 1963–1991. Categories are 3-day periods. The texture-filled bar represents species for which the deviation was significantly (p < 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; (c,f,i) Mudanjiang.
Forests 15 01712 g007
Figure 8. The optimum period (OP, solid bars) of temperature affected each phenophase most markedly. The texture-filled bar represents species for which the deviation was non-significantly (p > 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; and (c,f,i) Mudanjiang.
Figure 8. The optimum period (OP, solid bars) of temperature affected each phenophase most markedly. The texture-filled bar represents species for which the deviation was non-significantly (p > 0.05) different from zero. (a,d,g) Beijing; (b,e,h) Harbin; and (c,f,i) Mudanjiang.
Forests 15 01712 g008
Figure 9. Responses of FLD, LCD, and FFD of 3 stations to OP temperature. (left) Pearson’s correlation coefficients of FLD, LCD, and FFD against corresponding OP temperature. (right) Regression slopes of FLD, LCD, and FFD against corresponding OP temperature. Hollow geometries indicate that the correlation coefficient or regression slope was non-significant (p > 0.05). (a,d) Beijing; (b,e) Harbin; (c,f) Mudanjiang.
Figure 9. Responses of FLD, LCD, and FFD of 3 stations to OP temperature. (left) Pearson’s correlation coefficients of FLD, LCD, and FFD against corresponding OP temperature. (right) Regression slopes of FLD, LCD, and FFD against corresponding OP temperature. Hollow geometries indicate that the correlation coefficient or regression slope was non-significant (p > 0.05). (a,d) Beijing; (b,e) Harbin; (c,f) Mudanjiang.
Forests 15 01712 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cao, L.; Lin, S.; Liu, W.; Gao, C.; Bai, W.; Zhu, M.; Hao, Y.; Hao, X.; Dai, J. Shifts in Plant Phenology and Its Responses to Climate Warming in Three Temperate Cities of China during 1963–2020. Forests 2024, 15, 1712. https://doi.org/10.3390/f15101712

AMA Style

Cao L, Lin S, Liu W, Gao C, Bai W, Zhu M, Hao Y, Hao X, Dai J. Shifts in Plant Phenology and Its Responses to Climate Warming in Three Temperate Cities of China during 1963–2020. Forests. 2024; 15(10):1712. https://doi.org/10.3390/f15101712

Chicago/Turabian Style

Cao, Lijuan, Shaozhi Lin, Wei Liu, Chengxi Gao, Wenrui Bai, Mengyao Zhu, Yulong Hao, Xingming Hao, and Junhu Dai. 2024. "Shifts in Plant Phenology and Its Responses to Climate Warming in Three Temperate Cities of China during 1963–2020" Forests 15, no. 10: 1712. https://doi.org/10.3390/f15101712

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop