Ecological Archives E096-111-A1

Madelon Lohbeck, Lourens Poorter, Miguel Martínez-Ramos, and Frans Bongers. 2015. Biomass is the main driver of changes in ecosystem process rates during tropical forest succession. Ecology 96:1242–1252. http://dx.doi.org/10.1890/14-0472.1

Appendix A. A map of the study region, schematic overview of the variables used in the alternative models, and additional statistical details on the structural equation models.

Table A1. Ecosystem processes used in this study, their relation to ecosystem services (the benefits to people) and some key references.

Ecosystem process

Ecosystem service

Literature

Biomass productivity

Climate regulation (carbon uptake rates), wood provisioning

Clark et al. 2001; Janisch and Harmon 2002 

Litter productivity

Soil retention, decomposition, soil nutrients, net primary production

Ewel 1976; Ostertag et al 2008

Actual (in situ) litter decomposition

Soil quality, soil biotic activity, litter quality, soil nutrients, soil fertility, productivity, nutrient cycling

Swift 2004; Ewel 1976; Cornelissen 1996; Wardle 2004

Potential (ex situ) litter decomposition

Litter quality, soil nutrients, soil fertility, productivity, nutrient cycling

Swift 2004; Ewel 1976; Cornelissen 1996; Wardle 2004

 

Table A2. Results of the top five alternative models for each of the ecosystem processes. Given are the functional diversity component included in each model (functional richness, functional evenness, functional divergence or none), the community-weighted mean component included (CWM of Specific Leaf Area [SLA], Leaf Dry Matter Content [LDMC], litter Nitrogen [N], litter Phosporus [P] or none in case of decomposition, community-weighted mean of Specific Leaf Area [SLA], leaf Nitrogen [N]and Wood Density [WD] or none in case of productivity, see also Fig. A1). Given are the degrees of freedom, test statistic (minimum function chi-square), comparative fit index (CFI), model p-value, model Akaike Informaton Criterion (AIC), ∆ AIC compared to the best fitting model, standardized root mean square residual (SRMR), root mean square error of approximation (RMSEA) and the variance of the ecosystem process that was explained by the model (R²).

 

FD    

CWM

DF

χ2

CFI

p-value

AIC

∆AIC

SRMR

RSMEA

R2

Biomass productivity

 

 

 

1

Divergence

SLA

4

4.92

0.97

0.296

328.3

0

0.059

0.124

0.25

2

-

SLA

2

4.29

0.94

0.117

337.8

9.5

0.052

0.276

0.22

3

Richness

N

4

6.15

0.94

0.189

346.2

17.9

0.125

0.189

0.36

4

Divergence

N

4

6.83

0.91

0.145

354.0

25.7

0.148

0.217

0.33

5

* Only 4 models fit criteria

Litter productivity

 

 

 

1

Evenness

WD

4

6.89

0.94

0.142

165.5

0

0.078

0.219

0.77

2

Richness

WD

4

6.71

0.95

0.152

169.7

4.2

0.073

0.213

0.77

3

Divergence

WD

4

8.79

0.9

0.067

177.2

11.7

0.112

0.282

0.76

4

-

WD

2

5.15

0.93

0.076

186.6

21.1

0.037

0.324

0.76

5

Evenness

-

2

5.49

0.93

0.064

196.9

31.4

0.048

0.341

0.77

Actual litter decomposition

 

 

 

1

Evenness

Litter N

4

9.14

0.83

0.058

281.6

0

0.071

0.293

0.24

2*

* Only 1 model fits criteria

Potential litter decomposition

 

 

 

1

Divergence

LDMC

4

8.87

0.88

0.064

130.4

0

0.096

0.285

0.61

2

-

LDMC

2

2.67

0.98

0.263

142.5

12.1

0.043

0.15

0.46

3

Evenness

Litter P

4

8.81

0.88

0.066

165.7

35.3

0.092

0.283

0.59

4

Richness

-

2

2.99

0.97

0.224

182.4

52

0.046

0.182

0.41

5

Richness

Litter P

4

6.07

0.95

0.194

186.1

55.7

0.096

