Biomass and dominance of conservative species drive above-ground biomass productivity in a mediterranean-type forest of Chile

Ariel Isaías Ayma-Romay and Horacio E.Bown

Abstract

Keywords: Above-ground biomass productivity,Diversity,Niche complementarity,Mass-ratio,Mediterranean,Soil resources

Background

The aboveground biomass productivity (AGBP) of the world’s forests supports ecosystem functioning and key ecosystem services to human well-being (Chapin et al.1997). At small-scales and under natural conditions AGBP is driven by vegetation quantity, species diversity,particular dominant species and soil resources availability through four so far controversial hypothesis, recently compared (Lohbeck et al. 2015; Finegan et al. 2015;Prado-Junior et al. 2016; Kamruzzaman et al. 2017; van der Sande et al. 2017a), and namely:

a) The hypothesis of vegetation quantity proposes that the amount of vegetation is the main driver of changes in the rates of different ecosystem processes (Lohbeck et al.2015).Consistent with this hypothesis, some studies show that greater aboveground biomass, tree density, basal area and leaf area increase AGBP (Paquette and Messier 2011;Reich 2012; Prado-Junior et al.2016). However, greater vegetation quantity also might reduce AGBP, due to competition for resources and mortality (Finegan et al.2015). Thereby, forests with smaller biomass have shown be more productive than primeval forests (van der Sande et al. 2017a;Ojha and Dimov 2017).

b) The hypothesis of niche complementarity predicts that greater plant species diversity increases the rate of ecosystem processes maximizing resources use(water, nutrients), since different species use different resources, or the same resources, but with spatial and temporal variations(Tilman 1997).Then, it is expected that forest AGBP increase with a greater number of tree species (richness) and evenness (Zhang et al.2012),as well as with higher functional diversity (Paquette and Messier 2011).However, other studies show neutral or negative effects of these variables on AGBP (Finegan et al.2015; Prado-Junior et al. 2016; van der Sande et al. 2017a).

c) The hypothesis of mass-ratio proposes that the extent to which species affect ecosystem processes are proportional to their contribution to the total plant biomass (Grime 1998). Then, highly productive dominant species(species with greater biomass)with particular traits may drive forest AGBP (Ruiz-Jaen and Potvin 2010; Finegan et al. 2015; Prado-Junior et al. 2016).However, some studies showed that dominance had little effect on forest AGBP,compared with the effect of the species diversity and vegetation quantity (Lohbeck 2015; van der Sande et al. 2017a).

d) The belowground resource availability (here coined as a hypothesis) predicts that soil resource availability is a key factor limiting ecosystem processes (Chapin et al. 2014; van der Sande et al.2017b). Thereby,greater soil nitrogen (Campo and Vázquez-Yánes 2004)and/or soil water contents may increase AGBP (Sala et al. 1988; Poorter et al.2017).However, low levels of resources available in forest soils do not always translate into reduced AGBP (Prado-Junior et al. 2016), as tree species develop adaptations allowing them to grow in conditions limited by resources (Wright et al. 2011).

The relative importance of these mechanisms changes with the environmental context and forest composition,and their influences are not always clear (Díaz et al.2007). Therefore, one of the greatest challenges in current forest management and climate change mitigation depends on correctly identifying and manipulating the drivers to maintain or even enhance forest productivity (van der Sande et al. 2017a). However, there is uncertainty about the effects of these mechanisms working simultaneously on forest productivity.

Mediterranean-type climate ecosystems occupy 5% of the terrestrial surface of the world while holding 20% of the global plant diversity (Rundel et al. 2016). The matorral is a typical mediterranean-type forest of central Chile composed by evergreen sclerophyllous and deciduous malacophyllous plant species. The matorral is distributed across highly variable environmental conditions(topography, altitude) and exposed to permanent human disturbances (deforestation, fire, and livestock raising)(Rundel et al. 2016), - which in turn have modified the vegetation quantity (Bown et al. 2014), soil resources(Mooney and Dunn 1970), diversity and composition(Rundel 1981). The consequences of these changes in aboveground biomass productivity (AGBP) have been scarcely evaluated (Fuentes et al. 1995).

