Phenology of Five Shrub Communities along an Elevation Gradient in the Qilian Mountains, China

Abstract

Phenology indicates the response of ecosystem dynamics to climate change. Shrubs are an important component of alpine forests, and play a key role in forest ecosystem function, especially in semiarid alpine regions. In 2015, we observed the dynamics of phenology in five shrub communities along an elevation gradient (2600–3300 m) in the Qilian Mountains. Our results showed that the length of the growing season decreased from 153 days for Caragana tangutica to 150 days for Berberis diaphana, 149 days for Potentilla fruticosa L., 144 days for Caragana jubata (Pall.) Poir., and 134 days for Salix gilashanica C. Wang et P. Y. Fu. The accumulated temperature of the five shrub communities during the growing season varied from 1735.4 °C for C. tangutica to 1051.3 °C for C. jubata. The beginning of the growing season was earlier at lower, than at higher, elevations, while the end of the growing season was later at lower, than at higher, elevations. Elevation and aspect were two important environmental factors that affected shrub phenology. In our study, low temperature, coinciding with the occurrence of early frost, particularly in higher elevations, was the key factor in promoting end-of-season shrub growth cessation.

1. Introduction

Vegetation phenology is one of the most sensitive indicators of climate change [1,2,3]. Temperate forests represent one of the major types of vegetation at middle to high latitudes in the northern hemisphere, and globally act as a large carbon sink [4,5]. Thus, a better understanding of the phenological variables in temperate forests will improve the accuracy of vegetation models and estimates of regional carbon fluxes [6,7,8,9]. Although phenological changes may vary with species and geographic locations, numerous observations have shown that the advances in spring are followed by delays in autumn phenology, with a corresponding lengthening of the plant growing season [10,11,12,13]. These shifts, which are correlated with changes in temperature, are ecological responses to anthropogenic global change [14,15,16,17].

Northwestern China has experienced a rapid temperature increase of 0.26 °C per decade during the past 50 years, higher than the national average increase of 0.14 °C per decade [18]. Climate anomalies in the area resulted in profound changes in subalpine shrubs, including alterations in the temporal niche of phenophases and the dynamic interaction with ambient conditions [19]. A study of the relationships between shrub phenology and meteorological factors will help us understand the characteristics and driving mechanisms of the responses of shrub phenology to climate change under different climatic conditions [20].

Studies concerning the response of plant phenology to global climate change can be divided into two categories. Some used the Normalized Difference Vegetation Index (NDVI) to distinguish different vegetation types for simulation, while others used averages of the annual, seasonal, and monthly temperatures, or the maximum and minimum temperatures, combined with observations of plant phenology [19,21,22,23,24,25,26]. However, phenological studies should also emphasize measurements of leaf area indices, or radiation above and below the canopy to determine the relationship between environmental factors and leaf phenology. In addition, litter fall and leaf fall observations should also be given priority [27]. Soil moisture may also affect the phenology of woody cover [3]. In water-limited areas, variations of precipitation are expected to modify phenological patterns of species, while higher temperatures may lead to an advanced spring phenology [28,29].

Accurate estimates of canopy phenology are critical for quantifying carbon and water exchange between forests and the atmosphere, and predicting the response of forest ecosystems to climate change [26,30]. Previous studies of phenology in the Qilian Mountains in Northwestern China were largely based on remote sensing data, while field observations of phenology were scarce due to limited accessibility [18,19,31,32]. Although field measurements are time-consuming and laborious, their accuracy is unmatched by remote sensing observations. In this field study, we focused on the temperate subalpine shrublands in a semiarid alpine catchment in the Qilian Mountains. The phenophase variation in five typical shrubs along an elevation gradient was observed during the 2015 growing season. The objectives of this study were to (1) investigate the variability in phenology of shrubs along an elevation gradient; and (2) determine the relationships between temperature, topography, and phenophase changes.

2. Materials and Methods

2.1. Study Site