0.186

0.60

 

Table A3. Test statistics for the single best model for each studied ecosystem process (see also Fig. 2). Given are the unstandardized path coefficients (path-coefficient), the related standard error, Z-value, P value and standardized path coefficients (Std estimate) (as found in Fig. 2). The significant paths are highlighted in bold. For the overall fit of the model (Chi-square, AIC) see Table A2. Age refers to the time (in years) elapsed after agricultural field abandonment, F-evenness to functional evenness, F-divergence to functional divergence, SLA to community-weighted mean specific leaf area, WD to community-weighted mean wood density, LDMC to the community-weighted mean leaf dry matter content and Litter N to the nitrogen content of the litter.

Regressions

 Path-coefficient

Std error

Z-value

p-value

Std estimate

Biomass productivity

Biomass ~ age

1.293

0.156

8.311

0.000

0.906

F-divergence ~ biomass

0.005

0.003

1.592

0.111

0.380

SLA ~ biomass

-0.074

0.021

-3.456

0.001

-0.666

Biomass productivity ~ Biomass
F-divergence
SLA

 

0.201
-3.413
1.339

 

0.092
5.738
0.792

 

2.186
-0.595
1.691

 

0.029
0.552
0.091

 

0.686
-0.144
0.508

Litter productivity

Biomass ~ age

1.402

0.177

7.904

0.000

0.898 

F-evenness ~ biomass

-0.002

0.002

-1.243

0.214

-0.306

WD ~ biomass

0.001

0.001

0.859

0.390

0.217

Litter productivity ~
Biomass
F-evenness
WD

 

0.037
0.555
-0.058

 

0.006
0.718
1.030

 

6.736
0.773
-0.056

 

0.000
0.439
0.955

 

0.903
0.101
-0.007

Actual decomposition

Biomass ~ age

1.402

0.177

7.904

0.000

0.898 

F-evenness ~ biomass

-0.001

0.001

-1.029

0.303

-0.257

Litter N ~ biomass

-0.062

0.058

-1.073

0.283

-0.267

Litter decomposition ~ Biomass
F-evenness
Litter N

 

0.018
2.172
0.020

 

0.009
1.738
0.037

 

2.057
1.250
0.555

 

0.040
0.211
0.579

 

0.495
0.291
0.129

Potential decomposition

Biomass ~ age

1.402

0.177

7.904

0.000

.898

F-divergence ~ biomass

-0.001

0.003

-0.378

0.706

-0.097

LDMC ~ biomass

0.001

0.001

2.167

0.030

0.488

Litter decomposition ~ Biomass
F-divergence
LDMC

 

-0.005
0.411
-1.518

 

0.002
0.184
0.767

 

-2.200
2.232
-1.980

 

0.028
0.026
0.048

 

-0.406
0.360
-0.364

 

Table A4. Results of the correlation analyses based on the 15 successional plots with the variables that are used in the structural equation models, as well as species richness for comparison. This is presented for (a) the ecosystem processes biomass and litter productivity and for (b) the ecosystem processes actual and potential litter decomposition (because the variables used for the different set of processes differ slightly). Variables were grouped in successional indicators, diversity indicators, weighted mean and ecosystem process indicators. Sric refers to species richness, Fric to functional richness, Feve to functional evenness, Fdiv to functional divergence, SLA to community-weighted mean specific leaf area, WD to community-weighted mean wood density, N to community-weighted mean of leaf nitrogen content, LDMC to the community-weighted mean leaf dry matter content, Litter N to the nitrogen content of the litter and Litter P to the phoshorous content of the litter. Pearson pairwise correlation coefficients are given (* P < 0.05; ** P < 0.01, *** P < 0.001).

 

a. Correlation table for ecosystem processes biomass and litter productivity

 

Successional indicators

Diversity

Weighted mean

EP

 

Age

Biomass1

Biomass2

Sric

Fric

Feve

Fdiv

SLA

WD

N

Biomass product

Biomass1

 0.91***

 

 

 

 

 

 

 

 

 

 

Biomass2

0.90***

 0.99***

 

 

 

 

 

 

 

 

 