In line with the hypothesis of vegetation quantity,AGBP in the matorral migth increase with higher tree density, as it reduces environmental stresses, such as high radiation, high soil evaporation and soil erosion(Mooney and Dunn 1970). Besides, according to the hypothesis of soil resources availability, AGBP in the matorral might increase with increasing soil nitrogen(Kummerow et al. 1982) and soil water contents (Rodríguez-Soalleiro et al. 2018; Meza et al. 2018), since soil nutrients and water deficiencies are common factors limiting the growth of mediterranean-type forests(Mooney and Dunn 1970).

The effect of niche complementarity on Chilean matorral AGBP has not been evaluated before. However,the AGBP might increase with high species and functional diversity, as in previous studies carried out in the Mediterranean Basin (Vilà et al. 2003, 2007, 2013; Ruiz-Benito et al. 2014). Niche complementarity in the matorral might be expected, since deciduous malacophyllous and evergreen sclerophyllous species can use the soil water from different soil depths and grow in dissimilar seasons of the year, promoting spatial and temporal niche differentiation (Giliberto and Stay 1978;Montenegro et al. 1982a; Canadell and Zedler 1995).However, in line with the mass-ratio mechanism, the AGBP of the matorral might also increase in communities dominated by evergreen sclerophyllous species with a conservative strategy (sensu Reich 2012), such as was observed in the Mediterranean Basin (Ruiz-Benito et al.2014). These conservative species are characterized by longer leaf life span and higher plant height (Canadell and Zedler 1995). Such attributes may allow them to uptake carbon continuously over the year (Mooney 1977), with a greater competitive advantage for light acquisition compared to deciduous malacophyllous species (Armesto et al. 1994).

The aim of this study was to evaluate how alternative/complementary mechanisms based on the previously described hypotheses of vegetation quantity, niche complementarity, mass-ratio, and soil resources availability explain the AGBP of the Chilean matorral.

Methods

Study area

The study area is located in the National Reserve“Roblería de Los Cobres de Loncha” (NRRCL) (34°7′36′′ S;71°11′18′′ W), Commune of Alhué, Metropolitan Region, Chile. Cold, wet winters (May-September) and dry, warm summers (October-April) characterize this part of Central Chile as a mediterranean region (Hajek and di Castri 1975). This region has a temperate and sub-humid climate, with an average annual air temperature of 14.7°C and average annual precipitation of 474 mm(Amigo and Ramírez 1998).This typical Chilean matorral has high variability in composition and structure due to high environmental heterogeneity and historic human disturbances (Rundel et al. 2016). The vegetation in the NRRCL is described as “sclerophyllous woodlands of the Coastal mountain range”. This type of vegetation has been disturbed mainly by fire, cattle grazing and firewood extraction,which have formed different types of vegetation in a range from evergreen sclerophyllous woodlands to deciduous malacophyllous scrublands(Gajardo 1994). These types are distributed in medium altitudes (300 to 1200 m a.s.l.) in moderate slopes. The soil is derived from granitic and colluvial origin. The soil textures are loams, with a presence of quartz and gravel(Bown et al. 2014).

Fig.1 Location of 25 m×25 m plots distributed across evergreen sclerophyllous woodlands to malacophyllous scrublands in the National Reserve of Robleria de los Cobres de Loncha(Metropolitan Region,Chile). Legend: Types of matorral:Ca-Qs:Cryptocarya alba and Quillaja saponaria woodlands;b) Qs-Lc:Q. saponaria and Lithraea caustica shrublands;and,c)Acacia caven and Trevoa quinquenervia scrublands

In 2010 twelve plots of 25 m×25 m were established(Fig. 1). These plots were distributed in different types of matorral with broad variation in structure, species diversity, composition and soil resources availability across sclerophyllous woodlands and malacophyllous shrublands. Historic human disturbance is the key factor explaining such differences across the matorral(Armesto and Martínez 1978; Rundel 1981). The plots were distributed in a small altitudinal range from 300 to 800 m within a similar lower northern aspect slope extending for 13 km. The effect of the altitude upon the AGBP of the matorral was not evaluated into this study, because we assumed that differences in vegetation and soil resources were not explained by altitude. The vegetation of these plots was not subject to burning or felling in the last 22 years, since 1996, the creation date of the NRRCL. A more detailed description of the soil and vegetation of the plots can be found in Bown et al.(2014).

