Keywords Teak Seed provenance Growth characteristics Stress-wave velocity Pilodyn penetration. Introduction

Transcription

1 J Wood Sci (2013) 59: DOI /s NOTE Variation in tree growth characteristics, stress-wave velocity, and Pilodyn penetration of 24-year-old teak (Tectona grandis) trees originating in 21 seed provenances planted in Indonesia Fanny Hidayati Futoshi Ishiguri Kazuya Iizuka Kazuko Makino Yuya Takashima Sri Danarto Widaryanti Wahyu Winarni Denny Irawati Mohammad Na iem Shinso Yokota Received: 14 May 2013 / Accepted: 30 August 2013 / Published online: 22 September 2013 Ó The Japan Wood Research Society 2013 Abstract Growth characteristics [stem diameter (D), tree height (H), and bole volume (V)], stress-wave velocity (SWV), and Pilodyn penetration (P) were measured for 21 seed provenances of 24-year-old teak trees planted in Indonesia to characterize variation in tree growth characteristics, SWV, and P among seed provenances. Broadsense heritability and correlations between the measured characteristics were also determined. Significant differences for all measured characteristics were observed among provenances, indicating that these characteristics are genetically controlled. Broad-sense heritabilities of growth characteristics, SWV, and P are moderate values. These results indicate potential for improving growth characteristics and wood properties of teak trees with the help of breeding programs. Highly significant positive correlations were observed among the growth characteristics, suggesting that they are closely related. In contrast, no significant correlations were observed between the growth characteristics and SWV, indicating that they are Part of this work was presented at the 63rd Annual Meeting of the Japan Wood Research Society in Morioka, March F. Hidayati United Graduate School of Agricultural Science, Tokyo University of Agriculture and Technology, Fuchu, Tokyo , Japan F. Hidayati F. Ishiguri (&) K. Iizuka K. Makino Y. Takashima S. Yokota Faculty of Agriculture, Utsunomiya University, Utsunomiya , Japan ishiguri@cc.utsunomiya-u.ac.jp F. Hidayati S. Danarto W. W. Winarni D. Irawati M. Na iem Faculty of Forestry, Gadjah Mada University, Yogyakarta 55281, Indonesia independent. We conclude that mechanical properties are also important criteria for selecting plus trees in tree breeding programs. Principal component analysis revealed that seed provenances from Indonesia (Bangilan, Deling, and Randublatung) and India (Malabar and Central Province) have high scores of growth characteristics and SWV. Keywords Teak Seed provenance Growth characteristics Stress-wave velocity Pilodyn penetration Introduction Provenance trials of growth characteristics in teak (Tectona grandis) trees have been conducted with the purpose of improving the productivity and rotation period of plantations by selecting seed or sources of cutting. The results have revealed that growth characteristics and stem quality including stem form differ among seed provenances [1, 2]. However, these studies did not include wood quality traits. In our previous research, growth characteristics [stem diameter (D), tree height (H), and bole volume (V)], stresswave velocity (SWV), and Pilodyn penetration (P) were significantly different among 15 clones planted at two different sites in Indonesia [3]. Based on these results, we proposed that the wood properties of teak trees can be improved by implementing appropriate tree breeding programs [3]. However, in the test sites, it is very difficult to get the samples by felling a large number of trees. The stress-wave method is a well-known non-destructive method for evaluating wood quality, especially mechanical properties. There is a significant, positive relationship between the SWV of standing trees and the Young s modulus of logs or lumbers [4 7]. In addition, P is also a