The study was conducted at five shrubland sites along an elevation gradient in a typical catchment of the Qilian Mountains (100°17′–100°18′ E, 38°32′–38°33′ N). The area of the catchment is 2.85 km² at an elevation from 2500–3796 m. The catchment is characterized by a semiarid climate with an annual average temperature of 1.6 °C, and annual average precipitation in 2015 of 353.1 mm at the base of the mountains. From 2013–2015, seasonal precipitation was uneven, and precipitation from May to September accounted for almost 85% of that for the year. In 2015, the duration of annual sunshine reached 1892.6 h, and annual evaporation was 1081.7 mm. Spatial variability in topography, precipitation, and temperature led to significant differences in soil properties and vegetation communities. Along the elevation gradient, alpine meadow soil exhibited the highest accumulation of organic matter, followed by grey cinnamon soil, chestnut soil, and subalpine shrub meadow soil. Organic carbon and total nitrogen decreased with soil depth. From lower elevation to higher elevation, the vegetation types are mountain, forest grasslands, subalpine shrub meadows, and alpine meadows, respectively. The corresponding soil types are chestnut soil, grey cinnamon soil, subalpine shrub meadow soil, and alpine meadow soil, respectively [33]. Forests, dominated by Qinghai spruce (Picea crassifolia), are mainly distributed on the shady and semi-shady slopes at 2500–3300 m. Shrublands mainly occupy footslope positions of shady slopes and areas above the tree line. The main shrub species include Caragana tangutica, Berberis diaphana, and Potentilla fruticosa L. widely distributed at elevations <2900 m, and Caragana jubata (Pall.) Poir., and Salix gilashanica C. Wang et P. Y. Fu found above 3100 m. The five shrubs are woody perennials. The main herbaceous plants are Polygonum viviparum L., Carex atrata L., and Stipa capillata Linn. Study area with sample plots (numbers 1–5) and main shrub species was showed in Figure 1. Sample plots information was showed in

Table 1. Sample plot information.

Shrub Species	Soil Depth (cm)	Soil Types	Elevation (m)	Slope Degree (°)	Slope Aspect	Canopy Cover (%)
Caragana. tangutica	60	Chestnut soil	2600	26	SW	50
Berberisdiaphana	50	Chestnut soil	2600	30	W	65
Potentilla Fruticosa L.	60	Meadow soil	2900	31	E	80
CaraganaJubata (Pall.) Poir.	50	Meadow soil	3300	38	NE	65
Salix Gilashanica C. Wang et P. Y. Fu	50	Meadow soil	3300	36	NE	75

.

2.2. Phenological Observations

This work was guided by the “Observation Methodology for Long-term Forest Ecosystem Research” developed by the National Standards of the People’s Republic of China (GB/T 33027-2016). We selected five species for our study: C. tangutica, B. diaphana, P. fruticusa, C. jubata, and S. gilashanica, growing along an elevation gradient at 2600 m to 3300 m in 2015. We chose a growing season as our aim to investigate the effects of topography on forests under climate change. Three quadrats (10 m × 10 m) were established for each shrub species. We randomly chose five individual shrubs exhibiting normal growth, flowering, and fruiting in each quadrat to observe changes in phenophases as per GB/T 33027-2016. A total of 15 individual shrubs of each species were tracked in the entire elevation gradient. Phenological observations began at budding date, and ended when the leaves started to fall, spanning a growing season. At the beginning and end of the growing season, phenological observations were made daily. During the growing season, we adjusted the observation intervals because of faster growth of shrubs in the middle growing season. In addition, air temperature data were obtained from the nearby automated meteorological stations at 2600, 2900, and 3300 m.

2.3. Sample Plot Investigation

Three quadrats per species were established to investigate shrub height, and basal and crown diameters. The elevation, latitude, and longitude of the quadrat were obtained using a global positioning system. We used a compass to measure the slope aspect, and a declinometer to measure the slope degree. The soil profile was excavated to determine soil type and depth. The line transect method [20] was used to measure percent cover of five shrub species in this study.

2.4. Data Processing

The methodology to calculate the duration of each phenological stage included transforming the date and time to Julian date (1 January was day 1; 31 December was day 365) [34], and we used Julian days for dates of budding, leaf opening, flowering, fruiting, leaf discoloration, and leaf falling. The length and T_acc of these phenophases were calculated based on the formulas below [13,34]:

$$d_i=b-a$$

(1)

$$a{t}_i={\sum }_{j=a}^b{t}_j$$

(2)

where d_i was the length of phenophase i; a was the starting day of phenophase i; b was the end day of phenophase i; at_i was the T_acc of phenophase i; and t_j was the average temperature of phenophase i on day j. We used Microsoft Excel 2010 (Microsoft Corporation, Redmond, Washington, WA, USA) for data processing and mapping, and SPSS 21.0 software (SPSS Inc., IBM, Chicago, IL, USA) for statistical measures.

3. Results

3.1. Growth Parameters of Five Shrub Species

