Biomass-dominant species shape the productivity-diversity relationship in two temperate forests

Abstract

Key message

A negative productivity-diversity relationship was determined for biomass-dominant species at the community level. This study thus supports the hypothesis in which the effects of individual species on the productivity-diversity relationships at the community level are related to their biomass density, an important functional trait.

Context

The productivity-diversity relationships have been extensively studied in various forest ecosystems, but key mechanisms underlying the productivity-diversity relationships still remain controversial.

Aims

The objective of this study is to explore the productivity-diversity relationships at the community level, and to investigate the roles of individual species in shaping the community-level relationships between productivity and diversity under different forest types.

Methods

The study was conducted in two fully stem-mapped temperate mixed forest plots in Northeastern China: a natural secondary forest plot, and an old-growth forest plot. An individual-based study framework was used to estimate the productivity-diversity relationships at both species and community levels. A homogeneous Thomas point process was used to evaluate the significance of productivity-diversity relationship deviating from the neutral.

Results

At the species level, most of the studied species exhibit neutral productivity-diversity relationship in both forest plots. The percentage of species showing negative productivity-diversity relationship approaches linearly a peak value for very close neighborhoods (the secondary forest plot: r = 3 m, 38%; the old-growth forest plot: r = 4 m, 42%), and then decreases gradually with increasing spatial scale. Interestingly, only a few species displayed positive productivity-diversity relationship within their neighborhoods. Dominant species mainly exhibit negative productivity-diversity relationship while tree species with lower importance values exhibit neutral productivity-diversity relationship in both forests. At the community level, a consistent pattern of productivity-diversity relationship was observed in both forests, where tree productivity is significantly negatively associated with local species richness. Four biomass-dominant species (Juglans mandshurica Maxim., Acer mono Maxim.,Ulmus macrocarpa Hance and Acer mandshuricum Maxim.) determined a negative productivity-diversity relationship at the community level in the secondary forest plot, but only one species (Juglans mandshurica) in the old-growth forest plot.

Conclusion

The productivity-diversity relationship is closely related to the dominance of individual species at the species level. Moreover, this analysis is the first to report the roles of biomass-dominant species in shaping the productivity-diversity relationship at the community level.

1 Introduction

The influence of species diversity on ecosystem functioning, especially on biomass productivity, has become a major focus of ecological research (Loreau et al. 2001; Hooper et al. 2005). There are no universal relationships between productivity and diversity (hereinafter referred to as PDR). Positive relationships between productivity and diversity appear to be rather common (Schwartz et al. 2000; Hooper et al. 2005; Srivastava and Vellend 2005), whereas U-shaped, negative, and non-significant are also observed (Wardle et al. 1997; Wardle and Lagerström 2008; Ricklefs 1977; Tilman and Pacala 1993; Grace 1999; Mittelbach et al. 2001; Gillman and Wright 2006; Adler et al. 2011; Šímová et al. 2013; Gadow et al. 2016). The relationships between productivity and diversity are usually scale dependent. They are sometimes positive at the global/regional scale, but are often neutral, unimodal, or negative at the local spatial scales (Harrison et al. 2006). The mechanisms shaping the PDR have also received much attention. For example, species diversity may enhance productivity because of facilitation effects and/or niche complementarity among species, and because of the greater probability of including productive species in species-rich communities (Hector 1998; Tilman et al. 1997; Loreau and Hector 2001). But the ecological mechanisms which cause or contribute to a negative PDR remain controversial.

Species diversity affects ecosystem functioning (Tilman 1999; Hooper and Vitousek 1997; Symstad et al. 1998; Lohbeck et al. 2012), but there is no compelling evidence to show that ecosystem functions are consistently and causally connected with species diversity (Grime 1997). Dominant species were found to exhibit an asymmetrical influence on ecosystem functioning (Wardle 1999; Smith and Knapp 2003; Vile et al. 2006). For example, productivity is considered to be driven by a few dominant species that make up a large fraction of the plant biomass (Grime 1997). Ecosystem functioning is determined by functional traits of dominant species such as their relative biomass (Grime 1998; Wardle et al. 2000; Roscher et al. 2012). Numerous studies have demonstrated that ecosystem processes (Lepš et al. 1982) and ecosystem productivity (MacGillivray and Grime 1995; Tilman et al. 1997) are primarily driven by the functional traits of the dominant plants but are independent of species diversity.