2 J Wood Sci (2013) 59: Table 1 Mean values of growth characteristics, stress-wave velocity (SWV), and Pilodyn penetration (P) in each seed provenance Code Seed provenance Number of trees D (cm) H (m) V (m 3 ) SWV (km/s) P (mm) 1 Cepu, Java, Indonesia Bangilan, Java, Indonesia Muna, Southeast Sulawesi, Indonesia Ngliron, Java, Indonesia Margasari, Java, Indonesia Ponorogo, Java, Indonesia Gundih, Java, Indonesia Deling, Java, Indonesia Java, Indonesia Blora, Java, Indonesia Pati, Java, Indonesia Randublatung, Java, Indonesia Soe, Nusa Tenggara, Indonesia Malabar, India Central Province, India Godavari, India Thailand Kay, Vietnam Kouse, Vietnam Kouai, Vietnam Myanmar Seed provenance of code 9 is unknown D stem diameter, H tree height, V bole volume, SWV stress-wave velocity, P Pilodyn penetration useful non-destructive test to evaluate basic density [6, 8, 9]. In Indonesia, seedlings for establishing the commercial teak plantations are mostly obtained by sexual propagation. In general, seed sources in seed production areas have been selected from plantations based on growth characteristics and stem quality form. The collected seeds are used for establishing commercial teak plantation in Java, Indonesia. Therefore, it is desirable to genetically improve seed quality to establish more productive teak plantations in Indonesia. In the present study, to characterize genetic variations in growth characteristics and wood properties in the teak trees, tree growth characteristics (D, H, and V), SWV, and P were investigated for trees originating in 21 seed provenances planted in Indonesia. Materials and methods The experimental site of this study was located in the Education Forest of Wanagama, Gadjah Mada University, Yogyakarta, Indonesia ( S E). The provenance trial site for teak was established in February The environmental conditions of the provenance trial site were as follows: average temperature, 27.7 C; annual precipitation, 1954 mm/year; relative humidity, %; altitude, 214 m above sea level; soil, Mediterranean; and topography, flat. Seeds were collected from 21 different provenances in Indonesia and foreign countries (Table 1). The seedlings were initially planted at m spacing. No fertilization or thinning treatments were applied. A total of 155 of 24-year-old trees were used. D was measured at 1.3 m above the ground using diameter tape (F10-02DM, KDS). H was measured using an altimeter (Vertex IV, Haglöf), and V was calculated from the equation of Monteuuis et al. [1]. SWV of the stem was measured using a commercial, hand-held stress-wave timer (FAKOPP Microsecond Timer, FAKOPP Enterprise) according to the method previously described [3, 5]. Start and stop sensors were set at 150 and 50 cm from ground level, respectively. The start sensor was hit with a small hammer to create a stress wave. After the stress wave was received by the stop sensor, the stress-wave propagation time between the two sensors was recorded. Eleven measurements of stress-wave propagation time were obtained for each tree, and the mean value was calculated for each tree. SWV was calculated by dividing the distance between sensors (100 cm) by the averaged stress-wave propagation time. P was measured using a Pilodyn tester (strength of

3 514 J Wood Sci (2013) 59: spring, 6 J; diameter of pin, 2.5 cm, Proceq) at 1.3 m above ground level at three positions for each tree [3, 7, 8]. The bark at the measuring positions was removed before measuring P. The mean value of three P measurements was recorded for each tree. Analyses of variance (ANOVA) were performed to evaluate the differences in growth characteristics, SWV, and P among seed provenances. The following ANOVA model was used for each character: Y ij ¼ l þ S i þ C j þ e ij ð1þ where Y ij denotes the response of the jth seed provenance in the ith replication, l the overall mean, S i the site effect at the ith replication, C j the genetic effect of the jth seed provenance, and e ij is the error within Y ij. Broad-sense heritabilities of growth characteristics, SWV, and P were estimated using the following formula proposed by Falconer [10] to assess the magnitude of genetic effects: H 2 ¼ r 2 c = r2 c þ r2 e ð2þ where H 2 denotes the broad-sense heritability, r 2 c the variance component of seed provenance, and r 2 e denotes the variance component of the environment factor. Principal component analysis (PCA) was performed using R software ( to select seed provenances with good performance for both growth characteristics and wood properties. Results and discussion The mean values of D, H, V, SWV, and P in each seed provenance were cm, m, m 3, km/s, and mm, respectively (Table 1). Table 2 shows the result of ANOVA and broad-sense heritability of the teak trees from 21 seed provenances. Danarto and Hardiyanto [11] reported that 142 teak families originating in Java and East Nusa Tenggara, Indonesia, and planted in East Java indicated significant differences in D. Chaix et al. [2] reported that D, H, and mortality rate demonstrated significant differences among 41 different seed origins planted in Taliwas, Malaysia. Significant differences in growth characteristics (D, H, V, mortality rate, and fork height) were also observed among teak trees originating from 42 different genetic sources comprising 26 open-pollinated families from seedlings planted in Sabah, Malaysia [1]. In addition, we previously observed significant differences in D, H, and V among 15 clones of teak trees planted at two different sites in Indonesia [3]. The dynamic Young s modulus of teak trees planted in Costa Rica significantly differed among 20 clones [12]. Furthermore, SWV and P significantly differed among 15 clones of teak trees planted in Indonesia [3]. In the present study, the results for growth characteristics, SWV, and P obtained from trees originating in 21 different seed provenances were consistent with those obtained in previous studies. Therefore, it is considered that growth characteristics, SWV, and P are genetically controlled in teak. Danarto and Hardiyanto [11] reported the narrow-sense heritability of 0.23 for D in 142 families of 12-year-old teak trees planted in East Java, Indonesia. Monteuuis et al. [1] reported narrow-sense heritabilities of 0.24, 0.51, and 0.34 for D, H, and V, respectively, in 8.8-year-old trees from 42 different genetic sources comprising 26 openpollinated families from seedlings planted in Sabah, Malaysia. In addition, the narrow-sense heritabilities of D, H, and V in 3.5-year-old teak trees from 61 different families planted in Australia were 0.31, 0.22, and 0.29, respectively [13]. Chaix et al. [2] reported narrow-sense heritabilities of 0.46 and 0.76 for D and H, respectively, in 8.7-year-old teak trees from 26 clonal seed orchard families planted in Taliwas, Malaysia. In our previous study, the broad-sense heritabilities of D, H, V, SWV, and P in 15 clones of 12-year-old teak trees planted at two different sites in Indonesia were , , , , and , respectively [3]. Moya and Marín [12] reported the broad-sense heritability of the dynamic modulus of 0.34 in 10-year-old teak trees planted in Costa Rica. Moderate values for the broad-sense heritability of the measured characteristics in the present study indicate potential for improving growth characteristics and wood properties in teak with the help of breeding programs. Highly positive significant correlations were observed among D, H, and V in the present study (Table 3). Highly positive correlation coefficients of D and H have also been previously reported for teak trees [1 3, 13]. In addition, Table 2 Results of analysis of variance (ANOVA) and broadsense heritability of measured characteristics SD standard deviation, D stem diameter, H tree height, V bole volume, SWV stress-wave velocity, P Pilodyn penetration a Significance at the 1 % level Characteristic Mean SD Significance among seed provenances Variance of seed provenances Variance of environment D (cm) a H (m) a V (m 3 ) a SWV (km/s) a P (mm) a Heritability

4 J Wood Sci (2013) 59: Table 3 Correlation coefficients between growth characteristics, SWV, and P Characteristics D H V SWV P D H b V b b SWV ns ns ns P a ns ns ns Correlation coefficients were determined using mean values of each provenance D stem diameter, H tree height, V bole volume, SWV stress-wave velocity, P Pilodyn penetration, ns no significance a Significance at the 5 % level Table 4 Loading values from principal component analysis (PCA) of stem diameter (D), tree height (H), SWV, and P Characteristics PC 1 PC 2 D H SWV P D stem diameter, H tree height, SWV stress-wave velocity, P Pilodyn penetration, PC principal component b Significance at the 1 % level D and H have been reported to be positively correlated with V [1, 3]. These results indicate that growth characteristics are closely related to one another in teak trees. Furthermore, D is among the suitable criteria in teak breeding programs for selecting plus trees with high wood yield. No significant correlations were observed between growth characteristics and SWV in teak trees (Table 3), suggesting that SWV in teak is independent of growth characteristics. There is a significant, positive relationship between the SWV of standing trees and the Young s modulus of logs or lumbers [4 7]. Therefore, we concluded that mechanical properties such as Young s modulus are also important criteria for selecting the plus trees in teak breeding programs. In our previous study, we observed relatively high positive significant correlation between growth characteristics and P in teak trees planted at two different sites in Indonesia [3]. In the present study, a moderately significant positive correlation was observed between D and P (Table 3). It is well known that P is negatively correlated with wood density [6, 8, 9]. Therefore, it is considered that trees with good D may have lower wood density. Rao and Shashikala [14] reported that a negative significant correlation was found between growth rate (D) and basic density of Thithimathi clones of teak planted in India. However, further studies are required to clarify the relationship between D and P in teak tree. PCA is a multivariate technique that analyzes a data table in which observations are described by several intercorrelated quantitative dependent variables. Its goal is to extract the important information from the table, to represent it as a set of new orthogonal variables called principal components (PCs), and to display the pattern of similarity of the observations and of the variables as points in map. The correlation between a component and a variable estimates the information they share. In the PCA framework, Fig. 1 Principal component analysis (PCA) scores for 21 seed provenances of Tectona grandis. Closed circles and numbers indicate seed provenances. PC principal component this correlation is called a loading value [15]. As shown in Table 4, the absolute values of loading values for PC 1 were high in D and H, and loading values for PC 2 were high in SWV, indicating that growth characteristics and SWV contribute to PC 1 and PC 2, respectively. SWV of stems has been reported to be positively correlated with the Young s modulus of wood [4 7]. Thus, the Young s modulus contributes to PC 2. The seed provenances were well separated as shown in Fig. 1, indicating that PCA is a suitable analysis method for selecting the seed provenances with good performance in both growth characteristics and wood properties. In the present study, as per the results of PCA, some seed provenances demonstrated high scores for growth characteristics and Young s modulus, such as those from Indonesia [Bangilan (No. 2), Deling (No. 8), and Randublatung (No. 12)] and India [Malabar (No. 14) and Central Province (No. 15)]. Relatively high absolute value of loading value in PC 2 was observed for P, indicating that P also contributes to PC 2. P is closely related to the wood

5 516 J Wood Sci (2013) 59: density of stem [6, 8, 9]. Thus, seed provenances from Bangilan (No. 2) and Blora (No. 10) demonstrated good values for D and H, high values for SWV which are positively correlated with Young s modulus of wood, and low values for P (high basic density). Conclusion Tree growth characteristics (D, H, and V), SWV, and P were investigated in teak trees originating in 21 seed provenances planted in Indonesia, with the aim of characterizing the genetic variation in growth characteristics and wood properties. Growth characteristics (D, H, and V), SWV, and P varied among the seed provenances. The broad-sense heritabilities of the measured characteristics were moderate. Highly positive correlations were observed between the growth characteristics, but no significant correlations were observed between growth characteristics and SWV. Based on these results, we conclude that the growth characteristics and wood properties of teak plantations can be improved by implementing appropriate tree breeding programs in seed production areas. PCA results revealed that seed provenances from Indonesia (Bangilan, Deling, and Randublatung) and India (Malabar and Central Province) have high scores of growth characteristics and SWV. Acknowledgments The authors thank the Education Forest of Wanagama, Faculty of Forestry, Gadjah Mada University, Indonesia for permission to evaluate growth characteristics, SWV, and P of sample trees. References 1. Monteuuis O, Goh DKS, Garcia C, Alloysius D, Gidiman J, Bacilieri R, Chaix G (2011) Genetic variation of growth and tree quality traits among 42 diverse genetic origins of Tectona grandis planted under humid tropical conditions in Sabah, East Malaysia. Tree Genet Genomes 7: Chaix G, Monteuuis O, Garcia C, Alloysius D, Gidiman J, Bacilieri R, Goh DKS (2011) Genetic variation in major phenotypic traits among diverse genetic origins of teak (Tectona grandis L. f) planted in Taliwas, Sabah, East Malaysia. Ann For Sci 68: Hidayati F, Ishiguri F, Iizuka K, Makino K, Tanabe J, Marsoem SN, Na iem M, Yokota S, Yoshizawa N (2013) Growth characteristics, stress-wave velocity, and Pilodyn penetration of 15 clones of 12-year-old Tectona grandis trees planted at two different sites in Indonesia. J Wood Sci 59: Wang X, Ross RJ, McClellan M, Barbour RJ, Erickson JR, Forsman JW, McGinnis GD (2001) Nondestructive evaluation of standing trees with a stress wave method. Wood Fiber Sci 33: Ishiguri F, Eizawa J, Saito Y, Iizuka K, Yokota S, Priadi D, Sumiasri N, Yoshizawa N (2007) Variation in the wood properties of Paraserianthes falcataria planted in Indonesia. IAWA J 28: Wu SJ, Xu JM, Li GY, Risto V, Lu ZH, Li BQ, Wang W (2011) Estimation on basic density and modulus of elasticity of eucalypt clones in Southern China using non-destructive methods. J Trop For Sci 23: Ishiguri F, Matsui R, Iizuka K, Yokota S, Yoshizawa N (2008) Prediction of the mechanical properties of lumber by stress-wave velocity and Pilodyn penetration of 36-year-old Japanese larch trees. Holz Roh Werkst 66: Wu SJ, Xu JM, Li GY, Risto V, Lu ZH, Li BQ, Wang W (2010) Use of the Pilodyn for assessing wood properties in standing trees of Eucalyptus clones. J For Res 21: Raymond CA, MacDonald AC (1998) Where to shoot your Pilodyn: within tree variation in basic density in plantation Eucalyptus globulus and E. nitens in Tasmania. New For 15: Falconer DS (1989) Introduction to quantitative genetics, 3rd edn. Wiley, New York 11. Danarto S, Hardiyanto EB (2000) Results of the progeny test of teak at 12 years of age at Jember, East Java. In: Proceeding of third regional seminar on teak. Potentials and opportunities in marketing and trade of plantation teak: challenge for the new millennium. Yogyakarta, pp Moya R, Marín JD (2011) Grouping of Tectona grandis (L. f.) clones using wood color and stiffness. New For 42: Callister AN, Collins SL (2008) Genetic parameter estimates in a clonally replicated progeny test of teak (Tectona grandis Linn. f.). Tree Genet Genomes 4: Rao RV, Shashikala S (2005) Assessment of growth rate, basic density and heartwood content in selected teak clones of CSO, Thithimathi, Karnataka. In: Proceeding of the international conference on quality timber products of teak from sustainable forest management. Peechi, pp Abdi H, Williams LJ (2010) Principal component analysis. Wires Comput Stat 2: