SITE QUALITY CURVES FOR AFRICAN MAHOGANY PLANTATIONS IN BRAZIL

CURVAS DE CLASSIFICAÇÃO DE SÍTIO PARA PLANTIOS DE MOGNO AFRICANO NO BRASIL

ABSTRACT

Site quality estimation is an important tool in forest management since it is useful for modeling growth and yield for even-aged stands. Data from African mahogany (Khaya ivorensis A. Chev.) Brazilian plantations were used to develop a model to predict dominant height growth, comparing dynamic base-age invariant site index models with the guide curve method (static models). For the evaluation of the candidate models qualitative and quantitative criteria were used. We also verified the stability of the candidate models, preferring a model providing fewer site class changes when predicting site index from different ages. The Lundqvist-Korf function fitted with the guide curve method proved to be effective and accurate for site classification and dominant height predictions of African mahogany stands. The range of observed site index, at a reference age of 15, was between 17 and 33 meters.

Keywords: Forest management; Nonlinear equations; Dominant height; Site index

RESUMO

A estimativa da qualidade de sítio é uma importante ferramenta para o manejo florestal, uma vez que auxilia na modelagem do crescimento e produção de florestas equiâneas. Dados de plantios brasileiros de mogno africano (Khaya ivorensis A. Chev.) foram utilizados para avaliar modelos de crescimento em altura dominante comparando modelos dinâmicos invariantes com a idade padrão e modelos estáticos (método da curva guia). Para a avaliação dos modelos foram utilizados critérios quantitativos e qualitativos. Também foi verificada a estabilidade dos modelos candidatos, sendo preferido o modelo que resultou em menos mudanças de classes na predição da mesma parcela em diferentes idades. O modelo de Lundqvist-Korf com o método da curva guia se mostrou eficiente, sendo recomendado para a classificação de sítios e para a predição de altura dominante em povoamentos de mogno africano. A amplitude dos índices de sítio, na idade de referência de 15 anos, foi de 17 a 33 metros.

Palavras chave: Manejo florestal; Equações não lineares Altura dominante; Qualidade de sítio

INTRODUCTION

African mahogany, that includes the species Khaya ivorensis A. Chev., is a tree species that belongs to the botanical family Meliaceae and is well known in the international wood market (PINHEIRO et al., 2011). The recent cultivation of the species outside its native location is relatively new, starting in Brazil in 1976 at Pará state with five seedlings (FALESI; BAENA, 1999). These specimens yielded seeds which were used to start up plantations in other regions of the country (FRANÇA et al., 2016). The definition of an accurate prediction of site productivity is fundamental to predict timber yields and for meaningful simulation studies (VANCLAY, 1992; BRAVO; MONTERO, 2001). This concern is not new and the first land classification system of productivity is attributed to the roman Cato in 234-139 B.C. (TESCH, 1981) and several methods for evaluating forest site productivity have been studied (BURKHART; TOMÉ, 2012).

Many are the models based on the dominant height-age relationship applied to site index curves construction (SCOLFORO; MACHADO, 1988; DIÉGUEZ-ARANDA et al., 2006, DAVID et al., 2015), and those are the basis of most traditional stand-level forest management regimes (RENNOLLS, 1995). Among the techniques available to fit site index curves, most are based on three general methods (CLUTTER et al., 1983): guide curve; parameter prediction, and difference equation methods. Although there are many techniques available for site classification no, research has been published until now for African mahogany plantations in Brazil or in other parts of the world.

The guide curve represents the height development for the average site index in the data. Heights at all ages for all other site classes are typically assumed to be proportional to that of the guide curve (BURKHART; TOMÉ, 2012).

The difference equation approach can be applied to any height-age equation to produce families of anamorphic or polymorphic curves (BURKHART; TOMÉ, 2012). The height (hdom₂) at a future age must be expressed as a function of future age (t₂), current age (t₁) and current height (hdom₁), that is: hdom₂=(t₂, hdom₁, t₁). The formulation of a growth function as a difference equation can be obtained by substituting one parameter by its value computed from the initial condition (t₁, hdom₁). The substitution of the asymptotic parameter in the growth model produces anamorphic curves with multiple asymptotes and the substitution of any other parameter produces polymorphic curves with a common asymptote (PALAHÍ et al., 2004; TAMARIT-URIAS et al., 2014; ERCANLI et al., 2014, TEWARI et al., 2014).

