Sapling Structure and Regeneration Strategy in 18 ... - Semantic Scholar

Annals of Botany 96: 313–321, 2005 doi:10.1093/aob/mci179, available online at www.aob.oupjournals.org

Sapling Structure and Regeneration Strategy in 18 Shorea Species Co-occurring in a Tropical Rainforest M A S A H I R O A I B A 1,* and T O H R U N A K A S H I Z U K A 2 1

2

Center for Ecological Research, Kyoto University, Hirano, Otsu 520-2113, Japan and Research Institute for Humanity and Nature, Takashima 335, Kamigyo, Kyoto 602-0878, Japan

Received: 1 December 2004 Returned for revision: 15 February 2005 Accepted: 27 April 2005 Published electronically: 8 June 2005


Background and Aims Inevitable trade-offs in structure may be a basis for differentiation in plant strategies. Juvenile trees in different functional groups are characterized by specific suites of structural traits such as crown architecture and biomass distribution. The relationship between juvenile tree structure and function was tested to find out if it is robust among functionally and taxonomically similar species of the genus Shorea that coexist sympatrically in a tropical rain forest in Borneo.
Methods The sapling structures of 18 species were compared for standardized dry masses of 5 and 30 g. Pairwise simple correlation and multiple correlation patterns among structural traits of juveniles (01–15 m in height) of 18 Shorea species were examined using Pearson’s correlation and principal component analysis (PCA), respectively. The correlation was then tested between the PCA results and three indices of shade tolerance: the net photosynthetic rate, the wood density of mature trees and seed size.
Key Results The structural variation in saplings of the genus Shorea was as large as that found in sets of species with much more diverse origins. The PCA showed that both crown architecture and allocation to leaves are major sources of variation in the structures of the 18 species investigated. Of these two axes, allocation to leaves was significantly correlated with wood density and showed a limited correlation with photosynthetic rate, whereas crown architecture was significantly correlated to seed size.
Conclusions Overall, the results suggest that an allocation trade-off between leaves and other organs, which co-varied with wood density and to a certain extent with photosynthetic capacity, accounts for the difference in shade tolerance among congeneric, functionally similar species. In contrast, the relationship between the architecture and regeneration strategy differed from the pattern found between functional groups, and the function of crown architecture was ambiguous. Key words: Allocation, architecture, Borneo, dipterocarps, regeneration strategy, sapling, shade tolerance, seedling, Shorea.

INTRODUCTION Juvenile trees have suites of structural traits that differ markedly among species. This variety in traits such as crown architecture, biomass allocation and leaf characters is assumed to reflect their diversified regeneration strategies, i.e. regeneration niche (sensu Grubb, 1977), and hence to promote their coexistence (e.g. Horn, 1971). From this perspective, many ecologists have attempted to elucidate the ecological functions of the structure of tree seedlings and saplings. As a result, two inevitable trade-offs in structure are now known to be tightly correlated with potential growth speed and shade tolerance, i.e. the ability to survive in the shaded understorey. One is the architectural trade-off between a less-branched narrow crown supported by a slender trunk (small crown) and a well-branched wide crown on a thick trunk (large crown; King, 1990; Kohyama, 1987, 1991). A small crown could be advantageous for rapid height growth under relatively bright conditions, such as in gaps, whereas a large crown may increase light capture in shaded environments. The second is a trade-off relationship between a robust structure and a structure advantageous for carbon gain, which is mainly related to allocation patterns to organs and leaf structure. Despite early predictions that shade-tolerant * For correspondence. E-mail [email protected]

plants should maximize light interception (Givnish, 1988), they seem to have thicker lamina, smaller leaf area ratio and greater allocation to roots than light-demanding species, which minimize biomass loss rather than maximize carbon gain (Walters et al., 1993a; Kitajima, 1994; Cao and Ohkubo, 1998; Walters and Reich, 1999, 2000). Although these two trade-off axes have been well investigated, some critical questions remain unanswered. First, do these trade-offs in structural traits contribute to the differentiation of plant strategies even among species within a single taxonomic or functional group? Many previous studies compared contrasting functional groups and distantly related taxa, although species within a single taxon or functional group can also show considerable differences in both their structure and their shade tolerance (Ashton, 1995; Ackerly and Donoghue, 1998; Davies, 1998; Montgomery and Chazdon, 2002). Especially in tropical rain forests, dozens of species within a single large genus coexist sympatrically and show various structural characters. As competition is thought to be more severe among these phylogenetically and ecologically related species, whether niche differentiation exists and how it occurs is a more critical question for these species. Secondly, how are the two axes, i.e. the architectural trade-off and the biomass partitioning and leaf structure trade-off, related to each other? In some studies that compared species from different

ª The Author 2005. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For Permissions, please email: [email protected]

Aiba and Nakashizuka — Sapling Structure and Function in Shorea

314

