Document Detail

Higher levels of multiple ecosystem services are found in forests with more tree species.
Jump to Full Text
MedLine Citation:
PMID:  23299890     Owner:  NLM     Status:  MEDLINE    
Forests are of major importance to human society, contributing several crucial ecosystem services. Biodiversity is suggested to positively influence multiple services but evidence from natural systems at scales relevant to management is scarce. Here, across a scale of 400,000 km(2), we report that tree species richness in production forests shows positive to positively hump-shaped relationships with multiple ecosystem services. These include production of tree biomass, soil carbon storage, berry production and game production potential. For example, biomass production was approximately 50% greater with five than with one tree species. In addition, we show positive relationships between tree species richness and proxies for other biodiversity components. Importantly, no single tree species was able to promote all services, and some services were negatively correlated to each other. Management of production forests will therefore benefit from considering multiple tree species to sustain the full range of benefits that the society obtains from forests.
Lars Gamfeldt; Tord Snäll; Robert Bagchi; Micael Jonsson; Lena Gustafsson; Petter Kjellander; María C Ruiz-Jaen; Mats Fröberg; Johan Stendahl; Christopher D Philipson; Grzegorz Mikusiński; Erik Andersson; Bertil Westerlund; Henrik Andrén; Fredrik Moberg; Jon Moen; Jan Bengtsson
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't    
Journal Detail:
Title:  Nature communications     Volume:  4     ISSN:  2041-1723     ISO Abbreviation:  Nat Commun     Publication Date:  2013  
Date Detail:
Created Date:  2013-01-09     Completed Date:  2013-06-20     Revised Date:  2013-07-11    
Medline Journal Info:
Nlm Unique ID:  101528555     Medline TA:  Nat Commun     Country:  England    
Other Details:
Languages:  eng     Pagination:  1340     Citation Subset:  IM    
Department of Ecology, Swedish University of Agricultural Sciences, Box 7044, SE-75007 Uppsala, Sweden.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Models, Biological
Species Specificity
Trees / physiology*

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Nat Commun
Journal ID (iso-abbrev): Nat Commun
ISSN: 2041-1723
Publisher: Nature Pub. Group
Article Information
Download PDF
Copyright © 2013, Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.
Received Day: 07 Month: 06 Year: 2012
Accepted Day: 26 Month: 11 Year: 2012
Electronic publication date: Day: 08 Month: 01 Year: 2013
Volume: 4First Page: 1340 Last Page:
PubMed Id: 23299890
ID: 3562447
Publisher Item Identifier: ncomms2328
DOI: 10.1038/ncomms2328

Higher levels of multiple ecosystem services are found in forests with more tree species
Lars Gamfeldta12
Tord Snäll1
Robert Bagchi3
Micael Jonsson4
Lena Gustafsson1
Petter Kjellander5
María C. Ruiz-Jaen6
Mats Fröberg78
Johan Stendahl8
Christopher D. Philipson9
Grzegorz Mikusiński5
Erik Andersson1011
Bertil Westerlund12
Henrik Andrén5
Fredrik Moberg11
Jon Moen4
Jan Bengtsson1
1Department of Ecology, Swedish University of Agricultural Sciences, Box 7044, SE-75007 Uppsala, Sweden
2Department of Biological and Environmental Sciences, University of Gothenburg, Box 461, SE-40530 Gothenburg, Sweden
3Department of Biological and Biomedical Sciences, Durham University, South Road, Durham DH1 3LE, UK
4Department of Ecology and Environmental Science, Umeå University, SE-90187 Umeå, Sweden
5Grimsö Wildlife Research Station, Department of Ecology, Swedish University of Agricultural Sciences, SE-730 91 Riddarhyttan, Sweden
6Environmental Change Institute, South Parks Road, Oxford OX1 3QY, UK
7Department of Aquatic Sciences and Assessment, Swedish University of Agricultural Sciences, Box 7082, SE-75007 Uppsala, Sweden
8Department of Soil and Environment, Swedish University of Agricultural Sciences, Box 7014, SE-75007 Uppsala, Sweden
9Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Winterhurerstrasse 190, 8057 Zurich, Switzerland
10Southern Swedish Forest Research Centre, Swedish University of Agricultural Sciences, Box 49, SE-23053 Alnarp, Sweden
11Stockholm Resilience Centre, Stockholm University, SE-10691 Stockholm, Sweden
12Department of Forest Resource Management, Swedish University of Agricultural Sciences, Box 7001, SE-90183 Umeå, Sweden