The roles of dominant species in shaping PDR are often neglected. The presence of a dominant species may increase the primary production of an ecosystem, because it has higher ecological resource-use efficiency than subdominant individuals (Gong et al. 2011). Dominant species accounting for a relatively small percentage of all species may contribute considerably to total ecosystem functioning. Sala et al. (1996) have shown that the removal of the dominant species, which accounted for 40% of total biomass, had a greater effect on ecosystem functioning than the removal of a constant amount of biomass from each individual species in equal proportions. If most of a particular ecosystem functioning is performed by a few species, and when these species are also dominant, then it may be possible to obtain a negative relationship between ecosystem productivity and species richness (Creed et al. 2009).

Huston et al. (2000) found that there was no consistent effect of species diversity on productivity. Thus, it is necessary to establish fresh approaches to understand the multivariate links between productivity and diversity. On the basis of this need, an individual-based PDR framework was developed to explore relationships between tree productivity and local diversity (Wang et al. 2016). Limited seed dispersal determines the distributions of trees in a forest (Seidler and Plotkin 2006). Conspecifics occurring in distinct patches usually show positive spatial autocorrelation, which has been thought to be a fundamental ecological process influencing local diversity and ecosystem functioning (Pacala and Deutschman 1995; Kühn 2007). A number of theoretical studies have supported the hypothesis that species diversity is dispersal limited (Turnbull et al. 2000; Loreau et al. 2003; Cadotte 2006). The influence of dispersal limitations substantially alters the functions of natural ecosystems (Mouquet and Loreau 2002; Stein et al. 2008). Dispersal limitations contribute to a unimodal PDR not only in a theoretical study (Pärtel and Zobel 2007) but also in natural ecosystems (Stein et al. 2008).However, the spatial autocorrelation of species distributions has been largely neglected in many previous studies. In order to simulate the effects of spatial autocorrelation of species distribution on local diversity and ecosystem functions, a Thomas cluster process (Shen et al. 2009) should be incorporated into the individual-based PDR framework.

This study examines the structuring mechanisms of the community-level PDR in a natural secondary forest and an old-growth forest. We assume that the effects of individual species on the community-level PDR are closely related with an important functional trait, i.e., biomass density. Furthermore, we expect that community-level PDR is consistent in the two forests, but that the shaping roles of individual species should be different. Accordingly, the objective of this study is to explore the pattern of community-level PDR and to investigate the roles of individual species in shaping the community-level PDR in the two forest communities.

2 Materials and methods

2.1 Data collection

The observations for this study were assessed in two large permanent field plots in Northeastern China. The forest type of the study area is a typical temperate conifer and broadleaved mixed forest, in which the average rainfall is 606 mm/year, the average monthly temperature is − 18.6 °C in January and 21.7 °C in July. More details about forest distributions and climate in the Jilin area have been reported by Li et al. (2001). Two experimental field plots were established in the summer of 2010 and remeasured in the summer of 2015: a secondary forest plot and an old-growth forest plot. In each plot, all woody plants with a diameter at breast height (dbh) exceeding 1 cm were tagged and mapped; their species were identified and their dbhs, tree heights, crown widths, and heights to the first live branch were recorded. The two plots have similar woody plant species (Supplementary material Appendix 1 Table A1).

The observational study in the secondary forest plot covers an area of 21.84 ha (420 m × 520 m), and is located at 43° 58.207′ -43° 58.558′ N and 127° 44.092′ -127° 44.541′ E. The plot represents a natural secondary forest. The original vegetation in the secondary forest plot is the same as that in the old-growth forest plot, which emerged after heavy harvesting disturbance several decades ago. The topography in the secondary forest plot is flat with elevations ranging from 468 to 505 m above sea level. In terms of the importance value, the secondary forest plot is dominated by Juglans mandshurica, Acer mandshuricum, Acer mono, Fraxinus mandschurica Rupr., Ulmus davidiana Plabch. var. japonica (Rehd.) Nakai, Ulmus macrocarpa and Betula costata Trautv. (Supplementary material Appendix 2 Fig. A1). However, the top eight species in biomass density are Juglans mandshurica, Fraxinus mandschurica, Pinus koraiensis Sieb. et Zucc., Acer mono, Ulmus macrocarpa, Ulmus davidiana var. japonica, Acer mandshuricum, and Ulmus laciniata (Trautv.) Mayr. (Fig. 1a).These species with higher biomass density are hereafter referred to as biomass-dominant species. Our analysis in the secondary forest plot includes 24 woody species.