The objective of this study was to develop a site index model for Khaya ivorensis A. Chev. stands in Brazil comparing the difference equation approach with the guide curve methodology.

MATERIAL AND METHODS

Data were gathered in plantations located in three Brazilian regions (Southeast, Midwest and North) comprising 150 permanent sample plots with different ages and remeasurement times, usually around one year interval (Table 1 and Figure 1).

Table 1
Characteristics of the sample plots of African mahogany plantations used in the study.

Region	Variable	Minimum	Mean	Maximum	SD
Southeast	t	1.1	3.7	7.2	1.6
(n=125)	hdom	1.7	10.9	23.0	4.9
	nm	2	4	6	2
Midwest	t	4.9	5.7	6.7	0.7
(n=20)	hdom	11.4	16.5	20.3	2.3
	nm	2	2	2	0
North	t	14.0	19.2	39.0	9.9
(n=5)	hdom	23.7	29.1	38.7	4.7
	nm	2	2	2	0

• Where: n = number of plots; t = age (years); hdom = dominant height (m); nm = number of measurements; SD = standard deviation.

FIGURE 1
Observed dominant height development of Khaya ivorensis A. Chev. stands in Brazil

Forest management and genetic material of the plantations across Brazil are similar, planted with a wide spacing (4x5 to 6x6 for plots with less than 8 years and 12x12 m for older plots), intensive fertilization regimes and in some cases irrigation (areas with precipitation less than 800 mm per year). The plot area varies with the spacing used (800 to 4,400 m²). The oldest plot (in the North region) is composed of three trees planted in 1976 (about 20 meters from each other), so the dominant height was computed as the average of those trees. It was decided to keep these trees in the data base because they are the oldest and highest African mahogany trees in Brazil, thus allowing a more realistic estimation of the asymptote parameter of the models.

Given the low density of the plantations (the majority of the plots have less than 300 trees per hectare), the traditional definition of dominant height (Assmann’s mean height of the 100 thickest trees per hectare) would require height measurement of an excessive number of trees per plot. Also, some of plots have less than 100 trees per hectare, making Assmann’s dominant height definition inviable. To avoid these problems, we defined the dominant height as being the mean height of the 30 thickest trees per hectare. Similar strategies have been adopted by Paulo et al. (2011) for Quercus suber stands in Portugal and Danquah (2012) for African mahogany species planted in Ghana, selecting 25 and 40 thickest trees per hectare, respectively, to calculate the average of dominant height.

Several candidate functions to estimate the dominant height (Table 2) were selected, a total of 434 height-age pairs were used to fit the dynamic functions and 584 observations to fit the base functions (static) for the guide curve method. All candidate dynamic functions have the property of being time invariant, ensuring that projections using different starting ages but with equal final ages are equal, as described by Palahí et al. (2004).

Table 2
Candidate models for dominant height estimative in function of age for Khaya ivorensis A. Chev. stands in Brazil.