taxa and functional groups, species with a small crown showed more investment to the trunk and less to leaves than species with a large crown (Kohyama, 1987; Kohyama and Hotta, 1990; Takahashi et al., 2001), suggesting further synthesis of structural traits on a single trade-off axis. However, all of these studies of sapling architecture were conducted on the basis of height or above-ground mass, without evaluating the allocation to roots. Hence, the relationship between the two axes, both of which are related to shade tolerance, remains unclear. Diversification of congeneric tree species is a particular feature of Asian tropical forests, and structural variation seems to be one significant mechanism that allows the stable coexistence of these species (Thomas, 1996, 2004). One such diversified genus, Shorea (Dipterocarpaceae), dominates the canopy and the emergent layers of Bornean dipterocarp forests and many species coexist sympatrically there. All species in the genus Shorea are considered non-pioneer species with a restricted distribution in primary or relatively undisturbed forests, although they may show different levels of shade tolerance (Turner, 1990; Ashton, 1995; Still, 1996; Barker et al., 1997). These conditions provide an excellent opportunity to compare the sapling structure of many species within a single genus and with similar life histories. Despite the merits of using this genus, little is known about the shade tolerance of these species. Therefore, co-variation was tested between structure and three indices of shade tolerance: the maximum net photosynthetic rate, the wood density of mature trees and seed size. The maximum net photosynthetic rate can be considered to be an index of the potential for carbon assimilation in high-light environments (Walters et al., 1993a), and wood density is negatively correlated with the innate growth speed (e.g. Muller-Landau, 2004). Therefore, these two traits should be negatively and positively correlated with shade tolerance, respectively, and inversely correlated with the potential growth rate. While these two indices represent shade tolerance throughout the life history of the tree, seed size correlates with the temporal shade tolerance just after germination, and the correlation diminishes within a few years (Poorter and Rose, 2005). Such trends are often observed even among ecologically similar species (Reich et al., 1995; Still, 1996; Davies, 1998; Suzuki, 1999) and, by using these indices, subjective, qualitative measures of shade tolerance were avoided. In this study, the biomass allocation and plant architecture of juveniles (01–15 m tall) were analysed on a whole-plant basis, while considering ontogeny, in 18 co-occurring climax species of the genus Shorea in a Bornean tropical rainforest. The specific questions addressed were (a) which suite of morphological traits co-varies with the indices of shade tolerance and (b) how do the two aspects of plant structure, i.e. architecture and allocation, correlate with each other? MATERIALS AND METHODS Study site and species

This study was conducted at Lambir Hills National Park, Sarawak, Malaysia (4 20 N, 113 500 E, 150 m a.s.l.).

Consistently high temperatures, with no distinct dry season, characterize the climate of this region. The mean annual temperature and annual precipitation are approx. 26  C and 2700 mm (Nakagawa et al., 2000). The park is 6949 ha in area and is entirely covered with a primary lowland mixed dipterocarp forest on nutrient-poor sandy clay or clay-rich soils (Yamakura et al., 1995). Trees of the genus Shorea (Dipterocarpaceae) are dominant in the canopy and the emergent layer of the lowland mixed dipterocarp forest of West Malesia (Ashton, 1991). This genus is highly diversified in Borneo, and 129 species have been recorded on the island (Newman et al., 1996). Of the approx. 70 species of the genus Shorea that occur in Lambir, 18 species, for which enough saplings of a comparable size (01–15 m tall) could be collected, were analysed. Table 3 provides details of the species analysed. The species were sampled from six taxonomic sections, which are thought to differ in terms of shade tolerance (cf. Still, 1996), i.e. species of section Shorea are considered highly shade tolerant, while some species of sections Brachyptera, Mutica, Ovalis and Pachyparpae seem to be relatively light demanding. Most of the selected species were emergent trees (maximum height >50 m), although two lower species, S. biawak and S. macroptera, which reached up to 30 m were included. Biomass partitioning and other structural traits

In August 2003, approx. 20 sound saplings of various sizes (01–15 m in height) of each species were sampled to determine their allometric relationships. All sampled individuals grew under closed canopies. For S. ochracea, only 18 individuals were sampled because of its rarity. Also sampled were 10–15 saplings of each species excluding S. ochracea, in November 2003. After the trunk height (the vertical distance from the ground to the highest apex), crown width in two directions at right angles (including the widest width) and crown depth (the vertical distance from the lowest to the highest leaves) had been measured, the entire plant body was harvested. After harvest, hemispherical photographs were taken at 1-m height to measure canopy openness for three randomly selected individuals per species. Canopy openness was consistently low (59 6 07 %, mean 6 s.d.), indicating that the light environment under the closed canopy was homogenous. After recording the stem diameter at ground height in two directions at right angles, the sampled saplings were divided into three parts: root, stem and leaf lamina. In the latter sampling period, the stem was further divided into the main trunk and branches. Though S. ochracea was sampled only in the first sampling period, any individual of the species sampled lacked branches. Petioles were treated as stems in the first sampling and as the trunk or branches in the latter. For each sapling, the total leaf area and maximum individual leaf area were measured using the software LIA32 for Windows 95 (ver. 0371b1, K. Yamamoto) after scanning the images. All parts of the saplings were weighed after being oven-dried at 60  C for at least 5 d. The crown projection area and crown volume of each individual were approximated as an ellipse and an ellipsoid of revolution,

Aiba and Nakashizuka — Sapling Structure and Function in Shorea T A B L E 1. List of definitions and units of dimensions Parameter

Definition

Unit

TM RM SM BM LM LA Lmax CA CV H D SV Amax WD

Total dry mass Root dry mass Stem dry mass Branch dry mass Leaf lamina dry mass Total leaf area Maximum individual leaf area Crown projection area Crown volume Plant height Diameter of stem at ground height Seed volume Net photosynthetic rate per leaf area Wood density of matured tree

g g g g g 102 cm2 102 cm2 102 cm2 103 cm3 cm mm cm3 mmol CO2 m2 s1 g cm3

respectively. All parameters measured in this study are summarized in Table 1 with abbreviations and units. Indices of shade tolerance

The maximum net photosynthetic rate (Amax) of saplings under the closed canopy was quantified in the field using a portable infrared gas analyser (LI6400; LI-COR Incorporated, Lincoln, NB, USA) in November 2003. Although it is better to measure Amax under high light conditions to use as an index of shade tolerance, Amax was measured under shaded conditions to reduce the variation in site conditions where the plants grow, given that the light conditions of the shaded understorey are far more homogenous than that of a naturally generated high-light site, i.e. gaps of various size and age. All measurements were conducted in the morning (0700–1000 h). During all measurements, the airflow rate through the chamber was 500 ml min1, and the air temperature was maintained at 27  C. The CO2 concentration of the air was kept at 450 mg kg1, which is a typical value in the forest understorey, where CO2 from the canopy accumulates, and the humidity was sufficiently high (approx. 80 %). Light was provided artificially with a light bulb. Three to four approx. 1-m-tall saplings of each species that grew under a closed canopy were selected and the youngest fully expanded leaf on each sapling used for the measurements. After acclimation at a photon flux density of 400 mmol m2 s1, the photon flux density was gradually increased with sufficient acclimation until saturation of photosynthesis occurred. The saturated net photosynthetic rate was recorded as the maximum net photosynthetic rate. The wood density (WD) of mature trees was obtained from literatures on tropical timber (Sudou, 1970; PROSEA, 1994). For some species, wood density data were unavailable (S. acuta, S. andulensis, S. biawak, S. crassa, S. laxa and S. pilosa), so the average value for the timber group to which the species belonged was applied. A timber group is an arbitrary classification of dipterocarps used by foresters for convenience of trading, and the variation in WD within each group is not very large. Seed volumes (SV) were calculated as an ellipsoid of revolution from the nut length and width data obtained from literatures on dipterocarps (Ashton, 1964; Sarawak Forest Department, 1968).

315

Data analysis

The allometric relationships between the parts of an individual plant were used to analyse the structural properties and growth patterns of saplings. All parameters except branch mass (BM) were log-transformed before regression. For the interspecific comparison of regression lines, standardized major axis estimations and ANCOVA-like comparisons were applied using (S)MATR (Version 1; Falster et al., 2003). Standard regression methods were used to predict the dependent variables. Total dry mass (TM) was adopted as the independent variable in all regressions, and root dry mass (RM), stem dry mass (SM), leaf dry mass (LM), total leaf area (LA), projected crown area (CA), crown volume (CV), height (H) and stem diameter at ground height (D) were selected as the dependent variables to obtain eight regression lines. Quadratics were applied to the relationship of Lmax (maximum individual leaf area) to TM in the nonlinear regression, because in some species the increase in leaf size diminished with growth. To estimate BM, a nonlinear regression was applied following the general model of Klinkhamer et al. (1992), Y = a(X – c)b, because small saplings lacked branches. The relation of Lmax on TM and BM on TM was excluded from the ANCOVA-like test by (S)MATR, as there was no way to statistically compare such a complicated regression. Significant differences were found not only in the intercepts but also in the slopes; therefore, for further analyses, estimated values were obtained instead of the direct comparison of the regression parameters. Each sapling characteristic that was treated as a dependent variable was estimated against a TM of 5 and 30 g, which approximated the lower and upper limits of the common range of TM for the 18 species. Although regression will sometimes fail to properly predict values at the extreme of the range, these two values are actually not at the extreme of regression for most species. Even for the species for which prediction was done at the extreme, the high correlation coefficients indicated that the predicted values were reliable. Pairwise simple correlations and principal component analysis (PCA) were performed on the log-transformed estimated values to clarify the correlations between any pair of structural traits and the constellations among traits by uniting parameters into some principal components (PC). To incorporate BM, which included zero values, into the analysis, 01 was added to all values before log-transformation. Although other values such as 001 were tested, the resultant pattern was robust. In the PCA, factors with an absolute loading larger than 0468 were considered, which corresponds to a percentage of risk 80 % of the total variance was explained by three principal components. The factor loadings of the first and second principal component are

summarized as scatter plots in Fig. 2. At 5 g, the first axis, whose contribution was 34 %, had significant positive correlations with SM and H and negative correlations with LM, LA and CA. This axis summarizes a trade-off relationship between investment in leaves and investment in the trunk (and consequently, attained height). The second axis, which explained 30 % of the variance, was positively correlated with CV, BM, H and CA and negatively correlated with D, Lmax and RM. This axis contrasts species with wellbranched, large crowns consisting of small leaves supported by a thin trunk with species with small crowns consisting of large leaves supported by a thick trunk. The third axis, which contributed 17 %, showed a strong negative correlation with RM and a weak positive correlation with SM and Lmax (data not shown). This axis mainly represents variation in allocation to the roots. At 30 g, PC1 was relatively dominant, explaining 39 % of the total variance. The axis was correlated positively with CV, CA, BM and LA and negatively with Lmax, D and RM. This axis is very similar to the second axis at 5 g and summarizes the trade-off relationship between a large crown consisting of small leaves and a small crown consisting of large leaves. PC2, which accounted for 26 % of the total variance, had positive correlations with LM, LA and D. The interspecific variation in the investment in leaves is summarized in this axis. In addition, PC3, whose contribution was 22 %, had positive correlations with SM and H and a negative correlation with RM, which represented a trade-off relationship between the attained height and the allocation to roots at this ontogenetic stage. Correlation between structural characters and indices of shade tolerance

The Amax of saplings (approx. 1 m high) growing in nongap sites ranged from 26 (S. exelliptica) to 46 mmol m2 s1 (S. amplexicaulis) (Table 3b). Wood densities of these species obtained from the literatures ranged from 050 (S. acuta and S. cf. anculensis) to 097 (S. scrobiculata). These two indices of shade tolerance showed a significant negative correlation (r = 071, P < 0001; data not shown), suggesting their validity as an index of shade tolerance through the tree life history. Furthermore, these values generally support the empirical knowledge that some species of section Brachyptera, Mutica, Ovalis and Pachycarpae are relatively light-demanding and species of section Shorea are shade tolerant. Seed size, which represents a different aspect of the regeneration strategy, was not strongly correlated with the other two indices. The results of the correlation test between the PCA score of each species at both stages and three indices of regeneration strategy are shown in Table 5. At 5 g, only the correlation between PC1 and TD was marginally positive (r = 062, P = 006 after the sequential Bonferroni correction). At 30 g, PC1 had a significant negative correlation with seed size (r = 063, P < 005). The negative correlation between PC2 and TD was also significant at this stage (r = 067, P < 005), and the positive correlation between PC2 and Amax was relatively strong, although it was not significant (r = 054, P = 011). The positive correlation

317

Aiba and Nakashizuka — Sapling Structure and Function in Shorea

T A B L E 3. Structural characteristics of the saplings of 18 Shorea species growing in the shaded forest understorey (a) Size = 5 g Species

Section

RM

SM

BM

LM

LA

Lmax

CA

CV

H

D

S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S.

Mutica Anthoshorea Pachycarpae Brachyptera Pachycarpae Shorea Shorea Mutica Shorea Shorea Richetioides Mutica Anthoshorea Ovalis Mutica Pachycarpae Shorea Brachyptera

1.3 1.7 1.4 1.2 1.3 1.5 1.1 1.0 1.2 1.4 1.4 1.6 1.3 1.1 1.3 1.1 1.6 1.1

1.8 1.7 1.7 1.6 1.3 2.4 2.5 1.9 2.2 2.0 1.5 1.5 1.5 2.0 1.8 2.2 1.5 1.8

0.0 0.0 0.0 0.3 0.0 0.0 0.0 0.2 0.0 0.0 0.0 0.0 0.0 0.0 0.3 0.0 0.3 0.0

1.7 1.5 1.4 1.7 2.0 1.1 1.1 1.7 1.3 1.4 1.8 1.7 2.0 1.5 1.4 1.4 1.4 1.9

2.1 1.7 1.9 2.4 2.7 1.9 1.6 2.6 1.8 2.0 1.9 2.1 3.7 2.3 2.0 2.0 2.3 3.4

0.49 0.65 0.53 0.28 0.71 0.40 0.34 0.27 0.33 0.54 0.69 0.43 0.74 0.39 0.15 0.52 0.25 0.71

4.6 3.6 3.9 5.2 5.0 3.5 3.7 6.5 3.7 4.3 3.1 5.0 5.9 6.3 4.7 4.1 6.5 6.6

1.5 1.1 1.4 3.7 5.7 2.3 2.5 7.9 2.9 2.4 2.2 1.8 2.0 1.9 4.5 1.9 4.2 7.0

37 30 37 37 40 47 56 56 57 40 35 29 22 35 51 48 35 43

4.5 4.8 4.6 4.4 4.5 4.1 4.1 4.1 4.1 4.5 4.4 4.7 4.7 4.8 3.9 5.1 3.8 4.9

acuta agami amplexicaulis cf. andulensis beccariana biawak crassa curtisii exelliptica falciferoides ssp. glaucescens laxa macroptera ssp. baillonii ochracea ovalis ssp. sarawakensis ovata pilosa scrobiculata smithiana

(b) Size = 30 g Species S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S. S.

acuta agami amplexicaulis cf. andulensis beccariana biawak crassa curtisii exelliptica falciferoides laxa macroptera ochracea ovalis ovata pilosa scrobiculata smithiana

Amax

WD

SV

RM

SM

BM

LM

LA

Lmax

CA

CV

H

D

4.2 3.6 4.6 3.7 4.4 3.2 2.9 4.1 2.6 3.9 3.4 4.0 3.6 3.7 3.4 4.1 3.3 3.8

0.50 0.67 0.55 0.50 0.66 0.95 0.95 0.65 0.96 0.89 0.74 0.57 0.59 0.59 0.76 0.54 0.97 0.51

6.3 2.4 12.1 0.4 16.4 0.7 5.2 0.5 0.8 1.8 11.4 1.4 0.7 1.9 0.4 3.4 0.3 3.1

6.1 10.1 6.5 6.2 5.7 8.5 5.8 5.1 6.3 7.1 7.6 7.4 10.1 5.5 6.6 5.9 8.2 4.8

13.6 12.0 14.1 12.5 17.5 13.2 16.3 11.7 14.1 13.6 13.3 12.6 11.0 13.0 12.9 16.5 10.9 14.6

1.3 1.2 0.0 2.2 0.0 2.1 2.1 3.1 1.8 0.0 0.2 1.1 0.0 1.5 2.1 0.5 3.5 0.4

8.9 7.7 8.3 8.5 8.7 5.2 6.2 9.4 6.9 7.2 7.6 8.0 8.4 9.0 6.9 7.1 7.4 10.3

8.6 7.8 10.0 11.9 10.5 8.7 7.3 12.2 9.4 10.1 7.2 8.7 11.2 13.5 10.7 8.9 12.2 16.5

1.42 1.36 1.87 0.60 1.98 0.59 0.92 0.36 0.74 1.86 1.43 1.22 2.48 0.97 0.27 1.02 0.41 1.78

14 12 15 24 14 18 15 31 20 16 11 18 17 28 26 15 31 25

20 14 23 53 31 36 23 97 36 29 21 31 18 31 63 26 75 81

86 73 93 91 111 92 118 116 108 95 105 83 57 78 112 110 88 105

8.1 8.4 7.7 7.6 8.3 7.3 7.3 7.4 7.5 8.2 7.3 8.5 9.0 8.7 7.2 8.8 6.9 8.8

The maximum net photosynthetic rate of the saplings, seed volumes and wood density are also shown. For the definitions and units, see Table 1.

between PC3 and seed size was also strong, although it was only marginally significant (r = 061, P = 005). DISCUSSION Allometry and structure of 18 Shorea species

In six of the eight allometric relationships tested, significant differences in slopes were found among the species (Table 2). This indicates that the order of species for certain traits such as RM, SM, CA, CV and H per TM can shift with growth. Spearman’s rank correlation coefficients for these characters at two ontogenetic stages were 083, 037, 069, 076 and 083, respectively, indicating that a certain level of ontogenetic shift exists among these species. Such ontogenetic drift in plant structure is a general phenomenon

in juvenile trees as well as in herbs (Poorter and Pothmann, 1992; Walters et al., 1993b; Poorter, 1999; Walters and Reich, 1999). Predicted values of the morphological traits (Table 3) were comparable to values previously reported for saplings in tropical rain forests (Kohyama and Hotta, 1990; Poorter, 2001), temperate evergreen forests (Kohyama, 1987; Lusk, 2002), and temperate deciduous forests (Cao and Ohkubo, 1998; DeLucia et al., 1998; Takahashi et al., 2001). However, the variations among Shorea species in this study were unexpectedly high, given that all of the previous studies had compared taxonomically or ecologically diverse species. For example, Lusk (2002) compared eight species belonging to various functional groups (e.g. shrubs and tall trees, pioneer and climax species) and reported leaf mass ratios ranging from 011 to 037 and leaf area ratios ranging from 308 to 567 cm2 g1.

Aiba and Nakashizuka — Sapling Structure and Function in Shorea

318

T A B L E 4. Pearson’s correlation coefficients for every pair of ten morphological traits RM RM SM BM LM LA Lmax CA CV H D

SM 0.37

0.55 0.13 0.42 0.42 0.16 0.35 0.49 0.68* 0.05

0.30 0.11 0.22 0.30 0.41 0.19 0.61 0.18

BM 0.07 0.19 0.24 0.00 0.79*** 0.54 0.45 0.20 0.50

LM

LA

0.15 0.78*** 0.05

0.31 0.44 0.12 0.76**

0.59 0.31 0.22 0.20 0.17 0.50

Lmax 0.00 0.01 0.58 0.21 0.19

0.12 0.77*** 0.69* 0.08 0.24

0.67 0.65 0.38 0.70**

CA

CV

H

D

0.33 0.31 0.44 0.49 0.77*** 0.08

0.48 0.09 0.54 0.28 0.49 0.17 0.55

0.54 0.68* 0.18 0.54 0.37 0.21 0.18 0.47

0.07 0.14 0.62 0.40 0.29 0.71** 0.05 0.40 0.46

0.86*** 0.09 0.24

0.44 0.36

0.46

The results for 5-g plants are shown above the diagonal, and the results for 30-g plants are shown below the diagonal. ***P < 001, **P < 005, *P < 01 after a sequential Bonferroni correction.

101

A

BM (g)

TM (g)

100·50

100·25

B

100

R2 = 0·61 100·00 0·0 10 101·6

100·2 LM (g)

R2 = 0·62 10–1 –0·8 10

100·4

101·0

C

100·0 Lmax (102 cm2)

100·8

D

D (mm)

CA (102 cm2)

R2 = 0·59

101·3

101·0 0·8 10

101·1 LA (102 cm2)

101·4

R2 = 0·49 100·9

100·8 –0·8 10

100 Lmax (102 cm2)

100·8

F I G . 1. Representative correlations between structural traits at each ontogenetic stage: (A) trunk mass plotted against leaf mass at 5 g; (B) branch mass plotted against individual leaf size at 30 g; (C) crown area plotted against leaf area at 30 g; (D) stem diameter at ground height plotted against individual leaf size at 30 g. Note that all axes are log-scaled.

These ranges were only slightly wider for leaf mass ratio (017–034, re-calculated from Table 3), and even narrower for leaf area ratio (240–550 cm2 g1) than the ranges of values shown in the present study, although only non-pioneer species within a single genus were selected. In contrast, the variation in the leaf area index (total leaf area per crown projection area), which represents the extent of self-shading of crown architecture, was relatively small compared with the value for neotropical species reported by Sterck et al. (1999). The leaf area index in the 18 species in the present study ranged from 039 to 074, while Sterck et al. (1999) found a range from 028 to 095 in comparing

only four species. This may be partly because they selected species with compound leaves, which were excluded from the present study. Multiple correlation of structural traits

An allocation trade-off between SM and LM was clear at 5 g (Table 4 and Fig. 1), which suggests that the 18 species are diverse along the axis representing the functional tradeoff between height growth per mass and assimilation potential at this ontogenetic stage. However, this trend no longer existed at 30 g, and the weak correlation among allocation

319

Aiba and Nakashizuka — Sapling Structure and Function in Shorea 5g

30 g

PC2 1 CV

PC2 1 WL

WB

LA

D H

CA

Lmax

0·5