There were significant differences in height and canopy diameter of the five shrub species (p = 0.05) (

Table 2. Basic parameters of five shrubs species.

Shrub Species	Elevation (m)	Aspect	Shrub Height (cm)	Basal Diameter (cm)	West-East Crown Diameter (cm)	North-South Crown Diameter (cm)
C. tangutica	2600	Sunny	180.2 ± 5.8	1.3 ± 0.0	178.7 ± 9.4	183.9 ± 9.2
B. diaphana	2600	Semi-sunny	154.8 ± 4.8	1.2 ± 0.1	144.0 ± 5.5	156.9 ± 5.4
P. fruticosa	2900	Semi-sunny	45.5 ± 1.8	0.5 ± 0.0	44.5 ± 3.3	42.3 ± 3.0
C. jubata	3300	Semi-sunny	92.3 ± 2.2	1.2 ± 0.3	75.4 ± 3.7	89.2 ± 3.3
S. gilashanica	3300	Sunny	119.4 ± 3.5	1.1 ± 0.3	92.3 ± 5.0	90.5 ± 4.3

Shrub height, basal diameter, West-East crown diameter, and North-South crown diameter presented in Table 2 are expressed as mean ± S.D. C. tangutica: n = 26, P. fruticosa: n = 19, P. fruticosa: n = 57, C. jubata: n = 48, and S. gilashanica: n = 45.

). The branch diameter of P. fruticosa was significantly lower than that of the other four shrub species. Shrubs at lower elevations were taller than those at high elevations. In addition, the shrubs growing on the sunny slope were significantly taller than shrubs growing on semi-sunny slopes at the same elevation (Table 2).

3.2. Duration of Phenological Stages of the Five Shrub Species

All observed shrub species generally started to bud from April to May (

Table 3. Start and end dates of specific phenophases and their duration for five shrub species in 2015.

	Shrub Species	Budding		Leaf Opening		Flowering		Fruit Ripening		Leaf Discoloration		Leaf Fall
		Start	End	Start	End	Start	End	Start	End	Start	End	Start	End
Date (Month-Day)	C. tangutica	4–26	5–12	5–13	6–11	6–12	7–15	7–16	9–1	9–2	9–15	9–16	9–25
	B.diaphana	4–28	5–5	5–6	5–18	5–19	6–28	6–29	8–25	8–26	9–15	9–16	9–24
	P. fruticosa	5–1	5–15	5–16	6–30	7–1	7–31	8–1	8–14	8–15	9–15	9–16	9–26
	C. jubata	5–5	5–21	5–22	7–1	7–2	8–16	8–17	8–25	8–26	9–1	9–2	9–15
	S. gilashanica	5–5	5–21	5–22	6–24	6–25	8–14	8–15	8–25	8–26	9–14	9–15	9–25
Duration (days)	C. tangutica	17		30		34		48		14		10
	B. diaphana	8		13		41		58		21		9
	P. fruticosa	15		46		31		14		32		11
	C. jubata	17		41		46		9		7		14
	S. gilashanica	17		34		51		11		20		11

). C. tangutica at 2600 m budded about eight days earlier than C. jubata and S. gilashanica at 3300 m. B. diaphana at 2600 m was the earliest to open leaf, bear fruit, and change leaf color. The length of the growing season in the five shrub species decreased in the order of: C. tangutica (153 days), B. diaphana (150 days), P. fruticosa (149 days), C. jubata (144 days), and S. gilashanica (134 days). The average length of the growing season was 146 days, which was positively correlated to the average temperature of the growing season at the study site (R² = 0.8973) (Figure 2).

Interspecific differences and slope aspect may have led to the shorter growing season in C. jubata than that in S. gilashanica. Similarly, the length of the growing season in B. diaphana was less than that in C. tangutica. Five shrub species began to change color in late August, and leaves started falling in the middle of September as the lowest temperature was a key driving factor leading to a dormant state based on our observations.

3.3. Variability in Temperature and T_acc

3.3.1. Changes in Temperature

Daily average and minimum temperatures decreased with increasing elevation (Figure 3). Both temperatures peaked in late August, and started to drop in early September (Figure 3) (Supplementary Materials).

C. tangutica and B. diaphana began to change leaf color when the minimum temperature reached 9.5 °C, and their leaves began to fall when the minimum was 4.4 °C. P. fruticosa began to change leaf color when the minimum was 5.6 °C, and its leaf fall started when the minimum was 2.7 °C. C. jubata and S. gilashanica began to change leaf color when the minimum was 5.5 °C. Leaf fall in C. jubata began when minimum T was 2.3 °C, and in S. gilashanica, when it was −2.9 °C. The first frost occurred on 18 September. At this time, the minimum temperatures were 2.8 °C, 1.1 °C, and −1.2 °C at 2600 m, 2900 m, and 3300 m, respectively.

3.3.2. Variability in T_acc in 2015

The start of growing season varied from the 26 April in C. tangutica, to the 5 May for C. jubata and S. gilashanica. The end of growing season varied from 15 September in C. jubata, to 26 September for P. fruticosa (Table 3).

T_acc began to increase at the end of April, when the shrubs began to bud break. Therefore, T_acc is an important factor promoting shrub growth. In late September, T_acc no longer increased, and growth cessation began and shrubs started to become dormant. At this time, low temperature was an important factor promoting shrub dormancy. Fruit of all five species started to ripen in early August as the T_acc reached the maximum (Figure 4).

3.3.3. T_acc of Different Phenophases

Shrub phenophases changed from one to the next only when temperature accumulated to specific levels (

Table 4. T_acc of individual phenophases and their duration for each shrub.

	Shrub Species	Tacc Phenophase
		Budding		Leaf Opening		Flowering		Fruit Ripening		Leaf Color Change		Leaf Fall		Growing Season
		Start	End	Start	End	Start	End	Start	End	Start	End	Start	End
Tacc at start and end (°C)	C. tangutica	123.2	232.1	244.2	520.1	530.5	950.9	963.8	1659.8	1674.5	1781.1	1789.9	1858.6	1735.4
	B. diaphana	144.7	195.0	202.3	299.2	312.3	729.7	742.6	1555.5	1570.3	1781.1	1789.9	1854.2	1709.5
	P. fruticosa	119.6	193.0	201.2	603.3	615.4	1000.5	1018.8	1187.0	1199.7	1494.9	1501.9	1559.3	1439.7
	C. jubata	75.5	151.3	154.6	423.7	434.5	916.5	919.8	986.4	997.2	1062.4	1073.0	1126.8	1051.3
	S. gilashanica	75.5	151.3	154.6	351.0	359.3	893.2	903.6	986.4	997.2	1120.5	1126.8	1163.7	1088.2
Tacc for the duration of phenophase (°C)	C. tangutica	108.9		275.9		420.4		696.0		106.6		68.7		-
	B. diaphana	50.3		96.9		417.4		812.9		210.8		64.3		-
	P. fruticosa	73.4		402.1		385.1		168.2		295.2		57.4		-
	C. jubata	75.8		269.1		482.0		66.6		65.2		53.8		-
	S. gilashanica	75.8		196.4		533.9		82.8		123.3		36.9		-

). C. jubata and S. gilashanica began to bud at lower T_acc than C. tangutica and B. diaphana. Specifically, budbreak in C. jubata and S. gilashanica occurred at 75.5 °C, while that in B. diaphana occurred at 144.7 °C. C. tangutica began to leaf fall at higher T_acc than C. jubata. Specifically, leaf fall in C. tangutica occurred at 1858.6 °C, while that in C. jubata occurred at 1126.8 °C.

T_acc for both, the beginning and end of each phenophase in high-elevation shrubs was lower than that in low-elevation shrubs. T_acc on the sunny slope was higher than that on the semi-sunny slope. During the growing period at 2600 m, T_acc within each phenophase increased before fruit ripening, and decreased after that. A similar pattern was observed in shrubs at 2900 m and 3300 m (Table 4) (Supplementary Materials).

4. Discussion

4.1. Phenophases of Shrub Species

The beginning of the growing season in the Qilian Mountains fell between the 120th and 160th day of the year, that is, between early May and early June [31]. In general, vegetation in the Qilian Mountains showed a trend of greening from the southeast to the northwest [32]. This trend may be associated with a decrease in precipitation from the southeast to the northwest. The results of our study showed that the average growing season of typical shrubs in the Qilian Mountains started on the 122nd day, and ended on the 267th day, with the average length of the growing season of 146 days. Zhao et al. (2015) showed that the growing season in the Qilian Mountains started on the 139th day, and ended on the 250th day, with the average length of the growing season spanning 112 days [31]. Temperature and species differences mainly determined the differences in growth and phenology of the five shrub species in our study.

4.2. Effects of Temperature on Shrub Phenology