Matorral types

To have a better description of the species diversity and composition of the matorral across the 12 plots, we classified the matorral into three types: a) Cryptocarya alba(Lauraceae) and Quillaja saponaria (Quillajaceae) woodlands (Ca-Qs, six plots); b) Q. saponaria and Lithraea caustica (Anacardaceae) shrublands (Qs-Lc, three plots);and, c) Acacia caven (Fabaceae) and Trevoa quinquenervia (Rhamnaceae) scrublands (Ac-Tq, three plots)(Fig.2).This classification was performed according to the floristic dissimilarity among plots through cluster analysis.The floristic dissimilarity was calculated on the basis of the Bray and Curtis index using the relative proportion of basal area among species (Additional file 1: Online Resource 1). The cluster was carried out through the hierarchical average method. Data analysis was performed in R version 3.4.0(R Core Team 2017).We used the package“vegan”for this analysis(Oksanen et al.2017).

Abiotic and biotic variables measurements

In line with each hypothesis, we calculated in total 19 abiotic and biotic variables at the plot level (Table 1).

Vegetation quantity measurements

In 2010, tree density (TD), tree basal area (BA) and leaf area index of shrubs and trees (LAI) were measured as surrogates of vegetation quantity. TD and BA were quantified with trees having one or more stems alive with a DBH ≥5 cm (this criterion included species with multistemmed individuals). The leaf area was estimated with an AccuPAR/LAI ceptometer (Decagon Devices,Inc.) at 30 cm above the soil level (Table 1).

Functional traits measurements

Fig.2 Types of matorral distributed in the National Reserve“Roblería de Los Cobres de Loncha”(Central Chile)

Table 1 Descriptive statistics of all biotic and abiotic independent variables used to predict aboveground biomass productivity of the Chilean matorral across the 12 plots

As in previous studies, we selected the plant height (H)and leaf life span (LLS) as functional traits reflecting different strategies to use carbon/water among plant species (Paquette and Messier 2011; Ruiz-Benito et al.2014). The H was measured as the shortest distance between the upper crown and the ground level,expressed in meters (Pérez-Harguindeguy et al. 2013).The mean value of H for each species was obtained measuring trees DBH ≥5 cm of DBH across the 12 plots covering all species (292 trees were measured).The H was estimated with a clinometer PM 5/360(SUUNTO, Finland). The LLS was obtained from the literature (Montenegro et al. 1982b; Specht 1988).This functional trait was defined as the time period during which an individual leaf is alive and physiologically active, and expressed in months (Pérez-Harguindeguy et al. 2013).

Measuring species and functional diversity

As in previous studies, species richness (S), Shannon-Wiener index (H′) (Zhang et al. 2012) and the mean species pairwise dissimilarity (MPD) (de Bello et al.2016) were calculated as surrogates of niche complementarity (Table 1). S represented the number of different species within a plot. H′ was calculated following the equation:

where H′ is a metric of community diversity, by taking into account the number (S) and the distribution of the species relative abundance (pi) - which was calculated through the ratio between the basal area of each species and the total basal area per plot.

MPD was computed following the equation (de Bello et al. 2016):

where, S is the species richness in a community and, dijis the distance between pairwise species i and j. The piand pjare the relative abundances of species i and j, respectively. The dijwas calculated with the Gower distance using species’ traits, such as plant height (H) and leaf life spans (LLS). The relative species abundances were based on the basal area (Additional file 1: Online Resource 1).

Measuring species and functional dominance

In line with Conti and Diaz (2013) the relative abundance (pi) of two late-successional sclerophyllous (C.alba and P. boldus) and two early-successional sclerophyllous (L. caustica and Q. saponaria) (Armesto et al.1994) were quantified with the relative basal area of the corresponding species. These four species were those most dominant and frequent across the plots (Additional file 1:Online Resource 1).

Additionally, the community-weighted mean trait value (CWM) was calculated using the species mean values of plant height (H) and leaf life span (LLS) and the relative abundances through the following equation(Garnier et al. 2004):

