Important role of precipitation in controlling a more uniform spring phenology in the Qinba Mountains, China

Under global warming, the gradual pattern of spring phenology along elevation gradients (EG) has significantly changed. However, current knowledge on the phenomenon of a more uniform spring phenology is mainly focused on the effect of temperature and neglected precipitation. This study aimed to determine whether a more uniform spring phenology occurs along EG in the Qinba Mountains (QB) and explore the effect of precipitation on this pattern. We used Savitzky-Golay (S-G) filtering to extract the start of season (SOS) of the forest from the MODIS Enhanced Vegetation Index (EVI) during 2001-2018 and determined the main drivers of the SOS patterns along EG by partial correlation analyses. The SOS showed a more uniform trend along EG in the QB with a rate of 0.26 ± 0.01 days 100 m^-1 per decade during 2001-2018, but there were differences around 2011. A delayed SOS at low elevations was possibly due to the reduced spring precipitation (SP) and spring temperature (ST) between 2001 and 2011. Additionally, an advanced SOS at high elevations may have been caused by the increased SP and reduced winter temperature (WT). These divergent trends contributed to a significant uniform trend of SOS with a rate of 0.85 ± 0.02 days 100 m^-1 per decade. Since 2011, significantly higher SP (especially at low elevations) and rising ST advanced the SOS, and the SOS at lower altitudes was more advanced than at higher altitudes, resulting in greater SOS differences along EG (0.54 ± 0.02 days 100 m^-1 per decade). The SP determined the direction of the uniform trend in SOS by controlling the SOS patterns at low elevations. A more uniform SOS may have important effects on local ecosystem stability. Our findings could provide a theoretical basis for establishing ecological restoration measures in areas experiencing similar trends.

1 Introduction

Vegetation phenology is a natural phenomenon with an annual cycle that is formed by long-term adaptation to seasonal environmental changes (Liu et al., 2014; Ganjurjav et al., 2020; Cheng et al., 2021). Such phenomena often display a clear gradual pattern with increasing elevation. For example, the dates of leaf unfolding or senescence in many places show gradual postponement or advancement (Piao et al., 2011; Wu et al., 2021). This gradual variation in phenological characteristics presents a fascinating natural landscape (Gao et al., 2019). Moreover, this gradual pattern of vegetation phenology along elevation gradients (EG) plays a key role in maintaining the stability of ecosystem structure, such as carbon and nitrogen cycling, species distribution, climate feedback, and ecosystem service functions (Cong et al., 2012; Tao et al., 2018a; Shen et al., 2020; Sun et al., 2022).

However, a under climate warming, significant changes are occurring in phenological characteristics and their interactions (Shen et al., 2014; Wolf et al., 2017; Zhang et al., 2021a). For example, the reduction in chilling units due to warming could offset the increase in forcing units, and the negative impact of a higher temperatures on the start of the season (SOS) could be counterbalanced by higher precipitation (Li et al., 2020a; Wang et al., 2021). This interaction leads to a constant or delayed a SOS in some areas (Wolkovich et al., 2012; Meng et al., 2019). There is growing concern that this progressive pattern of elevation-induced phenological shifts may be changing. For instance, Chen et al. (2018) found that the spring phenology is becoming more uniform at different elevations in Europe. Temperature is generally considered the primary control of spring phenology (Piao et al., 2006; Li et al., 2016; Tao et al., 2018b). Specifically, winter warming may reduce chilling exposure at low elevations and increase spring forcing accumulation for leaf unfolding; the low temperatures at high elevations and relative increases in effective chilling accumulation may reduce the forcing requirement (Fu et al., 2015a; Asse et al., 2018; Vandvik et al., 2018). These divergent trends of leaf unfolding between high and low elevations contribute to a more uniform spring phenology.

However, the impact of temperature on the SOS is a nonlinear process, and warming (cooling) in winter will offset the advanced (delayed) SOS caused by warming (cooling) in spring to some extent (Cong et al., 2017; Piao et al., 2019; Ettinger et al., 2020; Fu et al., 2020). The SOS patterns along EG may not be fully explained temperature alone. Furthermore, recent studies suggest that precipitation may play a key role in spring phenology (Yuan et al., 2020; Sun et al., 2021; Gong et al., 2022). Precipitation somewhat determines the light and heat use efficiency of vegetation and then affects the spring phenology. More importantly, shifts in temporal trends of precipitation may directly alter the intensity of water stress on vegetation growth and change the sensitivity of vegetation to precipitation (Li et al., 2021b; Henry et al., 2022), especially in mountainous areas with high precipitation variability. Therefore, investigation of how the temporal trends in precipitation interact with SOS patterns is urgently needed. At present, relatively few studies have examined the more uniform spring phenology phenomenon. Moreover, these studies were mainly based on the assumption that temperature plays the dominant role (Chen et al., 2018; Vitasse et al., 2018; Dai et al., 2021), neglecting the effect of the temporal trends in precipitation on spring phenology. Therefore, the impact of the temporal trends in precipitation on a more uniform SOS along EG needs to be examined in depth.

As a north-south transition zone and a large-scale east-west ecological corridor in China, the Qinba Mountains (QB) has been a hot spot for ecological change research because of its high geographic complexity, biodiversity, and climate sensitivity (Xia et al., 2019; Zhang and Liang, 2020). With the warm-dry climate, the water stress on vegetation growth may be further enhanced, which may affect the progressive pattern of spring phenology along EG in the QB. To answer this question, it is necessary to deeply explore the role of precipitation intensity (water stress) in controlling the spring phenology patterns along EG in the QB. In this study, we compared the differences in the temporal trends in SOS along EG over the period 2001-2018 based on the MOD13Q1 Enhanced Vegetation Index (EVI) dataset. Partial correlation analysis was applied to identify the main controlling factors influencing the SOS patterns along EG. We aimed to answer the following questions: (1) whether there is a more uniform SOS along EG in the QB; (2) how do temperature and precipitation control the SOS patterns along EG; (3) what is the effect of temporal trends of precipitation on the SOS patterns along EG. A more uniform SOS may cause species to migrate along EG to adapt to environmental changes. This may impact species distribution and compromise the serviceability of mountain ecosystems (Inouye et al., 2000; Lenoir et al., 2008; Wolf et al., 2017). Therefore, addressing these issues is important for understanding and predicting vegetation patterns and their ecosystem functions under climate change.

2 Material and methods

2.1 Study area

The QB, which runs through Central China between 102°54′~112°40′E and 30°50′~34°59′N, is an important climatic and geographical boundary between northern and southern China (Figure 1). The entire region consists of three parts, namely, the Qinling Mountains, the Daba Mountains and the Jianghan Valley (Liu et al., 2016). The terrain clearly undulates, with an average annual precipitation of 700-1500 mm and an average annual temperature of 12-16 °C (Zhang and Liang, 2020). As the area is located in the transitional zone between the warm temperate zone and northern subtropical zone, vegetation types are clearly differentiated along EG, and their response to climate change is more sensitive (Deng et al., 2019; Qi et al., 2021).