Ecosystems provide a wealth of benefits to human society, and the provision of such ecosystem services depends fundamentally on functions performed by organisms1. This has led scientists to enquire how the diversity and composition of communities may regulate ecosystem functions2, 3. A large body of evidence has established that species diversity promotes ecosystem functions under experimental conditions4, 5. There are, however, many exceptions to the positive diversity–function relationship6. In addition, most experiments have been conducted at limited temporal or spatial scales (but see refs 7,8). It is thus uncertain if conclusions based on results from these studies can be extended to the scales relevant to policy makers.

Thus far, studies that manipulated biodiversity and monitored the consequences mainly focused on ecosystem functions, but because functions have complex links to ecosystem services1, 9, the policy implications of these studies are unclear5. There are relatively few large-scale observational studies explicitly aimed at quantifying the importance of biodiversity, and these have mostly analysed single-ecosystem functions or services10, 11, 12. Furthermore, while there are multiple functions that regulate services from ecosystems, few studies have investigated the role of biodiversity for multiple ecosystem functions jointly (but see, refs 13,14,15,16), and none of these have focused on services per se. In a recent comprehensive review, a majority of the included ecosystem services were related to biodiversity in the direction expected from predictions3. However, for many of the studied services, the evidence for beneficial effects of biodiversity was mixed, or there were not enough data for a thorough evaluation3. Studies that explicitly investigate the link between biodiversity and multiple ecosystem services across regional scales are thus urgently needed if research is to be informative for management and policy, such as the Intergovernmental Platform on Biodiversity and Ecosystem Services17.

We examined the relationships between multiple ecosystem services and both tree species richness and tree biomass in boreal and temperate forest. These extensive biomes are of high global importance1, constituting around 27% (8.7 million km2) and 16% (5.2 million km2) of the world's forests, respectively18. The majority of studies on the importance of tree diversity have focused on productivity, and there is evidence for both non-significant19 and positive effects20. We studied production forests, that is, forests that are subjected to a silvicultural system of harvesting and planting or natural regeneration from seed trees. Production forest is the dominant system in the sampled region, which covered an area of 400,000 km2 and a span of 13.7 degrees of latitude. These forests are relatively species poor. The maximum number of tree species in any one sampling plot in our study was 10 (Supplementary Fig. S1), and only 1.5% of the plots hosted more than five species. Six ecosystem services1 were included: tree biomass production, topsoil carbon storage, berry production, game production potential, understory plant species richness and dead wood occurrence. Production of tree biomass is a major global industry, contributing about US$ 400 billion to the annual gross world product in terms of wood1. Forests and their soils store about 45% of the terrestrial carbon, and act as a crucial sink for anthropogenic carbon emissions21. Both production of berries and game are of large monetary and recreational value in many countries1. Dead wood occurrence represents a supporting function for other ecosystem services, as it affects ecological processes such as nutrient cycling and soil formation1, 22, and it also indicates several other components of biodiversity. Understory plant species richness differs somewhat from the other five ecosystem services, as it is not linked to biological processes or biological production, but instead can be classified as a cultural ecosystem service that provides aesthetic, educational and recreational value to human beings1. Alternatively, and with recent terminology, understory plant species richness may be viewed as having a value of its own23. Safeguarding biodiversity at all levels of organization is also a central goal of the Convention on Biological Diversity.

Relationships between tree species richness and the services

The relationships between tree species richness and tree biomass production, berry production (of bilberries) and game production potential were positively hump-shaped (Fig. 1a), taking into account important environmental variables, such as climate, soil nutrients and forest age (Methods). On average, biomass production at the mean forest age was 54% greater with five tree species than with only one species (Fig. 1a). Corresponding percentages for berry production and game production potential were 45% and 20%, respectively. Soil carbon storage and understory plant species richness increased with tree species richness (Fig. 1b). Soil carbon storage was 11% greater, and understory plant species richness 31% greater, with five than with one tree species. The occurrence of dead wood increased at high tree species richness (Fig. 1f). An alternative way of evaluating how productivity changes with tree species richness is to quantify the probability that a given species richness would have higher productivity than a monoculture. In our analyses, tree biomass production is 41% more likely to be higher on plots with five species than on plots with one species (Supplementary Fig. S2 for details, additional statistics and corresponding percentages for the other five services). Excluding the effect of the biomass of individual tree species from the models did not change the positive relationships between the services and tree species (Supplementary Fig. S3).

Individual tree species and services

None of the six ecosystem services was consistently explained by the biomass of the same tree species (Fig. 2). For example, there were strong positive relationships between the biomass of spruce and tree biomass production and occurrence of dead wood, while pine biomass was positively related to bilberry production, and birch biomass to soil carbon storage. Beech was not positively related to any service.

Relationships between ecosystem services

Some of the services were positively related to each other (Fig. 3), but we also found some notable trade-offs. First, there were trade-offs between tree biomass production and occurrence of dead wood, game production potential and bilberry production. Second, dead wood occurrence and game production potential were negatively correlated.