Model	Base function (integral form)	Parameter related to site	Dynamic function
Lk	Lundqvist-Korf	k	$hdo{m}_2=A{\left(\frac{hdo{m}_1}{A}\right)}^{{(t_1/t_2)}^d}$
LA	$hdom=A{e}^{-k/t^d}$	A	$hdo{m}_2=hdo{m}_1{e}^{-k\left(\left(\frac{1}{t_1^d}\right)-\left(\frac{1}{t_2^d}\right)\right)}$
LG		A,k	$hdo{m}_2=e^{A-\left(k/t_2\right)+\left[\ln \left(hdom\right)-A+\left(k/t_1\right)\right].\left[\text{c}-\left(1/t_2\right)\right].\left[\text{c}-\left(1/t_1\right)\right]}$
Rk	Richards		$hdo{m}_2=A{\left\{1-{}^{t_2/t_1}\right\}}^{1/\left(1-m\right)}$
RA	$hdom=A{\left(1-e^{-kt}\right)}^{1/\left(1-m\right)}$	k	$hdo{m}_2=hdo{m}_1{\left(\frac{1-e^{-k{t}_2}}{1-e^{-k{t}_1}}\right)}^{1/\left(1-m\right)}$
Rm		A	$hdo{m}_2=A^{1-\left(\text{ln}\left(1-e^{-k{t}_2}\right)/\text{ln}\left(1-e^{-k{t}_1}\right)\right).hdo{m}_1{}^{\text{ln}(\left(1-e^{-k{t}_2}\right)/\text{ln}\left(1-e^{-k{t}_1}\right)}}$
MD	Hossfeld IV/ McDill-Amateis $hdom=A\frac{t^k}{A.c+t^k}$	m	$hdo{m}_2=\frac{A}{1-\left(1-\left(A/hdo{m}_1\right)\right).{\left(t_1/t_2\right)}^k}$

The Richards function has been widely utilized in forest growth modelling (PIENAAR; TURNBULL, 1973; AMARO et al., 1998; NUNES et al., 2011, RETSLAFF et al., 2015) as has the Lundqvist-Korf and McDill-Amateis functions (TOMÉ, 1988; CIESZEWSKI; BAILEY, 2000; FONTES et al., 2003, NUNES et al., 2011).

The algebraic difference approach (ADA) presents the structure of equations as base-age invariant, and the generalized algebraic difference approach (GADA) improves the traditional ADA, allowing more flexible dynamic equations which can be polymorphic and with multiple asymptotes (TAHAR et al., 2012). When only one parameter of the base model is related to a theoretical site quality measure, GADA is equivalent to ADA (NUNES et al., 2011).

The derivation of difference equation follows the methodology proposed by Amaro et al. (1998) for Richards and Lundqvist-Korf functions, consisting in equating subdefined ratios of base equations for the ADA method (BAILEY; CLUTTER, 1974).GADA (CIESZEWSKI; BAILEY, 2000) method was also used. The McDill-Amateis function (McDILL; AMATEIS, 1992) has the Hossfeld IV as integral form and has been used with success in several growth studies (WANG; PAYANDEH, 1996).

The parameter estimation was done using the nls function of the statistical environment R (R Core Team, 2015), performing a nonlinear regression analysis using the Gauss Newton algorithm. The model selection was based in qualitative and quantitative criteria (goodness-of-fit statistics). The root mean square error (RMSE) [1] in meters, the relative error in predictions (RE) [2] in percentage and the bias () [3] were calculated as follows, where: p is the number of parameters in the model, n is the number of observations used to fit the model, y_i is an observed dominant height, is the corresponding predicted dominant height and is the average of the observed dominant height.

$RMSE=\sqrt{\frac{{{\displaystyle \sum }}_{i=1}^n{\left(y_i-{\widehat{y}}_i\right)}^2}{n-p}}$ (1)

$RE\left(\%\right)=\frac{\sqrt{{{\displaystyle \sum }}_{i=1}^n{\left(y_i-{\widehat{y}}_i\right)}^2/n-p}}{\overline{y}}.100$ (2)

$\overline{e}=\frac{1}{n}{\displaystyle \sum }_{i=1}^n(y_i-{\widehat{y}}_i)$ (3)

For the qualitative criteria an evaluation of the different candidate models involved visual examination of the residuals against the estimated values and the fitted curves for different site index classes overlaid on the profile plots of the stands. Visual inspection is an essential point in selecting the most accurate model because curve profiles may differ considerably, even though goodness-of-fit statistics are similar (DIÉGUEZ-ARANDA et al., 2006). The models were chosen according to the goodness-of-fit, predictive ability, biological sense and compliance with the assumptions of nonlinear least square method (homoscedasticity, lack of autocorrelation and normality of residuals).

Due to the absence of continuous forest inventory data till the final rotation and the recent implementation of African mahogany culture in Brazil, an independent database for validation is not available, so the consistency of the models was analyzed by observing the stability of the site classification at the plot level, preferring classifications that were more stable, with fewer class changes for a single plot throughout the remeasurements. While this approach has been applied before (e.g. SCOLFORO, 1992; SELLE et al., 1994; MACHADO et al., 2011) it has been used to determine if the selected model produced an adequate classification using qualitative criteria. Here we propose a quantitative test to enable the comparison of the stability of the site classification generated by different models. This is done by quantifying the number of times each permanent plot changed site index over its life span. We then attributed weights for the number of changes with lower values for fewer changes (0 for no changes, 1 for one change and so on). The sum of the number of changes multiplied by this weight permitted the comparison of different classifications, identifying the model that presented the most stable classification as being the one with the lowest weighted sum.

RESULTS AND DISCUSSION

The model with the best performance using the guide curve method (integral form) was the Lundqvist-Korf. For the dynamic functions the models Lk, LG and MD had the best statistics (Table 3). The nonlinear regression analysis using the Richards dynamic function either did not converge (Rk) or didn’t have all the parameters significantly different from zero (RA), exception for model Rm with the m parameter as site specific.

Table 3
Goodness-of-fit statistics and parameter estimations of candidate functions, where ns indicates a non-significate parameter at 95% of probability.

Model	RMSE	RE	ē	Parameter	Estimate	SE	t-value
Guide curve method
Lundqvist-Korf	1.50	12.9	-0.008	A	57.9873	3.9614	14.64
				k	3.2358	0.0508	63.71
				d	0.5393	0.0266	20.27
Richards	1.58	13.6	-0.040	A	34.2267	0.9611	35.61
				k	0.1311	0.0080	16.32
				m	0.1276	0.0265	4.81
McDill-Amateis	1.54	13.3	-0.034	A	39.5919	1.4591	27.14
				k	1.2604	0.0332	37.97
				c	0.3195	0.0114	28.10
Difference equation method
Lk	0.98	7.71	0.020	A	35.6852	2.2150	16.11
				c	0.7301	0.0345	21.17
LA	1.13	8.89	0.057	k	2.8727	0.0756	38.02
				c	0.6776	0.0377	17.99
LG	0.94	7.36	0.014	A	3.3246	0.0166	200.66
				k	3.2222	0.0716	45.04
					1.0879	0.0134	81.36
Rk	Did not converge
RA	1.49	11.71	-0.203	k	0.0588	0.0694	0.847ns
				m	-0.0209	0.0330	-0.634ns
Rm	1.04	8.13	0.019	A	24.8564	0.7210	34.47
				k	0.24732	0.0149	16.53
MD	0.97	7.57	0.032	A	23.2717	0.5912	39.36
				k	1.7581	0.0398	44.16

The Lundqvist-Korf model had a rather high asymptote (A=58 m), more realistic values were found in equations Lk and Richards (A≈35 m), although Lemmens (2008) reported growth for K. ivorensis reaching 60 meters tall. Considering the goodness-of-fit statistics and the behavior of the asymptotic values of the model parameters we chose Richards, Lundqvist-Korf, McDill-Amateis for the guide curve method and Lk, LG and MD for the dynamic base-age invariant as the best candidate equations for site classification and discarded the others from further analysis.

Of the six selected models, the models fitted with the guide curve method presented a tendency towards heteroscedasticity (Figure 2). However, only the Lundqvist-Korf model did not present some tendency to underestimate values for smaller dominant height estimates (below 7 m) that could be observed for the other candidate models.

FIGURE 2
Standardized residuals plotted over the fi tted values for the best equations to estimate dominant height.

Before proceeding to the actual site classification, we plotted the site class limits of the best models to evaluate their behavior (Figure 3).

FIGURE 3
Comparison of height curves produced by the best models selected for site classifi cation of African mahogany plantations in Brazil

Although the values of the goodness of fit for functions Lk and MD were good, the curves showed unrealistic shapes. For both models even fixing the asymptote parameter still yielded curves with unrealistic classifications, but with better shape (Figure 4).

FIGURE 4
Height curves for equation Lk and MD with fixed asymptote parameter (A=45 and A=40, respectively).

Palahí et al. (2004) modelling site index for Pinus sylvestris in Spain also found that the McDill-Amateis model (MD) did not perform well resulting in poor estimates when extrapolated outside the age and site index range data, showing a similar curve shape as we found in our study to the highest site index before fixing the asymptote parameter.

For the final decision of the best model for site classification the classification stability (Table 5) was analyzed for the three more realistic models (Lundqvist-Korf and Richards using the guide curve method and LG using dynamic base-age invariant method).

Table 5
Classification stability of the best models selected for site classification.

	Stability of the site classification
Number of changes (weight attribute)	Lundqvist-Korf	Richards	LG
0 (x0)	50 (0)	42 (0)	43 (0)
1 (x1)	41 (41)	44 (44)	37 (37)
2 (x2)	40 (80)	33 (66)	34 (68)
3 (x3)	15 (45)	22 (66)	26 (78)
4 (x4)	4 (16)	9 (36)	9 (36)
5 (x5)	0 (0)	0 (0)	1 (5)
Total	150 (182)	150 (212)	150 (224)

The Lundqvist-Korf function was superior for the stability of the site classification. Model LG is a dynamic site index equation with 2 site-specific parameters (a representing the asymptote and b the inflexion point), providing two important and desirable attributes for site index classification, polymorphism and multiple asymptotes (CIESZEWSKI, 2002). In developing site index prediction systems, some desirable characteristics of site index models are polymorphism, multiple horizontal asymptote, one inflection point, and base-age invariance (ELFVING; KIVISTE, 1997; CIESZEWSKI; BAILEY, 2000).

Castro et al. (2015) evaluated the efficiency and reliability of two methods, guide curve and differences equations, concluding that the model developed by the difference equations method was more efficient to perform the prognosis of the studied Eucalyptus stand in Brazil. Arias-Rodil et al. (2015) comparing the age-dependent method with age-independent alternatives for estimating site index curves for Maritime Pine in Spain found that the age-independent equation performed best.

To check if our selected models presented autocorrelation, we used the plotted residuals versus lag residuals (Figure 5). The Lundqvist-Korf and Richards equations fitted with the guide curve method present a positive autocorrelation, confirmed by the significant values of the Durbin Watson test (d= 0.797 for Lundqvist-Korf and d= 0.815 for Richards). Equation LG did not show a strong temporal autocorrelation (d= 2.197) at P<0.05. The Q-Q plot of the standardized residuals for all equations showed that the residuals follow approximately the normal distribution. Due to the characteristics of the modelling data, the assumption of the independence of the error term is likely to be violated (PALAHÍ et al., 2004), however, hypothesis test may be quite unnecessary when the form of the model being fitted is limited, as the case of the traditional Lundqvist-Korf and Richards models.

FIGURE 5
Autocorrelation and Q-Q plot for the selected equation Lundqvist-Korf (a), Richards (b) and LG (c).

We are aware of the restrictions of our data, since few continuous forest inventories are available in Brazil and only three plots with measurements over the age of 10 compose the data base. Model selection was viewed as a compromise between statistical and biological considerations (reasonable values for prediction of dominant height for future ages as presented in Table 6), so the Lundqvist-Korf in its integral form is more indicated for site classification, but the polymorphic equation derived by the GADA approach from Lundqvist-Korf function (LG) resulted in a more adequate model for prediction of dominant height. Although advantages of dynamic site equations have been pointed out in several studies (AMARO et al., 1998; CALAMA et al., 2004; LOPEZ-SENESPLEDA et al., 2014) we decided to choose the more parsimonious model of Lundqvist-Korf using the guide curve method for site classification.

Table 6
Dominant height estimative limits in meters by age and site index classes for Khaya ivorensis A. Chev. in Brazil.

Age (years)	Site classes
	V (17)	IV (21)	III (25)	II (29)	I (33)
1	1.6 - 2.0	2.0 - 2.5	2.5 - 2.9	2.9 - 3.3	3.3 - 3.8
5	7.7 - 9.8	9.8 - 11.8	11.8 - 13.9	13.9 - 15.9	15.9 - 18.0
10	12.4 - 15.8	15.8 - 19.1	19.1 - 22.4	22.4 - 25.7	25.7 - 29.0
15	15.0 - 19.0	19.0 - 23.0	23.0 - 27.0	27.0 - 31.0	31.0 - 35.0
20	16.4 - 20.7	20.7 - 25.1	25.1 - 29.4	29.4 - 33.8	33.8 - 38.2
25	17.1 - 21.6	21.6 - 26.2	26.2 - 30.7	30.7 - 35.3	35.3 - 39.8
30	17.4 - 22.1	22.1 - 26.7	26.7 - 31.4	31.4 - 36.0	36.0 - 40.7
35	17.6 - 22.3	22.3 - 27.0	27.0 - 31.7	31.7 - 36.4	36.4 - 41.1

The site quality curves obtained in this study clearly show that Brazil is an excellent territory for African mahogany growth. The mean annual diameter increment for the different site quality classes (Figure 6) showed high values compared with other places that plant African mahogany, e.g. Heryati et al. (2011) comparing the growth of K. ivorensis A. Chev. in three different soil series in Malaysia reported values of mean annual increment for diameter from 2.32 to 2.88 cm year^-1. Erskine et al. (2005) conducted a trial with exotic species for five years in Australia, including K. nyasica, related a mean annual increment in diameter for Khaya of 1.58 cm.year^-1.

A more detailed study on quality of the proposed site classification should be carried out at the end of the first rotation of the plantations, with the aim of confirming whether the chosen model and reference age limits remain the most appropriate for the Khaya ivorensis A. Chev. plantations. Using the model selected on this study (Lundqvist-Korf integral form) for the five site indices, we calculated the limits for the dominant height by age and classes for the different site classes (Table 6).

FIGURE 6
Mean annual increment in diameter of plantation grown Khaya ivorensis A. Chev. in Brazil considering different site index classes obtained by the Lundqvist-Korf integral form model.

CONCLUSIONS

For site classification purposes, we found that the Lundqvist-Korf model, using the guide curve method, is the one with the most stable classification.

This model leads to the following equation for S prediction.

This was confirmed by the higher mean annual increment in diameter for the higher productivity classes.

ACKNOWLEDGMENTS

This work was conducted during a visiting scholar period at University of Lisbon, sponsored by the Capes Foundation within the Ministry of Education, Brazil (PDSE grant n. 99999.010833/2014-05).

REFERENCES

• AMARO, A.; REED, D.; TOMÉ, M.; THEMIDO, I. Modeling Dominant Height Growth: Eucalyptus Plantations in Portugal. Forest Science, v. 44, n. 1, p. 37-46, 1998.
• ARIAS-RODIL, M.; CRESCENTE CAMPOS, F.; BARRIO-ANTA, M.; DIÉGUEZ-ARANDA, U. Evaluation of age-independent methods of estimating site index and predicting height growth: a case study for maritime pine in Asturias (NW Spain). European Journal of Forest Research, v. 134, p. 223-233, 2015.
• BAILEY, R. L.; CLUTTER, J. L. Base-age invariant polymorphic site curves. Forest Science , v. 20, p. 155-159, 1974.
• BRAVO, F.; MONTERO, G. Site index estimation in Scots pine (Pinus sylvestris L.) stands in the high Ebro basin (northern Spain) using soil attributes. Forestry, v. 74, n. 4, p. 395-406, 2001.
• BURKHART, H. E.; TOMÉ, M. Modeling Forest Trees and Stands. Springer, 2012. 457p.
• CALAMA, R.; RÍO, M.; COQUILLAS, V.; CAÑELLAS, I.; MONTERO, G. Modelos de calidad de estación y de perfil de fuste para masas de Pinus uncinata Ram. en el Pirineo español. Invest Agrar: Sist Recur For, Fuera de serie, p.176-190, 2004.
• CASTRO, R. V. O.; CUNHA, A. B.; SILVA, L. V.; LEITE, H. G.; SILVA, A. A. L. Modelagem do crescimento e produção para um povoamento de Eucalyptus utilizando dois métodos para quantificação do índice de local. Scientia Forestalis, v. 43, n. 105, p. 83-90, 2015.
• CIESZEWSKI, C. J.; BAILEY, R. L. Generalized algebraic difference approach: theory based derivation of dynamic equations with polymorphism and variable asymptotes. Forest Science , v. 46, n. 1, p. 116-126, 2000.
• CIESZEWSKI, C. J. Comparing fixed and variable-base-age site equations having single versus multiple asymptotes. Forest Science , v. 48, p. 7-23, 2002.
• CLUTTER, J. L., FORTSON, J. C., PIENAAR, L. V., BRISTER, G. H., BAILEY, R. L. Timber management: a quantitative approach. John Wiley & Sons, 1983. 333p.
• DANQUAH, J. A. Restoration of degraded dry semideciduous forest ecosystems in Ghana: Effects of African mahogany species on soil chemistry, tree diversity and the application of leaf morphometrics for provisional seed zonation. 2012. 37f. Dissertation, University of Eastern Finland, Joensuu.
• DAVID, H. C.; PÉLLICO NETTO, S.; ARCE, J. E.; DALLA CORTE, A. P.; MARINHESKI FILHO, A.; BEHLING, A. Updating of Dominant Height Growth Modeling and site Index of Pinus taeda L. in southern Brazil. Australian Journal of Basic and Applied Sciences, v. 9, n. 2, p. 115-125, 2015.
• DIÉGUEZ-ARANDA, U.; GRANDAS-ARIAS, J. A.; ÁLVAREZ-GONZÁLEZ, J. G.; VON GADOW, K. Site Quality Curves for Birch Stands in North-Western Spain. Silva Fennica, v. 40, n. 4. p. 631-644, 2006.
• ELFVING, B.; KIVISTE, A. Construction of site index equations for Pinus sylvestris L. using permanent sample plots data in Sweden. Forest Ecology and Management, v. 98, p. 125-134, 1997.
• ERCANLI, I.; KAHRİMAN, A.; YAVUZ, H. Dynamic base-age invariant site index models based on generalized algebraic difference approach for mixed Scots pine (Pinus sylvestris L.) and Oriental beech (Fagus orientalis Lipsky) stands. Turkish Journal of Agriculture andForestry , v. 38, p. 134-147, 2014.
• ERSKINE, P.; LAMB, D.; BORSCHMANN, G. Growth performance and management of a mixed rainforest tree plantation. New Forest, v. 29, p. 117-134, 2005.
• FALESI, I. C.; BAENA, A. R. C. Mogno-africano em sistema silvopastoril com leguminosas e revestimento natural do solo. Embrapa Amazônia Oriental, Documentos 4, 1999. 52 p.
• FONTES, L.; TOMÉ, M.; COELHO, M. B.; WRIGHT, H.; SALES LUIS, J.; SAVILL, P. Modelling dominant height growth of Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) in Portugal. Forestry , v. 76, n. 5, p. 509-522, 2003.
• FRANÇA, T. S. F. A.; FRANÇA, F. J. N.; ARANGO, R. A.; WOODWARD, B. M.; ARANTES, M. D. C. Natural resistance of plantation grown African mahogany (Khaya ivorensisandKhaya senegalensis) from Brazil to wood-rot fungi and subterranean termites. International Biodeterioration & Biodegradation, v. 107, p. 88-91, 2016.
• HERYATI, Y.; BELAWAN, D.; ABDU, A.; MAHAT, M. N.; ABDUL-HAMID, H.; MAJID, N. M.; HASSAN, A.; HERIANSYAH, I. Growth Performance and Biomass Accumulation of a Khaya ivorensis Plantation in three Soil Series of Ultisols. American Journal of Agricultural and Biological Sciences, v.6, n. 1, p. 33-44, 2011.
• LEMMENS, R. H. M. J. 2008. Khaya ivorensis A. Chev. In: LOUPPE, D.; OTENG-AMOAKO, A. A.; BRINK, M. (Editors). PROTA (Plant Resources of Tropical Africa/Ressources végétales de l’Afrique tropicale), Wageningen, Netherlands. Available at: Available at: http://uses.plantnet-project.org/en/Khaya_ivorensis_(PROTA) Accessed in: 8 September 2016.
  » http://uses.plantnet-project.org/en/Khaya_ivorensis_(PROTA)
