Abstract
Soil eutrophication from atmospheric deposition and fertilization threatens biodiversity and the functioning of terrestrial ecosystems worldwide. Increases in soil nitrogen (N) and phosphorus (P) content can alter the biomass and structure of plant communities in grassland ecosystems; however, the impact of these changes on plant–pollinator interactions is not yet clear. In this study, we tested how changes in flowering plant diversity and composition due to N and P enrichment affected pollinator communities and pollination interactions. Our experiments, conducted in a Tibetan alpine grassland, included four fertilization treatments: N (10 g N m–2 year–1), P (5 g P m–2 year–1), a combination of N and P (N + P), and control. We found that changes in flowering plant composition and diversity under the N and P treatments did not alter the pollinator richness or abundance. The N and P treatments also had limited effects on the plant–pollinator interactions, including the interaction numbers, visit numbers, plant and pollinator species dissimilarity, plant–pollinator interaction dissimilarity, average number of pollinator species attracted by each plant species (vulnerability), and average number of plant species visited by each pollinator species (generality). However, the N + P treatment increased the species and interaction dissimilarity in flowering plant and pollinator communities and decreased the generality in plant–pollinator interactions. These data highlight that changes in flowering plants caused by N + P enrichment alter pollination interactions between flowering plants and pollinators. Owing to changes in flowering plant communities, the plant–pollinator interactions could be sensitive to the changing environment in alpine regions.
Introduction
The structure and function of biodiversity, such as the biomass and diversity of plants and the animals that depend on them for survival, are vital for the sustainability of ecosystems (). However, biodiversity loss caused by anthropogenic nutrient enrichment and climate change threatens the functions and services of terrestrial ecosystems, particularly grasslands (; ; ). For example, land-use change and environmental pollution have contributed to a decline in biodiversity worldwide (). To restore grassland services and functions, specific agricultural management strategies, such as livestock exclusion and chemical fertilizer application, are employed to improve grassland productivity and other ecological services, such as pollination (). Plant–pollinator interactions can change through an increase in dominance and a decrease in the diversity of plants in grasslands under nutrient enrichment through agricultural fertilization and atmospheric deposition (Zarzycki and Kopec, 2020; ). In addition, the number and diversity of wild pollinators in natural and agricultural ecosystems have declined because of climate change, habitat loss, pathogen transport, commercially managed pollinators, agrochemicals, and nutrient enrichment (; ; ; ), which alter the stability of pollination ecosystems. However, since recent reviews predict that N enrichment might disrupt or enhance individual plant–pollinator interactions in grassland communities (; ), it is not yet clear how the changes in plants due to multiple nutrient enrichment in soil alter plant–pollinator interaction networks.
The availability of nitrogen (N) is considered to be the most important nutritionally limiting factor for primary productivity in grasslands (; ; Xiao et al., 2020). Although the N supply is increasing worldwide (; ; Wang et al., 2017), N availability has decreased in many regions of the world due to increased temperatures and CO2 (; ), dramatically reducing the biodiversity of terrestrial ecosystems (; ). Nutrient co-limitation is common in the biomass of plant communities (; ). Phosphorus (P) enrichment can increase soil P bioavailability and accelerate N uptake by plants, which can increase plant biomass and change the structure of plant communities (; Xiao et al., 2020). This reshaping of the composition of plant communities through nutrient enrichment can change the animal and microbial diversity and abundance above and below ground (; ; Zi et al., 2022), in addition to potentially altering biodiversity at multiple trophic levels within the community’s food webs (; ; Xiao et al., 2020; ). Some studies on the responses of pollinator communities and plant–pollinator interaction networks to nutrient supply in grassland ecosystems have been conducted, but their results have been inconsistent (; ; ). For example, found that N enrichment exerted no effect on pollinator assemblages and the community network structure between plants and pollinators, but showed that fertilizer addition decreased plant species richness, floral abundance and bumblebee richness. Therefore, there is an urgent need to focus on the consequences of soil nutrient intake on pollinator assemblages and plant–pollinator interaction networks in other grassland ecosystems from agricultural activities and environmental changes.
The Tibetan grassland ecosystem covers over 60% of the Qinghai-Tibet Plateau; however, this region is experiencing increased atmospheric N deposition and local warming (; ). These grasslands are fragile ecosystems owing to their high altitude and low temperature and are sensitive to human activities and global changes (; ; ). Previous studies have shown that adding N and P can increase the net primary productivity of grasses, while predominantly reducing the biomass and abundance of legumes on the alpine grasslands (; ; ). Many recent studies have focused on the effects of N and P addition on plants and soil microorganisms in the alpine grasslands (Zong and Shi, 2019; ; ; Zi et al., 2022), showing that combined N and P enrichment decreased the diversity of plants but increased net primary productivity (; Xiao et al., 2020). However, because bottom-up effects driven by changes in flowering plant diversity and flower abundance could have a major impact on the plant–pollinator interactions, few studies have focused on how changes in plant community structure due to N and P enrichment affect the richness and abundance of pollinators in this alpine grassland ecosystem.
Here, we experimentally manipulated the addition of nutrients to investigate whether and how the addition of P and N affected flowering plant communities, pollinator assemblages, and plant–pollinator interactions (Figure 1). In addition, we used a network approach to describe the interaction networks between plants and pollinators, which is commonly used to study changes in plant–pollinator interactions in communities (; ; ). Previous research showed that nutrient addition reduced the diversity of plants in grasslands (; Xiao et al., 2020; ). Therefore, we hypothesized (Figure 1D) that the co-addition of N and P (1) would decrease the diversity of flowering plants and alter the composition of flowering plant communities (Hypothesis 1); (2) decrease the diversity of pollinators and change the structure of pollinator assemblages due to reduced legume biomass and relative abundance (Hypothesis 2); and (3) change plant–pollinator interactions across trophic levels due to changes in the relative abundances of flowering plants and pollinators (Hypothesis 3). The insights gained from this research will improve our understanding of how global changes, such as atmospheric nitrogen deposition and fertilization, can maintain and improve ecosystem biodiversity and the structure of interaction networks between plants and pollinators and make grassland ecosystems more sustainable.
FIGURE 1
Materials and methods
Study area
This study was conducted at the Haibei Alpine Grassland Ecosystem Research Station (37° 36′ N, 101° 12′ E, 3,250 m above sea level), Qinghai, China (Figure 1A). The research station has a typical plateau continental climate (short and cool summers, and long and cold winters), with a mean annual temperature and precipitation of –1.1°C and 480 mm, respectively. The precipitation (about 80%) is mainly concentrated in the growing season (May–September). The soil in this area is classified as Mat-Gryic Cambisols (WRB, 1998). In 2011, the soil content 10 cm underground in the station sampling site contained total carbon, 78.2 g kg–1; organic carbon, 63.1 g kg–1; nitrogen, 5.75 g kg–1; and phosphorus, 0.79 g kg–1 (Zhang et al., 2021). The plant community mainly consists of grasses and forbs (Figure 1B). The flowering plant species are mainly Gentiana straminea, Gentiana aristata, Morina kokonorica, Angelica nitida, and Dasiphora fruticosa (), which are primarily pollinated by insects such as honey bees, bumble bees, butterflies and flies (Zhang et al., 2015; Wang et al., 2021).
Nitrogen and phosphorus manipulations
From May 2011, four nutrient-addition experiments were conducted in an alpine grassland: N addition, P addition, a combination of N and P addition, and no nutrient addition (Figure 1C). The fertilizers were 100 kg of urea (CH4N2O) and 50 kg of triple superphosphate [Ca(H2PO4)2⋅CaHPO4] per hectare per year. Our experiment included 24 plots (4 treatments × 6 blocks) at 6 m × 6 m each (Figure 1C), which is comparable to other studies (; ). In the experimental plots, the legumes and forbs mainly include Oxytropis kansuensis, Gentiana straminea, A. nitida, Aster farreri, Tibetia himalaica, and Thermopsis lanceolata. To maintain the plant species’ similarity and reduce the impairment caused by different fertilization treatments on adjacent plots, the blocks were spaced 2 m apart, and the plots within each block were spaced 1 m apart. We divided the fertilizer into three portions and spread it evenly from June to August each year ().
Flower and pollinator survey
To examine the change in the diversity and flower abundance of flowering plants after the N and P addition treatments for eight consecutive years (2011–2018), we created four transects (1 m × 5 m) within each plot in 2019. From early July to late August 2019, we counted the number of flowering plant species and flowers for each plant four times per half month along the four transects of each plot.
To examine the diversity and number of pollinating insects, from early July to late August 2019, we monitored flowers for 30 min in each transect (1 m × 5 m) four times on clear days in the presence of strong wind (09: 00–19: 00). Using this observation method, each plant species was observed in each transect for 120 min. If an insect encountered the sexual components of the flowers (anthers or stigmas), we recorded a pollination visit without considering the effectiveness of the pollination (). As few pollinators visited the flowers in alpine meadows at low temperatures, we did not observe pollinator visits at night (). The pollinator abundance depended on the number of visits to flowering plants per plot for each pollinator species. With the help of taxonomic experts, we identified the species or higher-level groups of unknown visitors.
Pollination network construction and parameter calculation
To estimate the completeness of the sample, we divided the observed richness by the estimated abundance-based richness estimator, Chao1, using the vegan package (v2.5-7) with the R statistical software (v4.0.3, ).
To determine the effects of nutrient addition on the diversity and abundance of flowering plants and pollinators, we counted the flowering plant species, pollinator species, flowers, pollinator individuals, plant–pollinator unique interactions, and visits per plot.
To determine whether the core generalist flowering plant and pollinator species changed with different nutrient treatments, we used UCINET (v6.0) to calculate the eigenvector centrality score (network > centrality > eigenvector in UCINET) of each plant and pollinator species of different nutrient treatments (). In addition, when one plant or pollinator species with a high eigenvector centrality score participated in more than 5% of the visits and interacted with more than 25% of the taxa in each nutrient treatment, we delimited them as core generalists (). We then analyzed the core generalized plant and pollinator species to determine whether they had changed due to the nutrient addition treatments.
To determine the effects of nutrient addition on the dissimilarity between flowering plants and pollinators, we quantified the species and interaction dissimilarity (species or link differences in plant–pollinator interaction networks) based on the methods of . We ran all combinations of networks using the “betalinkr” function in the bipartite package (). According to , the dissimilarity of a pair of networks was divided into the dissimilarity due to differences in species composition (species dissimilarity) and dissimilarity due to interaction rewiring (interaction dissimilarity). We focused on the change in species and interaction composition involving species of plants and pollinators that were present for the control and nutrient-addition treatments. These two metrics quantify species and link changes in species interaction networks under different nutritional treatments (). A value of 0 indicates that the species/interaction compositions are identical, and a value of 1 indicates that the species/interaction compositions are completely different.
To determine the impact of nutrient addition on the changes in links between plants and pollinators, we calculated two quantitative network metrics, vulnerability and generality, using the methods of . We used the “networklevel” function in the bipartite package () to calculate the values of vulnerability and generality. The vulnerability and generality were weighted according to their marginal sums (). The vulnerability corresponds to the average number of pollinator species attracted by each plant species, and the generality corresponds to the average number of plant species visited by each pollinator species (). The changes in vulnerability and generality reflect the effects of nutrient addition on the distribution of links between plants and pollinators (). An increase in vulnerability indicates that pollinator communities are increasingly dependent on fewer plant species. Conversely, an increase in generality demonstrates that the range of plant species visited by pollinators has increased, or that generalized pollinator species have increased their abundance.
Data analysis
The data were checked for normality and homogeneity of variance before performing the ANOVA. If the assumptions of normality and variance were not met, the data were log-transformed. We then used two-way ANOVAs in the R basic package to determine the effects of nutrient enrichment on plants, pollinators, and their interactions. The N and P treatments were fixed factors, and the block was an error term. The variables included five flowering plant and pollinator variables (the flowering plant diversity, flower abundance, pollinator diversity, pollinator abundance, and species dissimilarity of plants and pollinators) and five pollination interaction variables (the interaction numbers, number of visits, interaction dissimilarity, vulnerability, and generality). In addition, three other statistical analyses (one-way ANOVA, the pairwise permutational multivariate analysis of variance, and piecewise structural equation modeling) were performed.
First, we used a one-way ANOVA to examine the effects of the nutrient-addition treatments (control, N, P, and N + P treatments) on each response variable (the flowering plant diversity, pollinator diversity, flower abundance, pollinator abundance, interaction numbers, number of visits, species dissimilarity, interaction dissimilarity, vulnerability, and generality).
Second, we applied the pairwise permutational multivariate analysis of variance (PERMANOVA) using the vegan package (v2.5-7) () to evaluate the effects of the addition of N and P on the plant and pollinator diversity. To illustrate the impacts of the N and P treatments on changes in the structure of plant and pollinator communities, we performed paired Bray–Curtis distance principal coordinate analysis (PCoA).
Finally, we applied piecewise structural equation modeling in the piecewise SEM package (v2.1.2) () to clarify the direct and indirect effects of the addition of N and P on changes in plant–pollinator interactions (Supplementary Figure 1). The hypothesized direct effects of N and P enrichment, and the indirect effects mediated by plant and pollinator community (i.e., plant diversity, flower abundance, pollinator diversity, and pollinator abundance) on vulnerability and generality are included in Supplementary Figure 1. The piecewise SEM included several linear mixed-effects models, and the block was a random effect. The full piecewise SEM comprised the effects of N and P’s effects on four mediators of plant and pollinator communities (plant diversity, flower abundance, pollinator diversity, and pollinator abundance), as well as the direct and indirect effects on the vulnerability and generality of plants and pollinators. To simplify the model, we did not consider the interactions between the N and P treatments.
Results
In total, we found that 20 flowering plant species (Supplementary Table 1) were visited by 54 pollinator species (Supplementary Table 2) that participated in 4954 individual interactions (Figure 2). Based on the sample completeness analysis, these data represented about 64% of the pollinator species for each treatment.
FIGURE 2
Structural changes from N and P addition to flowering plant communities
The N treatment decreased the flowering plant diversity, and the N + P treatment increased the negative effect on the flowering plant diversity (Figure 2 and Supplementary Figure 2A). In particular, the N treatment reduced the legume diversity (Figure 2, e.g., O. kansuensis, Oxytropis ochrocephala, and Oxytropis qinghaiensis). By contrast, the N + P treatment further reduced the diversity of rare species (Figure 2, e.g., Ranunculus membranaceus, Euphrasia regelii, and Galium verum). The abundance of flowers was similar to that for the control and N or P treatments (Supplementary Figure 2B). However, the N + P treatment increased the abundance of flowers (Supplementary Figure 2B), primarily by increasing the number of flowers of Descurainia sophia.
The N + P treatment altered the numbers and identities of the core generalist plants (Supplementary Figure 3A). For example, four plant species (Morina kokonorica, A. nitida, Anaphalis lacteal, and Aster farreri) were core generalists with the control treatment. By contrast, only two species (A. nitida and A. lacteal) were core generalists with the P treatment (Supplementary Figure 3A). Additionally, Potentilla saundersiana dominated the N treatment group as the core plant species at the expense of A. nitida. To the detriment of M. kokonorica and A. lacteal, D. sophia became a core plant species in the N + P treatment (Supplementary Figure 3A). PERMANOVA analyses revealed that the N + P treatment changed the structure of the flowering plant communities (Figure 3A).
FIGURE 3
Effects of structural changes in flowering plants on pollinator communities
Flies were the most abundant pollinators (Figure 2). The N and P treatments did not affect the diversity or abundance of pollinators (Supplementary Figures 2C,D). The N, P, and N + P treatments did not change the numbers and identities of the core pollinators (Supplementary Figure 3B), except that the N treatment decreased the diversity of the bees and the N + P treatment increased the diversity of the butterflies (Figure 2). None of the nutrient additions changed the structure of the pollinator communities (Figure 3B).
The impacts of nutrient additions on plant–pollinator interactions
The N, P, and N + P treatments did not change the number of unique interactions between plants and pollinators (Figure 4A) or the total number of visits to the pollination networks (Figure 4B). The N or P treatment alone did not change the species and interaction dissimilarity. However, compared to the control treatment, the N + P treatment increased the species and interaction dissimilarity in the pollination networks (Figures 4C,D). None of the treatments altered the vulnerability (Figure 4E). However, the N + P treatment decreased the generality (Figure 4F).
FIGURE 4
The full piecewise SEM (Supplementary Figure 4 and Supplementary Table 3; χ2 = 1.629, df = 6, P = 0.95, and AICc = –108.163) and the final model are qualitatively similar (Figure 5 and Supplementary Table 3; χ2 = 7.819, df = 14, P = 0.899, AICc = –123.88, and ΔAICc = 15.72), which explained about 75% of the variance of the generality (marginal R2 = 0.75). The final model revealed that the change in flowering plant species caused by the N and P treatments indirectly affected the generality through changing the flower abundance and vulnerability (Figure 5 and Supplementary Figure 4). The N or P treatment alone did not affect the richness and abundance of the pollinators and vulnerability (all P > 0.05, Figure 5). However, the N treatment directly decreased the flowering plant diversity (standardized path coefficient β = –0.809, P < 0.0001) and increased the generality (β = 0.323, P = 0.03). In addition, the P treatment directly increased the abundance of flowers (β = 0.481, P = 0.004) and deceased the generality (β = –0.654, P = 0.001).
FIGURE 5
Discussion
Examining how eutrophication due to climate change and anthropogenic activities affects communities of plants, animals, and species interactions is vital for understanding how biological communities respond to global changes (). In the current study, we investigated the effects of N and P addition on plant–pollinator interactions in an alpine grassland. Our findings revealed that N and P addition changed the structure of flowering plant communities. Although pollinator communities were resilient to nutrient addition, the plant–pollinator interactions were unstable under the N + P treatment in alpine grassland ecosystems. These results suggest that N- and P-fertilization-induced changes in flowering plant composition have a minor effect on pollinator communities. However, the plant–pollinator interactions can be sensitive to nutrient addition through agricultural activities and environmental changes.
Structural changes in flowering plant communities under N and P addition
Soil eutrophication due to N and P addition can increase the biomass of grasses and decrease the biomass and diversity of flowering plant species, such as forb and legumes, in nutrient-limited grasslands (; ). For example, fast-growing grasses can become taller, denser, and more competitive for sunlight under N enrichment. As a result, the biomass and diversity of nitrophobous forb species are often reduced and they can become locally extinct because of their small sizes and slow growth rates (; ). The results of the present study support our first hypothesis that N enrichment reduces the diversity of flowering plants, primarily that of legumes (Figure 2). However, our findings did not agree with those of , who found that the richness of legumes did not change under N enrichment in alpine grassland (). A possible explanation is that the diversity of legumes deceased and disappeared over time, suggesting that the time and intensity of N addition played a vital role in their local extinction. For example, the effect of fertilization was not fully visible within only three years () because the perennial legumes (e.g., O. kansuensis, O. ochrocephala, and O. qinghaiensis) can only change their distribution over generations, so they can become locally extinct under long-term fertilization.
Additionally, this research revealed that P addition further enhanced the negative effect of N enrichment on flowering plant species in alpine grasslands (Figure 2). For example, R. membranaceus, E. regelii, and G. verum disappeared after eight years of N + P enrichment. N enrichment affects the N and P balance in the soil, changing the situation from N limitation to P limitation (; ; Zhan et al., 2017). As P addition can also reduce the limitations for crucial nutrient resources (), this could consequently alter the competitive dynamics of plants, typically leading to a reduction in legumes and endangered forbs (; ). For example, N and P co-addition increased flower production because of the increased number of D. sophia plants (Supplementary Figure 2B), as N enrichment can significantly increase the height, biomass, and seed yield of D. sophia (). Furthermore, the changes in flowering plant diversity and flower abundance led to a further reorganization of the flowering plant community structure (Figure 3A), indicating that flowering plant communities are unstable under long-term P and N enrichment in alpine grassland ecosystems. However, future experiments with more time scales and fertilization intensities are still needed to examine the effects of nutrient enrichment on the structure of plant communities.
Weak effects of structural changes in flowering plants on pollinator communities
Contrary to our second hypothesis, our findings revealed that N and P addition did not alter the pollinator diversity and abundance and the core pollinator species (Supplementary Figures 2, 3). also found that N addition did not affect the richness of pollinator communities. A possible explanation might be that most of the remaining flowering plant species, such as Asteraceae, Ranunculaceae, and Rosaceae, had an open morphology with radially symmetrical flowers in the nutrient-addition plots (Figure 2 and Supplementary Table 1), which could attract many generalized pollinator species to visit these flowering plants (; Zhang et al., 2015; Wang et al., 2021). For example, we found that the flowers of D. sophia and Angelica nitida could attract more than ten pollinator species. Additionally, our results showed that the main pollinator species were flies, bees, and butterflies (Figure 2), which are generalized pollinators that can visit many flowering plants in alpine grasslands (Wang et al., 2021). Therefore, we could not determine any changes in pollinator diversity and abundance due to P and N addition (; ).
It is important to note that, similarly to in other studies (Supplementary Table 5), the scale of the experimental treatments for the plants was relatively small, which would probably have not affected the populations of pollinators but, rather, simply their behavior in terms of which plots they decided to forage in. However, our results showed that N and P co-addition decreased the richness of bumble bees but increased the diversity of butterflies (Figure 2). Changes in the flowering plant species composition can negate the ability of pollinators to choose different resource supplies (; ; ). Our results showed that bumble bees mainly visited the legumes, suggesting that they have specialized diets, which may explain why they are more susceptible to declines than flies and butterflies. For example, researchers have suggested that significant spatial changes in flowering plant species due to N enrichment would decrease the diversity of pollinators, such as bees and butterflies (). Due to the mismatch in the scales at which pollinators and plants respond to the nutrient treatments, large-spatial-scale studies on the effects of different nutrient enrichments on the populations of pollinators will be investigated in the future.
The impacts of nutrient additions on plant–pollinator interactions
Contrary to the results of previous studies, found no effects of N addition on the identity and frequency of plant–pollinator interactions. Our results demonstrate that pollinators increasingly rely on fewer flowering plant species and that pollinators visit a reduced range of plant species. Thus, plant–pollinator interactions are sensitive to changes in flowering plant composition and flower abundance, and the variation in pollinator behavior among different plants that were affected by the N, P, and N + P enrichment treatments. These results support our third hypothesis that nutrient enrichment can change plant–pollinator interactions.
Plant–pollinator interaction networks always exhibit nested properties and are centered around a core of generalized plant and pollinator species (; ). Thus, changes in the composition of the core generalized plant species and floral abundance can directly affect the structure of the pollination network, such as nestedness (i.e., interactions between the most generalist plants and animals create a dense core that includes other specialized plants and animals in its community). The nested structure of the network can reduce the extent to which species and interactions are affected by perturbation (; ), indicating that pollination networks are highly tolerant of plant extinctions due to their nestedness structure.
Furthermore, the addition of N or P alone did not change the core generalized flowering plants, which may not change the nestedness of plant–pollinator interaction networks, because of the generalized characteristics of pollination networks (Waser et al., 1996; ) and the rapid turnover of interactions due to the loss and gain of non-core plant and pollinator species (; ). We also cannot rule out the possibility that the experimental plots were too small to attract different pollinators or that the pollinators are not specific for a given flowering plant species, as mentioned above. Although previous studies revealed that N and P enrichment reduces the biomass and abundance of legumes and the overall diversity of plant species in alpine grasslands (; ), the loss of non-core species did not change the structure of the pollination networks (Weiner et al., 2011). Thus, changes in the composition and flower abundance of the core generalized flowering plant species due to the addition of N and P can influence the link distribution of the pollination networks.
Conclusion
We tested the importance of nutrient enrichment for changes in pollination interactions by combining network theory and pollination ecology. Our study revealed that changes in flowering plant diversity and composition, both of which are reduced by nitrogen and phosphorus enrichment, alter plant–pollinator interactions in a Tibetan alpine grassland. Future experiments in manipulating nutritional resources, and the composition and abundance of plant and pollinator species on larger temporal and spatial scales will provide important insights into how nutrient enrichment affects the response of pollinator assemblages.
Statements
Data availability statement
The original contributions presented in this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
L-LW, Y-PY, Y-WD, and J-SH conceived the ideas and designed methodology. L-LW, FR, CZ, Z-HZ, and X-JH collected the data. L-LW analyzed the data. L-LW, Y-PY, and Y-WD led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.
Funding
This work was supported by the Second Tibetan Scientific Expedition and Research (STEP) program (grant number: 2019QZKK0502), the Strategic Priority Research Program of Chinese Academy of Sciences (grant number: XDA20050000), China Postdoctoral Science Foundation (grant number: 2021M693235), Postdoctoral Directional Training Foundation of Yunnan Province to L-LW, National Natural Science Foundation of China (grant number: 31901171), and the Project of Qinghai Science & Technology Department (2018-ZJ-967Q to FR).
Acknowledgments
We thank L. Zhang, J. F. Wu, J. J. Gao, P.J. Zhang, and X. G. Zhang for the fieldwork and the Haibei Alpine Meadow Ecosystem Research Station staff for their assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.964109/full#supplementary-material
References
1
AlarcónR.WaserN. M.OllertonJ. (2008). Year-to-year variation in the topology of a plant-pollinator interaction network.Oikos1171796–1807. 10.1111/j.0030-1299.2008.16987.x
2
BascompteJ.JordanoP.MeliánC. J.OlesenJ. M. (2003). The nested assembly of plant-animal mutualistic networks.Proc. Natl. Acad. Sci. U. S. A.1009383–9387. 10.1073/pnas.1633576100
3
BaudeM.KuninW. E.BoatmanN. D.ConyersS.DaviesN.GillespieM. A. K.et al (2016). Historical nectar assessment reveals the fall and rise of floral resources in Britain.Nature53085–88. 10.1038/nature16532
4
BersierL. F.Banašek-RichterC.CattinM. F. (2002). Quantitative descriptors of food-web matrices.Ecology832394–2407. 10.2307/3071801
5
BurkleL.IrwinR. (2009). The importance of interannual variation and bottom-up nitrogen enrichment for plant-pollinator networks.Oikos1181816–1829. 10.1111/j.1600-0706.2009.17740.x
6
CaraDonnaP. J.PetryW. K.BrennanR. M.CunninghamJ. L.BronsteinJ. L.WaserN. M.et al (2017). Interaction rewiring and the rapid turnover of plant-pollinator networks.Ecol. Lett.20385–394. 10.1111/ele.12740
7
CarvalheiroL. G.BiesmeijerJ. C.FranzénM.Aguirre-GutiérrezJ.GaribaldiL. A.HelmA.et al (2019). Soil eutrophication shaped the composition of pollinator assemblages during the past century.Ecography43209–221. 10.1111/ecog.04656
8
CeulemansT.MerckxR.HensM.HonnayO. (2013). Plant species loss from European semi-natural grasslands following nutrient enrichment - is it nitrogen or is it phosphorus?Global Ecol. Biogeogr.2273–82. 10.1111/j.1466-8238.2012.00771.x
9
ChenH.ZhuQ. A.PengC. H.WuN.WangY. F.FangX. Q.et al (2013). The impacts of climate change and human activities on biogeochemical cycles on the Qinghai-Tibetan Plateau.Global Change Biol.192940–2955. 10.1111/gcb.12277
10
ClarkC. M.ClelandE. E.CollinsS. L.FargioneJ. E.GoughL.GrossK. L.et al (2007). Environmental and plant community determinants of species loss following nitrogen enrichment.Ecol. Lett.10596–607. 10.1111/j.1461-0248.2007.01053.x
11
CrawleyM. J.JohnstonA. E.SilvertownJ.DoddM.de MazancourtC.HeardM. S.et al (2005). Determinants of species richness in the park grass experiment.Am. Nat.165179–192. 10.1086/427270
12
DavidT. I.StorkeyJ.StevensC. J. (2019). Understanding how changing soil nitrogen affects plant-pollinator interactions.Arthropod Plant Interact.13671–684. 10.1007/s11829-019-09714-y
13
DongJ.CuiX.NiuH.ZhangJ.ZhuC.LiL.et al (2022). Effects of nitrogen addition on plants and microorganisms under high phosphorus addition level in the Tibetan alpine steppe.Front. Plant Sci.13:894365. 10.3389/fpls.2022.894365
14
DormannC. F.FruendJ.GruberB. (2021). Bipartite: Visualizing Bipartite Networks and Calculating Some (Ecological) Indices. R package version 2.15. Available online at: https://cran.r-project.org/web/packages/bipartite/(accessed April 20, 2022).
15
DuanY. W.ZhangT. F.LiuJ. Q. (2007). Interannual fluctuations in floral longevity, pollinator visitation and pollination limitation of an alpine plant (Gentiana straminea Maxim., Gentianaceae) at two altitudes in the Qinghai-Tibetan Plateau.Plant Syst. Evol.267255–265. 10.1007/s00606-007-0553-y
16
ElserJ. J.BrackenM. E. S.ClelandE. E.GrunerD. S.HarpoleW. S.HillebrandH.et al (2007). Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems.Ecol. Lett.101135–1142. 10.1111/j.1461-0248.2007.01113.x
17
FangQ.HuangS. Q. (2016). Plant-pollinator interactions in a biodiverse meadow are rather stable and tight for 3 consecutive years.Integr. Zool.11199–206. 10.1111/1749-4877.12190
18
GallowayJ. N.TownsendA. R.ErismanJ. W.BekundaM.CaiZ.FreneyJ. R.et al (2008). Transformation of the nitrogen cycle: recent trends, questions, and potential solutions.Science320889–892. 10.1126/science.1136674
19
GaoE. L.WangY. X.BiC.Kaiser-BunburyC. N.ZhaoZ. G. (2021). Restoration of degraded alpine meadows improves pollination network robustness and function in the Tibetan Plateau.Front. Ecol. Evol.9:632961. 10.3389/fevo.2021.632961
20
GhazoulJ. (2015). Qualifying pollinator decline evidence.Science348981–982. 10.1126/science.348.6238.981-b
21
HarpoleW. S.NgaiJ. T.ClelandE. E.SeabloomE. W.BorerE. T.BrackenM. E. S.et al (2011). Nutrient co-limitation of primary producer communities.Ecol. Lett.14852–862. 10.1111/j.1461-0248.2011.01651.x
22
HarpoleW. S.SullivanL. L.LindE. M.FirnJ.AdlerP. B.BorerE. T.et al (2016). Addition of multiple limiting resources reduces grassland diversity.Nature53793–96. 10.1038/nature19324
23
HautierY.NiklausP. A.HectorA. (2009). Competition for light causes plant biodiversity loss after eutrophication.Science324636–638. 10.1126/science.1169640
24
HeD.XiangX.HeJ.-S.WangC.CaoG.AdamsJ.et al (2016). Composition of the soil fungal community is more sensitive to phosphorus than nitrogen addition in the alpine meadow on the Qinghai-Tibetan Plateau.Biol. Fertility Soils521059–1072. 10.1007/s00374-016-1142-4
25
HooperD. U.AdairE. C.CardinaleB. J.ByrnesJ. E. K.HungateB. A.MatulichK. L.et al (2012). A global synthesis reveals biodiversity loss as a major driver of ecosystem change.Nature486105–108. 10.1038/nature11118
26
IPBES (2016). The Assessment Report on Pollinators, Pollination and Food Production. Secretariat of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services.Bonn: IPBES.
27
IPBES (2018). The IPBES Assessment Report on Land Degradation and Restoration. Secretariat of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services.Bonn: IPBES
28
IsbellF.ReichP. B.TilmanD.HobbieS. E.PolaskyS.BinderS. (2013). Nutrient enrichment, biodiversity loss, and consequent declines in ecosystem productivity.Proc. Natl. Acad. Sci. U. S. A.11011911–11916. 10.1073/pnas.1310880110
29
JohnsonS. D.SteinerK. E. (2000). Generalization versus specialization in plant pollination systems.Trends Ecol. Evol.15140–143. 10.1016/s0169-5347(99)01811-x
30
JordanoP.BascompteJ.OlesenJ. M. (2006). The Ecological Consequences of Complex Topology and Nested Structure in Pollination Webs.Chicago: University of Chicago Press.
31
Kaiser-BunburyC. N.MougalJ.WhittingtonA. E.ValentinT.GabrielR.OlesenJ. M.et al (2017). Ecosystem restoration strengthens pollination network resilience and function.Nature542223–227. 10.1038/nature21071
32
LefcheckJ. S. (2016). PIECEWISESEM: piecewise structural equation modelling in R for ecology, evolution, and systematics.Methods Ecol. Evol.7573–579. 10.1111/2041-210x.12512
33
LiuX.MaZ.CadotteM. W.ChenF.HeJ. S.ZhouS. (2019). Warming affects foliar fungal diseases more than precipitation in a Tibetan alpine meadow.New Phytol.2211574–1584. 10.1111/nph.15460
34
LiuX.ZhangL.HuangM. J.ZhouS. R. (2021). Plant diversity promotes soil fungal pathogen richness under fertilization in an alpine meadow.J. Plant Ecol.14323–336. 10.1093/jpe/rtaa099
35
LiuX.ZhangY.HanW.TangA.ShenJ.CuiZ.et al (2013). Enhanced nitrogen deposition over China.Nature494459–462. 10.1038/nature11917
36
LuoR.FanJ.WangW.LuoJ.KuzyakovY.HeJ. S.et al (2019). Nitrogen and phosphorus enrichment accelerates soil organic carbon loss in alpine grassland on the Qinghai-Tibetan Plateau.Sci. Total Environ.650303–312. 10.1016/j.scitotenv.2018.09.038
37
MasonR. E.CraineJ. M.LanyN. K.JonardM.OllingerS. V.GroffmanP. M.et al (2022). Evidence, causes, and consequences of declining nitrogen availability in terrestrial ecosystems.Science376:eabh3767. 10.1126/science.abh3767
38
McCannK. S. (2007). Protecting biostructure.Nature44629–29. 10.1038/446029a
39
MemmottJ. (1999). The structure of a plant-pollinator food web.Ecol. Lett.2276–280. 10.1046/j.1461-0248.1999.00087.x
40
Mokhtassi-BidgoliA.AghaAlikhaniM.Nassiri-MahallatiM.ZandE.Gonzalez-AndujarJ. L.AzariA. (2013). Agronomic performance, seed quality and nitrogen uptake of Descurainia sophia in response to different nitrogen rates and water regimes.Ind. Crop. Prod.44583–592. 10.1016/j.indcrop.2012.09.018
41
OksanenJ.BlanchetF. G.FriendlyM.KindtR.LegendreP.McGlinnD.et al (2020). Vegan: Community Ecology Package. R package version 2.5-7. Available online at: https://CRAN.R-project.org/package=vegan(accessed April 17, 2022).
42
OlffH.AertsR.BobbinkR.CornelissenJ. H. C.ErismanJ. W.GallowayJ. N.et al (2022). Explanations for nitrogen decline.Science3761169–1170. 10.1126/science.abq7575
43
OllertonJ. (2017). Pollinator diversity: distribution, ecological function, and conservation.Annu. Rev. Ecol. Evol. Syst.48353–376. 10.1146/annurev-ecolsys-110316-022919
44
PawarS. (2014). Why are plant-pollinator networks nested?Science345383–383. 10.1126/science.1256466
45
PoisotT.CanardE.MouillotD.MouquetN.GravelD. (2012). The dissimilarity of species interaction networks.Ecol. Lett.151353–1361. 10.1111/ele.12002
46
PottsS. G.BiesmeijerJ. C.KremenC.NeumannP.SchweigerO.KuninW. E. (2010). Global pollinator declines: trends, impacts and drivers.Trends Ecol. Evol.25345–353. 10.1016/j.tree.2010.01.007
47
R Core Team (2020). R: a Language and Environment for Statistical Computing. R version 4.0.3. Vienna: R Foundation for Statistical Computing.
48
RenF.SongW.ChenL.MiZ.ZhangZ.ZhuW.et al (2017). Phosphorus does not alleviate the negative effect of nitrogen enrichment on legume performance in an alpine grassland.J. Plant Ecol.10822–830. 10.1093/jpe/rtw089
49
RenF.YangX.ZhouH.ZhuW.ZhangZ.ChenL.et al (2016). Contrasting effects of nitrogen and phosphorus addition on soil respiration in an alpine grassland on the Qinghai-Tibetan Plateau.Sci. Rep.6:34786. 10.1038/srep34786
50
RogerN.MoermanR.CarvalheiroL. G.Aguirre-GuitierrezJ.JacquemartA. L.KleijnD.et al (2017). Impact of pollen resources drift on common bumblebees in NW Europe.Global Change Biol.2368–76. 10.1111/gcb.13373
51
StevensC. J.DavidT. I.StorkeyJ. (2018). Atmospheric nitrogen deposition in terrestrial ecosystems: its impact on plant communities and consequences across trophic levels.Funct. Ecol.321757–1769. 10.1111/1365-2435.13063
52
StorkeyJ.MacdonaldA. J.PoultonP. R.ScottT.KöehlerI. H.SchnyderH.et al (2015). Grassland biodiversity bounces back from long-term nitrogen addition.Nature528401–404. 10.1038/nature16444
53
SudingK. N.CollinsS. L.GoughL.ClarkC.ClelandE. E.GrossK. L.et al (2005). Functional- and abundance-based mechanisms explain diversity loss due to N fertilization.Proc. Natl. Acad. Sci. U. S. A.1024387–4392. 10.1073/pnas.0408648102
54
TylianakisJ. M.DidhamR. K.BascompteJ.WardleD. A. (2008). Global change and species interactions in terrestrial ecosystems.Ecol. Lett.111351–1363. 10.1111/j.1461-0248.2008.01250.x
55
Villa-GalavizE.SmartS. M.ClareE. L.WardS. E.MemmottJ. (2021). Differential effects of fertilisers on pollination and parasitoid interaction networks.J. Anim. Ecol.90404–414. 10.1111/1365-2656.13373
56
WangL. L.YangY. P.DuanY. W. (2021). Pollinator individual-based networks reveal the specialized plant-pollinator mutualism in two biodiverse communities.Ecol. Evol.1117509–17518. 10.1002/ece3.8384
57
WangR.GollD.BalkanskiY.HauglustaineD.BoucherO.CiaisP.et al (2017). Global forest carbon uptake due to nitrogen and phosphorus deposition from 1850 to 2100.Global Change Biol.234854–4872. 10.1111/gcb.13766
58
WaserN. M.ChittkaL.PriceM. V.WilliamsN. M.OllertonJ. (1996). Generalization in pollination systems, and why it matters.Ecology771043–1060. 10.2307/2265575
59
WeinerC. N.WernerM.LinsenmairK. E.BluethgenN. (2011). Land use intensity in grasslands: changes in biodiversity, species composition and specialisation in flower visitor networks.Basic Appl. Ecol.12292–299. 10.1016/j.baae.2010.08.006
60
WRB (1998). World Reference Base for Soil Resources.Rome: FAO.
61
XiaoH.WangB.LuS.ChenD.WuY.ZhuY.et al (2020). Soil acidification reduces the effects of short-term nutrient enrichment on plant and soil biota and their interactions in grasslands.Global Change Biol.264626–4637. 10.1111/gcb.15167
62
ZarzyckiJ.KopecM. (2020). The scheme of nutrient addition affects vegetation composition and plant species richness in different ways: results from a long-term grasslands experiment.Agric. Ecosyst. Environ.291:106789. 10.1016/j.agee.2019.106789
63
ZhanS.WangY.ZhuZ.LiW.BaiY. (2017). Nitrogen enrichment alters plant N: P stoichiometry and intensifies phosphorus limitation in a steppe ecosystem.Environ. Exp. Bot.13421–32. 10.1016/j.envexpbot.2016.10.014
64
ZhangC.WangL.-L.YangY.-P.DuanY.-W. (2015). Flower evolution of alpine forbs in the open top chambers (OTCs) from the Qinghai-Tibet Plateau.Sci. Rep.5:10254. 10.1038/srep10254
65
ZhangZ. H.WangG. S.WangH.QiQ.YangY. F.HeJ. S. (2021). Warming and drought increase but wetness reduces the net sink of CH4 in alpine meadow on the Tibetan Plateau.Appl. Soil Ecol.167:10406110.1016/j.apsoil.2021.104061
66
ZiH.HuL.WangC. (2022). Differentiate responses of soil microbial community and enzyme activities to nitrogen and phosphorus addition rates in an alpine meadow.Front. Plant Sci.13:829381. 10.3389/fpls.2022.829381
67
ZongN.ShiP. L. (2019). Enhanced community production rather than structure improvement under nitrogen and phosphorus addition in severely degraded alpine meadows.Sustainability11:2023. 10.3390/su11072023
Summary
Keywords
climate change, grassland ecosystems, biodiversity loss, nutrition addition, pollination network, Qinghai-Tibet Plateau
Citation
Wang L-L, Ren F, Zhang C, Huang X-J, Zhang Z-H, He J-S, Yang Y-P and Duan Y-W (2022) The effects of changes in flowering plant composition caused by nitrogen and phosphorus enrichment on plant–pollinator interactions in a Tibetan alpine grassland. Front. Plant Sci. 13:964109. doi: 10.3389/fpls.2022.964109
Received
08 June 2022
Accepted
05 July 2022
Published
25 July 2022
Volume
13 - 2022
Edited by
Janusz Zwiazek, University of Alberta, Canada
Reviewed by
Jinniu Wang, Chengdu Institute of Biology (CAS), China; Jianshuang Wu, Institute of Environment and Sustainable Development in Agriculture (CAAS), China
Updates
Copyright
© 2022 Wang, Ren, Zhang, Huang, Zhang, He, Yang and Duan.
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-Wen Duan, duanyw@mail.kib.ac.cnYong-Ping Yang, yangyp@mail.kib.ac.cnJin-Sheng He, jshe@pku.edu.cn
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.