The pros and cons of species mixtures for productivity and other ecosystem functions have been discussed at length since the early 19th century2, 24, 25. Not until recently, however, have scientists begun to explicitly investigate how species diversity might be important for the simultaneous provision of multiple functions or services13, 14, 15, 26, 27. Our results from boreal and temperate production forests show that the relationships between tree species richness and multiple ecosystem services were positive to positively hump-shaped, and that all services attained higher levels with five tree species than with one species. Although the relatively high level of tree biomass production with five compared with one tree species may seem both impressive and surprising, we note that similar effect sizes have been found previously11. Another observational study, in which climatic and environmental conditions were controlled for, also found a strong and positive effect of tree diversity on productivity12. Because controlled field studies on mixed-stand effects have revealed both non-significant and positive effects of diversity19, 20, future research should aim to understand and reconcile results from controlled experiments and observational studies. There are, unfortunately, no previous studies on the relationships between tree species richness and the other five services that allow for comparisons with our findings.

We found little support for a general importance of the biomass of particular species in explaining all services (Fig. 2). The finding that different species were most strongly related to different services indicates that monoculture practices will lead to reduced provision of at least some of the services we considered. For example, among all tree species considered, only pine was strongly and positively related to bilberry production (see also ref. 28) and birch showed the strongest positive relationship with soil carbon storage. As such, our results suggest the importance of tree species mixtures for the continued provisioning of ecosystem services from the 2-billion hectares of forest in the world currently managed as production forests or used for multiple purposes (55% of all forests)29.

The results also indicate trade-offs between some services (Fig. 3). The trade-off between tree biomass production and dead wood occurrence is probably an effect of stand age; tree production generally decreases with age, whereas the opposite is true for occurrence of dead wood30. The negative correlation between dead wood occurrence and game production potential is most likely also partly explained by forest stand age; herbivore forage plants are more common in young forests, whereas dead wood more commonly occurs in old forests. Negative relationships between tree biomass production and the production of both bilberry and food for game are likely consequences of light limitation and lower precipitation through-fall at the ground level in closed productive forests31. Unlike previous findings32, we thus found that some ecosystem services may come at the cost of others. Identifying such trade-offs is highly relevant for policy, enabling users of forest ecosystems to understand and balance the pros and cons of different management scenarios33. Although the trade-offs we found imply that it will be difficult to maximize all ecosystem services simultaneously at the stand scale, the positive relationships between tree diversity and individual services (Fig. 1) suggest that adjacently located monocultures would not optimize the provision of ecosystem services at the landscape scale. Instead, adjacent stands, each with multiple species but in different combinations, might be the best way to provide multiple ecosystem services at the landscape scale.

The negligible effects of the biomass of single tree species on the relationship between tree species richness and ecosystem services in our study suggest complementarity among tree species, for example, resource partitioning, often referred to as ‘complementarity effects' (or effects of species richness per se) in experimental studies5. Interactions between plants and associated microbiota, and higher resistance of a diverse community to disease, insects or pathogens may also be factors explaining the observed relationships between some of the services and tree species richness34.

Additional environmental conditions may directly or indirectly explain some of the unexplained variation in the services34, 35. For example, nitrogen availability (not measured in the Swedish National Forest Inventory) may influence forest productivity. However, the ratio of total carbon to total nitrogen in the topsoil that we used is a good indicator of nitrogen availability36. Also, non-measured environmental variables that are correlated with altitude or latitude could be important; for example, plots in the mountain region (in northwest Sweden) generally have fewer tree species (Supplementary Fig. S1). However, we accounted for temperature and humidity, which are the most important climate variables and correlated with altitude and latitude. Moreover, maximum tree species richness was relatively evenly distributed across the sampling region (Supplementary Fig. S1).

While our analyses were on the regional scale, there certainly may be differences at the local scale in the strength of the relationships between tree species richness and ecosystem services. All services may thus not be maximized similarly within landscapes across the whole sample region. Furthermore, we focused on the current distribution and composition of forests, but future changes in the environment may affect ecosystem services both directly and indirectly37. Global change may thus result in society not having all the management options currently available. Scenarios of changes in forest species composition and richness due to global change must be considered in analyses of future provisioning of ecosystem services38.