where, S=number of species in the community; ri=value of the functional trait of the ithspecies; and pi=relative abundance of the ithspecies, calculated as the relative basal area of the species.

Soil resources measurements

Spatial differences of soil resource availability were measured through the soil volumetric water content (VWC),total nitrogen (N), carbon-nitrogen ratio (C/N), soil inorganic nitrogen in the form of ammonium and nitrate(SIN) and soil texture (sandy, lime, clay) (Table 1). The VWC in each plot was obtained in winter and spring of 2009 and summer and autumn of 2010. In each season three soil samples were taken per plot using a soil corer(a metal cylinder of 5 cm in diameter and 30 cm in length). The VWC (m3H2O per m3soil) was calculated as the gravimetric water content (?g) multiplied by the soil bulk density (ρb) according to the following equation(Silva et al. 2000):

The soil bulk density was calculated by the wax paraffin method with aggregate samples taken from 10 to 40 cm deep(Gómez-Giraldo 2013).

The soil N and total soil carbon (C) were obtained for each plot taking three sample from 0 to 30 cm in soil depth the year 2012.The C/N ratio was used as a proxy of the nitrogen released by microbial decomposition. Higher C/N ratios (low nitrogen) reduces mineralization rates,since micro-organism requieres nitrogen for growth,while lower C/N ratios increase the mineralization proccess due to nitrogen surpassing micro-organism requirements for growth,and thus nitrogen become avaible to plants(Chapin et al.2014).For this reason,C/N ratios have an inverse relation with forest productivity (Yamakura and Sahunalu 1990). Total soil C and N were determined using a CNS analyzer at the Soils and Foliage Analyses Laboratory,Universidad Católica de Chile,Santiago,Chile.

SIN was obtained at the plot level by adding ammonium (NH+4) and nitrate (NO-3) chemical extractions.Four soil samples taken at 0-30 cm in soil depth were randomly cored from each plot, placing 10 g of fresh soil from each sample in 50 mL 2 mol·L-1KCl. Soil extracts were shaken for two hours to elucidate ions from the soil to the KCl solution. Then, soil extracts were left to settle for half an hour before they were filtered using Whatman N° 1 filter paper. Ammonium and nitrate concentration were determined using the Kjeldahl steam distillation method (Bremner 1965). The SIN was sampled in winter and spring of 2009 and summer and autumn of 2010. The measurements were averaged for the year. Our measurement of SIN was a surrogate of nitrogen availability to plants released by microbial activity in a year (Gray and Schlesinger 1981).

Soil texture data was obtained from a previous research carried out in the same 12 plots used in this study during 2012 (Nangari-Pi?a 2012). That study estimated silt, clay and sand fractions from 0 to 50 cm in depth.Soil texture was included in the analysis as it is a key component driving soil resource availability (Noy-Meir 1973). For instance, sandy soils may increase AGBP in ecosystems limited by water, enhancing water infiltration, water storage to greater soil depths avoiding evaporation (Sala et al. 1988), and facilitating root growth to deeper layers (van der Sande et al. 2017a).

Estimating forest above-ground biomass productivity

At the plot level, above-ground biomass productivity(AGBP) was calculated as the difference in aboveground biomass (AGB) of all the living trees in the summer of 2010 (to) versus those alive in the spring of 2017(tf) following the equation proposed by Prado-Junior et al. (2016):

AGBP=Above-ground biomass productivity of surviving trees between tfand to(kg·m-2·year-1),

AGBo=Above-ground biomass of trees alive in to(kg·m-2),

AGBf=Above-ground biomass of trees alive in tf(kg·m-2),

tf-to=Interval of time(years).

AGB was estimated adding the biomass of all the trees alive with a DBH ≥5 cm. The lack of species-specific models to estimate AGB in the Chilean matorral led us to use the following general allometric equation (Chave et al. 2005):

AGB=Aboveground biomass (kg),

ρ=Stem wood density (g·cm-3),

DBH=Diameter at breast height(cm).

The ρ of the nine dominant species was obtained from branch stem samples collected from 10 mature trees by species during the spring of 2016,following standardized protocols (Pérez-Harguindeguy et al. 2013) (Additional file 1:Online Resource 2).