Sric

0.69**

0.75**

0.75**

 

 

 

 

 

 

 

 

Fric

 0.53*

 0.56*

 0.56*

0.79***

 

 

 

 

 

 

 

Feve

-0.18

-0.31

-0.31

-0.21

-0.38

 

 

 

 

 

 

Fdiv

0.32

0.38

0.38

0.35

0.26

-0.05

 

 

 

 

 

SLA

-0.76**

-0.67**

-0.63*

-0.26

0.01

0.04

-0.12

 

 

 

 

WD

0.18

0.24

0.22

0.07

-0.13

0.19

-0.33

-0.42

 

 

 

N

0.44

0.39

0.41

0.25

0.26

0.03

 0.71**

-0.22

-0.27

 

 

Biomass product

0.21

0.3

0.39

0.18

0.06

-0.22

0.06

0.07

-0.08

0.16

 

Litter product

 0.90***

 0.87***

0.87***

0.67**

 0.57*

-0.18

0.29

-0.56*

0.21

0.45

0.28

1 Biomass used to in the model for biomass productivity; 2 Biomass used to in the model for litter productivity

 

b. Correlation table for ecosystem processes actual and potential litter decomposition

 

Successional ind

Diversity

Weighted mean

EP

 

Age

Biomass

Sric

Fric

Feve

Fdiv

SLA

LDMC

Litter N

Litter P

Actual decomp

Biomass

0.90***

 

 

 

 

 

 

 

 

 

 

Sric

0.69**

0.75**

 

 

 

 

 

 

 

 

 

Fric

 0.73**

 0.65**

0.70**

 

 

 

 

 

 

 

 

Feve

0.02

-0.26

-0.09

-0.18

 

 

 

 

 

 

 

Fdiv

-0.35

-0.1

-0.07

-0.44

0.11

 

 

 

 

 

 

SLA

-0.76**

-0.63*

-0.26

-0.28

-0.24

0.23

 

 

 

 

 

LDMC

0.31

0.49

0.45

0.35

-0.05

0.14

-0.12

 

 

 

 

litterN

-0.24

-0.27

-0.33

-0.24

-0.03

-0.06

0.38

0.01

 

 

 

litterP

-0.37

-0.26

-0.14

-0.45

-0.12

0.39

0.22

-0.16

-0.33

 

 

Actual decomp

0.3

0.39

0.39

0.29

0.16

0.09

-0.07

0.49

-0.01

-0.12

 

Potential decomp

-0.60*

-0.64*

-0.41

-0.48

0.02

0.36

 0.54*

-0.52*

0.24

 0.58*

-0.38

 

Table A5. Eigenvector scores of community functional properties on the first and second principal components (PC1 and PC2). Principal components analysis (PCA) was used to separate plots based on their community functional properties. Separate PCAs are constructed for functional diversity components (functional richness, functional evenness, functional divergence) and community-weighted means as well as for productivity processes and decomposition processes as a different set of traits is used. Values in parentheses indicate variance accounted for by each axis. The PCAs were used to construct the compound variables (PC1 and PC2) that were used in the structural equation models of Fig. A5. Explanations of the abbreviations can be found in the description of Table A4.

 

Functional diversity

Community-weighted mean

Biomass and litter productivity

 

PC1 (49%)

PC2 (32%)

 

PC1 (76%)

PC2 (20%)

Fric

0.69

-0.04

SLA

0.57

0.61

Feve

-0.58

0.57

WD

-0.52

0.79

Fdiv

0.43

0.82

N

0.64

0.10

Actual and potential decomposition

 

PC1 (51%)

PC2 (31%)

 

PC1 (35%)

PC2 (33%)

Fric

0.67

-0.17

litterN

-0.77

0.087

Feve

-0.37

-0.92

litterP

0.33

-0.67

Fdiv

-0.64

0.35

SLA

-0.54

-0.53

 

 

 

LDMC

0.01

0.51

 

FigA1