The majority of previous efforts to disentangle the effects of tree species richness on ecosystem functions or services in forests have focused on productivity of woody biomass11, 12, 39, 40. Mixed forest stands can indeed be more productive than monocultures20, and several studies have found large-scale positive relationships between tree diversity and wood production11, 12, 41. Positive relationships between biodiversity and individual ecosystem functions or services have also been observed in aquatic realms, such as the open ocean42, the deep sea43, coral reefs44 and lakes45. Our study, which accounts for important environmental conditions, shows consistent positive relationships between tree species richness (contrasting plots with five and one tree species) and multiple ecosystem services. It also highlights the importance of conserving a variation of tree species, to safeguard a future potential of high levels of multiple ecosystem services. As such, our research expands on previous work16, 46, none of which considers multiple services. Our study revealed complementarity among species; no single species can sustain multiple services at high levels simultaneously. This is particularly pertinent at a time when demands on ecosystems to provide humans with several important services are increasing47 at the same time as global change and human expansion is threatening the integrity of ecosystems, their biodiversity and functioning1, 3, 44, 48. Our results strongly support the integration of production and conservation in natural resource management, and align with the ‘ecosystem approach', the primary framework for action under the Convention of Biodiversity49. We show that moving towards multi-species management can better realize the full potential of several economically, ecologically and culturally valuable ecosystem services.

Description of the data

We used a nationwide forest data set from the Swedish National Forest Inventory and the Swedish Survey of Forest Soils and Vegetation. The inventory uses a randomly planned regular sampling grid50, and includes around 4,500 permanent tracts with each tract being surveyed once every 5 years. The tracts, which are rectangular in shape and are of different dimensions in different parts of the country, consists of 8 (in the north) to 4 (in the south) circular sample plots (each plot 314 m2). From the Inventory database, we extracted data from 1999 to 2002. We used only plots on ‘productive forest' (average production of standing volume, stem volume over bark >1 m3 ha−1 year−1), which had not been harvested, cleared or thinned during the previous 5 years before the survey. The previous time period of 5 years was also used for calculating biomass production. To be included in our analyses, the plots had to be located on only forest land, for example, not including any river, road, grassland and so on. We excluded plots where the biomass production of any of the focal tree species was greater than the 99% quantile (which is the standard procedure in analyses of data from the Swedish National Forest Inventory). This procedure excluded plots with unrealistically high values of production. Our selection resulted in a data set with a total of 4,335 sample plots distributed across 1,401 tracts. We included explanatory variables hypothesized to explain the variation in the six modelled ecosystem services, and squared variables and interactions, whenever biologically reasonable (see descriptions below, Fig. 4 and Supplementary Table S1).

Statistical modelling

To model the six ecosystem services as a function of tree species richness, we used hierarchical Bayesian-generalized linear models51. The Bayesian approach is convenient for fitting complex hierarchical models with a formal mathematical treatment of the natural variability and data uncertainty, for example, when the recording of the explanatory variables is incomplete.

For each ecosystem service, we fitted the following models: the main model, on which our results are based unless otherwise mentioned (Figs 1 and 2); a model aimed to investigate how the tree species biomasses affected the relationship between tree species richness and the services, specifically the main model but excluding the tree biomass variables (Supplementary Fig S3); five models aimed to investigate the pairwise relationships between the ecosystem services, excluding all other explanatory variables (Fig. 3). In the following, we exemplify the modelling using understory plant species richness as the response variable. Corresponding models were fitted for the other five ecosystem services.

We assumed that the understory plant species richness followed an overdispersed Poisson distribution with mean, μp,t (Supplementary Table S2). We modelled this mean, transformed by the ln link function (Supplementary Table S2), g(), on plot p=1, …, P⩽8, on tract t=1, …, T⩽1,394 as

[Formula ID: eq1]

where, Xp,t is a design matrix of n plot-level explanatory variables, βn is a vector of associated effect-size parameters, and Pp,t is the regional species pool (see description of understory plant species richness). Explanatory variables were centred and standardized. As we included the log of the regional species pool, we in fact modelled the proportion of the regional species pool. The parameters ɛp,t are overdispersion contributions, where ɛp,tN(0, σ), which means normally distributed with mean=0 and s.d.=σ. For normally distributed ecosystem services (Supplementary Table S2), ɛp,t are residual contributions. The tract-level ‘intercept' parameters αt were modelled further as

[Formula ID: eq2]

where σα is the s.d., and

[Formula ID: eq3]

where Zt is a design matrix of m=1, …, M⩽5 tract-level explanatory variables, which were centred and standardized, ρm is a vector of associated effect-size parameters, and γ is an ‘intercept' parameter.

We built the main model based on the Deviance Information Criterion (DIC52), on the posterior distribution of the effect-size parameters (βn in equation (1) and ρm in equation (3)), and on knowledge of the biological system studied. We first assessed the predictive power of each explanatory variable (Supplementary Table S1) based on the DIC and on the posterior distributions of βn and ρm. Next, we fitted a multiple model containing the retained explanatory variables. Finally, we simplified this multiple model by excluding and again including earlier excluded variables in a stepwise procedure. Parameter estimates for the final main model for each ecosystem service are presented in Fig. 4 and in Supplementary Table S3. For some ecosystem services, there were a few models with similar DIC (difference in DIC<2). However, the relationships between the response variables and (i) tree species richness, and (ii) the tree species biomasses differed only marginally between these models.

