Abstract
Changes in spring and autumn phenology and thus growing season length (GSL) pose great challenges in accurately predicting terrestrial primary productivity. However, how spring and autumn phenology in response to land-use change and nitrogen deposition and underlying mechanisms remain unclear. This study was conducted to explore the GSL and its components [i.e., the beginning of growing season and ending of growing season (EGS)] in response to mowing and nitrogen addition in a temperate steppe on the Mongolia Plateau during 2 years with hydrologically contrasting condition [dry (2014) vs. wet (2015)]. Our results demonstrated that mowing advanced the BGS only by 3.83 days, while nitrogen addition advanced and delayed the BGS and EGS by 2.85 and 3.31 days, respectively, and thus prolonged the GSL by 6.16 days across the two growing seasons from 2014 to 2015. When analyzed by each year, nitrogen addition lengthened the GSL in the dry year (2014), whereas it shortened the GSL in the wet year (2015). Further analyses revealed that the contrasting impacts of nitrogen on the GSL were attributed to monthly precipitation regimes and plant growth rate indicated by the maximum of normalized difference vegetation index (NDVmax). Moreover, changes in the GSL and its two components had divergent impacts on community productivity. The findings highlight the critical role of precipitation regimes in regulating the responses of spring and autumn phenology to nutrient enrichment and suggest that the relationships of ecosystem productivity with spring and autumn phenology largely depend on interannual precipitation fluctuations under future increased nitrogen deposition scenarios.
Introduction
Changes in the growing season timing and growing season length (GSL) can regulate biosphere-atmosphere interactions, with consequent carbon (Piao et al., 2007, 2008; Xia et al., 2015) and water cycling (White et al., 1999; Lian et al., 2020; Cheng et al., 2021). The lengthened growing season resulting from advanced beginning of growing season (BGS) and delayed ending of growing season (EGS) has increased primary productivity of terrestrial ecosystems (Piao et al., 2007; Dragoni et al., 2011; Cheng et al., 2021). In addition, carbon loss induced by earlier autumn phenology could counteract the carbon uptake associated with earlier spring phenology, leading to a net carbon loss from terrestrial ecosystems (Piao et al., 2008). Therefore, quantifying the changes in the spring and autumn phenology is critical for the accurate prediction of terrestrial ecosystem carbon balance.
It has been revealed that plant growth could balance from multiple resources (such as nutrient availability) and change allocation to maximize acquisition of the most limiting resources (Bloom et al., 1985; Chapin et al., 1990). This may suggest that elevated anthropogenic nitrogen deposition (Peñuelas et al., 2013; Yu et al., 2019) can alter vegetation activity and thus phenology (Piao et al., 2019, 2020; Vitasse et al., 2021) by providing nitrogen availability. Recent studies have demonstrated that variations in nitrogen or phosphorus availability could change the spring and autumn phenology (Yang et al., 2016; Yin et al., 2016; Fu et al., 2019; Vitasse et al., 2021), as well as the GSL (Wang and Tang, 2019). For example, increased soil nitrogen availability could supplement nutrient deficiencies and thus stimulate plant growth under low temperature in early autumn (McCormack et al., 2014; Delpierre et al., 2016; Yin et al., 2016), which can delay the EGS. In addition, nitrogen addition could also decrease the cellular maturation rate (Kalliokoski et al., 2013; Cuny et al., 2015), and consequently postpone the autumn phenology (Wingler et al., 2006). In nitrogen-limited ecosystems, species usually allocate more resources to growth (LeBauer and Treseder, 2008) under increased nitrogen availability (Xiang et al., 2016), which can delay the reproductive stages (Cleland et al., 2006), thus delaying the senescence date (Wang and Tang, 2019). Nevertheless, whether the effect of nitrogen enrichment on plant phenology regulated by other factors under land-use and climate change scenarios remains largely unknown.
As a widespread land-use practice, mowing has been applied in managed ecosystems (e.g., grasslands; Liu et al., 2017; Zhang et al., 2017) to maintain plant diversity and production. Mowing could affect the microenvironment (such as light, temperature, and moisture) and thus plant growth by removing aboveground litter (Collins et al., 1998; Huhta et al., 2001; Liu et al., 2017). In the early growing season, mowing can elevate soil temperature by increasing light availability and thus accumulated temperature, which is critical for driving plant phenology (Fu et al., 2015; Piao et al., 2019), and consequently advance the BGS. In the late growing season, mowed grasslands may have lower soil water availability due to increased evaporation associated with less litter cover, which can accelerate the EGS (Liu et al., 2016a,b). In addition, it has been proved that mowing may affect reproductive phenology of early flowering species on the Tibetan Plateau (Liu et al., 2017). Therefore, mowing could have the potential to regulate the responses of spring and autumn phenology to nutrient enrichment. However, the direct field experimental evidence and underlying mechanisms remain limited. Given the critical role of soil water availability in mediating plant growth and phenology in grasslands (Körner, 2015; Quan et al., 2019; Zhou et al., 2019, 2022), precipitation also can have the potential regulation on the response of spring and autumn phenology to mowing, nutrient enrichment, and their interactions.
Considering the above knowledge gaps, this study was conducted to explore the effects of nitrogen addition and mowing (annually) on the GSL and its two components (i.e., BGS and EGS) over two contrasting hydrologically growing seasons from 2014 (dry) to 2015 (wet) in a temperate grassland on the Mongolian Plateau. The specific questions we addressed in this study included: (1) How do nitrogen addition and mowing affect the GSL and its two components? (2) How do environmental factors mediate the responses of GSL and its two components to nitrogen addition and mowing?
Materials and Methods
Study Site
This study was located in a semiarid steppe in Duolun Restoration Ecology Research Station, Duolun County (42°02' N, 116°07' E, 1324 m a.s.l.), Inner Mongolia, China. Long-term (1961–2018) mean annual temperature and precipitation were 2.1°C and 385.5 mm, respectively (China Meteorological Data Sharing Service System). The long-term mean annual potential evaporation at this experimental site is 1,748 mm (Wang et al., 2021). The sandy soil is classified as Haplic Calcisol (Food and Agriculture Organization of the United Nations), with 62.75 ± 0.04% sand, 20.30 ± 0.01% silt, and 16.95 ± 0.01% clay. The study site has been fenced since 2001 to exclude cattle and sheep grazing. The dominant plant species in this temperate steppe were Stipa krylovii, Artemisia frigida, and Agropyron cristatum (Miao et al., 2020; Liu et al., 2021b).
Experimental Design
A factorial design with four treatments (control, mowing, nitrogen addition, and mowing plus nitrogen addition) and five replications for each treatment was employed in this experiment (Liu et al., 2018; Wang et al., 2020, 2021). There were 20 plots (each was 4 × 4 m2) arranged by four rows and five columns, and the buffer zone was 1 m between any two plots. To reduce the edge effect, a 0.5 m buffer zone to the edge of each plot was also designed, thus all the measurements were conducted in the 3 × 3 m2 core zone. The mowing treatment was conducted in the late August of each year since 2012. All the plants were mowed 5 cm aboveground to simulated hay harvesting, a widely land-use type in many grasslands (Luo et al., 2001; Niu et al., 2010; Liu et al., 2017; Du et al., 2018). The nitrogen treatment (10 g m−2 year−1) followed the range of airborne nutrient deposition observed in Northern China (Liu et al., 2021a). Nitrogen (N) was added in each nitrogen addition plot using the form of NH4NO3 at a rate of 100 kg N ha−1 in early May of each year since 2013.
Soil Temperature Measurement
To investigate the effects of temperature on spring and autumn phenology, we measured soil temperature at a depth of 10 cm that was recorded every 2 h using DS 1923 iButton (Maxim Integrated, San Jose, CA, USA) during the growing season (from May to October in each year). Due to not being waterproof, the sensors were sealed with balloons, which had been demonstrated to be an effective way in other systems to avoid direct exposure to precipitation (Lutterschmidt et al., 2006; Kearney et al., 2011). Then, the wrapped sensors were put into soil at a depth of 10 cm.
Normalized Difference Vegetation Index Measurement and Data Fitting
We calculated normalized difference vegetation index (NDVI) to understand the community development under different treatments and also obtain the spring (BGS) and autumn phenology (EGS), and thus, the GSL (the difference in day of year of the BGS and EGS). Spectral reflectance of a 1 m × 1 m permanent subplot in each plot was measured at 5- to 7-day intervals with a Tetracam Agricultural Digital Camera (ADC, Tetracam Inc., Chatsworth, CA, USA) at cloud-free noon during each growing season from May to October. The camera was held 1 m above each subplot with an iron shelf during each measurement. NDVI was calculated as Then, we fitted NDVI data with a 5-parameter Weibull function using Sigmaplot 14.0 software (Systat Software, Inc., San Jose, CA, USA). BGS, EGS, GSL, as well as the maximum value of NDVI, were extracted using the methods described by Xia et al. (2015).
Plant Sampling
To examine whether the changes in spring and autumn phenology affect the community productivity, we measured the aboveground net primary productivity (ANPP). One 0.5 m × 0.5 m quadrat was established in each experimental plot in each of the late growing season (May–October) of 2014 and 2015. Then, we clipped living aboveground biomass in each quadrat and separated into species level. Living aboveground biomass of each species was oven-dried at 65°C for 48 h and weighted to determine the dry mass (Guo et al., 2021).
Data Analysis
First, we calculated the monthly mean value of soil temperature and then averaged them from May to October as growing season mean values. Then, two-way ANOVAs were used to examine the effects of mowing and nitrogen addition and their interactions on the soil temperature, BGS, EGS, GSL, and the maximum of NDVImax, as well as the ANPP. The growing season mean values were used to calculate mowing and nitrogen effects on the above variables. Mowing effects were calculated as [(mowing – control)/control] in the unfertilized plots and [(mowing plus nitrogen addition – nitrogen addition)/nitrogen addition] in the fertilized plots. Nitrogen effects were calculated as [(nitrogen addition – control)/control] in the unmowed plots and [(mowing plus nitrogen addition – mowing)/mowing] in the mowed plots. Mowing and nitrogen effects were calculated in each year. The correlations among variables were explored by Pearson's correlation method. Significant differences were evaluated at the 0.05 probability level. Linear regressions were used to explore the relationships of the BGS, EGS, and GSL with soil temperature, precipitation, and ANPP. All analyses were conducted using SAS 8.0 (SAS Institute Inc., Cary, NC, USA). GraphPad Prism 9.0 (GraphPad Inc., San Diego, CA, USA) was used to plot the graphs.
Results
Precipitation Patterns and Variations of Soil Temperature Under Different Treatments
Growing season (May–October) precipitation in 2014 (305 mm) was 12% lower than the long-term mean (346 mm), whereas the precipitation of growing season in 2015 (359 mm) was 3% above the long-term mean. In addition, monthly precipitation fluctuated greater in 2015 than those in 2014 (Figures 1A,B), especially in August. Similarly, precipitation from June to September in 2014 (236 mm) was 21% lower than the long-term mean (299 mm), whereas precipitation from June to September in 2015 (307 mm) was 3% higher than the long-term mean (Figures 1A,B).
Figure 1
No interannual variation of soil temperature was observed over 2 years from 2014 to 2015. Mowing elevated soil temperature by 0.68°C (p = 0.04), whereas nitrogen addition had no impact on it over 2 years (Table 1). The effects of mowing on soil temperature did not change with year. There were no interactions between mowing and nitrogen addition on soil temperature (p > 0.05, Table 1). When analyzed by each year, both mowing and nitrogen addition had no impact on soil temperature in 2014 or 2015. When analyzed by different stages of growing season, mowing marginally increased soil temperature of the early growing season (May–June) by 1.31°C and 1.33°C (both p = 0.08, Supplementary Table S1) in 2014 and 2015, respectively, whereas it had no effects on that of the middle or late growing season in any of 2 years. When analyzed by each month during each growing season, mowing marginally elevated soil temperature only in May by 1.42°C (p = 0.08) and 1.38°C (Figures 1C,D, p = 0.07, Supplementary Table S2) in 2014 and 2015, respectively. Nitrogen addition had no impact on soil temperature in any of the months in 2 years (Supplementary Table S2).
Table 1
| Source of variation | Soil T | BGS | EGS | GSL | NDVImax |
|---|---|---|---|---|---|
| Year | 0.822 | <0.001 | <0.001 | <0.001 | <0.001 |
| Mowing | 0.043 | 0.036 | 0.273 | 0.387 | 0.928 |
| Nitrogen | 0.105 | 0.051 | <0.001 | <0.001 | <0.001 |
| Year * Mowing | 0.975 | 0.208 | 0.370 | 0.200 | 0.013 |
| Year * Nitrogen | 0.971 | <0.001 | <0.001 | <0.001 | <0.001 |
| Mowing * Nitrogen | 0.416 | <0.001 | <0.001 | <0.001 | 0.529 |
| Year * Mowing * Nitrogen | 0.988 | <0.001 | <0.001 | <0.001 | 0.132 |
Results (p-value) of three-way ANOVA on the effects of year, mowing, and nitrogen addition and their interactions on the soil temperature (Soil T), beginning (BGS), ending (EGS), and length of growing season (GSL), as well as the maximum of normalized difference vegetation index (NDVImax).
The bold numbers indicate the significance at p <0.05.
Effects of Mowing and Nitrogen Addition on the BGS, EGS, and GSL, as Well as the Maximum of NDVI
Significant interannual variations of BGS, EGS, and GSL, and the maximum of NDVI were found (all p < 0.001, Table 1). BGS and EGS of 2014 were 4.74 and 29.87 days earlier than those of 2015, respectively, leading to a shorter growing season in 2014 than that of 2015 (Figure 2A). In addition, the maximum of NDVI in 2014 (0.22 ± 0.00) was also lower than that of 2015 (0.52 ± 0.01, Figure 2B). Over 2 years, mowing advanced the BGS by 2.43 days (p = 0.04, Table 2), whereas it had no effects on the EGS or GSL. Nitrogen addition advanced the BGS and delayed the EGS by 2.26 days (p = 0.05) and 4.07 days (p < 0.001), and thus lengthened the growing season by 6.33 days (p < 0.001, Figure 2A). Mowing had no impact on the maximum of NDVI, whereas nitrogen addition increased it by 0.06 (p < 0.001, Table 1) over 2 years from 2014 to 2015. Moreover, the effects of nitrogen addition on the BGS, EGS, and GSL, as well as the maximum of NDVI significantly changed with year (all p < 0.001).
Figure 2
Table 2
| Year | Source of variation | BGS | EGS | GSL | NDVImax |
|---|---|---|---|---|---|
| 2014 | Mowing | 0.055 | 0.916 | 0.249 | 0.009 |
| Nitrogen | <0.001 | <0.001 | <0.001 | 0.045 | |
| Mowing*Nitrogen | <0.001 | <0.001 | <0.001 | 0.337 | |
| 2015 | Mowing | 0.421 | 0.003 | 0.602 | 0.171 |
| Nitrogen | <0.001 | <0.001 | <0.001 | <0.001 | |
| Mowing*Nitrogen | 0.505 | 0.396 | 0.371 | 0.237 |
Results (p-value) of two-way ANOVA on the effects of mowing, nitrogen addition, and their interactions on the BGS, EGS, GSL, and NDVImax in 2014 and 2015.
Refer to abbreviations in Table 1.
The bold numbers indicate the significance at p <0.05.
When analyzed by each year, mowing marginally advanced the BGS by 3.85 days (p = 0.06, Table 2), whereas it had no impact on the EGS or GSL, whereas nitrogen addition advanced the BGS and delayed the EGS by 16.54 and 10.20 days, respectively, and thus extended the growing season by 26.73 days in 2014 (Figure 2A, all p < 0.001, Table 2). Mowing only advanced the EGS by 1.73 days in 2015 (p < 0.01). In contrast to 2014, nitrogen addition delayed the BGS and advanced the EGS by 12.02 and 2.05 days, respectively, leading to a shortened growing season (Figure 2A, all p < 0.001, Table 2). Mowing decreased, whereas nitrogen addition elevated the maximum of NDVI by 0.02 (p = 0.04) and 0.01 (p < 0.01, Table 2), respectively, in 2014. Mowing had no impact on the maximum of NDVI (p > 0.05), while nitrogen addition increased it by 0.11 (p < 0.001) in 2015.
Relationships of the BGS, EGS, and GSL, as Well as the Maximum of NDVI With Soil Temperature and Precipitation in Different Stages of Growing Season
Our results showed that there were no relationships of the beginning, ending, and length of growing season, and the maximum of normalized difference index with soil temperature in any stages of growing season across the 2 years (Figure 3). When analyzed by each month, the results revealed positive relationships of EGS and the maximum of NDVI with soil temperature in August (Supplementary Figure S1). In contrast, negative dependences of ending and length of growing season, as well as the maximum of NDVI on soil temperature in October were observed (Supplementary Figure S1).
Figure 3
In addition, EGS showed negative dependence on early growing season precipitation, whereas positive dependence on middle and late growing season precipitation. Similar relationships of GSL and the maximum of NDVI with middle and late growing season precipitation were also found, respectively, over the 2 years from 2014 to 2015 (Figure 3). When analyzed by each month, there were no relationships of BGS with precipitation of any month (Supplementary Figure S2). EGS showed negative relationships with precipitation in May, June, and October, whereas positive relationships with precipitation in July, August, and September (Supplementary Figure S2). Similar patterns were also found in GSL and the maximum of NDVI.
Impacts of Changes in the BGS, EGS, and GSL, as Well as the Maximum of NDVI on Aboveground Net Primary Productivity
Positive dependences of ANPP on BGS (R2 = 0.10, p = 0.01), EGS (R2 = 0.65, P < 0.001), GSL (R2 = 0.24, p < 0.01), and the maximum of NDVI (R2 = 0.79, p < 0.001) were found over 2 years from 2014 to 2015 (Figures 4A–D). When analyzed by each year, there were no dependences of ANPP on the BGS, EGS, or GSL, or the maximum of NDVI in 2014 (Supplementary Figures S4a–d). In contrast, positive relationships of ANPP with the BGS (R2 = 0.41, p < 0.01) and the maximum of NDVI (R2 = 0.27, p = 0.02) were found in 2015 (Supplementary Figures S5a,d). Surprisingly, ANPP showed negative dependence on the GSL in 2015 (R2 = 0.31, p = 0.01, Supplementary Figure S5c).
Figure 4
Discussion
Effects of Mowing on GSL and Its Two Components
It has been well-documented that temperature plays a critical role in regulating spring and autumn phenology, and thus GSL in terrestrial ecosystems (Cleland et al., 2007; Peñuelas et al., 2009; Chuine et al., 2010). For example, spring temperature drives the onset of spring phenology in the northern hemisphere (Piao et al., 2007, 2015). In addition, higher temperatures could delay autumn phenology in temperate China (Piao et al., 2007; Liu et al., 2016a). As a result, increased temperature could lengthen the growing season by advancing the spring phenology and delaying the autumn phenology. In this study, although mowing increased soil temperature and advanced the BGS over the two growing seasons, no robust relationships of BGS with temperature were found (Figure 3, Supplementary Figure S1). In fact, recent studies have shown that the control role of temperature in mediating spring phenology is declining (Fu et al., 2015), which can partly support the weak relationships of GSL and its two components with soil temperature in this study. This observation also indicates that other factors (except for temperature) can mediate spring phenology under mowing.
Photoperiod has been demonstrated to affect plant phenology in terrestrial ecosystems (Körner and Basler, 2010; Flynn and Wolkovich, 2018). Longer photoperiod can advance the spring phenology of tree species (Chuine et al., 2010). In this study, by removing the standing litter, mowing could also increase the light availability for the short species (especially for Potentilla acaulis L., an early species in the study site) in the early growing season, during which the short species can have the potential to accelerate growth and thus advance the BGS across 2 years. In addition, leaf unfolding could be also regulated by the height-time hypothesis that short individuals can advance their leaf unfolding date (Sun and Frelich, 2011; Liu et al., 2021a). In this study, mowing indeed decreased plant height in the early growing season (unpublished data), which can consequently advance the spring phenology. Nevertheless, when analyzed by each year, the advancement effects of mowing were not significant, especially in 2015, which may be attributed to the insufficient soil water availability due to increased light availability and thus evaporation under mowing. Given the important role of soil water availability in regulating spring phenology in temperate steppes (Shen et al., 2015; Luo et al., 2021), the deficient available water in the soil may weaken the stimulated temperature and light availability associated with mowing in this study. In addition, the findings of mowing did not affect the EGS, suggesting that autumn phenology might not be sensitive to changes in microclimate factors associated with mowing in temperate grasslands. Moreover, given that changes in photoperiod could affect the response of plant phenology to temperature (Basler and Körner, 2014; Way and Montgomery, 2015), the increased light availability under mowing in the early growing season may also interact with temperature to affect spring phenology in the temperate grassland. Further manipulative evidence is still needed to support the assertions.
Effects of Nitrogen Addition on GSL and Its Two Components Depend on Precipitation
Nitrogen is one of the most important factors affecting plant growth and phenology (Jing et al., 2017; Liu et al., 2018; Vitasse et al., 2021), especially in nitrogen-limited ecosystems (such as grasslands). Nitrogen addition could stimulate plant growth and thus phenology (Fu et al., 2019; Wang and Tang, 2019). A meta-analysis has shown that nitrogen addition advances leaf senescence across all the biomes, including forest, grassland, cropland, and desert (Wang and Tang, 2019). Surprisingly, the effects of nitrogen addition on the GSL and its two components were opposite in the two hydrologically contrasting growing seasons in this study. In the dry growing season (2014), nitrogen addition advanced the BGS and delayed the EGS, leading to a lengthened growing season (Figure 2A). In contrast, nitrogen addition delayed the BGS and advanced the EGS, resulting in a shortened growing season (Figure 2A) in the wet growing season (2015). The contrasting findings are similar to those found in a previous study, which demonstrates the opposite phenological response in dry vs. wet years (Bao et al., 2021). Normalized difference vegetation index can be used an index of photosynthesis in ecosystems with low leaf area index or vegetation cover (Wohlfahrt et al., 2010; Del Grosso et al., 2018). In the dry growing season, plants had a lower growth rate indexed as the maximum of NDVI compared with those in the wet growing season (Figure 2B). The low growth rate of plants could reduce nitrogen use, which can result in retaining more nitrogen in the soil and subsequently used by plants when small precipitation events occur. In addition, precipitation in the early growing season is critical for driving the spring phenology (Shen et al., 2015; Ganjurjav et al., 2020; Wang et al., 2022). In this study, precipitation in May of 2014 was 40% greater than that of long-term mean (Supplementary Figure S3), which provides water requirement for plants to begin to grow in the early growing season (Shen et al., 2015; Luo et al., 2021). Along with the favorable precipitation condition, addition of nitrogen could further accelerate plant growth and consequently advance the BGS in 2014. In addition, relatively low values of the maximum of NDVI associated with lack of precipitation in the middle of the growing season could decrease the nitrogen use. Lower precipitation could also be favorable for retaining nitrogen into soil without being leached, which allows plants to use nitrogen more efficiently when precipitation events occur. In fact, small but steady precipitation indeed occurs in the late growing season of 2014, which stimulates the nitrogen effects on maintaining plant growth in the late growing season and thus delays the EGS, consequently lengthening the growing season in 2014. The relationships of the EGS and GSL with monthly precipitation (Figure 3, Supplementary Figure S2) support the above arguments.
In contrast, most of the precipitation in May occurred in the early stage of 2015 (Figure 1B, Supplementary Figure S3) and could not supply steady and consecutive water conditions in the middle and late May for plants to begin to grow. Therefore, the less precipitation in the early growing season might weaken the nitrogen effects on plant growth, which led to the delayed BGS in 2015 (Shen et al., 2015; Ganjurjav et al., 2020; Wang et al., 2022). In addition, although precipitation of 2015 was greater than that of the long-term mean, most of the precipitation occurred in June, July, and September (Figure 1B). On the one hand, high precipitation in the middle growing season (except for precipitation in August) significantly stimulated plant growth (indicated by the maximum of NDVI, Figure 2B) and thus for nitrogen use. On the other hand, high precipitation in the middle growing season could also increase nitrogen leaching (Brandt et al., 2010), which can lead to decreased soil-available nitrogen in the late growing season. Moreover, precipitation in the late growing season, especially in October, was low compared with that in 2014 (Figures 2A,B). The low soil-available nitrogen combined with low precipitation was not favorable for plant growth and resulted in accelerated leaf senescence in the late growing season (Liu et al., 2016a,b; Ren and Peichl, 2021). Thus, the delayed BGS and advanced EGS under nitrogen addition shortened the growing season in the wet growing season in 2015. Our findings indicate that precipitation plays a considerable role in regulating the nutrient effects on spring and autumn phenology in temperate grasslands. Given the diverse driving factors for plant phenology (Piao et al., 2019), the effects of other factors (except for precipitation) on spring and autumn phenology are still needed to be investigated in the future (Zhou et al., 2022).
Interactive Effects of Mowing and Nitrogen Addition on GSL and Its Two Components
To the best of our knowledge, our observations of interaction between mowing and nitrogen addition on spring and autumn phenology, as well as GSL, provide the first experimental evidence on the phenological responses under mowing and nitrogen enrichment in temperate steppes. However, the interactions between mowing and nitrogen addition were different in the two hydrologically contrasting years. The interactive effects of mowing and nitrogen addition on GSL and its two components observed in this study could be largely ascribed to those in 2014 (Tables 1, 2). We found that nitrogen addition had no impact on GSL and its two components without mowing, whereas advanced and delayed spring and autumn phenology with mowing, respectively, and thus extended GSL in 2014. Given that grassland in this study site is nitrogen limitation, nitrogen addition can stimulate plant growth and standing litter accumulation (Liu et al., 2018), which may decrease light availability and temperature accumulation and thus have negative effects on leaf unfolding of short species in the next growing season (Piao et al., 2015; Beil et al., 2021). Mowing can weaken the above negative effects by removing standing litter and thus providing light availability and temperature requirements for driving plant phenology (Flynn and Wolkovich, 2018; Piao et al., 2019), and consequently have substantial impacts on spring and autumn phenology, as well as GSL. Because of low precipitation in 2014, the accumulation of standing litter is less than that in 2013, which cannot be enough to cause the light and temperature limitation in 2015. As a consequence, no interactive effects of mowing and nitrogen addition on GSL and its two components were observed in 2015. These findings suggest that interactions between mowing and nutrient addition could also be mediated by precipitation regimes in different years. Given the changing precipitation regimes, including amount, frequency, intensity, and temporal distributions (IPCC, 2018), it is needed to conduct multifactor manipulative experiments to better understand the realistic response of vegetative phenology to global change.
Implications for Community Productivity in the Temperate Steppe
GSL and its components have critical roles in affecting ecosystem productivity and its interannual variation (Piao et al., 2007, 2008; Xia et al., 2015); however, their roles remain unclear in the grassland ecosystems, especially under land-use change and nitrogen deposition scenarios. In this study, ANPP in 2015 is greater than that in 2014; precipitation could be one of the limiting factors for the lower ANPP in 2014. In addition to precipitation, our observations of the relationships of ANPP with BGS, EGS, and GSL, as well as the maximum of NDVI (Figure 4), suggest that GSL and growth rate are also considerable factors in mediating ANPP in the temperate steppe. Nevertheless, the negative dependence of ANPP with GSL is not consistent with that reported by previous studies (White et al., 1999; Wu et al., 2012; Michaletz et al., 2014), which have shown that longer GSL generally increase ecosystem productivity. In fact, the negative relationship between ANPP and GSL can be affected by the nitrogen addition. The shortened growing season resulting from delaying the BGS and advancing the EGS combined by increasing maximum of NDVI in 2015 under nitrogen addition could explain the above findings. The positive effects of increasing the maximum of NDVI on ANPP (Supplementary Figure S5) under nitrogen addition can offset the negative effects of shortened growing season and consequently enhance ANPP in 2015. This indicates that the maximum of growth rate could determine the ecosystem productivity irrespective of shortened GSL. However, we did not observe the relationships of ANPP with GSL and its components in 2014, which indicate that low precipitation amount may affect the relationships among them. The findings highlight the jointly control roles of GSL and growth rate in mediating ecosystem productivity in temperate steppes under land-use change and increased nitrogen deposition scenarios.
Conclusion
Using a 2-year field manipulative experiment with mowing and nitrogen addition, we demonstrated that nitrogen addition showed divergent impacts on GSL and its two components (i.e., BGS and EGS), that is, nitrogen addition advanced the BGS and delayed the EGS in the dry year, whereas it delayed the BGS and advanced the EGS in the wet year. In addition, the effect of nitrogen addition on the maximum of NDVI was larger in the wet year than that in the dry year, indicating that the impacts of nitrogen on vegetation activity were dependent on precipitation regimes. The changes in the GSL and the maximum of NDVI had also diverse impacts on ANPP in years with different precipitation, suggesting that precipitation can enhance the dependences of ecosystem productivity on spring and autumn phenology under nitrogen deposition scenarios in the temperate steppe. The findings promote our understanding on the effects of land-use change and nitrogen deposition on vegetation activity and productivity among years with different precipitation regimes.
Funding
This study was financially supported by the National Natural Science Foundation of China (42107225, 31670477, 31200375) and the Postdoctoral Innovation and Practice Base of Anyang Institute of Technology (BSJ2020021, BHJ2021007).
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.
Statements
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/s.
Author contributions
ZZ, YL, and YM originally formulated the idea. KZ, WW, and JZ developed methodology. SC and HZ conducted fieldwork. ZZ, LZ, and YL generated data analyses, and wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors thank Juan Xuan and Jingwei Guo for their help with field work and sampling, and Xiang Zhao and Meiguang Jiang for their help managing the experiment.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.861794/full#supplementary-material
References
1
BaoF.XinZ. M.LiJ. Z.LiuM. H.CaoY. L.LuQ.et al. (2021). Effects of simulated enhancement of precipitation on the phenology of Nitraria tangutorum under extremely dry and wet years. Plants10, 1474. 10.3390/plants10071474
2
BaslerD.KörnerC. (2014). Photoperiod and temperature responses of bud swelling and bud burst in four temperate forest tree species. Tree Physiol.34, 377–388. 10.1093/treephys/tpu021
3
BeilI.KreylingJ.MeyerC.LemckeN.MalyshevA. V. (2021). Late to bed, late to rise–Warmer autumn temperatures delay spring phenology by delaying dormancy. Glob. Change Biol.27, 5806–5817. 10.1111/gcb.15858
4
BloomA. J.ChapinF. S.MooneyH. A. (1985). Resource limitation in plants-An economic analogy. Annu. Rev. Ecol. Syst.16, 363–392. 10.1146/annurev.es.16.110185.002051
5
BrandtL. A.KingJ. Y.HobbieS. E.MilchunasD. G.SinsabaughR. L. (2010). The role of photodegradation in surface litter decomposition across a grassland ecosystem precipitation gradient. Ecosystems13, 765–781. 10.1146/annurev.es.21.110190.002231
6
ChapinF. S.SchulzeE. D.MooneyH. A. (1990). The ecology and economics of storage in plants. Annu. Rev. Ecol. Syst.21, 423–447
7
ChengM.JinJ. X.JiangH. (2021). Strong impacts of autumn phenology on grassland ecosystem water use efficiency on the Tibetan Plateau. Sci. Total Environ.126, 107682. 10.1016/j.ecolind.2021.107682
8
ChuineI.MorinX.BugmannH. (2010). Warming, photoperiods, and tree phenology. Science329, 277–278. 10.1126/science.329.5989.277-e
9
ClelandE. E.ChiarielloN. R.LoarieS. R.MooneyH. A.FieldC. B. (2006). Diverse response of phenology to global changes in a grassland ecosystem. Proc. Natl. Acad. Sci U. S. A.103, 13740–13744. 10.1073/pnas.0600815103
10
ClelandE. E.ChuineI.MenzelA.MooneyH. A.SchwartzM. D. (2007). Shifting plant phenology in response to global change. Trends Ecol. Evol.22, 357–365. 10.1016/j.tree.2007.04.003
11
CollinsS. L.KnappA. K.BriggsJ. M.BlairJ. M.SteinauerE. M. (1998). Modulation of diversity by grazing and mowing in native tall grass prairie. Science280, 745–747. 10.1126/science.280.5364.745
12
CunyH. E.RathgeberC. B. K.FrankD.FontiP.MäkinenH.PrislanP.et al. (2015). Woody biomass production lags stem-grith increase by over one month in coniferous forests. Nat. Plants1, 15160. 10.1038/nplants.2015.160
13
Del GrossoS. J.PartonW. J.DernerJ. D.ChenM. S.TuckerC. J. (2018). Simple models to predict grassland ecosystem C exchange and actual evapotranspiration using NDVI and environmental variables. Agric. Forest Meteorol.249, 1–10. 10.1016/j.agrformet.2017.11.007
14
DelpierreN.VitasseY.ChuineI.GuillemotJ.BazotS.RutishauserT.et al. (2016). Temperate and boreal forest tree phenology: from organ-scale processes to terrestrial ecosystem models. Ann. For. Sci.73, 5–25. 10.1007/s13595-015-0477-6
15
DragoniD.SchmidH. P.WaysonC. A.PotterH.GrimmondC. S. B.RandolphJ. C. (2011). Evidence of increased net ecosystem productivity associated with a longer vegetated season in a deciduous forest in south-central Indiana, USA. Glob. Change Biol.17, 886–897. 10.1111/j.1365-2486.2010.02281.x
16
DuY.HanH. Y.WangY. F.ZhongM. X.HuiD. F.NiuS. L.et al. (2018). Plant functional groups regulate soil respiration responses to nitrogen and mowing over a decade. Funct. Ecol.32, 1117–1127. 10.1111/1365-2435.13045
17
FlynnD. F. B.WolkovichE. M. (2018). Temperature and photoperiod drive spring phenology across all species in a temperate forest community. N. Phytol.219, 1353–1362. 10.1111/nph.15232
18
FuY. S.ZhaoH. F.PiaoS. L.PeaucelleM.PengS. S.ZhouG. Y.et al. (2015). Declining global warming effects on the phenology of spring leaf unfolding. Nature526, 104–107. 10.1038/nature15402
19
FuY. S. H.PiaoS. L.DelpierreN.HaoF. H.HänninenH.GengX. J.et al. (2019). Nutrient availability alters the correlation between spring leaf-out and autumn leaf senescence dates. Tree Physiol.39, 1277–1284. 10.1093/treephys/tpz041
20
GanjurjavH.GornishE. S.HuG. Z.SchwartzM. W.WanY. F.LiY.et al. (2020). Warming and precipitation addition interact to affect plant spring phenology in alpine meadows on the central Qinghai-Tibetan Plateau. Agricult. Forest Meteorol.287, 107943. 10.1016/j.agrformet.2020.107943
21
GuoJ. W.ZhaoC. C.ZhangL. N.HanY. Y.CaoR.LiuY. Z.et al. (2021). Water table decline alters arthropod community structure by shifting plant communities and leaf nutrients in a Tibetan peatland. Sci. Total Environ.151944. 10.1016/j.scitotenv.2021.151944
22
HuhtaA.-P.RautioP.TuomiJ.LaineK. (2001). Restorative mowing on an abandoned semi-natural meadow: short-term and predicted long-term effects. J. Veget. Sci.12, 677–686. 10.2307/3236908
23
IPCC. (2018). Global warming of 1.5°C: An IPCC Special Report on the impacts of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty.
24
JingH.ZhouH. X.WangG. L.XueS.LiuG. B.DuanM. C. (2017). Nitrogen addition changes the stoichiometry and growth rate of different organs in Pinus tabuliformis seedlings. Front. Plant Sci.8, 1922. 10.3389/fpls.2017.01922
25
KalliokoskiT.MäkinenH.JyskeT.NöjdP.LinderS. (2013). Effects of nutrient optimization on intra-annual wood formation in Norway spruce. Tree Physiol.33, 1145–1155. 10.1093/treephys/tpt078
26
KearneyM.FergusonE.FumeiS.GallacherA.MitchellP.WoodfordR.et al. (2011). A cost-effective method of assessing thermal habitat quality for endotherms. Austral Ecol.36, 297–302. 10.1111/j.1442-9993.2010.02150.x
27
KörnerC.. (2015). Paradigm shift in plant growth control. Curr. Opin. Plant Biol.25, 107–114. 10.1016/j.pbi.2015.05.003
28
KörnerC.BaslerD. (2010). Phenology under global warming. Science327, 1461–1462. 10.1126/science.1186473
29
LeBauerD.TresederK. (2008). Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology89, 371–379. 10.1890/06-2057.1
30
LianX.PiaoS. L.LiL. Z.HuntingfordC.CiaisP.CescattiA.et al. (2020). Summer soil drying exacerbated by earlier spring greening of northern vegetation. Sci. Adv. 6, eaax0255. 10.1126/sciadv.aax0255
31
LiuQ.FuY. S.ZengZ. Z.HuangM. T.LiX. R.PiaoS. L. (2016a). Temperature, precipitation, and insolation effects on autumn vegetation phenology in temperate China. Glob. Change Biol.22, 644–655. 10.1111/gcb.13081
32
LiuQ.FuY. S.ZhuZ. C.LiuY. W.LiuZ.HuangM. T.et al. (2016b). Delayed autumn phenology in the Northern Hemisphere is related to change in both climate and spring phenology. Glob. Change Biol.22, 3702–3711. 10.1111/gcb.13311
33
LiuY. Z.LiG. Y.WuX. W.NiklasK. J.SunS. C. (2021a). Linkage between species traits and plant phenology in an alpine meadow. Oecologia195, 409–419. 10.1007/s00442-020-04846-y
34
LiuY. Z.MaG. G.ZanZ. M.ChenA. Q.MiaoY.WangD.et al. (2018). Effects of nitrogen addition and mowing on rodent damage in an Inner Mongolian steppe. Ecol. Evol.8, 3919–3926. 10.1002/ece3.3949
35
LiuY. Z.MiaoR. H.ChenA. Q.MiaoY.LiuY. J.WuX. W. (2017). Effects of nitrogen addition and mowing on reproductive phenology of three early-flowering forb species in a Tibetan alpine meadow. Ecol. Eng.99, 119–125. 10.1016/j.ecoleng.2016.11.033
36
LiuY. Z.ZhaoC. C.GuoJ. W.ZhangL. N.XuanJ.ChenA. Q.et al. (2021b). Short-term phosphorus addition augments the effects of nitrogen addition on soil respiration in a typical steppe. Sci. Total Environ.761, 143211. 10.1016/j.scitotenv.2020.143211
37
LuoM.MengF. H.SaC. L.DuanY. C.BaoY. H.LiuT.et al. (2021). Response of vegetation phenology to soil moisture dynamics in the Mongolian Plateau. Catena206, 105505. 10.1016/j.catena.2021.105505
38
LuoY. Q.WanS. Q.HuiD. F.WallaceL. L. (2001). Acclimatization of soil respiration to warming in a tall grass prairie. Nature413, 622–625. 10.1038/35098065
39
LutterschmidtD. I.LeMasterM. P.MasonR. T. (2006). Minimal overwintering temperatures of red-sided garter snakes (Thamnophis sirtalis parietalis): a possible cue for emergence?Can. J. Zool.84, 771–777. 10.1139/z06-043
40
McCormackM. L.AdamsT. S.SmithwickE. A.EissenstatD. M. (2014). Variability in root production, phenology, and turnover rate among 12 temperate tree species. Ecology95, 2224–2235. 10.1890/13-1942.1
41
MiaoY.LiuM. Z.XuanJ.XuW.WangS. L.MiaoR. H.et al. (2020). Effects of warming on soil respiration during the non-growing seasons in a semiarid temperate steppe. J. Plant Ecol.13, 288–294. 10.1093/jpe/rtaa013
42
MichaletzS. T.ChengD. L.KerkhoffA. J.EnquistB. J. (2014). Convergence of terrestrial plant production across global climate gradients. Nature512, 39–43. 10.1038/nature13470
43
NiuS. L.WuM. Y.HanY.XiaJ. Y.ZhangZ.YangH. J.et al. (2010). Nitrogen effects on net ecosystem exchange in a temperate steppe. Glob. Change Biol.16, 144–155. 10.1111/j.1365-2486.2009.01894.x
44
PeñuelasJ.PoulterB.SardansJ.CiaisP.van der VeldeM.BoppL.et al. (2013). Human-induced nitrogen-phosphorus imbalances alter natural and managed ecosystems across the globe. Nat. Commun.4, 2934. 10.1038/ncomms3934
45
PeñuelasJ.RutishauserT.FilellaI. (2009). Phenology feedbacks on climate change. Science324, 887–888. 10.1126/science.1173004
46
PiaoS. L.CiaisP.FriedlingsteinP.PeylinP.ReichsteinM.LuyssaertS.et al. (2008). Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature451, 49–52. 10.1038/nature06444
47
PiaoS. L.FriedlingsteinP.CiaisP.ViovyN.DemartyJ. (2007). Growing season extension and its impact on terrestrial carbon cycle in the Northern Hemisphere over the past 2 decades. Glob. Biogeochem. Cycles21, GB3018. 10.1029/2006GB002888
48
PiaoS. L.LiuQ.ChenA. P.JanssensI. A.FuY. S.DaiJ. H.et al. (2019). Plant phenology and global climate change: current progresses and challenges. Glob. Change Biol.25, 1922–1940. 10.1111/gcb.14619
49
PiaoS. L.TanJ. G.ChenA. P.FuY. S.CiaisP.LiuQ.et al. (2015). Leaf onset in the northern hemisphere triggered by daytime temperature. Nat. Commun.6, 6911. 10.1038/ncomms7911
50
PiaoS. L.WangX. H.ParkT.ChenC.LianX.HeY.et al. (2020). Characteristics, drivers and feedbacks of global greening. Nat. Rev. Earth Environ.1, 14–27. 10.1038/s43017-019-0001-x
51
QuanQ.TianD. S.LuoY. Q.ZhangF. Y.CrowtherT.ZhuK.et al. (2019). Water scaling on net ecosystem exchange in response to warming. Sci. Adv.5, eaav1131. 10.1126/sciadv.aav1131
52
RenS. L.PeichlM. (2021). Enhanced spatiotemporal heterogeneity and the climatic and biotic controls of autumn phenology in northern grasslands. Sci. Total Environ.788, 147806. 10.1016/j.scitotenv.2021.147806
53
ShenM. G.PiaoS. L.CongN.ZhangG. X.JanssensI. A. (2015). Precipitation impacts on vegetation spring phenology on the Tibetan Plateau. Glob. Change Biol.21, 3647–3656. 10.1111/gcb.12961
54
SunS. C.FrelichL. E. (2011). Flowering phenology and height growth pattern are associated with maximum plant height, relative growth rate and stem tissue mass density in herbaceous grassland species. J. Ecol.99, 991–1000. 10.1111/j.1365-2745.2011.01830.x
55
VitasseY.BaumgartenF.ZohnerC. M.KaewthongrachR.FuY. S.WaldeM. G.et al. (2021). Impact of microclimatic conditions and resource availability on spring and autumn phenology of temperate tree seedlings. N. Phytol.232, 537–550. 10.1111/nph.17606
56
WangC.TangY. (2019). Responses of plant phenology to nitrogen addition: a meta-analysis. Oikos128, 1243–1253. 10.1111/oik.06099
57
WangD.ChiZ. S.YueB. J.HuangX. D.ZhaoJ.SongH. Q.et al. (2020). Effects of mowing and nitrogen addition on the ecosystem C and N pools in a temperate steppe: a case study from northern China. Catena185, 104332. 10.1016/j.catena.2019.104332
58
WangD.HuangX. D.QiaoN.GengQ. L.LiuY. Z.SongH. Q.et al. (2021). Effects of mowing and fertilization on soil quality in a semiarid grassland of North China. Land Degrad. Dev.32, 1656–1666. 10.1002/ldr.3783
59
WangG. C.LuoZ. K.HuangY.XiaX. A.WeY. R.LinX. H.et al. (2022). Preseason heat requirements and days of precipitation jointly regulate plant phenological variations in Inner Mongolian grassland. Agricult. Forest Meteorol.314, 108783. 10.1016/j.agrformet.2021.108783
60
WayD. A.MontgomeryR. A. (2015). Photoperiod constraints on tree phenology, performance and migration in a warming world. Plant Cell Environ.38, 1725–1736. 10.1111/pce.12431
61
WhiteM. A.RunningS. W.ThorntonP. E. (1999). The impacts of growing-season length variability on carbon assimilation and evapotranspiration over 88 years in the eastern US deciduous forest. Int. J. Biometeorol.42, 139–145. 10.1007/s004840050097
62
WinglerA.PrudyS.MacLeanJ. A.PourtauN. (2006). The role of sugars in integrating environmental signals during the regulation of leaf senescence. J. Exp. Bot.57:391–399. 10.1093/jxb/eri279
63
WohlfahrtG.PilloniS.HörtnaglL.HammerleA. (2010). Estimating carbon dioxide fluxes from temperate mountain grasslands using broad-band vegetation indices. Biogeosciences7, 683–694. 10.5194/bg-7-683-2010
64
WuC. Y.GonsamoA.ChenJ. M.KurzW. A.PriceD. T.LafleurP. M.et al. (2012). Interannual and spatial impacts of phenological transitions, growing season length, and spring and autumn temperatures on carbon sequestration: a North America flux data synthesis. Glob. Planet. Change 92–93, 179–190. 10.1016/j.gloplacha.2012.05.021
65
XiaJ. Y.NiuS. L.CiaisP.JanssensI. A.ChenJ. Q.AmmannC.et al. (2015). Joint control of terrestrial gross primary productivity by plant phenology and physiology. Proc. Natl. Sci. Acad. U. S. A.112, 2788–2793. 10.1073/pnas.1413090112
66
XiangX.GibbonsS.HeH.-S.WangC.HeD.LiQ.et al. (2016). Rapid response of arbuscular mycorrhizal fungal communities to short-term fertilization in an alpine grassland on the Qinghai-Tibet Plateau. PeerJ.4, e2226. 10.7717/peerj.2226
67
YangN.ZavišićA.PenaR.PolleA. (2016). Phenology, photosynthesis, and phosphorus in European beech (Fagus sylvatica L.) in two forest soils with contrasting P contents. J. Plant Nutr. Soil Sci.179, 151–158. 10.1002/jpln.201500539
68
Yin T. Zheng L. Cao G. Song M. and Yu F. (2016). Species-specific phenological responses to long-term nitrogen fertilization in an alpine meadow. J. Plant Ecol.10, 301–309. 10.1093/jpe/rtw026
69
YuG. R.JiaY. L.HeN. P.ZhuJ. X.ChenZ.WangQ. F.et al. (2019). Stabilization of atmospheric nitrogen deposition in China over the past decade. Nat. Geosci.12, 424–429. 10.1038/s41561-019-0352-4
70
ZhangY. H.LoreauM.HeN. P.ZhangG. M.HanX. G. (2017). Mowing exacerbates the loss of ecosystem stability under nitrogen enrichment in a temperate grassland. Funct. Ecol.31, 1637–1646. 10.1111/1365-2435.12850
71
ZhouZ. X.LiY.SongJ.RuJ. Y.LeiL. J.ZhongM. X.et al. (2019). Growth controls over flowering phenology response to climate change in three temperate steppes along a precipitation gradient. Agricult. Forest Meteorol.274, 51–60. 10.1016/j.agrformet.2019.04.011
72
ZhouZ. X.ZhangK. P.SunZ. L.LiuY. Z.ZhangY. C.LeiL. J.et al. (2022). Lengthened flowering season under climate warming: evidence from manipulative experiments. Agricult. Forest Meteorol.312, 108713. 10.1016/j.agrformet.2021.108713
Summary
Keywords
grasslands, growing season, nitrogen, precipitation, phenology
Citation
Zhou Z, Zhang L, Liu Y, Zhang K, Wang W, Zhu J, Chai S, Zhang H and Miao Y (2022) Contrasting Effects of Nitrogen Addition on Vegetative Phenology in Dry and Wet Years in a Temperate Steppe on the Mongolian Plateau. Front. Plant Sci. 13:861794. doi: 10.3389/fpls.2022.861794
Received
25 January 2022
Accepted
21 March 2022
Published
25 April 2022
Volume
13 - 2022
Edited by
Veronica De Micco, University of Naples Federico II, Italy
Reviewed by
Maria Faticov, Université de Sherbrooke, Canada; Tao Zhang, Northeast Normal University, China
Updates
Copyright
© 2022 Zhou, Zhang, Liu, Zhang, Wang, Zhu, Chai, Zhang and Miao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuan Miao miaoyuan0921@126.com
†These authors have contributed equally to this work and share first authorship
This article was submitted to Functional Plant Ecology, a section of the journal Frontiers in Plant Science
Disclaimer
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.