Previous Article |
Table of Contents
| Next Article
Land Change Science Special Feature
LAND CHANGE SCIENCE SPECIAL FEATURE / BIOLOGICAL SCIENCES / BIOLOGICAL SCIENCES / RESEARCH ARTICLES / ECOLOGY / SUSTAINABILITY SCIENCE
Incorporating plant functional diversity effects in ecosystem service assessments





,
*Instituto Multidisciplinario de Biología Vegetal, Consejo Nacional de Investigaciones Científicas y Técnicas de Argentina, Universidad Nacional de Córdoba, Casilla de Correo 495, 5000 Córdoba, Argentina;
Laboratoire d'Ecologie Alpine, Unité Mixte de Recherche 5553, Centre National de la Recherche Scientifique, and Station Alpine Joseph Fourier, Unité Mixte de Services 2925, Centre National de la Recherche Scientifique, Université Joseph Fourier, BP 53, F-38042 Grenoble Cedex 09, France; and
Centro de Investigación Forestal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, Carretera de La Coruña km 7.5, Madrid 28040, Spain
Edited by B. L. Turner, II, Clark University, Worcester, MA, and approved September 17, 2007 (received for review May 2, 2007)
| Abstract |
|---|
|
|
|---|
biodiversity | land change | mass ratio hypothesis | plant functional traits
Conceptual models accounting for links between the functional trait values of local plant communities and EP are scattered in the literature, and the conceptual connections between them are not always clear. We articulate the most important of those models into a generic procedure and propose a small number of logical steps to reduce uncertainty in the prediction of EP and derived ES within the context of land cover change (Fig. 1). Reducing this uncertainty has important theoretical and applied implications. First, although there is increasing consensus on the fact that plant functional traits strongly affect EP and resulting ES (5, 6), very little is known about the relative role of different components of FD, such as the mean and frequency distribution of plant trait values (2, 7). Our generic method allows identification of cases in which EP can be satisfactorily predicted from different FD components and to quantify their relative importance. Second, the reduction of uncertainty will help identify the ES most vulnerable to biodiversity changes (1, 8, 9).
|
| FD at the Interface Between Global Change and EP |
|---|
|
|
|---|
| The Mass Ratio Hypothesis: A Cornerstone of FD–EP Relationships |
|---|
|
|
|---|
The mass ratio hypothesis is well supported by both theory and empirical evidence. There are conceptual models linking a relatively small number of plant (and especially leaf) traits with EP such as primary productivity, nutrient cycling, and trophic transfer to herbivores (3, 10, 14–20) that appear consistent across biomes, ecosystems, and floras (15–19). Links of local CWM of such plant traits with biogeochemistry-derived ES have been documented in several cases (21–25). After accounting for such effects, the distribution (e.g., FDvg) of plant trait values may play only a secondary role. Conversely, when the relationship between CWM and EP is poor, other components of FD, such as FDvg, might exert relatively strong effects.
| Results |
|---|
|
|
|---|
The procedure is organized in two stages comprising six steps. Stage I uses four steps to test sequentially for the effects of individual sets of factors. Within a given vegetation type, it starts by testing abiotic factors as drivers of EP (step 1). It then proceeds with community-level functional properties (CWM; step 2) where vegetation is simplified into an average trait value, strongly determined by the functional trait values of the more abundant species. In step 3, additional biotic complexity is taken into account through the degree of overlap in trait values within the community (FDvg; step 3). Finally, the procedure explores remaining idiosyncratic effects of particular species (step 4).
Although taken individually each step in stage I has already received considerable interest in the literature, we propose to go further and examine them in combination and to account for the relative effects of each component. Thus, stage II first combines all significant variables from steps 1–4 using a step-wise procedure (step 5). This bottom-up approach reflects the choice of being hypothesis-driven in our analysis. Finally, step 6 explores the possible discontinuous effects of FD or other variables, thereby accounting for remaining unexplained variance.
Stage I: Identifying Abiotic and Biotic Drivers of EP and ES. Step 1: Effects of abiotic factors. Abiotic variables often explain a significant amount of variation in EP (4, 19). This is especially the case when sampled land cover types span substantial variations in climate or soil, e.g., driven by topography, and/or when large spatial scales are considered. The first step in our procedure should thus check for a continuous relationship linking abiotic factors such as temperature, moisture, or soil chemical quality as independent variables against EP as a response variable (Fig. 1, step 1). For our approach, this could be seen as an initial null hypothesis, where biodiversity effects are of no detectable importance to the prediction of EP. There are many examples of abiotic effects on key EP. Water availability affects litter decomposition (14). Topography is an important determinant of soil stability (26). Primary productivity and nutrient cycling are strongly affected by soil fertility and/or water-holding capacity (27), abiotic factors that are in turn affected by land use (28). Many studies have also demonstrated strong legacies of past land use or disturbance regime on soils that carry over to current EP and ES (29, 30).
Step 2: Effects of community functional properties. In cases where the variations in EP cannot be satisfactorily explained by abiotic factors alone, one can assess whether FD exerts nonnegligible effects on EP. On that basis, we propose to start the process by testing for a continuous relationship based on the mass ratio hypothesis (Fig. 1, step 2). Model 2 (Fig. 1 legend) accounts for the relationship between EP and CWM for key plant functional traits. Traits are considered one at a time or may be combined to reflect trait syndromes. As a special case, the relative proportion of different growth forms (i.e., tussock grasses) can be used as CWM. Some examples of CWM strongly influencing EP and ES have been given in The Mass Ratio Hypothesis: A Cornerstone of FD–EP Relationships. Other examples include the often dramatic alteration of EP and ES through changes in CWM due to the expansion of invasive species with functional traits different from those of the native flora (31). CWM of single traits may be sufficient to explain the variance in EP as in specific leaf area and leaf tensile strength explaining specific productivity (22) and decomposition rate, respectively (14). However, it is common for multiple traits to affect a single EP (32). Here uncertainty in the prediction of EP from CWM is reduced by the inclusion of one or more traits (such as whole-plant and leaf traits or aboveground and root traits), but not by adding information on the distribution of trait values (next step).
Step 3: Effects of the distribution of trait values. In cases where abiotic factors or CWM of a single trait or multiple traits do not satisfactorily explain the variations in EP, the effects of the distribution of trait values (measured for example by FDvg) should be explored (Fig. 1, step 3). Examples of ecosystem effects attributed to differences in trait values among coexisting species include, for instance, enhanced biomass production and soil fertility by mixtures of species with and without nitrogen-fixing symbionts (23), with different leafing phenology (21), or with different rooting depths (33). However, most of these examples have been assessed as the effect of functional group richness rather than by quantifying FDvg. Another example is the alteration of litter decomposition and thus soil nutrient availability by litter mixing effects associated to the joint presence of fast-decomposing and slow-decomposing plants (34). Pollination of a particular species is influenced by the distribution of trait values in the surrounding community (9). Where the inclusion of FDvg significantly improves the predictive model (model 3), consideration of the distribution across species of values of one or more plant traits reduces uncertainty.
Step 4: Idiosyncratic effects of particular species. If EP cannot be satisfactorily explained by the weighted mean or distribution of the plant traits within land cover types, then it is possible that EP is associated in a significant and continuous way to the abundance of particular plant species (Fig. 1, step 4). Such effects could, among other reasons, be due to functional traits whose links with the EP of interest are not obvious from first principles or the literature and thus were not measured. For example, biomass production and nutrient retention in a grassland soil have been best explained by the colonization ability rather than by traits associated with nutrient use strategies (35). In the absence of data on colonization ability, the local abundance of certain species could assist in reducing uncertainty in predicting EP. Idiosyncratic effects could also be due to unknown species interactions with other plants, animals, or microorganisms (36). Idiosyncratic species effects can also be examined in an exploratory mode focusing on species known to strongly respond to specific global change drivers. Combined effects of several species may be tested after individual effects are demonstrated (32).
Stage II: Finding the Best Predictive Model. Step 5: Combining abiotic factors and FD components. Individual factors tested in steps 1–4 may provide significant predictions of EP, yet combining these models might further decrease the uncertainty in EP. It is therefore necessary to test various combinations of significant abiotic and biotic factors and to retain the final model that explains the largest amount of variance in the most parsimonious manner (step 5). This combination may lead to the rejection of some effects of FD components that are significant when tested individually (steps 2 and 3), i.e., because of collinearity with abiotic factor(s) affecting EP (step 1). In this case one can conclude that reducing the uncertainty in the prediction of EP does not strictly require additional information on plant traits. Conversely, FD effects may override abiotic effects on EP, e.g., when land use decisions lead to contrasting FD in otherwise homogeneous abiotic environments. Individual factors that do not significantly improve the explanatory power of the model can be considered as irrelevant for practical purposes.
There are only few examples in the literature where the effects of abiotic factors and several FD components have been explored simultaneously. As a notable exception, litter decomposition is often best explained by a combination of climatic variables and CWM of leaf traits such as nutrient concentration (14). As a result, the effects of nitrogen addition (e.g., by atmospheric deposition or through management) often depend on CWM of those leaf traits (37). Productivity is also the complex outcome of climate, management (e.g., fertilization), trait values of dominant species, and sometimes various effects of species richness or trait distribution (21, 23, 25). Most studies, however, especially the majority of those considering effects of species richness on productivity, have not conducted a systematic analysis of the different FD components.
Although our framework is flexible and can incorporate various procedures for model construction, we suggest building the model with those variables identified as significant in steps 1–4 using a step-wise process whereby complexity is progressively added to the model until a satisfactory level of predictive power is reached (model 5). The default option presented here implies a hierarchy of drivers of EP starting with abiotic factors, but reordering is possible depending on the particularities of the study at hand (e.g., spatial scale of interest).
Step 6: Discontinuous effect of CWM or trait distribution on EP. In all models discussed so far, we have assumed that single abiotic factors and/or FD components, or their combinations, explained EP. There can be cases, however, in which none of steps 1–5 yields a satisfactory model of EP. A plausible alternative explanation here could be that EP is strongly driven by a biotic or abiotic factor, but their relationship is discontinuous (model 6 in Fig. 1). For example, in some grasslands soil nitrogen losses depend in a discontinuous way on the C:N ratio of dominant grasses (38). Flammability of ecosystems also has a strongly discontinuous relationship with the abundance of flammable species (39). The presence of key species (step 4) can further result in discontinuous biotic effects on EP. For EP where continuous relationships prove not to be satisfactory, more intensive sampling of the gradient(s) or additional experimental approaches are needed to identify thresholds (T in model 6, Fig. 1 legend).
Case Study: FD and ES in Subalpine Grasslands. We used a previously published data set depicting changes in FD and EP in subalpine grasslands under different combinations of past and present land use (24, 40, 41). Key ES of importance to regional stakeholders, including farmers, other local residents, and tourists, were related to grassland EP (Table 1) (41). The data set was submitted to the six-step analysis to search for the best model explaining ES delivery from abiotic factors and FD components. Results on the significant factors identified via general linear models are presented in Table 1. Below we highlight three cases where different abiotic or FD factors explain variations in EP relevant to the ES identified in the land system.
|
Cases 2 and 3. Litter accumulation was identified by stakeholders as negatively affecting a range of ES including prevention of snow gliding, cultural heritage, aesthetic value, and provision of habitat for butterflies. In contrast to fodder production, abiotic factors did not influence EP associated with litter accumulation, which were best explained by CWM of leaf traits (leaf nitrogen and lignin contents; leaf tensile strength) and CWM of VH (in the case of total accumulated litter in the spring) (step 2). A marginally significant relationship was also found with the distribution of VH within communities (step 3), but for parsimony it was eliminated from the final model. Both for standing litter (the pool) and litter decomposability (the flux), FD components explained a very high proportion of the total variance. Similarly, CWM of leaf and root traits explained a very high proportion of the total variance in the maintenance of soil fertility (step 2). Hence, services associated with litter accumulation and soil fertility could be strongly related to average functional properties of the vegetation, providing support for the mass ratio hypothesis.
Case 4. Finally, sustained fodder production of green fodder through summer, which was associated with soil water content (SWC) in the summer period, illustrated the case of a complex model associating both abiotic and FD factors. SWC was explained by local soil conditions determining water-holding capacity and by an additional AGB term that accounted for direct biophysical effects of AGB on soil evaporation (water conservation) and plant transpiration (water demand). However, as shown in case 1, no significant model was found relating AGB to either local abiotic conditions or FD components. Effects of community-level means of leaf area (net positive effect representing reduction of soil evaporation) and root length (RL; negative effect representing water uptake) were significant and were included in the final model, with a very high proportion of the total variance explained. Hence, the prolonged delivery of green fodder depended on the combination of local abiotic factors, their effects on aboveground standing biomass, and average functional properties of the vegetation.
| Discussion |
|---|
|
|
|---|
Reducing Uncertainty in the Prediction of ES Delivery in Response to Environmental Change. Our case study illustrates the variety of ways in which EP underlying the delivery of ES can be related to abiotic factors and different components of FD. The quantitative models produced can also be applied to make projections of the future delivery of ES under contrasting land cover change scenarios (41). The links between CWM of plant traits and EP are strong and widespread, as shown by theory and accumulating empirical evidence under contrasting ecological and biogeographical contexts. In the specific example of the subalpine grasslands, when FD was involved, there was support for the mass ratio hypothesis, because community-level mean (CWM), rather than the distribution (quantified here by FDvg) of traits, ultimately explained EP, but other components of FD have been found to reduce uncertainty in predictions of EP (2, 8). Nevertheless, our case study also showed that the knowledge of FD does not always decrease uncertainty with respect to EP and that abiotic factors can sometimes be sufficient for practical purposes. Our objective here is therefore to highlight the fact that predicting ES from FD is a tractable yet complex task and that measuring several components of FD can be useful. In this context, the main strength of our method is to provide a generic procedure for selecting models and ranking FD components by their importance to local EP and ES in different land systems.
Although we have applied the proposed framework to plants for demonstration, we expect it to be applicable to other organisms and associated ES (e.g., insect diversity and pollination or crop protection; soil organism diversity and soil fertility). Some ES (e.g., pollination and nutrient cycling) will further require the combined consideration of FD across trophic levels (9, 32), short of which uncertainty in EP cannot be reduced from plant FD alone. This should make it possible to gain a more precise understanding of the extent to which biodiversity needs to be taken into account to predict changes in ES and human well-being under different scenarios.
It should be noted that, although it can rank different FD components in terms of explained variation in EP and ES, our method is not designed for the critical testing of different mechanisms underlying such FD effects. The mass ratio and idiosyncratic species effect models have been theoretically linked to the "selection" mechanism, whereas the model that focuses on FDvg, by stressing the presence of a variety of functional trait values, can be seen as linked to the "niche complementarity effect" (which also includes facilitation) (6, 7). Rather than providing a crucial test for these mechanisms, our method explores whether there are recurrent combinations of abiotic and FD factors that explain EP for different ES.
Practical Implications and Perspectives. To link land cover changes, whether through land use or other drivers, to changes in ES, our procedure requires data on abiotic and FD factors. For the latter, information is needed on both the traits and the relative abundance of species in the different land cover types under investigation. Recording species lists is often easier than measuring their relative abundance and trait values, yet our approach provides one more reason for compiling such databases on the basis of the rich yet scattered information available (2). Projections of EP and ES consequences of CWM shifts are currently well advanced (42), but projections for other factors, such as FDvg or discontinuous effects, are still in their infancy.
Adequate societal response to global change require quantitative assessments of the ES consequences of different land change alternatives (1). In this context, our approach offers an innovative way to develop quantitative predictive models of land cover change impacts on ES. By explicitly including FD effects on EP and ES in these models, our procedure also allows for a formal incorporation of biodiversity as a driving factor in the sensitivity of ES to global change.
| Methods |
|---|
|
|
|---|
EP. Biomass stocks, primary productivity (ANPP and SANPP, the logarithm of ANPP per unit of production time) (22), litter accumulation, and decomposition were quantified during the 2004 season following standardized protocols (22, 24). The leaf lignin:nitrogen ratio and the decomposability of litter (14) were determined for samples collected in July 2003. Components of the nitrogen cycle including total nitrogen pools, available nitrates (estimated with resin bags), and microbial nitrification and denitrification potentials were quantified throughout the 2004 and 2005 seasons (30). Finally, effects of vegetation on soil volumetric water content were assessed experimentally by a vegetation removal experiment in 2005 and repeated measurements through the growing season in vegetated versus unvegetated plots (40).
Quantification of FD. Plant FD within each grassland was quantified by using two metrics: CWM trait values calculated by using species' relative abundance (22) and FDvg quantified by using a modification of Rao's index that takes into account intraspecific trait variability (45). These two metrics were calculated for vegetative height (VH), specific leaf area, leaf dry matter content, LNC, leaf nitrogen to phosphorus ratio (N:P), leaf tensile strength, and flowering phenology. Root traits and aboveground lignin:nitrogen ratio were available only as aggregated values. Finally we combined specific leaf area, VH, and seed mass as a single FD index representative of species functional strategies (LHS, sensu ref. 46) by averaging FDvg across these three traits.
Statistical Analysis. We tested the relationships between EP and FD using general linear models (GenStat 9.0) in the six steps of our proposed procedure. For each of the six steps we proceeded preferentially based on expert knowledge about abiotic factors and traits relevant to the EP of interest.
Step 1. Effects of fertility (represented by NNI) and soil water-holding capacity were first tested individually, then in combination. Effects of disturbance date or intensity on EP were tested in ref. 24 but were not significant. We thus chose not to include them in this analysis.
Step 2. We used an ascending procedure (i) testing CWM for individual traits and then (ii) combining significant traits.
Step 3. We included FDvg of individual traits into previously selected models (significant terms only) for either traits already significant or additional traits.
Step 4. If no significant effects could be found in steps 2 and 3 we included the abundance of F. paniculata, the grass species that strongly dominates unmown grasslands (24). Previous analyses highlighted the effects of this species' abundance on EP and ES identified by local stakeholders (41).
Steps 1–4 were completed sequentially, then in step 5 all significant terms from steps 1–4 were combined in a step-wise ascending procedure. The most parsimonious model was selected by using the Akaike criterion.
Step 6. If none of the linear models from steps 1–4 (involving FD metrics, abiotic factors, and abundance of F. paniculata) were satisfactory, we reinspected data for possible discontinuities.
| ACKNOWLEDGMENTS. |
|---|
|
|
|---|
| Footnotes |
|---|
Abbreviations: AGB, aboveground biomass; ANPP, aboveground net primary productivity; SANPP, specific ANPP; CWM, community weighted mean value; EP, ecosystem property; ES, ecosystem service; FD, functional diversity; FDvg, functional divergence; LNC, leaf nitrogen content; NNI, nitrogen nutrition index; RL, root length; VH, vegetative height.
To whom correspondence should be addressed. E-mail: sdiaz{at}com.uncor.edu
Author contributions: S.D. and S.L. designed research; S.D., S.L., F.d.B., F.Q., K.G., and T.M.R. performed research; S.L. and K.G. analyzed data; and S.D., S.L., F.d.B., F.Q., K.G., and T.M.R. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
© 2007 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
Correction for Turner et al., Land Change Science Special Feature: The emergence of land change science for global environmental change and sustainability PNAS, February 19, 2008; 105(7): 2751 - 2751. [Full Text] [PDF] |
||||
![]() |
B. L. Turner II, E. F. Lambin, and A. Reenberg Land Change Science Special Feature: The emergence of land change science for global environmental change and sustainability PNAS, December 26, 2007; 104(52): 20666 - 20671. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||