For a description of parameters, symbols and prior distributions, see Supplementary Table S4. We assumed uninformative prior distributions for all parameters. The models were fitted using the softwares OpenBUGS 2.1 (ref. 53). For the models aimed at investigating the pairwise relationships between services, we excluded all plots with missing data. Moreover, bilberry production was also arcsine transformed when used as an explanatory variable. Below follows a detailed description of our response and explanatory variables.

Ecosystem services

Tree biomass production was estimated as the yearly change in tree biomass (kg m−2 year−1), calculated over a period of 5 years, and for all tree individuals higher than 1.3 m. For plots visited in the years 1999–2002, the baseline for biomass production was thus the years 1994–1997. Biomass was calculated with biomass functions54, 55 and was the sum of the biomass from the stem, twigs and branches, the stump and roots. For deciduous tree species, there is only a function for Betula spp., and this function was thus applied for all other deciduous tree species. The function for Pinus sylvestris was applied for Larix decidua and Pinus contorta. Even though this creates a slight tree species bias, it has minor effects on our production estimates, as we compared the biomass of the trees at two points in time. Sample size for tree production was 4,335 plots, and due to fallen trees, observed values could be negative (as was the case in 5 plots).

Soil carbon storage was measured as the amount of carbon (g m−2) in the topsoil, which consisted of either purely organic horizons, for example, mor layers (63%) or peat layers (21%), or less frequently of minerogenic A-horizons (16%). This is the part of the soil most affected by the current above-ground biota. To compensate for the conceptual difference in topsoil types, mean soil carbon stocks were set equal for purely organic soils (measured in the whole organic horizon down to a maximum depth of 30 cm) and minerogenic topsoils (measured in the top 10-cm horizon). The soil fraction <2 mm was analysed. Soil sampling was carried out on around 50% of the inventory plots, totalling 1,953 plots. Data were log transformed.

Bilberry production was measured as the percentage of each plot covered by bilberry, Vaccinium myrtillus. The cover of V. myrtillus is strongly correlated to annual production, and bilberry is one of the economically most important wild berry species in Northern Europe28. Bilberry sampling was carried out on around 50% of the inventory plots, totalling 2,127 plots. Data were converted to proportions and arcsine transformed, after correcting for the area where berries could not grow, for example, the area of stems and boulders.

Game production potential was calculated as the percentage cover of each plot that was covered by any of the 15 plant species important for large herbivorous mammals, for example, moose (Alces alces) and roe deer (Capreolus capreolus), see Supplementary Table S5 and ref. 56. The cover of these field layer species was measured in the same way as the bilberry cover. Not all plants are equally preferred, however. Based on known preferences and nutritional values57, we therefore weighted the least and most important plant species by 0.66 and 1, respectively (Supplementary Table S5). As the herbivores also browse on coppice, we included the cover of preferred species in this layer (Supplementary Table S5). Field layer and coppice were summed to yield a measure of total food for herbivores, which could amount to more than 100% because of multiple vegetation layers. Game production is the basis for both recreational hunting and meat, and was estimated to a value of US$ 460 million in Sweden as of 2006 (ref. 58). Sample size was 2,127 plots and data were log transformed after adding a constant of 0.01.

Understory plant species richness: The Swedish National Forest Inventory includes about 270 ground vegetation species, and of these, we selected those with forest as main habitat, in total 141 species (Supplementary Table S6). We modelled the proportion of the richness of these species in each plot in relation to the regional species richness of the same species, estimated as all of the 141 species occurring within a 600 × 600 km window centred on each plot. The regional species pool thus captured the change in species composition with latitude (as the longitudinal distance in the sampling region never exceeds 600 km). We found that the median number of plant species in the regional pool reached an asymptote at a distance of ∼300 km from each plot, meaning that further increasing the area of the regional species pool did not, on average, incorporate more species. Sample size was 4,327 plots.

We modelled the probability of dead wood occurrence using data on presence/absence of dead wood in each plot (minimum diameter for dead wood pieces to be counted equals 40 mm). We chose this variable, as the abundance of dead wood showed a highly skewed distribution in which small fragments were common and a few plots had very high abundances. Sample size was 4,335 plots. Dead wood is the essential habitat for many organisms: more than 50% of the red-listed forest species in Sweden depend on dead wood59.

The six ecosystem services included in our analyses are a subset of all possible services that forests provide. The results and conclusions should be viewed with this in mind.

Explanatory variables