Data analysis

We analyzed the effects of the independent variables on AGBP through linear regressions, having as reference the sequential procedure proposed by Díaz et al. (2007).A pairwise linear regression analysis was performed to test the effect of each independent variable on aboveground biomass productivity (AGBP). Prior to this analysis, linear and nonlinear relationships of each independent variable with AGBP were checked. Nonlinear (as logarithmic, potential and exponential relationships) were analyzed through simple linearized regressions using logarithmic transformations. If the nonlinear models were marginally better in predicting AGBP, then for purposes of parsimony, simpler linear models were used (Additional file 1:Online Resource 3). Then, the best predictors into each mechanism explaining the AGBP were identified through a backward elimination procedure. These best predictors were introduced in a global linear regression model and then we derived a set of models with all the possible combinations, after which were ranked based on Akaike’s information criterion for small samples (AICc) (Anderson 2008). The AICc differences (Δs) between the best model and any other model were used to rank them from best to worst. The Akaike’s weight (w) was derived for each model to describe the probability that a given model is the best in a scale from 0 to 1(Burnham and Anderson 2011).The best subset of models was selected using a confidence of 95% (those models whose cumulative Akaike weight,wi, ≤0.95) (Symonds and Moussalli 2011). Collinearity among predictors was assessed in the best subset of models with the variance inflation factor (VIF). Those models having a VIF ＜10 were maintained in the best subset (Zuur et al. 2010). Then, the best subset was averaged calculating a weighted averaged of the partial coefficients(β) and the relative importance (RI) of predictors (probability that a given predictor is involved in the best model).All the data analysis was performed in the R system for statistical computing version 3.4.0 (R Core Team 2017).The species and functional diversity indices were calculated using the “vegan” package (Oksanen 2011). The CWM was calculated using the “FD” package (Laliberte and Legendre 2010). The package “MuMIn” was used for calculating the AICc and multimodel inference approach(Barton 2018).

Results

Aboveground biomass productivity (AGBP) was on average of 0.14 kg·m-2·year-1. Generally, AGBP increased

from malacophyllous deciduous shrublands to evergreen

sclerophyllous woodlands ranging from 0 to 0.42 kg·m-22·year-1. The values of AGBP were positive in almost all the plots, but in two plots dominated by Q. saponaria and L. caustica the growth was nil. The analysis of simple pairwise regression suggested that AGBP significantly increased with increasing vegetation quantity (tree density, basal area, leaf area index) (Fig. 3a-c), species and functional dominance(relative abundance of C. alba and P. boldus, CWM of plant height and CWM of leaf life span) (Fig. 3d-g) and some soil attributes related to resources availability (sand %, volumentric water content, soil nitrogen content) (Fig. 3h-j). However, some variables related to species dominance (relative dominance of Q. saponaria and L. caustica) (Additional file 1:Online Resource 4), soil resources availability (SIN, C/N)(Fig. 3k-l) and all the variables related to niche complementarity (S, H′, MPD) (Fig. 3m-o) did not show any significant effect on AGBP. Further details about the simple linear regressions fitted can be found in Additional file 1: Online Resource 4.

The backward elimination procedure allowed identifying the best predictors driving AGBP into each subset of hypotheses. Such predictors were the following: tree density (hypothesis of vegetation quantity), relative abundance of C. alba (hypothesis of dominance), sand fraction(hypothesis of resources availability);and species richness and Shannon index (niche complementarity)(Additional file 1: Online Resource 5). These key predictors were also correlated to other variables linked to the same hypothesis (evaluated through the Pearson correlation coefficient). For instance, tree density was a variable closely related to BA (r=0.82; P ＜0.001) and LAI(r=0.82; P ＜0.01). The relative abundance of C. alba had a positive correlation with the relative abundance of P. boldus (r=0.67; P=0.02), CWM of H (r=0.73; P=0.01) and CWM of LLS (r=0.97; P ＜0.001). The sand %had a positive correlation with VWC (r=0.75; P=0.005)and N (r=0.72; P=0.01) and a negative association with the clay (r=-0.85; P ＜0.001) and silt fraction (r=-0.92;P ＜0.001). The ln (S) was correlated with the Shannon index (r=0.85; P ＜0.001) and MPD (r=0.94; P ＜0.001).See further details in Additional file 1:Online Resource 6.