0·5 CA

LA WT –1

CV

WT

WL

PC1 –0·5

0·5

1

PC1 –1

–0·5

0·5 H

–0·5 D

WR

WR

1 WB

–0·5

Lmax –1

–1

F I G . 2. Two-dimensional scatter diagrams of the factor loadings of the two major principal components in the two ontogenetic stages for ten structural traits in 18 Shorea species.

T A B L E 5. Pearson’s correlation coefficient between the three major principal component scores and the three indices of shade tolerance PC1 5g Amax WD SV 30 g Amax WD SV

PC2

0.52 0.62* 0.02

0.22 0.23 0.56

0.23 0.19 0.63**

0.58 0.67** 0.16

PC3

0.07 0.23 0.31 0.09 0.03 0.61*

**P < 005, *P < 01 after a sequential Bonferroni correction.

traits indicates that interspecific differences in allocation cannot be explained by any single trade-off axis at this stage. The strong correlation between LM and LA at 5 g means that LA is mainly determined by leaf mass ratio at this stage. Conversely, at the later stage, the difference in specific leaf area also seems important because the correlation is not as strong at 30 g. A positive correlation between LA and CA was found at both stages. As the incremental rate is larger for CA for interspecific comparisons (in case of 30 g, logCA = 11logLA+013), leaf area index decreases in large-crowned leafy species. A positive correlation between Lmax and D was also observed at both stages. This relationship seems analogous to Corner’s Rules, which predict a correlation between individual leaf size and twig diameter (e.g. Corner, 1949; Cornelissen, 1999; Westoby and Wright, 2003), although the functional or physiological meaning of this is not clear. At 30 g, the correlation between BM and Lmax was also strong. It is reasonable that species with smaller leaves need to allocate more to branches to increase leaf area, while avoiding self-shading. The correlation between CA and CV was very strong at 30 g, but not as

strong at 5 g. This is because at a younger stage, several large leaves arranged in a plane or small leaves arranged three-dimensionally can realize a large CA at that stage, while at a later stage, a large CA cannot be achieved without branches. The correlation between allocation (RM, SM, BM and LM) and architecture (CA, CV, H and D) was generally not as strong, suggesting that architecture is not as strictly determined by biomass allocation. For example, H was not dependent on BM, while it was partly dependent on SM, at both ontogenetic stages, although it has been hypothesized that the allocation to the lateral expansion of branches reduces the attained height for a given mass. That is, since the trade-offs between SM and other organs (LM at 5 g and RM at 30 g) were very strong, the trade-off between SM and BM was negligible. The PCA satisfactorily summarized the structural variation of these 18 species into three principal components at both stages (cumulative contribution >80 %). At 5 g, the trade-off between height growth and investment in leaves, which results from an allocation trade-off between SM and LM, and a gradient from species with a wellbranched, three-dimensional crown supported by a thin trunk to species with a crown consisting of several large leaves without branches supported by a thick trunk, explained similar portions of the total morphological variation. The relatively unimportant third axis mainly represents the independence of the allocation to roots from these two major axes. At 30 g, the first axis summarizes the tradeoff relationship between well-branched, large crowns consisting of small leaves and less-branched, small crowns consisting of large leaves. This axis is relatively independent of allocation traits other than BM, indicating that crown architecture (CA and CV) and allocation do not converge on a single axis, as has been reported in some studies that have compared species from contrasting groups. The other two axes account for the interspecific difference in allocation. The second axis represents the difference in the investment

320

Aiba and Nakashizuka — Sapling Structure and Function in Shorea

to leaves and the third axis shows the difference in SM and resultant height growth efficiency. Sapling structure and regeneration strategy

At 5 g, the correlation between the PCA score of any axis and the three indices of shade tolerance among 18 species was not very strong. Only the positive correlation between PC1 score and TD was marginally significant. This correlation means that the trade-off between height growth efficiency and carbon gain co-varies with growth speed or shade tolerance. Since light-demanding species have a larger leaf area than shade-tolerant species, the result that fast growers have a larger LM and LA, which enable them to grow fast and be competitive in high-light conditions, agrees with previous studies. Conversely, we found that slow growers had a larger investment in the trunk and greater stature per unit mass. With these traits, slow growers might escape from the shaded understorey in relatively short period. In addition to this ecological reason, their dense wood may physically enable them to be slim compared with fast growing, light-wood species. Moreover, at 30 g, the correlation between PC2 score and WD was significant, suggesting that fast growers are characterized by a larger LM and LA, even within a functional group. The correlation between PC2 score and Amax was also relatively strong, although not significant at this stage, further supporting this view. At this stage, LM and LA were in weak trade-off with RM (Fig. 1 and Table 4); hence, slow growers were characterized by a large RM and small LM and LA. This tendency is also reasonable, as the large allocation to roots, especially to non-structural carbon, which enables juveniles to temporarily tolerate a negative carbon balance arising from herbivory, physical damage or shading by competitors (Canham et al., 1999), is a strategy that is advantageous under low-light conditions (Kobe, 1997). The score of first axis at 30 g, which mainly summarizes the trade-off in crown architecture and leaf size, showed a significant negative correlation with seed size. As previous studies have indicated that shade-tolerant or ‘pessimistic’ species (species which can tolerate waiting longer for a gap to form) are characterized by well-branched, larger crowns, while light-demanding or ‘optimistic’ species (species supposing a shorter time before gap formation) have less-branched, small crowns (Kohyama 1987, 1991), it was expected that species with a higher wood density, lower maximum photosynthetic rate and larger seeds would have larger crowns and vice versa. However, the present results show that in the 18 species studied, crown architecture (and leaf size) is virtually independent of wood density and photosynthetic rate, and the trend for seed size was contrary to the prediction. Small-seeded species formed large crowns and large-seeded species formed small crowns. The present simple correlation analysis showed that H per unit mass is not affected by crown size or the investment in branches (Table 4), while leaf area index, an index of self shading, was dependent on crown size, as described above. This independence of H per unit mass from crown architecture suggests that the functional trade-off between height growth and light capturing may not exist in these 18 species.