Temperature is defined as the total accumulated average daily temperature (in °C) over 5 °C during the vegetation period. Humidity is measured as the difference in mm between the addition of water via precipitation and the loss of water due to transpiration and evaporation during the vegetation period. For nitrogen deposition, spatially distributed N deposition data in kg N ha−1 year−1 (wet and dry) were downloaded from the Swedish Meteorological and Hydrological Institute website. As temperature, humidity and nitrogen deposition are ‘modelled' variables, their resolution is not at the scale of local plots, and were thus modelled at the scale of tracts. pH of the topsoil was measured in water suspension. The ratio of total carbon to total nitrogen is the ratio (C:N) in the topsoil. Soil moisture was originally measured on a categorical scale of 1 (driest) to 5 (wettest). These categories correspond to a continuous scale of % volumetric soil water contents ranging from 0.15 to 0.35 (ref. 60). Peat was modelled as peat (1) or non-peat (0) soil. The biomass of the six most abundant or regionally important (‘focal') tree species in Sweden were included. These were: Norway spruce (Picea abies), Scots pine (Pinus sylvestris), birch (Betula spp.; B. pendula and B. pubescens), oak (Quercus robur), aspen (Populus tremula) and beech (Fagus sylvatica), measured in kg m−2. Tree species richness was our estimate of tree species diversity. It was calculated as the observed presence of 20 species in each plot, and ranged from 1 to 10. The tree species included our six focal species together with Acer platanoides, Acer pseudoplatanus, Alnus glutinosa, Alnus incana, Carpinus betulus, Fraxinus excelsior, Larix decidua, Pinus contorta, Prunus avium, Salix caprea, Sorbus aucuparia, Sorbus intermedia, Tilia cordata and Ulmus glabra. We chose species richness as our diversity estimate because it is a simple and intuitive measure. Species richness estimates are relatively easy to obtain and relate to the ability to manage ecosystems under limited information about species' relative distributions and traits. Forest stand age was estimated in each plot as the basal area weighted mean age from at least three trees.

Author contributions

Study idea by L.G., J.M. and J.B.; L.G., T.S., M.J., Le.G., P.K., J.M. and J.B. designed the research; L.G., T.S., R.B., Le.G., P.K., M.F., J.S., B.W., C.D.P. and J.M. extracted and prepared the data; L.G., T.S., R.B., M.C.R.-J., C.D.P., J.M. and J.B. analysed the data; All authors interpreted the results and wrote the paper.

Additional information

How to cite this article: Gamfeldt, L. et al. Higher levels of multiple ecosystem services are found in forests with more tree species. Nat. Commun. 4:1340 doi: 10.1038/ncomms2328 (2013).

Supplementary Material Supplementary information

Supplementary Figures S1-S3 and Supplementary Tables S1-S6

This work was conducted as part of the Future Forests research program in Sweden. L.G. was supported by grant 621-2009-5457 from the Swedish research council VR, and T.S. by grant 2005-933 from Formas. We thank David Wardle and Andy Hector for comments that helped improve an earlier version of the paper. Johan Mentzer helped drafting Fig. 2.