The combination of these key predictors (tree density,piC. alba, richness and sand) yielded 16 candidate linear regression models - which were ranked and compared through the AICc multi-model inference procedure.This analysis showed that there was not an exclusive best model explaining AGBP. The predictors from all mechanisms hypothesized were included in the best subset of linear regression models (six models) with a confidence ＞95%. All these models combining two predictors had acceptable collinearity (VIF ＜10). The best models had tree density (vegetation quantity hypothesis) and/or the relative dominance of C. alba (mass-ratio hypothesis) as common predictors having positive effect on AGBP. The models that included either S (niche complementarity) or sand (%) (resource availability)were less commons. Furthermore, S had negative effects on AGBP(Table 2).

Fig.3 Aboveground biomass productivity of a Chilean Mediterranean-type matorral explained by four ecological mechanisms.Legend:Mechanisms hypothesized are:vegetation quantity(a,b,c);mass-ratio(d,e,f,g);soil resources availability(h,i,j,k,l)and niche complementarity(m,n,o).Similar to Conti and Díaz(2013)we show the changes of AGBP as function of abiotic and biotic variables across different types of matorral(see Fig.2).The pairwise relationships of AGBP with pi of Q.saponaria,pi of L.caustica,silt(%)and clay(%)are shown in Additional file 1:Online Resource 4

The averaging model procedure showed that tree density was the strongest predictor increasing AGBP(β=0.76; Importance=0.69) followed by the relative abundance of C. alba (β=0.30; Importance=0.36). In contrast, species richness had a small effect decreasing AGBP (β=-0.15; Importance=0.44) and, the sand (%)had the smallest positive effect on AGBP (β=0.006; Importance=0.04) (Table 3).

The analysis showed that vegetation quantity and massratio were the key mechanisms favoring matorral AGBP.In contrast, soil resources availability and niche complementarity were unlikely to largely increase AGBP.

Discussion

The above-ground biomass productivity (AGBP) found in this study is within the range for mediterranean-type forests (0.12 to 0.65 kg·m-2·year-1) (Ehleringer and Mooney 1983; Chapin et al. 2014). However, some matorral plots showed lower or null AGBP. Such marginal response might be due to recurring drought events from 2010 to 2015 in central Chile(Garreaud et al.2017).

Vegetation quantity and conservative species drive AGBP

Tree above-ground biomass productivity (AGBP) increased mainly by vegetation quantity, and dominance of species rather than by soil resources availability and niche complementarity, from deciduous malacophyllous shrublands to evergreen sclerophyllous woodlands. This finding agrees with studies in tropical and temperate forests showing that vegetation quantity (Paquette and Messier 2011; Lohbeck et al. 2015; Ojha and Dimov 2017) and dominance of particular species were more important for AGBP than other mechanisms (Finegan et al. 2015; Prado-Junior et al. 2016). Furthermore, our findings seem to contradict those from mediterraneantype forests of Spain (Troumbis and Memtsas 2000; Vilà et al. 2003) and California (Rundel 1982), where species diversity and nutrients are key drivers of productivity,respectively.

A greater vegetation quantity may increase AGBP by increasing the efficiency for capturing solar energy (Reich 2012) and reducing stressful micro-environmental conditions. Previous studies suggested that higher tree density in the Mediterranean Basin - type matorral had a positive relationship with carbon uptake (Vilà et al. 2013). It is possible that higher tree density and leaf area increase shading, which in turn, reduce the evaporation from the soil, and mitochondrial respiration to sustain growth and repair damages due to high radiative stresses (Holmgren et al.1997).

The species dominance may increase AGBP because some particular tree species, such as C. alba and P. boldus (evergreen sclerophyllous) increased the CWM of H and the CWM of LLS (all these variables were correlated) favoring the accumulation of greater vegetation quantity and longer periods intercepting light and assimilating carbon in the canopy of the matorral (Ehleringer and Mooney 1977; Rundel 1982). In contrast,deciduous malacophyllous have a smaller plant height and shorter period of growth,which might reduce AGBP(Mooney 1977;Montenegro 1982).