T_acc is one of the most important factors affecting plant growth and, therefore, is of great significance in ecology and forest management [34]. Temperature is the main driver of phenological dynamics for many plants; in many cases, higher temperature is likely to accelerate plant development and lead to earlier switching to the next ontogenetic stage [8]. During our field observation, we found snow and frost at the beginning and end of the growing season, respectively. When snowfall occurred, vegetation growth slowed. Frost occurred in early September, and was followed by the onset of leaf drop. We speculated that the low temperature had a profound effect on phenological dynamics. A similar conclusion was drawn by Du et al. (2014) and He et al. (2015), who demonstrated the primary role of minimum temperature in controlling the start and the end of the growing season [18,19]. Gough et al. (2010) suggested that the end of season phenology may be less responsive to environmental variability than to genetics [5]. Furthermore, precipitation was coupled to a shift in maximum NDVI [28]. The restrictive role of minimum temperatures in the growth of subalpine plants requires further study as it is expected to be closely associated with the frequency of frost risk in spring and autumn [35,36].

Our study showed that the T_acc within each of the five phenophases in the five shrub communities first increased, then declined. The variability in T_acc was consistent with the length of the phenophase. Chang et al. (2012) also indicated that the variation in T_acc was largely caused by the length of the phenophase, while the temperature-induced variation was comparatively less. Chang et al. (2016) found that the length of phenophase could not be shortened or extended indefinitely [34], which was consistent with our study.

4.3. Effects of Topography on Shrub Phenology

4.3.1. Elevation

In this study, the growing season for C. tangutica and B. diaphana at 2600 m started earlier and ended later than that for C. jubata and S. gilashanica at 3300 m. As the elevation decreased, the beginning and end of the growing season showed a trend of advance and delay, respectively, and the length of the growing season increased. This was possibly associated with the limited heat at higher elevations.

The spatial patterns of vegetation phenology in the Dan Jiang Kou reservoir area indicated that the relatively longer length of season in the south resulted from an early start and late end of season. Regression models and correlation analysis indicated that elevation was moderately related to vegetation [32]. Further, a close connection was established between the growth of the main vegetation types and precipitation in the mid-Qilian Mountains, while temperature was the vital climatic factor for some special vegetation types, such as alpine meadows and coniferous forests, located in regions with high elevation or sufficient water resources [37]. Elevation may also be a key factor affecting vegetation phenology due to the different thermal conditions at different elevations.

4.3.2. Slope Aspect

In a mountainous ecosystem, slope aspect is also important topographic factors affecting water and heat dynamics [32]. Vegetation growth on shady slopes may be restricted by heat conditions, while increases in precipitation and warming may enhance vegetation productivity in such areas; the duration and flux intensity of solar radiation vary from sunny to shady slopes [38]. However, on sunny slopes in mountainous areas, increased evapotranspiration and an effect of soil water stress may be more significant than on shady slopes [3,39]. Yin et al. (2016) found that topographic variables, such as slope gradient and growing-season direct solar radiation may have minor influences on the growth-climate relationships [38]. Our study indicated that the growing season of C. tangutica on the sunny slope started earlier than that of B. diaphana on the semi-sunny slope. This result was consistent with previous results; however, the growing season of C. jubata on the sunny slope was shorter than that of S. gilashanica on the semi-sunny slope, which may be caused by interspecific differences.

Many other climatic factors need to be considered in any future studies; these may include potential changes in soil moisture and temperature, duration of insolation, effective solar radiation, and the depth of frozen soil. Moreover, identifying and quantifying climatic constraints on phenological dynamics in forests is a complex task that needs to be addressed at appropriate scales. This study is more or less descriptive in nature because of the limited sample sizes. Thus, we will choose greater sample sizes to study the relationship between phonology and climatic change in future.

5. Conclusions

In our study, the lowest temperature may be the environmental factor precipitating leaf discoloration and fall in shrubs, and the highest temperature may be the environmental factor leading to leaf expansion and fruit ripening. Shrub phenology will progress from the previous to the next phenophase as long as the T_acc reaches a certain level. The start date of high-elevation growth is relatively late, while the end date is relatively early. Low temperature at the beginning of the growing season may not be conducive to shrub bud break. Lower temperatures at the end of the growing season also promote shrub growth cessation, which may lead to a shortening of the growing season at high elevations. Frostbite may affect the growing season of shrubs and result in a shorter growing season at higher elevations. In addition, the slope aspect may be an important factor affecting plant phenology due to the difference in insolation and illumination among aspects.