Among the various plant traits, phenology is one of the most sensitive to environmental changes. Nitrogen (N) has been reported to be an important factor constraining plant growth, while N deposition is an increasing concern due to social and economic development. However, studies exploring influences of atmospheric N deposition on alpine plant phenology have been barely reported. We conducted a N addition gradient experiment in an alpine meadow on the plateau to examine plants' phenological response. The results suggest that N addition plays an important role in regulating plant phenology, but only for the dominant species; specifically, the budding time of Kobresia pygmaea advanced in response to N addition. The prolonged lengths of flower phase and reproductive period of K. pygmaea were mainly ascribed to an earlier start of the growing season and more resources being allocated to reproductive growth in the N-rich soil. The N addition caused no effect on the phenology of the accompanying species. There were no obvious phenology differences among each N addition treatment. The earlier budding time and the prolonged growing season of the dominant species in response to the N addition means a longer provision duration of food sources to local herders, which can facilitate local economic development.
Phenology is the study of periodic plant and animal life cycle events and how they are influenced by seasonal and interannual climate variations (Parmesan, 2006). Plant phenology, usually composed of plant budding, leaf unfolding, anthesis, fructification, and senescence, is one of the most sensitive traits in response to environmental change (Parmesan and Yohe, 2003; Cleland et al., 2007; Post and Inouye, 2008). In minimizing competition for pollinators and other resources, plants employ diversified phenology, which is an important mechanism for plants to maintain their coexistence in plant communities (Rathcke and Lacey, 1985; Martinkova et al., 2002). The timing of growth onset and senescence may also determine the growing season length, thus driving annual carbon uptake in terrestrial ecosystems (Cleland et al., 2007; Forrest and Miller-Rushing, 2010).
In recent decades, effects of warming on vegetation phenophase have been plentifully reported (Ahas et al., 2000; Sparks and Menzel, 2002; Piao et al., 2007; Zhang et al., 2013a). Besides temperature, other factors like soil moisture, nitrogen (N) availability, light, and elevated CO2 are all likely to affect plant phenology (Nord and Lynch, 2009). What is more, other environmental factors may confound the temperature sensitivity of plant phenology in natural ecosystems (Wolkovich et al., 2012).
N is a basic nutrient for plants. However, an excessive usage of N can lead to eutrophication, thus changing ecosystem structure and function. The consequences can be altered growth, altered organism chemistry composition, shifted phenology, changed nutrient cycling and, ultimately, alteration in plant community composition (Xia and Wan, 2008; Yang et al., 2012). It is estimated that the global land surface has received approximately doubled N input as a result of the worldwide application of artificial N fertilizers (Gruber and Galloway, 2008). Some researchers have reported advanced or delayed plant phenology due to N addition (Cleland et al., 2006; Power et al., 2006; Lupi et al., 2012; Xia and Wan, 2013). But related studies have been barely reported for ecologically fragile areas, such as the alpine meadows of the Tibetan Plateau (Smith et al., 2012).
The Tibetan Plateau is among the most sensitive regions to global changes because of its unique geographical environment, including high altitude, strong solar radiation and low air pressure (IPCC, 2007). During the past 30 years, the N deposition has climbed up from 1 kg N ha-1 to 5 kg N ha-1 (Liu et al., 2013). The Tibetan Plateau alpine meadow system is strongly restricted by N availability, partially due to a slow N mineralization. Exogenous N input can increase resource availability and lessen the nutrient restriction for the alpine ecosystem (Bowman et al., 2006). Previous studies have reported changed biomass allocation strategy and de creased species abundance due to an extra N input and increased N mineralization around the world (Luo et al., 2006; Niu et al., 2008; Bai et al., 2011). However, little is known about how plant phenology responds to N addition, especially in alpine environments. To date, most published findings about the phenology response on the Tibetan Plateau have been focused on effects of warming at a regional scale (Yu et al., 2010; Piao et al., 2011; Zhang et al., 2013a) or a site scale (Klein et al., 2014).
The traits of each species monitored in this study.
Plant phenology plays an important role in maintaining species coexistence (Miller-Rushing and Inouye, 2009; Forrest and Miller-Rushing, 2010). Studies on the responses of phenology to N deposition can shed light on its effects on community structure. To this end, we conducted a gradient fertilization experiment and monitored phenological phases of seven species from different functional types. Specifically, we are trying to answer the following two questions. (1) How does N addition affect the plant phenology of an alpine meadow on the Tibetan Plateau? (2) Will multiple-level N additions cause differences in phenology responses?
Material and Methods
The experiment was conducted in a typical alpine meadow (31°38.513′N, 92°0.921′E, 4585 m a.s.l.) on the Tibetan Plateau. The long-term (1955–2012) annual mean temperature and precipitation are -1.16 °C and 430 mm, respectively. The annual cloud-free solar radiation is about 2788 h. There is no absolute frost-free period throughout the year. The snow period or soil freezing period is from October through following May. The duration of the growing season is from June to August. The vegetation, typical of Tibetan alpine meadows, is dominated by Kobresia pygmaea, accompanied by Potentilla saundersiana, Youngia simulatrix, and P. cuneata, etc. The soil is classified as meadow soil with sandy loam. To exclude livestock grazing effects, the experimental site has been fenced since October 2011.
The experimental plots were established on a fairly flat terrain to minimize the precipitation leaching effects on added N. The experiment used a complete block design with four N treatment levels, including N addition at 0 (control), 2.5, 5, and 10 g N m-2 yr-1 levels (added as urea). Hereafter, the treatments were denoted as N0, N2.5, N5, and N10. The four treatments were randomly arranged within a block and replicated five times. Thus, twenty 6 × 6 m plots were arranged in a 4 × 5 matrix, with a 1 m buffer zone between adjacent plots and a 2 m buffer zone between adjacent blocks. The pelletized fertilizer was hand-broadcasted onto plots twice a year, each time with 50% of the total applied before rain events in the middle of June and July, respectively. The plots have been fertilized since 2012. As monitored, there was no rainstorm right after fertilizer application. Also, the flat terrain can minimize fertilizer leaching. To check the fertilizing efficiency, we measured ecosystem carbon exchange (net ecosystem exchange [NEE], ecosystem respiration [Reco], and gross primary productivity [GPP]) for each treatment.
CLIMATE CONDITIONS AND SOIL MOISTURE
We used precipitation and air temperature data from a meteorology station ~18 km away from the experimental site. Soil moisture was measured every seven days from June to September at 3.8, 7.6, and 12 cm depths in all 20 plots using the TDR 300 (Spectrum Technologies, Plainfield, Illinois, U.S.A.).
We monitored phenological phases of seven species, including five forbs (P. saundersiana, P. cuneata, P. bifurca, Y. simulatrix, and Saussurea stolickai), one graminoid specie (Festuca coelestis), and one sedge (K. pygmaea) (Table 1). K. pygmaea accounts for about 70% of the total community aboveground biomass and 90% of the total community coverage, respectively. They are all anemophilous flowers. The shallow-rooted K. pygmaea begins to blossom in June. It has a fibrous rooting system, which is mostly confined to the top 10 cm surface soil. Also, as a shallow-rooted plant, Potentilla begins to blossom in June. They have small taproots (generally smaller than 3 mm in diameter), which normally do not penetrate into soil deeper than 10 cm. The deep-rooted F. coelestis begins to blossom in July. It has a strong fibrous rooting system. Y. simulatrix and S. stolickai, as deep-rooted plants, begin to blossom in late July. These two species have strong taproots, which can reach 30 cm deep into soil (Lv et al., 2009; Dorji et al., 2013).
We randomly flagged 15 individuals of each species in total at the start of the growing season, with 3 individuals in each plot for each species, and monitored their phenological phases at approximately 4 day intervals. We used the scoring method to record their phenology (Price and Waser, 1998; Dunne et al., 2003). For forbs, we used six phenological stages: “0” for plant turning green, “0.5” for most culms being visible, “1” for unopened buds, “2” for open flowers, “3” for initiated fruit, “4” for enlarged fruit, and “5” for dehisced fruit. For grass and sedge, we used five phenological stages: “0” for plant turning green, “1” for appearance of spikelets, “2” for exerted anthers and styles from the spikelet florets, “3” for developed seed, and “4” for disarticulated seeds. Because most of P. bifurca, S. stolickai, and Y. simulatrix have no complete reproductive growth in our study area, we defined vegetative growth from 0 to 0.9, and leaf withering from 5.1 to 6.
It is difficult to continuously monitor the exact timing of each phenological stage. We obtained them by fitting the growth curvatures of the Richards growth equation against the day in Julian units (X) for each plot. Based on our monitoring, P. saundersiana, P. cuneata, F. coelestis, and K. pygmaea almost completed their reproductive growth. So we fitted the observation scores to the Richards growth equation directly (Richards, 1959; Sherry et al., 2007; Xia and Wan, 2013).
For the Richards growth equation:
Once the parameters of the Richards equation were obtained for each plant, budding, flowering, or fruiting timing can be calculated from the equation, and then also the plot-average. The budding and flowering timing was calculated as Y = 1 and Y = 2 for all species, respectively. The fruiting time was calculated as Y = 3 and Y = 2.5 for forbs and grass/sedge, respectively. The phenological duration equaled the phase Y = 4 for the forbs minus the phase Y = 1 (or phase Y = 3 minus phase Y = 0.5 for grass and sedge).
Few individuals of P. bifurca, Y. simulatrix, and S. stolickai started or completed their reproductive growth, so we fit the Richards growth equation model on the phenological scores for vegetative growth and leaf withering, separately. Using the Richards equation, we fitted the phenological scores to the two models. The phenological duration was obtained as Y = 5.1 from the second model minus Y = 0.5 from the first model.
Treatment differences were tested using the one-way analysis of variance (ANOVA), and the least significant difference (LSD) test was employed for the multiple factor comparisons. Soil moisture difference during each month was compared by the independent-sample-T-test. The fittings of the Richards equation were carried out in Matlab (Mathworks, Natick, Massachusetts, U.S.A.) with the contraction-expansion algorithm (Gu et al., 1998), and all statistical analyses were conducted using SPSS software (SPSS, Chicago, Illinois, U.S.A.).
CLIMATE CONDITION AND SOIL MOISTURE
In comparison with the long-term (1955–2012) mean (-1.16 °C), the mean air temperature in 2013 (-0.1 °C) was about 1 °C higher. Meanwhile, total precipitation in 2013 (478 mm) was 11.2% higher than the long-term mean (430 mm) (Fig. 1, part A). However, total precipitation over the entire growing season (June to August) in 2013 (274 mm) was 3.5% lower than the long-term mean (284 mm).
The seasonal dynamics of air temperature followed a unimodal pattern, with a peak value in summer (Fig. 1, part B). The mean maximum air temperature (16.8 °C) is almost four times the minimum (4.4 °C) over the growing season. Precipitation mostly occurred in July as rainfall and September as snowfall.
Soil moisture fluctuated with precipitation during the growing season (Fig. 2, part A). Soil moisture was significantly higher in July than in June and August at each depth (P < 0.01). Soil moisture at the 7.6 cm depth was slightly higher than in other depths. The multiple-level N additions had no distinct effects on soil moisture (Fig. 2, parts B, C, and D).
NITROGEN ADDITION EFFECTS ON ECOSYSTEM CARBON EXCHANGE
During peak growing seasons in 2013, NEE under N2.5, N5, and N10 treatment significantly increased by 64% (P = 0.011), 54% (P = 0.027), and 92% (P < 0.001), respectively. Different N additions also significantly increased GPP by 44% (P = 0.002), 42% (P = 0.003), and 71% (P < 0.001), respectively. Only the N10 treatment significantly increased ecosystem Reco by 49% (P = 0.045).
NITROGEN ADDITION EFFECTS ON FLOWER PHASE
Among the seven species observed, most of K. pygmaea, F. coelestis, P. saundersiana, and P. cuneata completed their reproductive growth and the N addition resulted in no obvious difference on the flowering time of the four species (Fig. 3, part A). Different levels of N additions failed to cause distinct flowering time responses from the four species. Two response patterns were detected in the flower phase length. The length of K. pygmaea was elongated by 0.7 d and the elongation was statistically significant (P < 0.05); the length of F. coelestis, P. saundersiana, and P. cuneata showed no responses (Fig. 3, part B).
NITROGEN ADDITION EFFECTS ON FRUIT PHASE
There was no significant difference of fruiting time and fruit phase length between the N additions and the control treatments (Fig. 4, part A). In the meantime, there were no significant difference of fruiting time and fruit phase length among each level of N addition treatment (Fig. 4, part B).
NITROGEN ADDITION EFFECTS ON GROWTH DURATION
The reproductive phase is composed of bud, flower, and fruit phases. Among the seven studied species, K. pygmaea, F. coelestis, P. saundersiana, and P. cuneata had reproductive phases, while most of P. bifurca, S. stolickai, and Y. simulatrix failed to start reproductive growth in 2013 (Fig. 5). Different levels of N additions failed to cause distinct growth duration responses from the seven species.
For the budding time or the start of vegetative phase of the seven species, only K. pygmaea showed significantly positive responses to N addition (Fig. 5, part A), whose budding time advanced by 5.9 d (P < 0.05).
In terms of ending the vegetative phase, only S. stolickai exhibited statistically different responses to N addition (Fig. 5, part A). It advanced the end of the vegetative growth by 3.3 d (P < 0.05). Other species showed no change in ending the vegetative growth in response to N addition.
Due to the changes in budding timing and reproductive growth termination (the start and the end of vegetative growth), the reproductive (vegetative) duration of the seven species responded differently to N addition (Fig. 5 part B,). Among the seven species, only K. pygmaea stretched the duration by 6.0 d in response to N addition (P < 0.05), while other species showed no change.
NITROGEN ADDITION EFFECTS ON SPECIES-LEVEL PHENOLOGY
As revealed in this study, plant species responded differently to the short-term N addition. The responses may depend on their morphological, physiological, and life-history traits (Kliber and Eckert, 2004). Specifically, the phenology of K .pygmaea responded to the N addition in a more obvious manner than other species.
We monitored phenological phases of seven species, including five forbs (P. saundersiana, P. cuneata, P. bifurca, Y. simulatrix, and S. stolickai), one grass (F. coelestis), and one sedge (K. pygmaea). Some other phenological studies have observed N-addition delayed flowering in grasses but accelerated flowering in forbs (Cleland et al., 2006; Smith et al., 2012; Xia and Wan, 2013). In our study, the phenology of the grass and forbs showed no response to N additions. The plant traits associated with resource acquisition, such as rooting depth and life history might mediate plant phenology responses to changing nutritional conditions. K. pygmaea, as a shallow-rooted and dominant species, could be more benefited by fertilizing-induced increases in soil N compared to other species. However, also as shallow-rooted plants, Potentilla exhibited no response to N additions. The reason might lie in their disadvantage in competing for nutrition resources compared to the dominant species (Miller and Bowman, 2003; Miller et al., 2007; Ashton et al., 2010). F. coelestis, Y. simulatrix, and S. stolickai, as deep-rooted and late flowering plants, might be slow in absorbing the added N.
On the other hand, the N addition treatment in this study has been conducted only for one year and the studied plants are all perennial. Their belowground reproductive growths can last one or more years, thus delaying an obvious aboveground response to altered environmental conditions (Mark, 1970; Aydelotte and Diggle, 1997; Smith et al., 2012).
A limited number of studies have reported how sedge phenology responds to N addition. For example, under fertilization, K. capillifolia flowered earlier or later (Bayaerta et al., 2010; Zhang et al., 2013b); K. setchwanensis shortened the flower phase length; K. macrantha showed no response (Han, 2012). The different responses from the similar plant functional types may be due to their distinct growth strategies under distinct environmental conditions.
As the dominant species, each phenological phase of K. pygmaea has exhibited a divergent response to N addition. First, K. pygmaea budded earlier in response to N addition. Most alpine plants rely on a conservative growth strategy (Chapin, 1980; Aerts and Chapin, 2000). But K. pygmaea, as a superior N competitor and shallow-rooted plant, may be able to respond quickly to increased upper-soil N availability and build their growth advantages. For the next phase, the flowering time of K. pygmaea showed no change. However, its flowering phase length was prolonged significantly. Its fruit phase exhibited no change either. In N-rich soils, plant species usually allocate more resources to reproductive growth (Tilman and Wedin, 1991). As a result, the entire reproductive period of K. pygmaea was significantly prolonged in N-added treatment. The longer growth duration of K. .pygmaea in N addition plots could be ascribed to an earlier start of growing season and a non-changed end of the reproductive phase.
The phenological phase change of K. pygmaea, as the dominant species in the alpine meadow ecosystem, has important implications for various aspects of studies, including ecosystem dynamics, plant demography and community interactions, trophic dynamics, and ecosystem energy balance. Furthermore, it has fundamental implications for pastoralists and wildlife in the region. In response to increased N, K. pygmaea becomes increasingly dominant in the community due to its more effective use of the added nutrition. In alpine ecosystems, species composition would change with altered productivity (Grabherr et al., 1994; Pauli et al., 2012). K. pygmaea can become more productive due to its greater dominance under nutrition addition treatment, which in turn would change the community structure and species composition. On the other hand, advanced phenology may make K. pygmaea increasingly susceptible to frost damage, which may reduce its production. A further increased N deposition is possible to have significant effects on plant species composition of the alpine meadows, which can cause impacts on pastoral livelihoods and the delivery of key ecosystem services from the Tibetan Plateau (Yu et al., 2012).
MULTIPLE-LEVEL NITROGEN ADDITION EFFECTS ON PLANT PHENOLOGY
Such plant traits as plant phenology, leaf properties, and biomass of different modules (inflorescences, leaves, stems, and rootstocks) have been conventionally investigated to identify global change-related treatment effects on plant development and growth allocation patterns (Franzaring et al., 2008). The majority of the previous N addition treatment studies have focused on such common indicators as biomass or leaf properties (Gough et al., 2000; Zhang et al., 2007; Niu et al., 2008; Xia and Wan, 2008). Only a few have reported on plant phenology (Lupi et al., 2012; Smith et al., 2012; Xia and Wan, 2013), let alone different levels of N supply effect on phenology (Franzaring et al., 2008).
In the present study, the effects of multiple-level N additions on phenology were not obvious. Similar studies have reported effects of different amounts of nutrition supply on plant development timing of Molinia caerulea from three origins (Franzaring et al., 2008). However, it is not clear to what extent the phenology responded to the multiple-level N additions. For one origin of M. caerulea, leaf senescence decreased with the amount of added nutrition, while it increased for another one. For the third origin, leaf senescence is composed of two opposite responses to the two levels of nutrition addition. The divergent plant responses to a gradient amount of added N may depend on plant compositions (Saarinen, 1998; Green and Galatowitsch, 2002; Xia and Wan, 2008).
Knowledge about environmental change effects caused by different levels of N addition might improve our understanding of the above-mentioned phenomena. In this study, the influence of different amounts of N addition on soil moisture was not obvious. Due to the harsh physical environment, the growing season is relatively short over the Tibetan Plateau (Liu et al., 2006; Ernakovich et al., 2014). At our research site, temperature and moisture conditions are suitable for vegetation growth only in July. So it is also possible that phenological responses may require multiple growing seasons to fully materialize.
Our results suggest that N plays an important role in regulating plant phenology of the dominant species. The addition of N significantly advanced budding time and prolonged the flower phase length and reproductive period duration of dominant species, whereas N addition produced no obvious effects on accompanying species. Different amounts of N addition resulted in no distinct phenological responses of each species. Due to the conservative growth strategy of alpine plants, their phenological responses may require multiple growing seasons to fully materialize. To further elucidate on and predict the phenological responses of organisms to nutritional condition modifications, long-term experiments in various ecosystems are required.
This research is supported by the 973 Program (2013CB956302) of the Ministry of Science and Technology of China, One Hundred Talents Program of Chinese Academy of Sciences, and National Natural Science Foundation of China (31300356, 41201055).
- R. Aerts , and F. S. Chapin , 2000: The mineral nutrition of wild plants revisited: a re-evaluation of processes and patterns. Advances in Ecological Research , 30: 1–67. Google Scholar
- R. Ahas , J. Jaagus , and A. Aasa , 2000: The phenological calendar of Estonia and its correlation with mean air temperature. International Journal of Biometeorology , 44: 159–166. Google Scholar
- I. W. Ashton , A. E. Miller , W. D. Bowman , and K. N. Suding , 2010: Niche complementarity due to plasticity in resource use: plant partitioning of chemical N forms. Ecology , 91: 3252–3260. Google Scholar
- A. Aydelotte , and P. Diggle , 1997: Analysis of developmental preformation in the alpine herb Caltha leptosepala (Ranunculaceae). American Journal of Botany , 84: 1646. Google Scholar
- J. Bai , X. Xu , G. Fu , M. Song , Y. He , and J. Jiang , 2011: Effects of temperature and nitrogen input on nitrogen mineralization in alpine soils on Tibetan Plateau. Agricultural Science and Technology , 12: 1909–1912. Google Scholar
- P. Bayaerta, Jia , X. Yang , and G. Z. Du , 2010: Response of dominating species flowering phenology to fertilization in Qinghai-Tibetan alpine meadow. Acta Prataculturae Sinica , 19: 233–239. Google Scholar
- W. D. Bowman , J. R. Gartner , K. Holland , and M. Wiedermann , 2006: Nitrogen critical loads for alpine vegetation and terrestrial ecosystem response: Are we there yet? Ecological Applications , 16: 1183–1193. Google Scholar
- F. S. Chapin , 1980: The mineral-nutrition of wild plants. Annual Review of Ecology and Systematics , 11: 233–260. Google Scholar
- E. E. Cleland , N. R. Chiariello , S. R. Loarie , H. A. Mooney , and C. B. Field , 2006: Diverse responses of phenology to global changes in a grassland ecosystem. Proceedings of the National Academy of Sciences , 103: 13740–13744. Google Scholar
- E. E. Cleland , I. Chuine , A. Menzel , H. A. Mooney , and M. D. Schwartz , 2007: Shifting plant phenology in response to global change. Trends in Ecology and Evolution , 22: 357–365. Google Scholar
- T. Dorji , O. Totland , S. R. Moe , K. A. Hopping , J. B. Pan , and J. A. Klein , 2013: Plant functional traits mediate reproductive phenology and success in response to experimental warming and snow addition in Tibet. Global Change Biology , 19: 459–472. Google Scholar
- J. A. Dunne , J. Harte , and K. J. Taylor , 2003: Subalpine meadow flowering phenology responses to climate change: integrating experimental and gradient methods. Ecological Monographs , 73: 69–86. Google Scholar
- J. G. Ernakovich , K. A. Hopping , A. B. Berdanier , R. T. Simpson , E. J. Kachergis , H. Steltzer , and M. D. Wallenstein , 2014: Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Global Change Biology , 20: 3256–3269. Google Scholar
- J. Forrest , and A. J. Miller-Rushing , 2010: Toward a synthetic understanding of the role of phenology in ecology and evolution. Philosophical Transactions of the Royal Society B-Biological Sciences , 365: 3101–3112. Google Scholar
- J. Franzaring , I. Holz , and A. Fangmeier , 2008: Different responses of Molinia caerulea plants from three origins to CO2 enrichment and nutrient supply. Acta Oecologica , 33: 176–187. Google Scholar
- L. Gough , C. W. Osenberg , K. L. Gross , and S. L. Collins , 2000: Fertilization effects on species density and primary productivity in herbaceous plant communities. Oikos , 89: 428–439. Google Scholar
- G. Grabherr , M. Gottfried , and H. Pauli , 1994: Climate effects on mountain plants. Nature , 369: 448. Google Scholar
- E. K. Green , and S. M. Galatowitsch , 2002: Effects of Phalaris arundinacea and nitrate-N addition on the establishment of wetland plant communities. Journal of Applied Ecology , 39: 134–144. Google Scholar
- N. Gruber , and J. N. Galloway , 2008: An Earth-system perspective of the global nitrogen cycle. Nature , 451: 293–296. Google Scholar
- S. Gu , D. Hui , and A. Bian , 1998: The contraction-expansion algorithm and its use in fitting nonlinear equations. Journal of Biomathematics , 13: 426–434. Google Scholar
- J. F. Han , 2012: Effects of Simulated Precipitation Changes and Increased N Deposition on the Phenology and Biomass of Plants in an Alpine Meadow. Master's thesis, College of Animal Science and Technology, Sichuan Agricultural University, Yaan, Sichuan, China, 23 pp. Google Scholar
- IPCC, 2007: The AR4 Synthesis Report. Cambridge: Cambridge University Press. Google Scholar
- J. A. Klein , K. A. Hopping , E. T. Yeh , Y. Nyima , R. B. Boone , and K. A. Galvin , 2014: Unexpected climate impacts on the Tibetan Plateau: local and scientific knowledge in findings of delayed summer. Global Environmental Change , 28: 141–152. Google Scholar
- A. Kliber , and C. G. Eckert , 2004: Sequential decline in allocation among flowers within inflorescences: proximate mechanisms and adaptive significance. Ecology , 85: 1675–1687. Google Scholar
- X. Liu , Y. Zhang , W. Han , A. Tang , J. Shen , Z. Cui , P. Vitousek , J. W. Erisman , K. Goulding , P. Christie , A. Fangmeier , and F. Zhang , 2013: Enhanced nitrogen deposition over China. Nature , 494: 459–462. Google Scholar
- X. D. Liu , Z. Y. Yin , X. M. Shao , and N. S. Qin , 2006: Temporal trends and variability of daily maximum and minimum, extreme temperature events, and growing season length over the eastern and central Tibetan Plateau during 1961–2003. Journal of Geophysical Research-Atmospheres , 11: 11–19. Google Scholar
- Y. J. Luo , G. L. Qin , and G. Z. Du , 2006: Importance of assemblage-level thinning: a field experiment in an alpine meadow on the Tibet plateau. Journal of Vegetation Science , 17: 417–424. Google Scholar
- C. Lupi , H. Morin , A. Deslauriers , S. Rossi , and D. Houle , 2012: Increasing nitrogen availability and soil temperature: effects on xylem phenology and anatomy of mature black spruce. Canadian Journal of Forest Research-Revue Canadienne De Recherche Forestiere , 42: 1277–1288. Google Scholar
- X. Lv , S. Kang , L. Zhu , Y. Zhang , and W. Han , 2009: Phenology characters of dominant plants in the Nam Co Basin and its response to climate, Tibet. Journal of Mountain Research , 27: 648–654. Google Scholar
- A. F. Mark , 1970: Floral initiation and development in New Zealand alpine plants. New Zealand Journal of Botany , 8: 67–75. Google Scholar
- J. Martinkova , P. Smilauer , and S. Mihulka , 2002: Phenological pattern of grassland species: relation to the ecological and morphological traits. Flora , 197: 290–302. Google Scholar
- A. E. Miller , and W. D. Bowman , 2003: Alpine plants show species-level differences in the uptake of organic and inorganic nitrogen. Plant and Soil , 250: 283–292. Google Scholar
- A. E. Miller , W. D. Bowman , and K. N. Suding , 2007: Plant uptake of inorganic and organic nitrogen: neighbor identity matters. Ecology , 88: 1832–1840. Google Scholar
- A. J. Miller-Rushing , and D. W. Inouye , 2009: Variation in the impact of climate change on flowering phenology and abundance: an examination of two pairs of closely related wildflower species. American Journal of Botany , 96: 1821–1829. Google Scholar
- K. C. Niu , Y. J. Luo , P. Choler , and G. Z. Du , 2008: The role of biomass allocation strategy in diversity loss due to fertilization. Basic and Applied Ecology , 9: 485–493. Google Scholar
- E. A. Nord , and J. P. Lynch , 2009: Plant phenology: a critical controller of soil resource acquisition. Journal of Experimental Botany , 60: 1927–1937. Google Scholar
- C. Parmesan , 2006: Ecological and evolutionary responses to recent climate change. Annual Review of Ecology Evolution and Systematics , 37: 637–669. Google Scholar
- C. Parmesan , and G. Yohe , 2003: A globally coherent fingerprint of climate change impacts across natural systems. Nature , 421: 37–42. Google Scholar
- H. Pauli , M. Gottfried , S. Dullinger , O. Abdaladze , M. Akhalkatsi , J. L. B. Alonso , G. Coldea , J. Dick , B. Erschbamer , R. F. Calzado , D. Ghosn , J. I. Holten , R. Kanka , G. Kazakis , J. Kollár , P. Larsson , P. Moiseev , D. Moiseev , U. Molau , J. M. Mesa , L. Nagy , G. Pelino , M. Puşcaş , G. Rossi , A. Stanisci , A. O. Syverhuset , J.-P. Theurillat , M. Tomaselli , P. Unterluggauer , L. Villar , P. Vittoz , and G. Grabherr , 2012: Recent plant diversity changes on Europe's mountain summits. Science , 336: 353–355. Google Scholar
- S. L. Piao , P. Friedlingstein , P. Ciais , N. Viovy , and J. Demarty , 2007: Growing season extension and its impact on terrestrial carbon cycle in the northern hemisphere over the past 2 decades. Global Biogeochemical Cycles , 21: 1–11. Google Scholar
- S. Piao , M. Cui , A. Chen , X. Wang , P. Ciais , J. Liu , and Y. Tang , 2011: Altitude and temperature dependence of change in the spring vegetation green-up date from 1982 to 2006 in the Qinghai-Xizang Plateau. Agricultural and Forest Meteorology , 151: 1599–1608. Google Scholar
- E. S. Post , and D. W. Inouye , 2008: Phenology: response, driver, and integrator. Ecology , 89: 319–320. Google Scholar
- S. A. Power , E. R. Green , C. G. Barker , J. N. B. Bell , and M. R. Ashmore , 2006: Ecosystem recovery: heathland response to a reduction in nitrogen deposition. Global Change Biology , 12: 1241–1252. Google Scholar
- M. V. Price , and N. M. Waser , 1998: Effects of experimental warming on plant reproductive phenology in a subalpine meadow. Ecology , 79: 1261–1271. Google Scholar
- B. Rathcke , and E. P. Lacey , 1985: Phenological patterns of terrestrial plants. Annual Review of Ecology and Systematics , 16: 179–214. Google Scholar
- F. J. Richards , 1959: A flexible growth function for empirical use. Journal of Experimental Botany , 10: 290–301. Google Scholar
- T. Saarinen , 1998: Internal C:N balance and biomass partitioning of Carex rostrata grown at three levels of nitrogen supply. Canadian Journal of Botany , 76: 762–768. Google Scholar
- R. A. Sherry , X. Zhou , S. Gu , J. A. Arnone , D. S. Schimel , P. S. Verburg , L. L. Wallace , and Y. Luo , 2007: Divergence of reproductive phenology under climate warming. Proceedings of the National Academy of Sciences , 104: 198–202. Google Scholar
- J. G. Smith , W. Sconiers , M. J. Spasojevic , I. W. Ashton , and K. N. Suding , 2012: Phenological changes in alpine plants in response to increased snowpack, temperature, and nitrogen. Arctic, Antarctic, and Alpine Research , 44: 135–142. Google Scholar
- T. H. Sparks , and A. Menzel , 2002: Observed changes in seasons: an overview. International Journal of Climatology , 22: 1715–1725. Google Scholar
- D. Tilman , and D. Wedin , 1991: Plant traits and resource reduction for five grasses growing on a nitrogen gradient. Ecology , 72: 685–700. Google Scholar
- E. M. Wolkovich , B. I. Cook , J. M. Allen , T. M. Crimmins , J. L. Betancourt , S. E. Travers , S. Pau , J. Regetz , T. J. Davies , N. J. B. Kraft , T. R. Ault , K. Bolmgren , S. J. Mazer , G. J. McCabe , B. J. McGill , C. Parmesan , N. Salamin , M. D. Schwartz , and E. E. Cleland , 2012: Warming experiments underpredict plant phenological responses to climate change. Nature , 485: 494–497. Google Scholar
- J. Y. Xia , and S. Q. Wan , 2008: Global response patterns of terrestrial plant species to nitrogen addition. New Phytologist , 179: 428–439. Google Scholar
- J. Y. Xia , and S. Q. Wan , 2013: Independent effects of warming and nitrogen addition on plant phenology in the Inner Mongolian steppe. Annals of Botany , 111: 1207–1217. Google Scholar
- H. Yang , L. Jiang , L. Li , A. Li , M. Wu , and S. Wan , 2012: Diversity-dependent stability under mowing and nutrient addition: evidence from a 7-year grassland experiment. Ecology Letters , 15: 619–626. Google Scholar
- C. Q. Yu , Y. J. Zhang , H. Claus , R. Zeng , X. Z. Zhang , and J. S. Wang , 2012: Ecological and environmental issues faced by a developing Tibet. Environmental Science and Technology , 46: 1979–1980. Google Scholar
- H. Y. Yu , E. Luedeling , and J. C. Xu , 2010: Winter and spring warming result in delayed spring phenology on the Tibetan Plateau. Proceedings of the National Academy of Sciences , 107: 22151–22156. Google Scholar
- G. Zhang , Y. Zhang , J. Dong , and X. Xiao , 2013a: Green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011. Proceedings of the National Academy of Sciences , 110: 4309–4314. Google Scholar
- L. H. Zhang , C. C. Song , D. X. Wang , and Y. Y. Wang , 2007: Effects of exogenous nitrogen on freshwater marsh plant growth and N2O fluxes in Sanjiang Plain, Northeast China. Atmospheric Environment , 41: 1080–1090. Google Scholar
- Z. L. Zhang , X. D. Liu , P. Jia , and G. Z. Du , 2013b: Effects of nitrogen addition on flowering phenology, functional traits and community structure of alpine meadow in the Eastern Qinghai-Tibet Plateau. Pratacultural Science , 30: 728–735. Google Scholar