Millenium Ecosystem Assessment. Ecosystems and Human Well-Being: Biodiversity Synthesis World Resources Institute. 100 (Island Press (Year: 2005).
Naeem S.,, Bunker D. E.,, Hector A.,, Loreau M.. , Perrings C. H.,Biodiversity, Ecosystem Functioning, & Human Wellbeing, an Ecological and Economic Perspective Oxford University Press, (Year: 2009).
Cardinale B. J.,et al. . Biodiversity loss and its impact on humanity. Nature486, 59–67 (Year: 2012).22678280
Schmid B.,et al. in Biodiversity, Ecosystem Functioning, & Human Wellbeing: An Ecological and Economic Perspective eds Naeem S. Oxford University Press (Year: 2009).
Cardinale B. J.,et al. . The functional role of producer diversity in ecosystems. Am. J. Bot.98, 572–592 (Year: 2011).21613148
Wardle D. A.,, Jonsson M.,Biodiversity effects in real ecosystems—a response to Duffy. Front. Ecol. Environ.8, 10–11 (Year: 2010).
Wardle D. A.,, Zackrisson O.,Effects of species and functional group loss on island ecosystem properties. Nature435, 806–810 (Year: 2005).15944702
Reich P. B.,et al. . Impacts of biodiversity loss escalate through time as redundancy fades. Science336, 589–592 (Year: 2012).22556253
Naeem S.,, Bunker D. E.,, Hector A.,, Loreau M.,, Perrings C. H.. in Biodiversity, Ecosystem Functioning, & Human Wellbeing : An Ecological and Economic Perspective (eds Naeem S Ch. 1, 3–13 Oxford University Press (Year: 2009).
Grace J. B.,et al. . Does species diversity limit productivity in natural grassland communities?Ecol. Lett.10, 680–689 (Year: 2007).17594423
Vila M.,et al. . Species richness and wood production: a positive association in Mediterranean forests. Ecol. Lett.10, 241–250 (Year: 2007).17305807
Paquette A.,, Messier C.,The effect of biodiversity on tree productivity: from temperate to boreal forests. Global Ecol. Biogeogr.20, 170–180 (Year: 2011).
Hector A.,, Bagchi R.,Biodiversity and ecosystem multifunctionality. Nature448, 188–190 (Year: 2007).17625564
Gamfeldt L.,, Hillebrand H.,, Jonsson P. R.,Multiple functions increase the importance of biodiversity for overall ecosystem functioning. Ecology89, 1223–1231 (Year: 2008).18543617
Zavaleta E. S.,, Pasari J. R.,, Hulvey K. B.,, Tilman G. D.,Sustaining multiple ecosystem functions in grassland communities requires higher biodiversity. Proc. Natl Acad. Sci. USA107, 1443–1446 (Year: 2010).20080690
Maestre F. T.,et al. . Plant species richness and ecosystem multifunctionality in global drylands. Science335, 214–218 (Year: 2012).22246775
Larigauderie A.,, Mooney H. A.,The Intergovernmental science-policy Platform on Biodiversity and Ecosystem Services: moving a step closer to an IPCC-like mechanism for biodiversity. Curr. Opin. Environ. Sustain.2, 9–14 (Year: 2010).
Hansen M. C.,, Stehman S. V.,, Potapov P. V.,Quantification of global gross forest cover loss. Proc. Natl Acad. Sci. USA107, 8650–8655 (Year: 2010).20421467
Frivold L. H.,, Frank J.,Growth of mixed birch-coniferous stands in relation to pure coniferous stands at similar sites in south-eastern Norway. Scand. J. Forest Res.17, 139–149 (Year: 2002).
Pretzsch H.,, Schutze G.,Transgressive overyielding in mixed compared with pure stands of Norway spruce and European beech in Central Europe: evidence on stand level and explanation on individual tree level. Eur. J. Forest Res.128, 183–204 (Year: 2009).
Bonan G. B.,Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science320, 1444–1449 (Year: 2008).18556546
Janisch J. E.,, Harmon M. E.,Successional changes in live and dead wood carbon stores: implications for net ecosystem productivity. Tree Physiol.22, 77–89 (Year: 2002).11830405
Mace G. M.,, Norris K.,, Fitter A. H.,Biodiversity and ecosystem services: a multilayered relationship. Trends Ecol. Evol.27, 19–26 (Year: 2012).21943703
Cotta H.v.,Anweisung zum Waldbau Arnoldische Buchhandlung (Year: 1828).
Pretzsch H., in Forest Diversity and Function Vol. Ecological Studies eds Scherer-Lorenzen M., Körner C., Schulze E. D. 41–64Springer-Verlag (Year: 2005).
Hooper D. U.,, Vitousek P. M.,Effects of plant composition and diversity on nutrient cycling. Ecol. Monogr.68, 121–149 (Year: 1998).
Duffy J. E.,, Richardson J. P.,, Canuel E. A.,Grazer diversity effects on ecosystem functioning in seagrass beds. Ecol. Lett.6, 637–645 (Year: 2003).
Miina J.,, Hotanen J. P.,, Salo K.,Modelling the abundance and temporal variation in the production of bilberry (Vaccinium myrtillus L.) in Finnish mineral soil forests. Silva. Fenn.43, 577–593 (Year: 2009).
Food and Agricultural Organization of the United Nations. in Global Forest Resources Assessment 2010: Main Report FAO Forestry Paper Rome (Year: 2010).
Swedish Forest Agency. Statistical Yearbook of Forestry 2012: Official Statistics of Sweden Jönköping (Year: 2012).
Anderson R. C.,, Loucks O. L.,, Swain A. M.,Herbaceous response to canopy cover light intensity and throughfall precipitation in coniferous forests. Ecology50, 255–263 (Year: 1969).
Nelson E.,et al. . Modeling multiple ecosystem services, biodiversity conservation, commodity production, and tradeoffs at landscape scales. Front. Ecol. Environ.7, 4–11 (Year: 2009).
Raudsepp-Hearne C.,, Peterson G. D.,, Bennett E. M.,Ecosystem service bundles for analyzing tradeoffs in diverse landscapes. Proc. Natl Acad. Sci. USA107, 5242–5247 (Year: 2010).20194739
Schnitzer S. A.,et al. . Soil microbes drive the classic plant diversity-productivity pattern. Ecology92, 296–303 (Year: 2010).21618909
Godbold J. A.,, Solan M.,Relative importance of biodiversity and the abiotic environment in mediating an ecosystem process. Mar. Ecol. Prog. Ser.396, 273–282 (Year: 2009).
Dise N. B.,, Rothwell J. J.,, Gauci V.,, van der Salm C.,, de Vries W.,Predicting dissolved inorganic nitrogen leaching in European forests using two independent databases. Sci. Total Environ.407, 1798–1808 (Year: 2009).19185335
Bulling M. T.,et al. . Marine biodiversity-ecosystem functions under uncertain environmental futures. Phil. Trans. R. Soc. B365, 2107–2116 (Year: 2010).20513718
Hickler T.,et al. . Projecting the future distribution of European potential natural vegetation zones with a generalized, tree species-based dynamic vegetation model. Global Ecol. Biogeogr.21, 50–63 (Year: 2012).
Scherer-Lorenzen M.,, Körner C.,, Schulze E. D.,Forest Diversity and Function Vol. 176, (Springer-Verlag (Year: 2005).
Morin X.,, Fahse L.,, Scherer-Lorenzen M.,, Bugmann H.,Tree species richness promotes productivity in temperate forests through strong complementarity between species. Ecol. Lett.14, 1211–1219 (Year: 2011).21955682
Caspersen J. P.,, Pacala S. W.,Successional diversity and forest ecosystem function. Ecol. Res.16, 895–903 (Year: 2001).
Worm B.,et al. . Impacts of biodiversity loss on ocean ecosystem services. Science314, 787–790 (Year: 2006).17082450
Danovaro R.,et al. . Exponential decline of deep-sea ecosystem functioning linked to benthic biodiversity loss. Curr. Biol.18, 1–8 (Year: 2008).18164201
Mora C.,et al. . Global human footprint on the linkage between biodiversity and ecosystem functioning in reef fishes. PLoS Biol.9, e1000606 (Year: 2011).21483714
Ptacnik R.,et al. . Diversity predicts stability and resource use efficiency in natural phytoplankton communities. Proc. Natl Acad. Sci. USA105, 5134–5138 (Year: 2008).18375765
Isbell F.,et al. . High plant diversity is needed to maintain ecosystem services. Nature477, 199–202 (Year: 2011).21832994
Foley J. A.,et al. . Solutions for a cultivated planet. Nature478, 337–342 (Year: 2011).21993620
Hooper D. U.,et al. . A global synthesis reveals biodiversity loss as a major driver of ecosystem change. Nature486, 105–108 (Year: 2012).22678289
Shepherd G.,The Ecosystem Approach: Learning from Experience IUCN (Year: 2008).
Axelsson A.-L.,et al. in National Forest Inventories—Pathways for Common Reporting eds Tomppo E., Gschwantner Th., Lawrence M., McRoberts R. E. 541–553Springer (Year: 2010).
Gelman A.,, Carlin J. B.,, Stern H. S.,, Rubin D. B.,Bayesian Data Analysis 2nd edn CRC Press (Year: 2004).
Spiegelhalter D. J.,, Best N. G.,, Carlin B. R.,, van der Linde A.,Bayesian measures of model complexity and fit. J. R. Stat. Soc. Ser. B Stat. Methodol.64, 583–616 (Year: 2002).
Thomas A.,, O'Hara R. B.,, Ligges U.,, Sturtz S.,Making BUGS open. R News6, 12–17 (Year: 2006).
Marklund L. G.,Biomass Functions for Pine, Spruce and Birch in Sweden Department of Forest Survey, SLU (Year: 1988).
Petersson H.,, Ståhl G.,Functions for Below-Ground Biomass of Pinus Sylvestris, Picea Abies, Betula Pendula and Betula Pubescens in Sweden Department of Forest Resource Management and Geomatics, SLU (Year: 2006).
Cederlund G.,, Ljungqvist H.,, Markgren G.,, Stålfelt F.,Foods of moose and roe deer at Grimsö in central Sweden—results of rumen content analysis. Swed. Wildlife Res.11, 167–247 (Year: 1980).
Hanley T. A.,A nutritional view of understanding and complexity in the problem of diet selection by deer (Cervidae). Oikos79, 209–218 (Year: 1997).
Boman M.,, Mattsson L.,, Ericsson G. R.,, Kriström B.,Moose hunting values in Sweden now and two decades ago: the Swedish hunters revisited. Environ. Resource Econ.50, 515–530 (Year: 2011).
National Atlas of Sweden. The Royal Swedish Academy of Agriculture and Forestry. Agriculture and Forestry in Sweden Since 1900 – a Cartographic Description232 (Stockholm (Year: 2011).
Warfvinge P.,, Sverdrup H.,Critical Loads of Acidity to Swedish Forest Soils. Methods, Data and Results104 (Lund (Year: 1995).

Article Categories:
  • Article

Previous Document:  Probing the electronic structure at semiconductor surfaces using charge transport in nanomembranes.
Next Document:  Kinship reduces alloparental care in cooperative cichlids where helpers pay-to-stay.