Fig. A1. Schematic representation of the two different facets of the vegetation (quality and quantity), how these translate into three testable hypotheses on alternative drivers of ecosystem process recovery and what indices are used to test each of these (see also the conceptual model in Fig. 1). EPs stands for ecosystem processes, CWM for community-weighted mean, WD for wood density, leaf N for leaf nitrogen content, SLA for specific leaf area, LDMC for leaf dry matter content, Litter P and Litter N for litter phosphorus and nitrogen content.


 

FigA2

Fig. A2. Variation of each studied ecosystem process as a function of the strongest explanatory variable (highest standardized path coefficient), which was aboveground standing biomass in all four ecosystem processes (see Table A2). (a) Biomass productivity, (b) litter productivity, (c) actual litter decomposition, (d) potential litter decomposition. The fitted line is the ecosystem process as a function of the main explanatory variable, keeping the other factors contributing to the ecosystem process constant at their mean value.


 

FigA3

Fig. A3. Image of Mexico indicating the study location in the state of Chiapas in the white circle. Small white circles are the locations of the 15 secondary forest plots (< 1–29 yr) surrounding the village of Loma Bonita in Marques de Comillas. The continuous forest area north of the river Lacantún represents the Montes Azules Biosphere Reserve. The yellow line is the border with Guatemala. Plots of similar age (< 2 years difference) were at least 800 m apart.


 

FigA4

Fig. A4. Best fitting structural equation models for the four ecosystem processes when using compound variables instead of individual indicators (see Fig. 2 where individual indicators were used). As compound variables we used the principal component axes (PC1 and PC2) when secondary forest plots were separated based on their functional diversity indices (FD_PC1 and FD_PC2) and based on their community-weighted means (CWM_PC1 and CWM_PC2). See table A5 for the eigenvectors scores of the community functional properties with PC1 and PC2, note that a slightly different set of traits was used for different processes.
(a) Biomass productivity (Biomass product) was significantly explained by biomass and together with FD(PC2) and CWM(PC1) explained 24% of the variation (df = 4), (b) Litter productivity (Litter product) was explained by biomass and together with FD(PC2) explained 76% (df = 2), (c) Actual decomposition (Actual decomp) was explained by biomass and together with FD(PC1) and CWM(PC1) explained 22% of the variation (df = 4) and (d) Litter potential decomposition (Potential decomp) was explained by biomass and explained 40% (df = 1). The most suitable structural equation model presented here follows the selection criteria described in the methods, and based on different combinations of PC1 and PC2 for functional diversity and community-weighted mean indicators. Some paths are left out from the figure (unlike in Fig. 2) because this reduced the AIC based on which the best model was selected. This difference suggests that individual indicators for FD and CWM (as in Fig. 2) capture more relevant information for explaining ecosystem processes than does the use of compound variables. Black arrows represent significant relations (P < 0.05), gray arrows nonsignificant ones, given are the standardized path-coefficients.


 

Literature cited

Clark, D. A., S. Brown, D. W. Kicklighter, J. Q. Chambers, J. R. Thomlinson, J. Ni, and E. A. Holland. 2001. Net primary production in tropical forests: an evaluation and synthesis of existing field data. Ecological Applications 11:371–384.

Cornelissen, J. H. C. 1996. An experimental comparison of leaf decomposition rates in a wide range of temperate plant species and types. Journal of Ecology 84:573–582.

Ewel, J. J. 1976. Litter fall and leaf decomposition in a tropical forest succession in Eastern Guatemala. Journal of Ecology 64:293–308.

Janisch, J. E., and M. E. Harmon. 2002. Successional changes in live and dead wood carbon stores: implications for net ecosystem productivity. Tree Physiology 22:77–89.

Ostertag, R., E. Marín-Spiotta, W. L. Silver, and J. Schulten. 2008. Litterfall and decomposition in relation to soil carbon pools along a secondary forest chronosequence in Puerto Rico. Ecosystems 11:701–714.

Swift, M. 2004. Biodiversity and ecosystem services in agricultural landscapes? Are we asking the right questions? Agriculture, Ecosystems & Environment 104:113–134.

Wardle, D. A. 2004. Ecological linkages between aboveground and belowground biota. Science 304:1629–1633.


[Back to E096-111]