Vegetation quantity and dominance of taller species with a conservative leaf strategy to use carbon(such as C.alba) were pivotal/complementary mechanisms drivingAGBP in our study and also in other ecosystems under harsh environmental conditions (Paquette and Messier 2011; Vilà et al. 2013; Prado-Junior et al. 2016). However,under more favorable environmental conditions, such as in humid tropical forests, AGBP was both related negatively to vegetation quantity and positively to taller dominant species with acquisitive strategies to use carbon (＜CWM LLS) (Finegan et al. 2015; van der Sande et al.2017b). Possibly, a greater accumulation of biomass and the dominance of particular species might be the best suit of attributes to grow in a community under harsh environmental stresses.

Table 2 Best subset of models explaining aboveground biomass productivity(AGBP)

Table 3 Averaging model procedure showing standardized coefficients(β) and relative importance (RI) of predictors influencing AGBP

Soil resources availability have less importance for productivity

In line with previous studies carried out in semiarid ecosystems, we found that the sand fraction, the soil water content and soil nitrogen had a significant positive effect on AGBP of the matorral (Kummerow et al. 1982; Sala et al. 1988). A greater sand fraction might increase AGBP due to greater water infiltration and less evaporation (Sala et al. 1988), longer periods for releasing soil nitrogen by micro-organisms (Schwinning and Sala 2004), and root growth to deeper soil layers (van der Sande et al. 2017a). However, this study was unable to demonstrate that these soil variables would better explain AGBP than vegetation quantity and species dominance. Some studies suggested that AGBP is scarcely driven by spatial differences in soil resources,since dominant species may exhibit particular adaptative traits allowing them to develop under semi-arid and infertile conditions (Prado-Junior et al. 2016). For instance, evergreen sclerophyllous species have developed strategies to conserve water by reducing stomatal and mesophyll conductance (Giliberto and Stay 1978; Canadell and Zeedler 1995). Besides, evergreen sclerophyllous species may develop deep roots to capture water throughout the year independently of the water content in the topsoil (Hoffmann and Kummerow 1978). Likewise, it has been shown that the growth of sclerophyllous species can be marginally influenced by soil nutrients availability (Rundel 1982;Espinoza et al.2017).Some malacophyllous deciduous species (Rhamnaceae and Fabaceae) are N-fixers(Rundel 1982), which might increase AGBP. However, in semi-arid ecosystems, nitrogen soil availability depends upon soil water (Schwinning and Sala 2004), which is scarce in matorrals dominated by these species.

Furthermore, the SIN and C/N ratio in our study were unrelated to productivity, possibly because the nitrogen released by microbial activity in mediterraneanecosystems is scarce, since stressful environmental conditions and recalcitrant soil organic matter are limiting decomposition and nutrient release (Gray and Schlesinger 1981).

Negative effect of niche complementarity on productivity

In contrast to evidence showed by previous studies(Zhang et al. 2012; Liang et al. 2016) this study shows that either species diversity (richness and Shannon index) or functional diversity might not increase the matorral AGBP. In some studies, species richness had a positive effect on forest AGBP in parallel to vegetation quantity and dominance predictors (Ruiz-Jaen and Potvin 2010; Paquette and Messier 2011; Vilà et al. 2013).However, in this study species richness, which is correlated with evenness and functional diversity, had consistently a negative effect on AGBP, as in some previous studies (Finegan et al. 2015; Prado-Junior et al. 2016). It is possible that species diversity is only the result of disturbances promoting the establishment of a higher number of species rather than the result of complementary resource use (Grime 2006). Thus, a moderate disturbance of the matorral could increase species diversity(Rundel 1981), but in parallel could affect the competitive ability and plant density of dominant species, which are the potential drivers of AGBP.

Limitations