Rather, also this trend seems similar to Corner’s Rules, i.e. small-leaved species tend to have small fruit and branch more frequently than large-leaved species. In a comparison of 17 Acer species, Ackerly and Donoghue (1998) found that the crown architecture of saplings was independent of Corner’s Rules through evolution. However, the strong correlation among sapling architecture, leaf size, and seed size found in the present study is thought to be rather general, because a tight relationship between leaf size and crown architecture has been reported (Kohyama, 1987; Cao and Ohkubo, 1998; King, 1998, 1999), as well as a correlation between leaf size and seed size (Cornelissen, 1999; Westoby and Wright, 2003). The strength and validity of the relationship may be size dependent; hence, the difference may be due to the inclusion of larger saplings in the study by Ackerly and Donoghue (1998). In either case, the ecological merit of such a correlation among crown size, leaf size and seed size remains enigmatic, primarily because Corner’s Rules are empirical rules. In addition, at 30 g, the positive correlation between PC3 score and seed size was marginally significant. This relationship can be understood as poorly dispersed large-seeded species require efficient height growth to escape shaded conditions. Overall, the present results suggest that the allocation trade-off between leaves and other organs, which co-varies with wood density and to a certain extent with photosynthetic capacity, accounts for the difference in regeneration strategy, even in congeneric non-pioneer species. On the other hand, the correlation between the crown architecture, which is relatively independent of biomass allocation, and shade tolerance seems to be weak in closely related species, though it is obvious that species of different functional groups are characterized by specific crown architecture (King, 1990). Further comparative studies of the structure of closely related species, as well as a discussion of the results with respect to growth dynamics, are required to clarify the generality of these trends and to understand how tree species with a variety of structural traits compete and coexist. ACKNOWLEDGEMENTS We thank Dr T. Ichie for his technical advice on measuring photosynthesis and Mr S. Tan for his help with the identification of species. Lourens Poorter and an anonymous reviewer provided helpful comments on earlier version of this manuscript. This study was partly supported by a Research Institute for Humanity and Nature project (P2-2) and a grant for Biodiversity Research for the 21st Century COE (A14). LITERATURE CITED Ackerly DD, Donoghue MJ. 1998. Leaf size, sapling allometry, and Corner’s rules: phylogeny and correlated evolution in maples (Acer). American Naturalist 152: 767791. Ashton MS. 1995. Seedling growth of cooccurring Shorea species in the simulated light environments of a rain-forest. Forest Ecology and Management 72: 112. Ashton PS. 1964. Manual of the dipterocarp trees of Brunei state. Oxford: Oxford University Press.