FIGURE 1

Figure 1 Location and elevation map of the Qinba Mountains.

2.2 Extraction of spring phenology

The SOS in the QB was extracted from the MOD13Q1 EVI dataset provided by the official NASA website (https://ladsweb.modaps.eosdis.nasa.gov/search/) for the period 2001-2018 with a spatial resolution of 250 m and a temporal resolution of 16 days. Due to the influence of the sensors themselves and other external factors, there were missing data or outliers, resulting in differences between the EVI time series curves and the real vegetation growth patterns (Qi et al., 2021). Therefore, we used Savitzky-Golay (S-G) filtering based on Timesat 3.3 software to eliminate noise that deviated from the normal growth trend line. By comparing the results of multiple calculations, we set the window size to 5, the envelope iterations to 2, and the adaptation strength to 8 to reconstruct the EVI time series of the QB from 2001 to 2018. Then, we defined the SOS as the date when the EVI increased to 20% of the seasonal amplitude (Wu et al., 2021). To reduce the uncertainty from the outliers, we excluded pixels with SOS earlier than the 30th day of the year or later the 180th day of the year.

2.3 Climate data

The Chinese meteorological forcing dataset (CMFD) is widely used for its continuous time coverage and consistent quality (http://data.tpdc.ac.cn/zh-hans/), with an accuracy between meteorological observation data and satellite remote sensing data, and a spatial resolution of 0.1° × 0.1° (Yang et al., 2010; He et al., 2020). Therefore, we selected temperature and precipitation datasets from the CMFD during 2000-2018 and analysed the temporal trends in winter temperature (WT, from the previous November to the beginning of greening in March), spring temperature (ST, from the beginning of greening in March to the end of greening in June), and spring precipitation (SP, from the beginning of greening in March to the end of greening in June) along EG in the QB at 100 m altitude intervals.

Considering the low spatial resolution of the CMFD dataset, the number of pixels in the QB is relatively small, and there may be errors based on pixel statistical analysis. To this end, we used a geographically weighted regression (GWR) model to downscale the WT, ST, and SP data to a spatial resolution of 0.01° (Lu et al., 2018; Lu et al., 2019). We selected the EVI and elevation as influencing factors and downscaled the WT, ST, and SP data by three levels with an intermediate resolution of 0.05° (Figure S1) (Zhang et al., 2020; Zhang and Cheng, 2020; Xu and Cheng, 2021). The values were calculated as follows:

$\begin{array}{ll}{y}_i={\beta }_0\left(u_i,v_i\right)+{\beta }_1\left(u_i,v_i\right)x_{i1}+{\beta }_2\left(u_i,v_i\right)x_{i2}+{ϵ}_i& (1)\end{array}$

Where (u_i,v_i) . denotes the spatial coordinates at a spatial position i, x_i1 and x_i2 . present the values of elevation and EVI, respectively. β₀, β₁, and β₂ . present the constant terms, the regression coefficients of elevation and the regression coefficient of EVI at raster i, espectively.

2.4 DEM and land use data

The digital elevation model (DEM) data were obtained from the 90 m resolution SRTM product provided by the Geospatial Data Cloud (https://www.gscloud.cn/). Land Use Data, derived from Globe Land Cover in 2010 (http://www.globallandcover.com/), were used to extract the distribution of the forest vegetation in the QB. The DEM and Land Use Data were resampled to 250 m to match the resolution of the EVI by using ArcMap 10.3.

2.5 Statistical analysis

The slope of SOS along EG for each year from 2001 to 2018 was used to compare the difference in SOS along EG among different years (Table S1). The smaller the slope was, the smaller the difference in SOS along EG. The quadratic curve fitting method based on Origin 2018 was used to detect abrupt changes in the slope of SOS along EG. The year corresponding to the point where the slope line of the curve intersected with the 0 slope line was considered an abrupt change point. Then, Theil-Sen trend analysis was used to analyse the temporal trends of SOS, WT, ST, and SP in the QB from 2001 to 2018 and before and after the abrupt years, and the significance of trends was tested by the Mann-Kendall (M-K) statistical test by using MATLAB R2016a (Sen, 2012; Liu et al., 2016). This value was calculated as follows:

$\begin{array}{ll}\beta =mean\left(\frac{x_j-x_i}{j-i}\right)\forall j>i& (2)\end{array}$

where β is the temporal trend of the SOS; j and i denote the time series; and x_j and x_i denote the SOS at times j and i, respectively. β > 0 indicates that the temporal trend of the SOS has a delayed trend, and β< 0 indicates an advanced trend.

Partial correlation analysis based on MATLAB R2016a was used to analyse the relationship between SOS and WT, ST, and SP. The degree of association between the two variables was measured by the partial correlation coefficient after excluding the effects of other control variables (Shen et al., 2020). The climate factors (WT, ST, and SP) were used as independent variables, and the SOS was the dependent variable. Statistical significance was determined at a level of P< 0.05 based on a two-tailed t-test. Additionally, the temporal trend and partial correlation coefficient patterns in SOS and the climate factors patterns along EG were analysed at 100 m altitude intervals. Considering the rarity of forest pixels at high elevations, regions with fewer than 100 pixels along EG were excluded, and only the regions below 3100 m were analysed in the QB.

3 Results

3.1 Spatial distribution of SOS from 2001 to 2018

The mean SOS from 2001 to 2018 was early in the valleys and late on the mountains. Among them, 89.71% of the mean SOS were concentrated around 73-105 days. Along EG, the mean SOS was significantly (P< 0.001) delayed by 1.7 days per 100 m increase (Figure 2A). The temporal trends of delayed SOS (48.96% of total pixels) were mainly distributed in the valley and at marginal low elevations, and the magnitude of delay was concentrated around 0-0.8 days·a^-1 (Figure 2B). Significant delayed (P< 0.05) were mainly distributed in the northeastern low elevation regions of the QB (Figure S2A). The advanced SOS (51.04%) was mainly distributed at high elevations, and the magnitude of advance was concentrated in 0-0.8 days·a^-1 (Figure 2B). Significant advanced (P< 0.05) were mainly distributed in the western, northern, and southern alpine regions (Figure S2A). Along EG, the temporal trends in delayed SOS gradually decreased with increasing elevation (below 1000 m), and the advanced SOS gradually increased (above 1000 m). These divergent trends in SOS between high and low elevations showed a more uniform trend of SOS along EG in the QB.

FIGURE 2

Figure 2 Spatial distribution of mean the start of season (SOS) (A) and temporal trends in SOS (B) during 2001-2018. The shaded area represents the standard deviation.

3.2 Temporal trends in the SOS before and after 2011

The change in slope of SOS along EG from 2001 to 2018 decreased by 0.26 ± 0.01 days 100 m^-1 per decade (P = 0.06), the differences in SOS along EG decreased continuously, and the SOS showed a more uniform trend along EG (Figure 3B). Furthermore, 2011 was identified as a year of abrupt change according to the change in the SOS slope along EG from 2001 to 2018 based on the quadratic curve fitting method (Figure 3A). The change in the SOS slope along EG significantly (P = 0.001) decreased by 0.85 ± 0.02 days 100 m^-1 per decade between 2001 and 2011. Since 2011, the change was reversed by 0.54 ± 0.02 days 100 m^-1 per decade.

FIGURE 3

Figure 3 Slope change and curve slope of SOS along EG from 2001 to 2018. (A) The binomial curve of the slope of SOS along EG. The red line represents the binomial curve, the blue line represents the slope of the binomial curve. (B) The change in slope of SOS along EG from 2001 to 2018. The black linerepresents the change during 2001-2018, the red line represents 2001-2011, and the blue line represents 2011-2018.

The temporal trends in SOS before 2011 mainly showed a delayed trend (67.54%) that was distributed in the regions below 1500 m, and the magnitude decreased with increasing elevation (Figures 4A–D). Significant delayed (P< 0.05) were mainly distributed at low elevations in the eastern and southern regions (Figure S2B). While the advanced SOS (32.46%) was mainly distributed at high elevations above 1500 m, and the magnitude was increased (Figures 4A–D). Significant advanced (P< 0.05) were mainly distributed in the western, northern, and southern alpine regions (Figure S2B). These divergent trends showed a more uniform trend of SOS along EG in the QB.

FIGURE 4

Figure 4 Spatial distribution of temporal trends in SOS during 2001-2011 (A) and 2011-2018 (B). (C) The proportion of temporal trends in SOS. The redcolumns represent the proportion of SOS during 2001-2011,the blue columns represent the proportion of SOS during2011-2018. (D) The slope of trends in SOS along EG. The redline represent 2001-2011, the blue line represents 2011-2018.

However, the temporal trends in SOS after 2011 mainly showed an advanced trend (73.38%) at all elevations, and the magnitude decreased by 0.06 days·a ^-1 for each 100 m increase (Figures 4B–D). Significant advanced (P< 0.05) were mainly distributed in the eastern and southern regions (Figure S2C). The advanced SOS (53.97%), which changed from a delayed SOS before 2011, was mainly distributed at low elevations (Figure S3), and the magnitude was stronger than that at high elevations. The difference in SOS gradually widened along EG. Shifts in the temporal trend of SOS at low elevations determined the direction of the uniform trend in SOS along EG.

3.3 Temporal trends in WT, ST and SP around 2011

The results of the trend analysis showed that the temporal trends in WT before 2011 mainly showed a cooling trend (79.82%), which was mainly distributed in marginal low elevation regions (Figure 5A, C). Significant decreased (P< 0.05) were mainly distributed in the eastern high elevation regions (Figure S4A). Along EG, the temporal trends in WT decreased by 0.02°C·(10a) ^-1 for each 100 m increase (Figure 5D). After 2011, the temporal trends in WT mainly showed a warming trend (87.47%), which was mainly distributed in the marginal low elevation areas (Figures 5B, C). Significant increased (P< 0.05) were mainly distributed in the eastern and southern high elevation regions (Figure S4D). Along EG, the temporal trends in WT increased by 0.05°C·(10a) ^-1 for each 100 m increase (Figure 5D).

FIGURE 5

Figure 5 Temporal trends in winter temperature (WT) before and after 2011. (A) The temporal trends in WT during 2001-2011. (B) The temporal trends in WT during 2011-2018. (C) The proportion of temporal trends in WT. The red columns represent the proportion of WT during 2001-2011, the blue columns represent the proportion of WT during 2011-2018. (D) The slope of trends in WT along EG. The red line represents 2001-2011, the blue line represents 2011-2018.

The temporal trends in ST before 2011 mainly showed a cooling trend (60.26%), which was mainly distributed in high elevation regions (Figures 6A, C). Significant decreased (P< 0.05) were mainly distributed in the northern and eastern high elevation regions (Figure S4B). While the warming ST (39.74%) was concentrated in the marginal low elevation regions and the central Daba Mountains (Figure 6A). Along EG, the temporal trends in ST decreased by 0.02°C·(10a) ^-1 for each 100 m increase, and the temporal trends in ST showed a slight increase below 800 m and then cooled (Figure 6D). After 2011, the temporal trends in ST mainly showed a warming trend (73.85%), which was mainly distributed in the central and western regions. While the cooling ST (26.15%) was concentrated in the eastern low elevation regions (Figures 6B, C), with a significant decreased (P< 0.05) in the northeastern regions (Figure S4E). Along EG, the temporal trends in ST increased by 0.09°C·(10a) ^-1 for each 100 m increase, and the temporal trends in ST showed a decrease below 500 m and then warmed (Figure 6D).

FIGURE 6

Figure 6 Temporal trends in spring temperature (ST) before and after 2011. (A) The temporal trends in ST during 2001-2011. (B) The temporal trends in ST during 2011-2018. (C) The proportion of temporal trends in ST. The red columns represent the proportion of ST during 2001-2011, the blue columns represent the proportion of ST during 2011-2018. (D) The slope of trends in ST along EG. The red line represents 2001-2011, the blue line represents 2011-2018.

The temporal trends in SP before 2011 increased on average by 0.45 mm·a^-1. Among them, the increased regions (59.4%) were mainly distributed at high elevations, while the decreased regions (40.6%) were distributed at marginally low elevations. The temporal trends in SP were of low amplitudes and mainly concentrated around 0-1 mm·a^-1 (Figures 7A, C). Along EG, the temporal trends in SP increased by 0.18 mm·(10a)^-1 for each 100 m increase (Figure 7D). After 2011, the temporal trends in SP increased by 4.41 mm·a^-1 on average, which was approximately 10 times higher than that before 2011. Among them, the increased regions (96.45%) were mainly concentrated around 2-8 mm·a^-1 (Figures 7B, C), with a significant increased (P< 0.05) in the eastern low elevation regions (Figure S4F). Along EG, the temporal trends in SP decreased by 0.29 mm·(10a)^-1 for each 100 m increase (Figure 7D).

FIGURE 7

Figure 7 Temporal trends in spring precipitation (SP) before and after 2011. (A) The temporal trends in SP during 2001-2011. (B) The temporal trends in SP during 2011-2018. (C) The proportion of temporal trends in SP. The red columns represent the proportion of SP during 2001-2011, the blue columns represent the proportion of SP during 2011-2018. (D) The slope of trends in SP along EG. The red line represents 2001-2011, the blue line represents 2011-2018.

3.4 Relationship between SOS and its potential drivers

Based on the results from the partial correlation analysis, we found that the SOS and WT were mainly negatively partially correlated (76.18%) in the QB before 2011, with a significant negative partial correlation (P< 0.05) in the central and eastern low elevation regions (Figure S5A), while the regions with a positive partial correlation (20.18%) were concentrated in the western and southern high elevation areas (Figure 8A). Along EG, the negative partial correlation coefficient gradually weakened with increasing elevation (below 2100 m), and the positive correlation coefficient continued to increase (above 2100 m) (Figure 8D). The SOS before 2011 showed a partial negative correlation with ST (75.36%), with a significant negative partial correlation (P< 0.05) in the central and eastern low elevation regions (Figure S5B), while the regions with a positive partial correlation coefficient (24.64%) were concentrated in the western and southern high elevation areas (Figure 8B). Along EG, the negative partial correlation coefficient significantly weakened (P< 0.001) with increasing elevation (Figure 8E). The SOS showed a partial negative correlation with SP (57.91%), with a significant negative partial correlation (P< 0.05) in the eastern low elevation regions and the western and southern high elevation regions (Figure S5C), while the regions with a positive partial correlation coefficient (42.09%) were distributed within all elevation gradients (Figure 8C). Along EG, the negative partial correlation coefficient did not change much (P = 0.28) with increasing elevation (Figure 8F).

FIGURE 8

Figure 8 Partial correlation coefficient between SOS and its potential drivers during 2001-2011. (A–C) are correlations of SOS with WT, ST, and SP, respectively. (D–F) are the distributions of coefficients along EG corresponding to (A–C), respectively. The shaded area indicates the 95% confidence interval.

For 2011-2018, the SOS and WT were mainly positively partially correlated (63.81%), with a significant positive partial correlation (P< 0.05) in the eastern low elevation regions and the central and western regions (Figure S5D), while the regions with a positive partial correlation (36.19%) were concentrated in the northeastern and southern high elevation areas (Figure 9A). Along EG, the positive partial correlation coefficient significantly increased (P< 0.001) with increasing elevation (Figure 9D). The SOS showed a partial negative correlation with ST (74%), with a significant negative partial correlation (P< 0.05) in the eastern low elevation regions and the central and western regions (Figure S5E), while the regions with a positive partial correlation coefficient (26%) were concentrated in the northeastern and southern high elevation areas (Figure 9B). Along EG, the negative partial correlation coefficient significantly increased (P< 0.001) with increasing elevation (Figure 9E). The SOS showed a partial negative correlation with SP (67.52%), with a significant negative partial correlation (P< 0.05) in the eastern low elevation regions and the western and southern high elevation regions (Figure S5F), while the regions with a positive partial correlation coefficient (32.48%) were sporadically distributed within all elevation gradients (Figure 9C). Along EG, the negative partial correlation coefficient significantly weakened (P = 0.03) with increasing elevation (Figure 9F).

FIGURE 9

Figure 9 Partial correlation coefficient between SOS and its potential drivers during 2011-2018. (A–C) are correlations of SOS with WT, ST, and SP, respectively. (D–F) are the distributions of coefficients along EG corresponding to (A–C), respectively. The shaded area indicates the 95% confidence interval.

4 Discussion

4.1 SOS response to driving factors

In this study, we found a decreased WT and a slightly increased ST below 500 m in the QB from 2001 to 2011, which somewhat provides more sufficient accumulation of chilling and forcing for the onset of spring phenology (Fu et al., 2015a; Ettinger et al., 2020; Pan et al., 2022). However, the SOS in this region was delayed (Figure 4D). Temperature alone did not sufficiently explain the SOS patterns in this region. Moreover, the reduced SP in this region (Figure 7D) somewhat limited the utility of water and heat conditions and offset the advanced SOS that would have been caused by the decreased WT and increased ST (Jewaria et al., 2021; Zheng et al., 2021), resulting in a delay in SOS. Hence, the decreased SP was the predominant controlling factor for the delayed SOS in this region. Similarly, at elevations of 500-1500 m, the decreased WT and increased SP somewhat facilitated the accumulation of chilling and the water demand for vegetation growth (Shen et al., 2015; Ettinger et al., 2020), which would have advanced the SOS. However, this advance would have likely been counteracted by reduced forcing as a result of the reduced ST (Figure 6D), which better explained why the SOS in this region was delayed (Piao et al., 2019; Ganjurjav et al., 2020). The decreased ST may have a stronger impact on SOS relative to the decreased WT and the increased SP. Meanwhile, at elevations above 1500 m, the decreased ST reduced the accumulation of forcing and delayed the SOS. Moreover, the decreased WT and the increased SP in this region (Figures 5D, 7D) somewhat increased the accumulation of chilling and satisfied the water demand for the advanced SOS (Lin et al., 2022; Wang et al., 2022). The SOS response to reduced ST may have been less than that to the reduced WT and increased SP, resulting in the advanced SOS in this region. This showed that the joint control of the reduced WT and increased SP, rather than the ST, was the driving factor for the advanced SOS.

For 2011-2018, the increased WT and the reduced ST at elevations below 500 m were not conducive for more chilling and forcing accumulation (Yu et al., 2010; Wang et al., 2020a), which somewhat delayed the SOS. However, this delay was counteracted by the increased water demand as a result of the increased SP (Figure 7D) and led to the advance of SOS (Shen et al., 2015). The advanced SOS below 500 m was influenced by the increased SP rather than temperature. Similarly, at elevations above 500 m, the increased WT may reduce chilling accumulation and delay the SOS (Wang et al., 2020b; Jewaria et al., 2021). However, this delay was offset by the increased water demand and forcing accumulation (Figures 6D, 7D), resulting in the advanced SOS (Piao et al., 2019; Li et al., 2021a). The effects of increased ST and SP on SOS were stronger than those of increased WT. The present results suggest that the SP also plays a crucial role in controlling the patterns of SOS along EG, and analyses conducted from the perspective of temperature alone may not fully elucidate the intrinsic mechanism of these patterns in SOS.

4.2 Effect of precipitation in controlling a more uniform SOS

Our research demonstrated that the magnitude of SOS response to SP gradually decreased and the response to temperature gradually increased with increasing elevation, which was consistent with previous studies in the QB (Liu et al., 2016; Chen et al., 2019). The SOS at high elevations, where the SP showed a consistent increase from 2001 to 2018 (Figure 7D), was less sensitive to SP. In other words, sufficient SP at high elevations resulted in less water stress on vegetation growth (Jennifer and Florence, 2018; Gupta et al., 2020), and larger amounts of SP would not advance SOS. Along with this was the potential for greater temperature sensitivity of SOS to maximise thermal benefits (Shen et al., 2015; Fu et al., 2021). This hypothesis was further supported by the weaker advance of SOS in the context of stronger precipitation after 2011 and the stronger partial correlation between SOS and temperature (Figures 4D, 9). In contrast, at low elevations, where the SP showed a decrease from 2001 to 2011 (Figure 7D), maximised water usage led to the stronger sensitivity of SOS to SP (Shen et al., 2015; Chen et al., 2019). Stronger water stress led to a delayed SOS even under better temperature conditions (Figures 5D, 6D). In addition, we speculated that the SOS would advance if more rainfall occurred after stronger water stress, even if the temperatures were less optimal. This hypothesis was confirmed by the significant advance of SOS at low elevations under the significant increase in SP after 2011 (Figure 4D, Figure 7D). This result showed that the temporal trends in SP played a crucial role in controlling the patterns of SOS at low elevations, which determined the direction of the uniform trend in SOS along EG.

Recent research has predicted that by the end of this century, the climate in east-central China will have a continuous warm-dry trend (Ma et al., 2019; Zhang et al., 2022), which may enhance water evaporation at low elevations in the QB. In addition, the structural overshoot due to the warm winter may further increase the water stress during the growing season (Bastos et al., 2020; Zhang et al., 2021b), thus strengthening the uniform trend of SOS along EG in the QB. Therefore, the sensitivity of SOS to precipitation may be further enhanced, and future changes in the spatiotemporal distribution of precipitation, rather than temperature, are likely to have a stronger control on the direction of the uniform trend in SOS. Since a more uniform SOS along EG may compromise the stability and serviceability of mountain ecosystems (Inouye et al., 2000; Lenoir et al., 2008; Wolf et al., 2017), it is crucial to improve the water-use efficiency of vegetation at low elevations. Previous studies have shown that a higher species diversity could notably enhance drought resistance (Grossiord, 2020; Liu et al., 2022), and different tree species have different strategies and abilities to cope with water stress (Anderegg et al., 2018; Li et al., 2020b). In the new round of ecological restoration projects, conversion of the current monoculture to mixed-species tree plantations and the planting of resilient tree species with a high water-use efficiency could relieve potential water stress in the future.

4.3 Comparison with other studies

In this study, we confirmed that the SOS extracted by the MOD13Q1 EVI can be used to accurately trace the spring phenology in the QB. This conclusion was predominantly based on the strong consistency with previous studies in related areas that were based on different data sources (Table S2). In addition, our findings demonstrated that the SOS showed a more uniform trend along EG between 2001 and 2018 in the QB. These divergent trends of SOS between low and high elevations led to a uniform trend of SOS along EG in the QB, which was consistent with the patterns of SOS in the Alps (Vitasse et al., 2018). However, the driving mechanisms of the more uniform SOS in the QB were different from those reported in the Alps. Vitasse et al. (2018), based on the hypothesis that temperature plays the dominant role, found that the reduced chilling accumulation at low elevations caused by the warming WT moderated the magnitude of SOS advance compared to high elevations. The SOS was more advanced at higher altitudes than at lower altitudes, resulting in a more uniform trend of SOS along EG. In this study, we added the effects of SP and found that the SP, rather than temperature, determined the direction of the uniform trend in SOS by controlling the SOS patterns at low elevations. These differences may have been due to potential mechanisms playing diverse roles across different areas.

In addition, the results demonstrated negative relationships between ST, SP and SOS during 2001-2018 and a positive relationship between WT and SOS between 2011 and 2018, which was consistent with most previous phenological studies (Piao et al., 2006; Zheng et al., 2021). However, the WT was negatively correlated with SOS below 2100 m during 2001-2011 and the cooling WT delayed the SOS in the QB. The simultaneous reduction of ST between 500 and 2100 m and the SP below 500 m may have counteracted the advance of SOS caused by decreased WT (Cong et al., 2017; Asse et al., 2018), possibly explaining why the SOS was delayed under a cooling WT.

In contrast to our results, Gao et al. (2019) found that there were no prevalent trends of elevational homogenization of SOS in most regions worldwide over the last 30 years. This discrepancy was likely due to the difference in the main drivers of vegetation growth at different spatial scales. For example, the stronger water stress in arid areas results in a greater sensitivity of SOS to precipitation, but the sensitivity of SOS to temperature and precipitation may change at a larger scale (Cleverly et al., 2016; Cong et al., 2021). In addition, the magnitudes and even the directions of the response to climate change are largely different among vegetation types (Lesica and Kittelson, 2010; Fu et al., 2015b). The proportion of vegetation types at different spatial scales varies, which may lead to differences in the variation patterns of SOS.

4.4 Uncertainty analysis

First, using different abrupt change tests may lead to different results. We used the M-K test, sliding t-test, and quadratic curve fitting test to detect abrupt changes in the slope of SOS along EG from 2001 to 2018 (Figure S6). We found that the abrupt change years were approximately 2003, 2008, and 2011 in the M-K test, sliding t test, and quadratic curve fitting test, respectively. Moreover, previous research has shown that the vegetation coverage reversed in 2010 in the QB during 2001-2014 (Liu et al., 2016). Chen et al. (2021) found that 2010 was the significant acceleration point for gross and net primary production variations in China during 2001-2018. Therefore, using the quadratic curve fitting test, we defined 2011 as an abrupt change year. Second, since the CMFD dataset was reanalysis data, the temperature and precipitation data already contained elevation information (He et al., 2020). Downscaling by selecting the EVI and elevation as influencing factors may overestimate the effect of elevation on temperature and precipitation, which somewhat impacts the downscaling results. Finally, different vegetation types have different responses to climate change. As the dividing line between north subtropical and temperate forests, vegetation types in the QB are diverse and have an obvious altitudinal differentiation (Deng et al., 2019; Qi et al., 2021). Therefore, the patterns of SOS along EG on vegetation types and the north and south slopes need to be further studied.

5 Conclusions

The variation patterns of SOS along EG in the QB from 2001 to 2018 were explored. We found that the SOS showed a more uniform trend along EG in the QB. Furthermore, this uniform trend of SOS along EG was not completely continuous and reversed around 2011. The ST and SP were, in general, negatively partially correlated with SOS, and WT was positively partially correlated with SOS (except for regions below 2100 m during 2001-2011). Before 2011, the decreased ST and SP led to a delayed SOS at low elevation and the increased SP and decreased WT led to an advanced SOS at high elevation. These opposite SOS trends at high and low elevations led to a more uniform SOS along EG. After 2011, increased SP and ST led to a much greater advance of SOS at low elevations than at high elevations, resulting in a gradual widening of the difference in SOS along EG. Moreover, the temporal trends in SP played a crucial role in controlling the SOS patterns at low elevations, which determined the direction of the uniform trend in SOS along EG. Our study deepens the understanding of the altitudinal sensitivity of SOS under climate change and provides a theoretical basis for regions that experience a uniform trend of SOS along EG to develop appropriate ecological measures to mitigate its adverse effects.

Data availability statement

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

Author contributions

JL and JZ conceived and designed this study. BL, DY, and JG provided guidance on research methods and content, WH and RT assisted with the data analysis, JL wrote the first draft of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (41461035), and major projects to protect and restore important ecosystems in Xinjiang (202005140014).

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.

Supplementary material

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

References

Anderegg, W., Konings, A., Trugman, A., Yu, K., Bowling, D., Gabbitas, R., et al. (2018). Hydraulic diversity of forests regulates ecosystem resilience during drought. Nature 561, 538–541. doi: 10.1038/s41586-018-0539-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Asse, D., Chuine, I., Vitasse, Y., Yoccoz, N. G., Delpierre, N., Badeau, V., et al. (2018). Warmer winters reduce the advance of tree spring phenology induced by warmer springs in the Alps. Agr For. Meteorol 252, 220–230. doi: 10.1016/j.agrformet.2018.01.030

CrossRef Full Text | Google Scholar

Bastos, A., Ciais, P., Friedlingstein, P., Sitch, S., Pongratz, J., Fan, L. (2020). Direct and seasonal legacy effects of the 2018 heat wave and drought on European ecosystem productivity. Sci. Adv. 6, 2724. doi: 10.1126/sciadv.aba2724

CrossRef Full Text | Google Scholar

Chen, Y., Feng, X., Tian, H., Wu, X., Gao, Z., Feng, Y., et al. (2021). Accelerated increase in vegetation carbon sequestration in China after 2010: A turning point resulting from climate and human interaction. Global Change Biol. 27, 5848–5864. doi: 10.1111/GCB.15854

CrossRef Full Text | Google Scholar

Cheng, W., Li, Z., 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

Chen, L., Huang, J., Ma, Q., Hänninen, H., Rossi, S., Piao, S., et al. (2018). Spring phenology at different altitudes is becoming more uniform under global warming in Europe. Global Change Biol. 24, 3969–3975. doi: 10.1111/gcb.14288

CrossRef Full Text | Google Scholar

Chen, C., Zhu, L., Tian, L., Li, X. (2019). Spatial-temporal changes in vegetation characteristics and climate in the qinling-daba mountains. Acta Ecol. Sin. 39, 3257–3266. doi: 10.5846/stxb201801300252

CrossRef Full Text | Google Scholar

Cleverly, J., Eamus, D., Restrepo, C. N., Chen, C., Maes, W., Li, L., et al. (2016). Soil moisture controls on phenology and productivity in a semi-arid critical zone. Sci. Total Environ. 568, 1227–1237. doi: 10.1016/j.scitotenv.2016.05.142

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, N., Huang, K., Zhang, Y. (2021). Unsynchronized driving mechanisms of spring and autumn phenology over northern hemisphere grasslands. Front. For Glob Change 3. doi: 10.3389/ffgc.2020.610162

CrossRef Full Text | Google Scholar

Cong, N., Piao, S., Chen, A., Wang, X., Lin, X., Chen, S., et al. (2012). Spring vegetation green-up date in China inferred from SPOT NDVI data: A multiple model analysis. Agr For. Meteorol 165, 104–113. doi: 10.1016/j.agrformet.2012.06.009

CrossRef Full Text | Google Scholar

Cong, N., Shen, M., Piao, S., Chen, X., An, S., Yang, W., et al. (2017). Little change in heat requirement for vegetation green-up on the Tibetan plateau over the warming period of 1998–2012. Agr For. Meteorol 232, 650–658. doi: 10.1016/j.agrformet.2016.10.021

CrossRef Full Text | Google Scholar

Dai, J., Zhu, M., Mao, W., Liu, R., Wang, H., Alatalo, J. M., et al. (2021). Divergent changes of the elevational synchronicity in vegetation spring phenology in north China from 2001 to 2017 in connection with variations in chilling. Int. J. Climatol 41, 6109–6121. doi: 10.1002/JOC.7170

CrossRef Full Text | Google Scholar

Deng, C., Bai, H., Gao, S., Zhao, T., Ma, X. (2019). Differences and variations in the elevation-dependent climatic growing season of the northern and southern slopes of the qinling mountains of China from 1985 to 2015. Theor. Appl. Climatol 137, 1159–1169. doi: 10.1007/s00704-018-2654-7

CrossRef Full Text | Google Scholar

Ettinger, A. K., Chamberlain, C. J., Morales-Castilla, I., Buonaiuto, D. M., Flynn, D. F. B., Savas, T., et al. (2020). Winter temperatures predominate in spring phenological responses to warming. Nat. Clim Change 10, 1137–1142. doi: 10.1038/s41558-020-00917-3

CrossRef Full Text | Google Scholar

Fu, Y., Li, X., Zhou, X., Geng, X., Guo, Y., Zhang, Y. (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., Piao, S., Vitasse, Y., Zhao, H., De Boeck, H. J., Liu, Q., et al. (2015b). Increased heat requirement for leaf flushing in temperate woody species over 1980-2012: Effects of chilling, precipitation and insolation. Glob Change Biol. 21, 2687–2697. doi: 10.1111/gcb.12863

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y., Zhou, X., Li, X., Zhang, Y., Geng, X., Hao, F., et al. (2021). Decreasing control of precipitation on grassland spring phenology in temperate China. Global Ecol. Biogeogr 30, 490–499. doi: 10.1111/geb.13234

CrossRef Full Text | Google Scholar

Ganjurjav, H., Gornish, E. S., Hu, G., Schwartz, M. W., Wan, Y., Li, Y., et al. (2020). Warming and precipitation addition interact to affect plant spring phenology in alpine meadows on the central qinghai-Tibetan plateau. Agr For. Meteorol 287, 107943. doi: 10.1016/j.agrformet.2020.107943

CrossRef Full Text | Google Scholar

Gao, M., Piao, S., Chen, A., Yang, H., Liu, Q., Fu, Y. H., et al. (2019). Divergent changes in the elevational gradient of vegetation activities over the last 30 years. Nat. Commun. 10. doi: 10.1038/s41467-019-11035-w

CrossRef Full Text | Google Scholar

Gong, F., Chen, X., Yuan, W., Su, Y., Yang, X., Liu, L., et al. (2022). Partitioning of three phenology rhythms in American tropical and subtropical forests using remotely sensed solar-induced chlorophyll fluorescence and field litterfall observations. Int. J. Appl. Earth Obs. Geoinf 107, 102698. doi: 10.1016/J.JAG.2022.102698

CrossRef Full Text | Google Scholar

Grossiord, C. (2020). Having the right neighbors: how tree species diversity modulates drought impacts on forests. New Phytol. 228, 42–49. doi: 10.1111/nph.15667

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, A., Rico-Medina, A., Caño-Delgado, A. I. (2020). The physiology of plant responses to drought. Science 368, 266–269. doi: 10.1126/science.aaz7614

PubMed Abstract | CrossRef Full Text | Google Scholar

Henry, E. H., Terando, A. J., Morris, W. F., Daniels, J. C., Haddad, N. M. (2022). Shifting precipitation regimes alter the phenology and population dynamics of low latitude ectotherms. Clim Change Ecol. 3, 2666–9005. doi: 10.1016/j.ecochg.2022.100051

CrossRef Full Text | Google Scholar

He, J., Yang, K., Tang, W., Lu, H., Qin, J., Chen, Y., et al. (2020). The first high-resolution meteorological forcing dataset for land process studies over China. Sci. Data 7. doi: 10.1038/s41597-020-0369-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Inouye, D. W., Barr, B., Armitage, K. B., Inouye, B. D. (2000). Climate change is affecting altitudinal migrants and hibernating species. Proc. Natl. Acad. Sci. U.S.A. 97, 1630–1633. doi: 10.1073/pnas.97.4.1630

PubMed Abstract | CrossRef Full Text | Google Scholar

Jennifer, A. B., Florence, A. V. (2018). Implications of plant functional traits and drought survival strategies for ecological restoration. J. Appl. Ecol. 55, 631–640. doi: 10.1111/1365-2664.12979

CrossRef Full Text | Google Scholar

Jewaria, P. K., Hänninen, H., Li, X., Bhalerao, R. P., Zhang, R. (2021). A hundred years after: endodormancy and the chilling requirement in subtropical trees. New Phytol. 231, 565–570. doi: 10.1111/nph.17382

PubMed Abstract | CrossRef Full Text | Google Scholar

Lenoir, J., Geígout, J. C., Marquet, P. A., de Ruffray, P., Brisse, H. (2008). A significant upward shift in plant species optimum elevation during the 20th century. Science 320, 1768–1771. doi: 10.1126/science.1156831

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesica, P., Kittelson, P. M. (2010). Precipitation and temperature are associated with advanced flowering phenology in a semi-arid grassland. J. Arid Environ. 74, 1013–1017. doi: 10.1016/j.jaridenv.2010.02.002

CrossRef Full Text | Google Scholar

Li, X., Fu, Y. H., Chen, S., Xiao, J., Yin, G., Li, X., et al. (2021b). Increasing importance of precipitation in spring phenology with decreasing latitudes in subtropical forest area in China. Agr For. Meteorol, 304–305. doi: 10.1016/j.agrformet.2021.108427

CrossRef Full Text | Google Scholar

Li, J., Huang, J., Tardif, J. C., Liang, H., Jiang, S., Zhu, H., et al. (2020a). Spatially heterogeneous responses of tree radial growth to recent El niño southern-oscillation variability across East Asia subtropical forests. Agr For. Meteorol 287, 107939. doi: 10.1016/j.agrformet.2020.107939

CrossRef Full Text | Google Scholar

Li, H., Jiang, J., Chen, B., Li, Y., Xu, Y., Shen, W. (2016). Pattern of NDVI-based vegetation greening along an altitudinal gradient in the eastern Himalayas and its response to global warming. Environ. Monit Assess. 188. doi: 10.1007/s10661-016-5196-4

CrossRef Full Text | Google Scholar

Lin, S., Wang, H., Ge, Q., Hu, Z. (2022). Effects of chilling on heat requirement of spring phenology vary between years. Agr For. Meteorol 312, 108718. doi: 10.1016/j.agrformet.2021.108718

CrossRef Full Text | Google Scholar

Li, X., Piao, S., Wang, K., Wang, X., Wang, T., Ciais, P., et al. (2020b). Temporal trade-off between gymnosperm resistance and resilience increases forest sensitivity to extreme drought. Nat. Ecol. Evol. 4, 1075–1083. doi: 10.1038/s41559-020-1217-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Tao, J., Cheng, B., Wu, Q., Peng, H. (2021a). Sensitivity of spring phenology to elevation in qinling mountains, China. Chin. J. Appl. Ecol. 32, 2089–2097. doi: 10.13287/j.1001-9332.202106.016

CrossRef Full Text | Google Scholar

Liu, L., Liu, L., Liang, L., Donnelly, A., Park, I., Schwartz, M. D. (2014). Effects of elevation on spring phenological sensitivity to temperature in Tibetan plateau grasslands. Chin. Sci. Bull. 59, 4856–4863. doi: 10.1007/s11434-014-0476-2

CrossRef Full Text | Google Scholar

Liu, D., Wang, T., Peñuelas, J., Piao, S. (2022). Drought resistance enhanced by tree species diversity in global forests. Nat. Geosci 15, 800–804. doi: 10.1038/s41561-022-01026-w

CrossRef Full Text | Google Scholar

Liu, X., Zhu, X., Pan, Y., Li, S., Ma, Y., Nie, J. (2016). Vegetation dynamics in qinling-daba mountains in relation to climate factors between 2000 and 2014. J. Geogr. Sci. 26, 45–58. doi: 10.1007/s11442-016-1253-8

CrossRef Full Text | Google Scholar

Lu, B., Brunsdon, C., Charlton, M., Harris, P. (2019). A response to 'A comment on geographically weighted regression with parameter-specific distance metrics'. Int. J. Geogr. Inf Sci. 33, 1300–1312. doi: 10.1080/13658816.2019.1585541

CrossRef Full Text | Google Scholar

Lu, B., Yang, W., Ge, Y., Harris, P. (2018). Improvements to the calibration of a geographically weighted regression with parameter-specific distance metrics and bandwidths. Comput. Environ. Urban 71, 41–57. doi: 10.1016/j.compenvurbsys.2018.03.012

CrossRef Full Text | Google Scholar

Ma, D., Deng, H., Yin, Y., Wu, S., Zheng, D. (2019). Sensitivity of arid/humid patterns in China to future climate change under a high-emissions scenario. J. Geogr. Sci. 29, 29–48. doi: 10.1007/s11442-019-1582-5

CrossRef Full Text | Google Scholar

Meng, F., Zhang, L., Zhang, Z., Jiang, L., Wang, Y., Duan, J., et al. (2019). Opposite effects of winter day and night temperature changes on early phenophases. Ecology 100, e02775. doi: 10.1002/ecy.2775

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, Y., Zeng, X., Chen, W., Tang, X., Dai, K., Du, Y., et al. (2022). Chilling rather than photoperiod controls budburst for gymnosperm species in subtropical China. J. Plant Ecol. 15, 100–110. doi: 10.1093/jpe/rtab076

CrossRef Full Text | Google Scholar

Piao, S., Cui, M., Chen, A., Wang, X., Ciais, P., Liu, J., et al. (2011). Altitude and temperature dependence of change in the spring vegetation green-up date from 1982 to 2006 in the qinghai-xizang plateau. Agr For. Meteorol 151, 1599–1608. doi: 10.1016/j.agrformet.2011.06.016

CrossRef Full Text | Google Scholar

Piao, S., Friedlingstein, P., Ciais, P., Zhou, L., Chen, A. (2006). Effect of climate and CO₂ changes on the greening of the northern hemisphere over the past two decades. Geophys Res. Lett. 33, 0. doi: 10.1029/2006gl028205

PubMed Abstract | 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. Global Change Biol. 25, 1922–1940. doi: 10.1111/gcb.14619

CrossRef Full Text | Google Scholar

Qi, G., Bai, H., Zhao, T., Meng, Q., Zhang, S. (2021). Sensitivity and areal differentiation of vegetation responses to hydrothermal dynamics on the northern and southern slopes of the qinling mountains in shaanxi province. J. Geogr. Sci. 31, 785–801. doi: 10.1007/S11442-021-1871-7

CrossRef Full Text | Google Scholar

Sen, P. K. (2012). Estimates of the regression coefficient based on kendall's tau. J. Am. Stat. 324, 1379–1389. doi: 10.1080/01621459.1968.10480934

CrossRef Full Text | Google Scholar

Shen, M., Jiang, N., Peng, D., Rao, Y., Huang, Y., Fu, Y. H., et al. (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., Piao, S., Cong, N., Zhang, G., Jassens, I. A. (2015). Precipitation impacts on vegetation spring phenology on the Tibetan plateau. Global Change Biol. 21, 3647–3656. doi: 10.1111/gcb.12961

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

Sun, H., Chen, Y., Xiong, J., Ye, C., Yong, Z., Wang, Y., et al. (2022). Relationships between climate change, phenology, edaphic factors, and net primary productivity across the Tibetan plateau. Int. J. Appl. Earth Obs. Geoinf 107, 102708. doi: 10.1016/J.JAG.2022.102708

CrossRef Full Text | Google Scholar

Sun, Y., Guan, Q., Wang, Q., Yang, L., Pan, N., Ma, Y., et al. (2021). Quantitative assessment of the impact of climatic factors on phenological changes in the qilian mountains, China. For. Ecol. Manag 499, 119594. doi: 10.1016/j.foreco.2021.119594

CrossRef Full Text | Google Scholar

Tao, Z., Wang, H., Dai, J., Alatalo, J., Ge, Q. (2018b). Modeling spatiotemporal variations in leaf coloring date of three tree species across China. Agr For. Meteorol 249, 310–318. doi: 10.1016/j.agrformet.2017.10.034

CrossRef Full Text | Google Scholar

Tao, J., Xu, T., Dong, J., Yu, X., Jiang, Y., Zhang, Y., et al. (2018a). Elevation-dependent effects of climate change on vegetation greenness in the high mountains of southwest China during 1982–2013. Int. J. Climatol 38, 2029–2038. doi: 10.1002/joc.5314

CrossRef Full Text | Google Scholar

Vandvik, V., Halbritter, A. H., Telford, R. J. (2018). Greening up the mountain. Proc. Natl. Acad. Sci. U.S.A. 115, 833–835. doi: 10.1073/pnas.1721285115

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitasse, Y., Signarbieux, C., 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

Wang, J., Fang, Y., Zou, X., Chai, G., Wang, Z. (2022). Diverse responses of vegetation phenology to changes in temperature and precipitation in northern China. Geocarto Int., 1–19. doi: 10.1080/10106049.2022.2068677

CrossRef Full Text | Google Scholar

Wang, H., Wang, H., Ge, Q., Dai, J. (2020a). The interactive effects of chilling, photoperiod, and forcing temperature on flowering phenology of temperate woody plants. Front. Plant Sci. 11, 443. doi: 10.3389/fpls.2020.00443

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Wu, C., Ciais, P., Peñuelas, J., Dai, J., Fu, Y., et al. (2020b). Overestimation of the effect of climatic warming on spring phenology due to misrepresentation of chilling. Nat. Commun. 11, 4945. doi: 10.1038/s41467-020-18743-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Xi, Z., He, X., Chen, S., Rossi, S., Smith, N. G., et al. (2021). Contrasting temporal variations in responses of leaf unfolding to daytime and nighttime warming. Glob Change Biol. 27, 5084–5093. doi: 10.1111/gcb.15777

CrossRef Full Text | Google Scholar

Wolf, A. A., Zavaleta, E. S., Selmants, P. C. (2017). Flowering phenology shifts in response to biodiversity loss. Proc. Natl. Acad. Sci. U.S.A. 114, 3463–3468. doi: 10.1073/pnas.1608357114

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolkovich, E., Cook, B., Allen, J., Crimmins, T., Betancourt, J., Travers, S., et al. (2012). Warming experiments underpredict plant phenological responses to climate change. Nature 485, 494–497. doi: 10.1038/nature11014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, L., Ma, X., Dou, X., Zhu, J., Zhao, C. (2021). Impacts of climate change on vegetation phenology and net primary productivity in arid central Asia. Sci. Total Environ. 796, 149055. doi: 10.1016/j.scitotenv.2021.149055

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, H., Qin, Y., Feng, G., Meng, Q., Cui, Y., Song, H., et al. (2019). Forest phenology dynamics to climate change and topography in a geographic and climate transition zone: The qinling mountains in central China. Forests 10, 1007. doi: 10.3390/f10111007

CrossRef Full Text | Google Scholar

Xu, S., Cheng, J. (2021). A new land surface temperature fusion strategy based on cumulative distribution function matching and multiresolution kalman filtering. Remote Sens Environ. 254, 112256. doi: 10.1016/j.rse.2020.112256

CrossRef Full Text | Google Scholar

Yang, K., He, J., Tang, W., Qin, J., Cheng, C. C. K. (2010). On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan plateau. Agr For. Meteorol 150, 38–46. doi: 10.1016/j.agrformet.2009.08.004

CrossRef Full Text | Google Scholar

Yuan, M., Wang, L., Lin, A., Liu, Z., Li, Q., Qu, S. (2020). Vegetation green up under the influence of daily minimum temperature and urbanization in the yellow river basin, China. Ecol. Indic 108, 105760. doi: 10.1016/j.ecolind.2019.105760

CrossRef Full Text | Google Scholar

Yu, H., Luedeling, E., 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

Zhang, Q., Cheng, J. (2020). An empirical algorithm for retrieving land surface temperature from AMSR-e data considering the comprehensive effects of environmental variables. Earth Space Sci. 7, 1–26. doi: 10.1029/2019EA001006

CrossRef Full Text | Google Scholar

Zhang, Y., Keenan, T. F., Zhou, S. (2021b). Exacerbated drought impacts on global ecosystems due to structural overshoot. Nat. Ecol. Evol. 5, 1490–1498. doi: 10.1038/s41559-021-01551-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Liang, C. (2020). Analysis of annual and seasonal precipitation variation in the qinba mountain area, China. Sci. Rep-Uk 10. doi: 10.1038/s41598-020-57743-y

CrossRef Full Text | Google Scholar

Zhang, C., Ren, Y., Cao, L., Wu, J., Zhang, S., Hu, C., et al. (2022). Characteristics of dry-wet climate change in China during the past 60 years and its trends projection. Atmosphere 13, 275. doi: 10.3390/atmos13020275

CrossRef Full Text | Google Scholar

Zhang, Q., Wang, N., Cheng, J., Xu, S. (2020). A stepwise downscaling method for generating high-resolution land surface temperature from AMSR-e data. IEEE J-Stars 13, 5669–5681. doi: 10.1109/JSTARS.2020.3022997

CrossRef Full Text | Google Scholar

Zhang, R., Wang, F., Zheng, J., Lin, J., Hänninen, H., Wu, J. (2021a). 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

Zheng, L., Shi, P., Song, M., Zhou, T., Zong, N., Zhang, X. (2021). Climate sensitivity of high altitude tree growth across the Hindu kush himalaya. For. Ecol. Manag 486, 118963. doi: 10.1016/j.foreco.2021.118963

CrossRef Full Text | Google Scholar