Overall, observational studies in natural ecosystems have inherent problems for testing causal relationships, as a result that abiotic and biotic variables while being correlated are influencing simultaneously ecological processes such as AGBP. However, multiple linear regression and model selection procedures might contribute to separate individual effects and assess such relationships (Díaz et al. 2007). Furthermore, we aware that our sample size is relatively small and we have collinearity among most of the independent variables, which might reduce the power of our inferences by increasing the probability to commit a type II error (i.e. failure to reject the null hypothesis when it is untrue) (Zuur 2010). However, despite these potential sources of error we were able to detect significant relationships among AGBP and most of the independent variables, which gave us chance to test all the hypotheses proposed. There were previous studies, with similar limitations, in which there was also possible to find significant effects of variables related to biodiversity on ecological process through linear regression models (Conti and Diaz et al. 2013; Lohbeck et al. 2015).

Besides, there might be a potential unexplained variation in our results due to geographic and climatic differences among plots, although as the spatial extent of our study was relatively small such variation was probably small and randomly allocated within the plots.

Conclusions

The findings of this study suggest that vegetation quantity and species dominance with greater plant height and a conservative strategy to use carbon are pivotal drivers for productivity. The importance of both mechanisms suggests that canopy quantity and quality may be key to sustain AGBP and linked ecosystem services in the mediterranean-type matorral, since these dominant species mediate greater biomass gains and assimilating carbon continuously during the year. In contrast, taxonomic species diversity, functional diversity and soil resources might be less important to drive AGBP, at least at this small scale of observation. Possibly, strong competition for resources avoids maximizing species diversity in the highly productive matorral,while topsoil resources are not a factor limiting productivity, due to sclerophyllous species uptake soil resources at greater soil depths.

Supplementary information

Supplementary information accompanies this paper at https://doi.org/10.1186/s40663-019-0205-z.

Additional file 1: Online Resource 1. Relative basal area of tree species across 12 plots in the National Reserve Roblería de los Cobres de Loncha (Chile). Online Resource 2.Wood density of nine tree species of the Chilean mediterranean-type matorral. Online Resource 3.Analyzing linear and nonlinear relationships of biotic and abiotic variables with AGBP.Online Resource 4.Simple pairwise linear regressions analyzing the effect of species diversity, vegetation quantity, species relative abundance and soil resources availability on above-ground biomass productivity of the mediterranean-type matorral of Chile. Online Resource 5.Backward elimination procedure for selecting variables explaining aboveground biomass through multiple linear regressions. Online Resources 6.Pairwise correlation matrix among independent variables.

Abbreviations

AGB: Above-ground biomass; AGBP:Above-ground biomass productivity;AICc: Akaike’s information criterion for small samples; BA: Basal area; C/N: Carbon-nitrogen ratio; CWM: Community Weight Median; CY: Clay;DBH: Diameter at breast height; ER:Evidence ratio; H: Plant height; H′: Shannon index; LAI: Leaf area index; LLS: Leaf life span; MPD:Mean Pairwise Distance; N: Total soil nitrogen;NRRCL: National Reserve Roblería de Los Cobres de Loncha; pi:Relative abundance of species; S:Species richness;SD: Sand; SIN: Soil inorganic nitrogen;ST: silt; TD: Tree density;VWC: Volumetric water content; w: Akaike’s weight; ρ: Stem wood density

Acknowledgments

We appreciate the collaboration of Roberto Cerda, park ranger,of the CONAF (National Forestry Corporation) for their help during the fieldwork in the National Reserve of “Roblería de Los Cobres de Loncha”. Besides, we thank Juan Pablo Fuentes(Universidad de Chile) for his collaboration recollecting and analyzing soil samples from all the plots of this study. Also,we thank Lucas Enrico (Universidad Nacional de Córdoba, Argentina) for his valuable suggestions for this paper.

Authors’ contributions

AA and HB read and approved the final manuscript. AA conceived the presented idea and developed the theory and analytical methods. AA and HB planned and carried out the measurements.AA and HB discussed the methodologies and results.AA wrote the manuscript with assistance from HB. Both authors read and approved the final manuscript.

Funding

Funding for this research was obtained from CONICyT (Comisión Nacional de Investigación Científica y Tecnológica) for the grant Fondecyt No 1150877. Additional funding was derived from the CONICyT doctoral grant No 21150802.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Received: 31 January 2019 Accepted: 15 October 2019