Aiba and Nakashizuka — Sapling Structure and Function in Shorea Ashton PS. 1991. Toward a regional classification of the humid tropics of Asia. Tropics 1: 112. Barker MG, Press MC, Brown ND. 1997. Photosynthetic characteristics of dipterocarp seedlings in three tropical rain forest light environments: a basis for niche partitioning? Oecologia 112: 453463. Canham CD, Kobe RK, Latty EF, Chazdon RL. 1999. Interspecific and intraspecific variation in tree seedling survival: effects of allocation to roots versus carbohydrate reserves. Oecologia 121: 111. Cao KF, Ohkubo T. 1998. Allometry, root/shoot ratio and root architecture in understory saplings of deciduous dicotyledonous trees in central Japan. Ecological Research 13: 217227. Cornelissen JHC. 1999. A triangular relationship between leaf size and seed size among woody species: allometry, ontogeny, ecology and taxonomy. Oecologia 118: 248255. Corner E. 1949. The Durian theory, or the origin of the modern tree. Annals of Botany 13: 368414. Davies SJ. 1998. Photosynthesis of nine pioneer Macaranga species from Borneo in relation to life history. Ecology 79: 22922308. DeLucia EH, Sipe TW, Herrick J, Maherali H. 1998. Sapling biomass allocation and growth in the understory of a deciduous hardwood forest. American Journal of Botany 85: 955963. Falster DS, Warton DI, Wright IJ. 2003. (S)MATR: standardised major axis tests and routines. Version 1.0 (http://www.bio.mq.edu.au/ ecology/SMATR). Givnish T. 1988. Adaptation to sun and shade: a whole-plant perspective. Australian Journal of Plant Physiology 15: 6392. Grubb P. 1977. The maintenance of species richness in plant communities: the importance of the regeneration niche. Biological Review 52: 107145. Horn HS. 1971. The adaptive geometry of trees. Princeton: Princeton University Press. King DA. 1990. Allometry of saplings and understorey trees of a Panamanian forest. Functional Ecology 4: 2732. King DA. 1998. Influence of leaf size on tree architecture: first branch height and crown dimensions in tropical rain forest trees. Trees – Structure and Function 12: 438445. King DA. 1999. Tree architecture in relation to leaf dimensions and tree stature in temperate and tropical rain forests. Journal of Ecology 87: 10121024. Kitajima K. 1994. Relative importance of photosynthetic traits and allocation patterns as correlates of seedling shade tolerance of 13 tropical trees. Oecologia 98: 419428. Klinkhamer PGL, Meelis E, Dejong TJ, Weiner J. 1992. On the analysis of size-dependent reproductive output in plants. Functional Ecology 6: 308316. Kobe RK. 1997. Carbohydrate allocation to storage as a basis of interspecific variation in sapling survivorship and growth. Oikos 80: 226233. Kohyama T. 1987. Significance of architecture and allometry in saplings. Functional Ecology 1: 399404. Kohyama T. 1991. A functional-model describing sapling growth under a tropical forest canopy. Functional Ecology 5: 8390. Kohyama T, Hotta M. 1990. Significance of allometry in tropical saplings. Functional Ecology 4: 515521. Lusk CH. 2002. Leaf area accumulation helps juvenile evergreen trees tolerate shade in a temperate rainforest. Oecologia 132: 188196. Montgomery RA, Chazdon RL. 2002. Light gradient partitioning by tropical tree seedlings in the absence of canopy gaps. Oecologia 131: 165174. Muller-Landau HC. 2004. Interspecific and inter-site variation in wood specific gravity of tropical trees. Biotropica 36: 2032. Nakagawa M, Tanaka K, Nakashizuka T, Ohkubo T, Kato T, Maeda T, et al. 2000. Impact of severe drought associated with the 1997–1998 El Nino in a tropical forest in Sarawak. Journal of Tropical Ecology 16: 355367.

321

Newman MF, Burgess PF, Whitmore TC. 1996. Manuals of dipterocarps: Borneo island light hard woods. Edinburgh/Jakarta: Royal Botanic Garden Edinburgh/Center for International Forestry Research. Poorter H, Pothmann P. 1992. Growth and carbon economy of a fastgrowing and a slow-growing grass species as dependent on ontogeny. New Phytologist 120: 159166. Poorter L. 1999. Growth responses of 15 rain-forest tree species to a light gradient: the relative importance of morphological and physiological traits. Functional Ecology 13: 396410. Poorter L. 2001. Light-dependent changes in biomass allocation and their importance for growth of rain forest tree species. Functional Ecology 15: 113123. Poorter L, Rose SA. 2005. Light-dependent changes in the relationship between seed mass and seedling traits: a meta-analysis for rain forest tree species. Oecologia 142: 378387. PROSEA 1994. Plant resources of South-East Asia. 5 (1) Timber trees: major commercial timbers. Bogor/Wageningen: Prosea Foundation/ Pudoc-DLO. Reich PB, Ellsworth DS, Uhl C. 1995. Leaf carbon and nutrient assimilation and conservation in species of differing successional status in an oligotrophic Amazonian forest. Functional Ecology 9: 6576. Sarawak Forest Department. 1968. A manual of the dipterocarp trees of Brunei state and of Sarawak. Kuching: Borneo Literature Bureau. Sterck FJ, Clark DB, Clark DA, Bongers F. 1999. Light fluctuations, crown traits, and response delays for tree saplings in a Costa Rican lowland rain forest. Journal of Tropical Ecology 15: 8395 Still MJ. 1996. Rates of mortality and growth in three groups of dipterocarp seedlings in Sabah, Malaysia. In: Swaine MD, eds. The ecology of tropical forest tree seedlings. Parthenon: UNESCO, 315–332. Sudou S. 1970. Tropical timbers Tokyo: Chikyu. Suzuki E. 1999. Diversity in specific gravity and water content of wood among Bornean tropical rainforest trees. Ecological Research 14: 211224. Takahashi K, Seino T, Kohyama T. 2001. Responses to canopy openings in architectural development of saplings in eight deciduous broadleaved tree species. Canadian Journal of Forest Research 31: 13361347. Thomas SC. 1996. Asymptotic height as a predictor of growth and allometric characteristics in Malaysian rain forest trees. American Journal of Botany 83: 556566. Thomas SC. 2003. Comparative biology of tropical trees: a perspective from Pasoh. In: Okuda T, Manokaran N, Matsumoto Y, Niiyama K, Thomas SC, Ashton, PS, eds. Pasoh: ecology of a lowland rain forest in southeast Asia. Tokyo: Springer, 181–208. Turner IM. 1990. The seedling survivorship and growth of 3 Shorea species in a Malaysian tropical rain-forest. Journal of Tropical Ecology 6: 469478. Walters MB, Reich PB. 1999. Low-light carbon balance and shade tolerance in the seedlings of woody plants: do winter deciduous and broadleaved evergreen species differ? New Phytologist 143: 143154. Walters MB. Reich PB. 2000. Seed size, nitrogen supply, and growth rate affect tree seedling survival in deep shade. Ecology 81: 18871901. Walters MB, Kruger EL, Reich PB. 1993a. Growth, biomass distribution and CO2 exchange of northern hardwood seedlings in high and low light—relationships with successional status and shade tolerance. Oecologia 94: 716. Walters MB, Kruger EL, Reich PB. 1993b. Relative growth-rate in relation to physiological and morphological traits for northern hardwood tree seedlings—species, light environment and ontogenic considerations. Oecologia 96: 219231. Westoby M, Wright IJ. 2003. The leaf size–twig size spectrum and its relationship to other important spectra of variation among species. Oecologia 135: 621628. Yamakura T, Kanzaki M, Itoh A, Ohkubo T, Ogino K, Chai EOK, et al. 1995. Topography of a large-scale research plot established within a tropical rain forest at Lambir, Sarawak. Tropics 5: 4156.