Climate change and human activities have caused a shift in vegetation composition and soil biogeochemical cycles of alpine wetlands on the Tibetan Plateau. The primary goal of this study was to test for associations between community-weighted mean (CWM) traits, functional diversity, and soil properties during wetland drying. We collected soil samples and investigated the aboveground vegetation in swamp, swamp meadow, and typical meadow environments. Four CWM trait values (specific leaf area is SLA, leaf dry matter content is LDMC, leaf area is LA, and mature plant height is MPH) for 42 common species were measured across the three habitats; three components of functional diversity (functional richness, functional evenness, and functional divergence) were also quantified at these sites. Our results showed that the drying of the wetland dramatically altered plant community and soil properties. There was a significant correlation between CWM of traits and soil properties, but not a significant correlation between functional diversity and soil properties. Our results further showed that CWM-LA, CWM-SLA, and CWM-LDMC had positive correlations with soil readily available nutrients (available nitrogen, AN; available phosphorus, AP), but negative correlations with total soil nutrients (soil organic carbon is SOC, total nitrogen is TN, and total phosphorus is TP). Our study demonstrated that simple, quantitative plant functional traits, but not functional diversity, are directly related to soil C and N properties, and they likely play an important role in plant–soil interactions. Our results also suggest that functional identity of species may be more important than functional diversity in influencing ecosystem processes during wetland drying.
About one-third of China's natural wetlands are situated on the Tibetan Plateau, according to a state-of-the-art remote-sensing based classification (Gong et al., 2010). The wetlands of Maqu on the Tibetan Plateau are one of the highest wetlands, and they play a very important role in water resource conservation and in regulating the river systems and regional climates across the entire Yellow River basin (Cai and Guo, 2007). However, these wetlands are also climate-sensitive and have been gradually shrinking and drying due to climate change and human activities, such as heavy grazing and land use changes (Cai and Guo, 2007; Peng et al., 2015; Wang et al., 2016). Hence, it is imperative that we examine the changes in vegetation and soil and the relationships between them during this process of wetland drying for prediction and conservation efforts.
Previous studies have suggested that there may be a secondary successional change from swamp to swamp meadow and mature meadow and that soil water availability may be the main driver of this successional process (Ma et al., 2014; He and Richards, 2015; Lin et al., 2015). At the same time, plant functional traits play an important role in predicting the patterns of species composition, community structure, and their responses to environment change, which has drawn substantial ecological interest (e.g., Wright et al., 2004; McGill et al., 2006; Pérez et al., 2014; Li et al., 2015a); however, it is not clear whether functional traits and functional diversity (i.e., the identity, abundance, and range of species in a given community) can reflect the changes in soil properties during wetland drying on the Tibetan Plateau. A growing body of evidence has shown that functional diversity is directly linked with ecosystem processes (McGill et al., 2006; Flynn et al., 2011). However, recent studies have shown that environmental disturbance (e.g., fertilization or grazing) does not influence functional diversity (Li et al., 2015a; Niu et al., 2015), and functional traits of species may be more important than functional diversity in influencing ecosystem processes in an alpine meadow community (Li et al., 2015a).
Since the relationships between the plant community and environmental factors are dynamic, certain traits may be used to predict the changes in environmental factors (Lavorel and Garnier, 2002; Suding et al., 2008). At the species scale, some leaf traits that respond to water or nutrient availability (e.g., leaf dry matter content is LDMC and leaf nitrogen content is LNC) also affect leaf palatability and litter decomposability (Diaz et al., 2004). Similarly, at the community scale, the community-level means of some functional traits (e.g., specific leaf area (SLA), leaf area (LA), LDMC or plant height) weighted by their relative abundances can respond to environmental change in grasslands and affect nutrient cycling and dynamics (Garnier et al., 2004; Quétier et al., 2007).
Brief descriptions for vegetation properties (species richness,
vegetation cover, aboveground biomass, and dominant species), grazing
intensity, and standing water status during the wetland drying. Swamp is S, swamp meadow is SM, and typical meadow is TM. Significant differences across treatments
within each variable were determined using Tukey's honest significant
difference test (
The “mass ratio hypothesis” proposed by Grime (1998) implies that the functioning of ecosystems is determined to a large extent by the trait values of the dominant species. This hypothesis was confirmed by some studies that account for ecosystem functions such as primary productivity (Garnier et al., 2004) and nitrification (Laughlin, 2011). Conversely, the “diversity hypothesis” postulates that the trait diversity within a community can affect ecosystem processes (Tilman, 1997). The high functional diversity may allow for a more complete use of resources among species, thereby increasing biomass production and nutrient retention in ecosystems (Loreau, 2000; Díaz and Cabido, 2001; Villéger et al., 2008; Mouchet et al., 2010). The plant functional traits selected therefore pertain to the processing of resources at the species level, with the underlying assumption that this would scale up to ecosystem functioning (Lavorel and Garnier, 2002). We chose four functional traits, known to affect nutrient cycles at the leaf, whole-plant, and ecosystem levels (Cornelissen et al., 1999; Lavorel and Garnier, 2002): SLA (the ratio of water-saturated leaf area to leaf dry mass), LDMC (the ratio of leaf dry mass to water-saturated fresh mass), LA, and mature plant height (MPH). SLA is an important variable in comparative plant ecology because it is closely related to relative growth rate (Wright et al., 2001) and leaf net assimilation rate (Shipley and Lechowicz, 2000); it is also a good predictor of plant response to resource availability (Grime, 1977). LDMC is tied to plant nutrient retention and water (Poorter and Garnier, 1999). LA has important consequences for the leaf energy and water balance (Cornelissen et al., 2003). MPH has been associated with competitive ability in herbaceous plant communities, with the general idea that light competition becomes more intense at high soil fertility when aboveground biomass or vegetation stature increases (Hautier et al., 2009). Fast-growing species from nutrient-rich habitats usually have high SLA, high LNC, and low LDMC, while opposite trends characterize species from nutrient-poor habitats (Diaz et al., 2004). These responses reflect a fundamental trade-off (leaf economics spectrum) between traits related to nutrient conservation and traits related to nutrient acquisition and turnover (Wright et al., 2004).
The location of the Tibetan Plateau in China and our study site (the Research Station of Alpine Meadow and Wetland Ecosystems of Lanzhou University) at the Tibetan Plateau.
In this study, we investigate whether soil properties are affected by plant functional traits and functional diversity during wetland drying. Ongoing extensification of human activities has affected the functional composition and structure of wetlands on the Tibetan Plateau (Cai and Guo, 2007; Peng et al., 2015). At the same time, the soil nutrient availability has also been modified (Robson et al., 2007), and this has been directly linked to community-level trait changes in response to environmental dynamics (Diaz et al., 2007). The aim of this study was to answer the following questions: (1) How do community composition and structure and soil nutrient characteristics change following wetland drying on the Tibetan Plateau? (2) Can the community-weighted mean (CWM) trait values and functional diversity be regarded as predictors of soil properties during wetland drying?
This study was conducted at the Research Station of Alpine Meadow and
Wetland Ecosystems of Lanzhou University (33
In each
meadow type, 25 sampling plots (10 m
Soil samples (0–15 cm depth) were collected from each quadrat in each
sampling plot of each habitat after the aboveground material was
harvested. In each quadrat, soil was randomly collected from three points
(3.8 cm in diameter) using a bucket auger and mixed into a single soil
sample. All of the soil samples were brought into the laboratory in airtight
plastic bags. All of the soil samples were air-dried and then filtered
through a 0.2 mm sieve, discarding the visible roots and other plant debris.
Soil pH was measured using a pH meter with a glass electrode (soil / KCl
ratio is
Following Perez-Harguindeguy et al. (2013), we measured four functional
traits (SLA, LDMC, LA, and MPH) of 42 common species for three habitats in
early September 2010. These species represented 85–95 % of the peak
standing biomass and 80–90 % of the vegetation cover of the total plant
community in the studied plots. We randomly sampled one individual and three
mature leaves at flowering time for each of the 42 species in each 6 m
From the vegetation harvest data, we calculated the species richness (Pielou, 1969), represented
by the number of species recorded in each quadrat. Following Garnier et al. (2004), CWM trait values for each trait were
calculated for every sample using species mean trait values and species
relative cover: CWM
Although various indices have been proposed to measure the functional diversity of a community, there is still no consensus on which are most suitable. Villéger et al. (2008) suggested that some of the indices of functional diversity are redundant and they recommended using three independent components of functional diversity: functional richness (FRic), functional evenness (FEve), and functional divergence (FDiv). In this study, we chose FRic, FEve, and FDiv to examine how different components of functional diversity responded to different habitat types. We used the FDiversity software program to calculate FRic, FEve, and FDiv after the traits were standardized to ensure equal contributions from each trait.
We used one-way analysis of variance (ANOVA) to test the effect of different habitat types on plant taxonomic diversity (species richness), the CWM traits, functional diversity (FRic, FEve, and FDiv), and soil properties. Correlations between the CWM traits and soil properties were tested using Pearson's correlation coefficients, and a principal components analysis (PCA) on standardized data was conducted to analyze the overall pattern of correlations between soil properties, functional diversity, and aggregated trait values. These statistical analyses were performed using the R 3.2.3 software (R Development Core Team, 2011).
Relationships between CWM traits and soil properties from a Pearson's correlation analysis.
CWM is community-weighted mean trait values, SLA is leaf area per unit
dry mass, LDMC is leaf dry matter content, LA is leaf area, and MPH is mature plant
height.
Wetland drying had obvious effects on vegetation composition and structure,
and our study demonstrated that the species richness, vegetation cover, and
aboveground biomass significantly increased with wetland drying (Table 1).
There were also significant differences in species composition among the
community types (Table 1). Many typical wetland plants gradually disappeared
from the vegetation after wetland drying occurred; these included
Community-weighted means of traits differed greatly between the three habitats. The CWM-LA, CWM-SLA, and CWM-LDMC significantly increased with wetland drying (Fig. 2a, b, and c); however, the CWM-MPH of swamp (Fig. 2d) was significantly higher than swamp meadow and typical meadow. Functional richness (Fig. 3a), functional evenness (Fig. 3b), and functional divergence (Fig. 3c) had no obvious changes with wetland drying. These results also showed that community-level trait responses (community-weighted means of traits and functional diversity) were completely different, and the single-trait-level responses may be more sensitive than multi-trait-level responses to the process of wetland drying.
The changes in CWM trait across three habitats.
CWM is
community-weighted mean trait values, SLA is leaf area per unit dry mass,
LDMC is leaf dry matter content, LA is leaf area, MPH is mature plant height,
S is
swamp, SM is swamp meadow, and TM is typical meadow. Significant differences
indicated by dissimilar letters above each bar were determined using Tukey's
honest significant difference (HSD) test (
Wetland drying had obviously different effects on different soil characteristics. Soil organic carbon (SOC, Fig. 4a), soil total nitrogen (TN, Fig. 4b), soil total phosphorus (TP, Fig. 4c), and soil moisture (SM, Fig. 4f) of the swamp were significantly greater than swamp meadow and the typical meadow. However, soil available nitrogen (AN, Fig. 4d), soil available phosphorus (AP, Fig. 4e), and soil pH (Fig. 4g) of the swamp were significantly lower than those of swamp meadow and the typical meadow. These results clearly demonstrated that the changes in soil readily available nutrients (AN and AP) and total soil nutrients (SOC, TP, and TP) may be habitat-dependent, and the soil water content may be the dominant driving factor in the wetland drying.
The changes in functional diversity (functional richness, functional
evenness, and functional divergence) across three habitats. Significant
differences indicated by dissimilar letters above each bar were determined
using Tukey's honest significant difference (HSD) test (
Simple correlations and a PCA combining data on soil properties, community functional diversity, and the CWM traits were conducted. The analyses revealed a significant correlation between soil properties and CWM of traits but not a significant correlation between soil properties and functional diversity (Table 2). Fig. 5 further shows that there were significantly positive correlations between SOC, TN, TP, and SM and CWM-MPH but significantly negative correlations between AN, AP, and pH and CWM-MPH (Table 2, Fig. 5). In contrast, there was a significantly negative correlation between SOC, TN, TP, and SM and CWM-SLA, and CWM-LDMC and CWM-LA but significantly positive correlations between AN, AP, and pH and CWM-SLA, CWM-LDMC, and CWM-LA (Table 2, Fig. 5). In addition, the FRic had a positive relationship with FDiv but a negative relationship with FEve (Fig. 5). Thus, these results showed that community responses of single functional traits may play a crucial role in influencing soil properties during wetland drying.
Our study demonstrated that species richness, vegetation cover, and aboveground biomass significantly increased with wetland drying and that there were also significant differences in species composition among the community types. Many typical wetland plants disappeared from the community after wetland drying occurred, and many other new species (adapted to a high-drought environment) emerged from the typical meadow.
The changes in soil properties across three habitats. S is swamp,
SM is
swamp meadow, and TM is typical meadow. Significant differences indicated by
dissimilar letters above each bar were determined using Tukey's honest
significant difference (HSD) test (
Previous studies have shown that SLA and LA are often closely positively
correlated with photosynthesis and transpiration rate (Reich et al., 1999;
Westoby et al., 2002), and plants with an exploitative strategy (e.g., high
SLA and transpiration rate) are often drought intolerant (Reich et al.,
1999; Diaz et al., 2004). Our results also demonstrated that some forb
species (e.g.,
In contrast, functional diversity (FRic, FEve, and FDiv; Fig. 3a ,b, and c) did not show any significant changes among the three habitats (Fig. 3). We speculated that these changes in functional diversity would primarily be determined by the balance between competitive exclusion and stabilizing niche differences (Chesson, 2000; Hille Ris Lambers et al., 2012). Functional diversity may decrease when stabilizing niche differences of species are smaller than are needed to overcome competitive exclusion, and functional diversity may increase when stabilizing niche differences are greater than competitive exclusion. In addition, the responses of functional diversity can be constrained by close functional linkages to any number of traits determining performance at the whole plant level (Reich, 2014). Trade-offs in the responses within a suite of traits can create approximately equal performance at the whole plant level and hence a low value for change in functional diversity (Shipley et al., 2006). These results are consistent with many recent observations that the response of functional diversity to fertilization and grazing are relatively stable (Niu et al., 2015; Li et al., 2015a), and they highlight the potential complexity of community responses to environmental change.
Principal component analysis combining data on soil properties, CWM trait, and functional diversity. Only the first two axes (PC1 and PC2), which account for 80 % of the total variation, are retained here. The first axis can account for 71 % of the total inertia, and it differentiates communities according to soil moisture: the communities of higher soil moisture showed higher SOC, TN, TP, and CWM-MPH, while the communities of lower soil moisture showed higher AN, AP, pH, CWM-LA, CWM-SLA, and CWM-LDMC. The second axis can account for 9 % of the total inertia, and it showed that functional richness (FRic) had a positive relationship with functional divergence (FDiv) but a negative relationship with functional evenness (FEve).
Our results showed that the soil of the swamp had significantly higher moisture, organic carbon, total phosphorus, and total nitrogen but lower pH, available nitrogen, and available phosphorus than the swamp meadow and typical meadow. These results indicate that the soil total nutrient content had a decreasing trend with wetland drying. The wetlands can accumulate a large amount of soil organic matter because of lower decomposition rates of litter due to permanent or temporary anaerobiosis associated with waterlogged conditions (Nieder and Benbi, 2008). However, the drying of wetlands can stimulate microbial activity and increase mineralization of soil organic matter (Denef et al., 2001), therefore leading to an increase in soil readily available nutrients, such as nitrogen and phosphorus. However, other studies have reported that a reduction in mineralization with drying could be caused by reduced microbial activity and mobility (Pulleman and Tietema, 1999). In swamp meadow and typical meadow, certain factors can contribute to reduced soil total nutrient content. First, trampling by livestock (yak and Tibetan sheep) can lead to compaction and changes in infiltration rates, bulk density, and decline of edaphon activity (Li et al., 2011b; Yu and Jia, 2014). Second, under the long-term pressure of grazing, some energy and nutrients are transferred to livestock (Li et al., 2011b; Lu et al., 2015). In addition, the available forms of nitrogen and phosphorus in the wetland are lower than in the typical meadow because they are very water-soluble and move rapidly in the wetland environment.
Simple correlations and a PCA analysis strongly support the idea that
fast-growing plant species can promote the rates of soil C, N, and P
cycling, which in turn results in high nutrient availability but relatively
low C and N sequestration. A striking result of the present study is that
aboveground plant traits can determine the changes in belowground soil
properties during wetland drying. There are two likely reasons to explain
this phenomenon. First, it is possible that high-quality litter results in
high decomposition rates in suitable environmental conditions, whereas
lower-quality litter might result in a more consistent level of microbial
activity due to the longer residence time of the litter in soil. For the
swamp, the accumulation of SOC, TN, and TP was greater than in swamp meadow
and typical meadow because primary production exceeds the slow
mineralization rates caused by permanent or temporary anaerobiosis
associated with waterlogged conditions (Nieder and Benbi, 2008). Second,
other factors that covaried with leaf and litter quality may have been more
important drivers of soil properties. For example, the results of our study
support the fact that greater litter diversity associated with species diversity (in
typical meadow) can stimulate a large and active soil microbial community,
which further results in higher rates of litter decomposition and higher
soil nutrient availability (Wardle et al., 2006). In addition, in this
study, fast-growing species (with higher SLA and LA) that produce
high-quality litter may exceed slow-growing species (with higher LDMC), which
produce a litter of lower quality, therefore yielding the higher rates of
litter decomposition in swamp meadow and typical meadow. Some ecosystem
model studies such as BIOME-BGC (Running and Hunt, 1993) or CENTURY (Parton
et al., 1993) also suggest that some key functional traits, such as SLA,
LDMC, LA, and MPH, can directly affect ecosystem properties. Additionally,
effects of plant functional traits on soil properties may also result in
positive or negative feedback effects, which influence plant community
dynamics and composition (McLaren, 2006). In this study, higher soil
available nutrients changed plant community structure and productivity to
favor graminoid species, which have greater competitive and colonization
abilities (e.g.,
Our results showed that functional diversity had no measurable effect on soil properties. Many studies have found that functional diversity does not influence litter decomposition and nutrient cycling (Barantal et al., 2011; Schindler and Gessner, 2009). However, these studies typically use only one metric of functional diversity. Functional diversity is complex, and single metrics may not capture its effects (Mason et al., 2005). In this study, we used three metrics of functional diversity (FRic, FEve, and FDiv), yet we still did not find evidence for a relationship between functional diversity and soil properties. We found little evidence to support the diversity hypothesis, which would predict that trait diversity would lead to a greater efficiency of nutrient use. These results also suggested that functional diversity (traits taken in aggregate) may not be directly related to ecosystem functioning (Li et al., 2015a). Zuo et al. (2016) directly measured the relationship between plant functional traits and the C and N storage in plant, litter, root, and soil along a restoration gradient of sandy grassland in northern China, and their results showed that ecosystem C and N pools are primarily associated with the traits of the most abundant species in communities, thereby also supporting the biomass ratio hypothesis. Overall, these results suggest that soil properties are likely to be predictable based on plant functional traits and/or litter quality, but not functional diversity. However, predicting soil properties from plant functional traits is likely to be a complex process.
Our study has demonstrated that simple, quantitative plant functional traits, but not functional diversity, could provide relevant information on key aspects of ecosystem functioning. We have also shown that the particular traits selected here, specific leaf area, leaf dry matter content, leaf area, and mature plant height, could be used to capture the functioning of plant species and communities. Our observations also promote a better understanding of plant–soil feedbacks related to plant functional traits in plant communities, but long-term monitoring is needed in future research, especially in different times of the year, because the cycling and sequestration of soil nutrients may be a “fast-out, slow-in” process with high variation due to regulation by environmental factors (Schulze, 2006).
The data are not publicly available due to copyright issues.
The authors declare that they have no conflict of interest.
This study was supported by the Natural Science Foundation of China (41601586, 41671289, 41230852), Natural Science Foundation of Shaanxi Province (2016JQ3008), State's Key Project of Research and Development Plan (2016YFC050704), Northwest A & F University (QN2013069), and the Key Cultivation Project of the Chinese Academy of Sciences – the promotion and management of ecosystem functions of restored vegetation on the Loess Plateau, China. Edited by: M. Oliva Reviewed by: M. Ma and one anonymous referee