• LOPEZ-SENESPLEDA, E.; BRAVO-OVIEDO, A.; ALONSO, R.; MONTERO, G. Modeling dominant height growth including site attributes in the GADA approach for Quercus faginea Lam. in Spain. Forest Systems, v. 23, n. 3, p. 494-499, 2014.
• MACHADO, S. A.; SOUZA, R. F.; JASKIU, E.; CAVALEIRO, R. Construction of site curves for native Mimosa scabrella stands in the metropolitan region of Curitiba. Cerne, v. 17, n. 4, p. 489-497, 2011.
• McDILL, M. E.; AMATEIS, R. L. Measuring forest site quality using the parameters of a dimensionally compatible height growth function. Forest Science , n. 38, p. 409-429, 1992.
• NUNES, L.; PATROCÍNIO, M.; TOMÉ, J.; TOMÉ, M. Modeling dominant height growth of maritime pine in Portugal using GADA methodology with parameters depending on soil and climate variables. Annals ofForest Science , v. 68, p. 311-323, 2011.
• PALAHÍ, M.; TOMÉ, M.; PUKKALA, T.; TRASOBARES, A.; MONTERO, G. Site index model for Pinus sylvestris in north-east Spain. Forest Ecology and Management , v. 187, p. 35-47, 2004.
• PAULO, J.; TOMÉ, J.; TOMÉ, M. Nonlinear fixed and random generalized height-diameter models for Portuguese cork oak stands. Annals ofForest Science , v. 68, n.2, p. 295-309, 2011.
• PIENAAR, L. V.; TURNBULL, K. J. The Chapman-Richards generalization of Von Bertalanffy’s growth model for basal area growth and yield in even-aged stands. Forest Science , v. 19, n. 1, p. 2-22, 1973.
• PINHEIRO, A. L.; COUTO, L.; PINHEIRO, D. T.; BRUNETTA, J. M. F. C. Ecologia, silvicultura e tecnologia de utilizações dos mognos-africanos (Khaya ssp.). Viçosa: Sociedade Brasileira de Agrossilvicultura, 2011. 102p.
• R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2015.
• RENNOLLS, K. Forest height growth modelling. Forest Ecology and Management , v. 71, p. 217-225, 1995.
• RETSLAFF, F. A. S.; FIGUEIREDO FILHO, A.; DIAS, A. N.; BERNETT, L. G.; FIGURA, M. A. Curvas de sítio e relações hipsométricas para Eucalyptus grandis na região dos Campos Gerais, Paraná. Cerne , v. 21, n. 2, p. 219-225, 2015.
• SCOLFORO, J. R. S.; MACHADO, S. A. Curvas de índice de sítio para plantações de Pinus taeda nos estados do Paraná e Santa Catarina. Floresta, v. 8, n. 1-2, p. 159-173, 1988.
• SCOLFORO, J. R. S. Curvas de índice de sítio para Pinus caribaea var. hondurensis IPEF, n. 45, p. 40-47, 1992.
• SELLE, G. L.; SCHNEIDER, P. R.; FINGER, C. A. G. Classificação de sítios para Pinus taeda L., através da altura dominante, para a região de Cambará do Sul, RS, Brasil. Ciência Floresta l , v. 4, n. 1, p. 77-95, 1994.
• TAHAR, S.; PALAHI, M.; GARCHI, S.; BONET, J. A.; AMMARI, Y.; PIQUE, M. Modeling Dominant Height Growth in Planted Pinus pinea Stands in Northwest of Tunisia. International Journal of Forestry Research, v. 1, p. 1-13, 2012.
• TAMARIT-URIAS, J. C.; DE LOS SANTOS-POSADAS, H. M.; ALDRETE, A.; VALDEZ-LAZALDE, J. R.; RAMÍREZ-MALDONADO, H.; GUERRA-DE LA CRUZ, V. Dynamic site index equations for Tectona grandis at Campeche, Mexico. Agrociencia, v. 48, p. 225-238, 2014.
• TESCH, S. D. The evolution of forestry yield determination and site classification. Forest Ecology and Management , v. 3, p. 169-182, 1981.
• TEWARI, V. P.; JUAN G. ÁLVAREZ-GONZÁLEZ, J. G.; GADOW, K. Dynamic base-age invariant site index models for Tectona grandis in peninsular India. Southern Forests, v. 76, n. 1, p. 21-27, 2014.
• TOMÉ, M. Modelação do crescimento da árvore individual em povoamentos de Eucalyptus globulus Labill. (1^a rotação). Região Centro de Portugal. 1988. 277 p. PhD thesis. Technical University of Lisbon, Lisbon.
• VANCLAY, J. K. Assessing site productivity in tropical moist forests: a review. Forest Ecology and Management , v. 54, p. 257-287, 1992.
• WANG, Y; PAYANDEH, B. Evaluating a dimensionally compatible height growth modelling approach. The Forestry Chronicle, v. 72, n. 6, p. 647-650, 1996.

Publication Dates

• Publication in this collection
  Oct-Dec 2016

History

• Received
  29 July 2016
• Accepted
  13 Sept 2016