Fig. 1

Species rankings ordered by (a) biomass density in 2010 and (b) the biomass productivity of individual species. Biomass productivity is the annual increment of aboveground woody biomass from 2010 to 2015. Abbreviated species names in the figure are shown in the Supplementary material Appendix 1 Table A1

The old-growth forest plot covers an observational study area of 30.0 ha (500 m × 600 m), which is located at 43° 57.928′ -43° 58.214′ N and 127° 45.287′ -127° 45.790′ E. The old-growth forest plot is far away from villages where human disturbance has been virtually unknown. The elevation of the old-growth forest plot ranges from 576 to 784 m above sea level and the topography involves a valley between two slopes. The top ranking species regarding the importance value are Acer mono, Ulmus laciniata, Carpinus cordata, Acer mandshuricum, Tilia amurensis, Acer barbinerve, and Pinus koraiensis (Supplementary material Appendix 2 Fig. A1). The top eight species in biomass density are Ulmus laciniata, Acer mono, Tilia amurensis Rupr., Pinus koraiensis, Betula costata, Juglans mandshurica, Phellodendron amurense Rupr., and Quercus mongolica Fisch. ex Ledeb (Fig. 1a). Again, our analysis in the old-growth forest plot includes 24 woody species.

Although the secondary forest plot has a species composition that is similar to the old-growth forest plot, the dominance, abundance, and productivity of specific species have changed after human disturbance. For example, the importance value of Fraxinus mandschurica in the secondary forest plot is greater than that in the old-growth forest (Supplementary material Appendix 2 Fig. A1).

2.2 Species-level relationships between tree productivity and diversity

We evaluated the species-level PDR through the individual-based PDR framework (Wang et al. 2016). The individual-based PDR framework represents a long-recognized tradeoff between competition and facilitation among species. It provides a powerful tool to identify the role of individual species in shaping the PDR. For the combination of tree productivity and species richness surrounding the trees of a focal species, the index of association is calculated as follows:

$$ \mathrm{AI}(r)=\frac{1}{N}\sum \limits_{i=1}^N\sqrt{S{R}_i\times {P}_i} $$

(1)

where AI(r) is the index of association; N is the number of individuals. SR_i is the species richness within a circular area with varying radius r around the ith tree; ten circular areas whose radius increase from 1 to 10 m by 1 m (r = 1, 2…10) are considered for each target tree; P_i is biomass productivity (in kg/year) of the ith tree. To eliminate the dimensional effect of the variable units, SR_i and P_i were standardized by the Min-Max normalization. When using a = πr², the AI(r) function can also be expressed in terms of a circular area a. To avoid the circular area of the focal species being located outside study area, we selected a buffer zone of 10 m along the border of the study area. Only individuals of the focal species within the inner area were used to calculate the P_i values. However, all individuals within the entire area were used to calculate the SR_i values.

There is a growing consensus that ecosystem processes are determined by functional traits of individual species rather than by the number of species (Chapin III et al. 2000; Díaz and Cabido 2001). By considering species richness, other important components of biodiversity are underestimated. Therefore, it is worth trying to calculate the index of association by including a functional diversity index instead of species richness to gain better insights on forest functioning in the future. However, in this study, we still follow the original framework of Wang et al. (2016) due to a lack of particular plant traits.

2.3 Community-level relationships between tree productivity and diversity

The individual-based PDR framework was also used to estimate the PDR at the community level. All trees within the entire research plot were regarded as focal individuals regardless of the species identity. The biomass productivity of each tree (P_i) and species richness surrounding each tree (SR_i) were included in Eq. (1) to calculate the index of association.

2.4 Estimating tree productivity

Tree biomass productivity, specifically annual aboveground biomass increments between 2010 and 2015, was calculated using the following equation:

$$ {P}_i=\frac{B{15}_i-B{10}_i}{5} $$

(2)

where P_i is biomass productivity of live tree i between the measurements in 2010 and 2015; B10 and B15 are the biomass estimates of tree i in 2010 and 2015, respectively. Biomass productivity was calculated using live trees at both measuring times, while recruits and dead trees were excluded in the calculation. Dry biomass of each individual tree was estimated using dbh and aboveground biomass equations, which were described in detail in a previous study (Wang et al. 2016).

2.5 The significance of PDR deviating from neutral

Ripley’s method was used to examine the spatial pattern of species distribution (Ripley 1976). The results show that all studied species exhibited a clumped distribution in the two forests (Supplementary material Appendix 2 Fig. A2 and Fig. A3). We assumed that the spatial autocorrelation of species distribution determine local diversity and tree growth. Thus, a Thomas cluster model was performed to examine the significance of PDR deviating from neutral by using AI as a test statistic in the individual-based PDR framework. The Thomas cluster process is a classic Cox process. In the Thomas cluster process, the parents are generated by a homogeneous Poisson process with intensity k and each parent location gives rise to a random number of offspring according to a Poisson distribution with mean μ offspring per parent (Diggle 1978; Thomas 1949).For each parent location y_i, the offspring points x_ij are independently and identically distributed, with a probability density of offspring locations estimated by an isotropic Gaussian density (Baddeley et al. 2016).

The confidence envelope was computed using 100 simulated distributions based on the homogeneous Thomas point process. If the observed AI is greater/less than that of the fifth largest/smallest value of the 100 realizations, then the observed AI is regarded to differ significantly from a neutral pattern at the 5% level (Fig.2). To satisfy reasonable sample sizes, the analysis was conducted only for the most abundant species comprising at least 30 individuals in each plot.

Fig. 2

Conceptual schematic diagram of individual-based PDR framework. The solid line is the index of association which is calculated from the observed dataset; the dashed lines are the confidence envelopes which can be calculated from the Monte Carlo simulation. (a) For a positive PDR, the observed line is located above the upper limit of the envelope; (b) For a neutral PDR, the observed line is within the envelope; (c) For a negative PDR, the observed line is below the lower limit of the envelope

2.6 The roles of individual species

In order to identify the roles of biomass-dominant species, an individual species or species group are randomly removed from the computation process of the community-level PDR in the sequence of the rankings shown in Fig.1a. Thus, it was possible to identify the relative roles of these tree species in shaping the community-level PDR by comparing the emerging patterns of the community-level PDR.

3 Results

3.1 Species-level relationships between tree productivity and diversity

At the individual species level, the departure of observed PDR from neutrals was tested using a homogeneous Thomas process in the two forests (Supplementary material Appendix 2 Fig. A4 and Fig. A5). The analysis revealed that most species exhibit a neutral PDR in both plots. The percentage of species showing a negative PDR approach a peak value linearly at very close neighborhoods (the secondary forest plot: r = 3 m, 38%; the old-growth forest plot: r = 4 m, 42%), and then decreases gradually with increasing scale.

Interestingly, very few species displayed a positive PDR within their neighborhoods (Fig. 3). Tree species with higher importance values (i.e., dominant species) mainly exhibit negative PDRs while those with lower importance values tend to exhibit neutral PDRs within the two forests (Supplementary material Appendix 2 Fig. A1).

Fig. 3

Proportion of species showing significant positive or negative PDR in the secondary forest plot (a) and the old-growth forest plot (b). The significant departure from a neutral PDR for focal species was evaluated using a homogeneous Thomas process

3.2 The effects of specific species on community-level PDR

The individual-based PDR framework was used to estimate the relationships between productivity and diversity at the community level. Our results show that there was a consistent community-level PDR pattern in the two forests, and that tree productivity is significantly negatively associated with tree species richness (Fig. 4).

Fig. 4

Relationships between tree productivity and species richness for all living trees greater than 5 cm in dbh within the secondary forest plot (a) and the old-growth forest plot (b). The dashed lines are the confidence intervals calculated by the homogeneous Thomas process. The solid lines are the observed AI values. The map in each figure shows the distribution density of all living trees greater than 5 cm in dbh

The rank-biomass distribution curves reflect the relative contribution of individual species to ecosystem functioning (Fig. 1a). The rank-biomass curve is much steeper in the secondary forest plot than that in the old-growth forest plot. The first seven species account for 89% of total biomass in the secondary forest plot. However, in the old-growth forest plot, the relative regular rank-biomass curve suggests that most dominant species contribute similarly to the total biomass.

At the community level, the PDR shifts from negative to neutral after removing the biomass-dominant species (Figs. 5 and 6). To identify individual species or species groups which may shape the community-level PDR, biomass-dominant species were sequentially removed during the computation process of the community-level PDR. Our results show that four species (Juglans mandshurica, Acer mono, Ulmus macrocarpa and Acer mandshuricum) jointly determine the negative community-level PDR in the secondary forest plot, but it was only one species (Juglans mandshurica) in the old-growth forest plot (Supplementary material Appendix 1Table A2 and Table A3).

Fig. 5

The effects of biomass-dominant species on community-level PDR in the secondary forest plot. An individual species or a group of individual species was removed during the computation process of the community-level PDR according to the species biomass density rank ordering

Fig. 6

The effects of biomass-dominant species on community-level PDR in the old-growth forest plot. An individual species or a group of individual species was removed during the computation process of the community-level PDR according to the species biomass density rank ordering

4 Discussion

The species-level PDR shows significant scale dependence in both study areas. The average spacing between trees, calculated by the square root of 10,000/N (where N is number of trees per ha; Nilson 2006), with a dbh greater than 5 cm is 3.4 m in the secondary forest plot and 3.9 m in the old-growth forest plot. Interestingly, up to a circle radius corresponding to that average spacing the percentage of species showing negative PDR achieves the maximum (Fig.3). The average spacing between trees is usually closely related to competition among trees. Competitive interactions will enlarge the average spacing between trees when weaker trees close to others are eliminated (Meyer et al. 2008). However, the percentages of species showing negative PDR are less than 30% in the majority of the studied scales in both forests. According to the results of the removal analysis, four dominant species showing negative PDR jointly determine the negative community-level PDR in the secondary forest plot, but only one productive species showed a neutral PDR in the old-growth forest plot. This result indicates that competitive effects are the fundamental driving forces influencing PDR’s patterns in the secondary forest, but neutral effects in the old-growth forest. This observation is consistent with the fact that the secondary forest plot is still characterized by plant competitive interactions, while the old-growth forest plot is a relative mature and stable forest community. It is worth noting that the two forest plots present different topographies and different altitudinal ranges. Climatic and other abiotic factors are likely to affect productivity in different ways. Therefore, it is hard to conclude the differences in PDR at the community level are completely resulted from the successional stages. The two large forest plots provide a wealth of evidence for assessing the PDR at the community level. However, due to the limited number of plots, the results may be considered exploratory rather than confirmatory.

The species-level PDR and the species dominance are intimately related. Dominant species mainly exhibit a negative PDR while tree species with lower importance values tend to be neutral PDR in both forests. The dominant species have more access to essential resources such as light, water, and soil nutrients. In addition, they usually exhibit a greater interspecific competitive ability (Sala et al. 1996), and may exert greater pressure on other species. Interspecific competition is common in natural communities (Connell 1983; Schoener 1983), and dominant species could negatively affect species richness by competitively excluding other species (Kunte 2008; Tilman 1980).The biomass productivity is mainly generated by a few dominant species in our study (Fig. 1b).

Our study revealed a consistent pattern of negative community-level PDR in the two forests which may further increase the superiority of a few dominant species (Grime 1973; Laughlin and Moore 2009). Dominant species play a decisive role in regulating species diversity and ecosystem productivity. For example, numerically dominant species in a forest community can decrease species diversity by influencing the establishment of new species (Gilbert et al. 2009; Crutsinger et al. 2010). Removing biomass from many species in proportion to their abundance has a smaller effect than removing biomass from only the dominant ones (Sala et al. 1996). We found that the community-level PDR largely depends on a particular subset of biomass-dominant species. In the secondary forest plot, four species presenting higher importance values (Juglans mandshurica, Acer mono, Ulmus macrocarpa, and Acer mandshuricum) jointly determine the negative community-level PDR, but only a single species with a relatively lower importance value (Juglans mandshurica) in the old-growth forest plot. Therefore, Juglans mandshurica seems to be a key species regarding its critical role in shaping the community-level PDR.

The findings of this study are not only interesting to ecologists, but also to foresters who manage multi-species natural forests. For example, structure-based forest management (Li et al. 2014) and near-natural forest management (Gadow et al. 2016; Puettmann et al. 2015) have been proposed to maintain species richness as well as forest productivity. These management practices focus on the cultivation and protection of target trees, but the dominance of tree species is often underestimated. In the future, forest management should take account of the importance of dominant species in shaping the relationships between productivity and diversity in a multi-species uneven-aged forest (Zhang et al. 2014). Similar studies are required to identify the specific competitive role of dominant species (Seifert et al. 2014) and thus improve the effectiveness of management in a natural ecosystem.