
 
ARTÍCULOS
Determination of stand structure in evenaged Oriental Beech Forests in Turkey Determinación de la estructura de rodales de Fagus orientalis coetáneos en Turquía
Aydin Kahriman ^{a*}, Lokman Altun ^{b}, Engin Güvendi ^{c} * Corresponding author: ^{a}Artvin Coruh University, Faculty of Forestry, Department of Forest Engineering, 08000, Artvin, Turkey, phone +904662151035, fax +904662151035, kaydin61@hotmail.com
SUMMARY This study aimed at developing a stand yield model for evenaged and pure Oriental beech (Fagus orientalis) stands in Sinop Province, located in the north of Turkey. The data used in this study were obtained from 238 sample plots. The sample plots are 24131 years of stand age range, 19.035.9 m of site index range, and 2.210.8 of stand density range (relative density index). With the generated equation systems, stand age, site index and density of the three main factors in Oriental beech stands, including the effects of individual, dual and triple factor interactions on the stands were numerically determined. As for evenaged Oriental beech stands, the main stand volume values increase depending on age, site index and density. The number of stems per hectare decreases as the stand grows older, and site index rises; the number, however, increases in parallel with stand density. The current annual volume increment reaches a peak at 3035 years at all site index groups and density levels, and begins to fall afterwards. Mean annual volume increment peaks up at 6575 years at all site index classes and density levels, before beginning to decline. The results were in arrangement with the acknowledged growth rules. Stand yield models were assessed as DensityVariable Yield Table used for estimation of the remaining, removal and total stands. Key words: stand model, oriental beech, densityvariable yield table, pure stand. RESUMEN Este estudio tuvo como objetivo desarrollar un modelo de rendimiento de rodales para bosques puros y uniforme de Fagus orientalis en la provincia de Sinop, norte de Turquía. Los datos se obtuvieron de 238 parcelas de muestreo de rodales con edad de 24131 años, índice de sitio de 19,035,9 m e índice de densidad relativo de 2,210,8. Con sistemas de ecuaciones generadas fueron determinados numéricamente la edad del rodal, el índice de sitio y la densidad de los rodales, incluyendo la determinación numérica de los efectos de las interacciones de los factores individuales, dobles y triples. En los rodales coetáneos, el aumento principal de los valores dependió de la edad, el índice de sitio y la densidad. La cantidad de tallos por hectárea disminuyó a medida que el rodal envejecía y el índice de sitio aumentaba; la cantidad, sin embargo, aumentó en paralelo con la densidad del rodal. El incremento de volumen anual actual alcanzó un pico a los 3035 años en todos los grupos de índice de sitio y de niveles de densidad, y comenzó a caer después. El incremento anual del volumen medio escala su pico a los 6575 años en todas las clases de índice de sitio y niveles de densidad, antes de comenzar a disminuir. Los resultados se encontraron de acuerdo con las normas reconocidas de crecimiento. El modelo de rendimiento de rodales se evaluó como tabla de rendimiento de variable de densidad utilizada para la estimación de los rodales restantes, eliminados y totales. Palabras clave: modelo de rodal, hayas orientales, tabla de rendimiento de variable de densidad, rodal puro.
INTRODUCTION Growth and yield models have contributed to recruitment, growth and mortality models. Growth models are divided into three groups according to the information available for predictions: Stand Model, SizeClass Model and IndividualTree Model (Vanclay 1994). Stand models provide forecast values of stand parameters such as number of trees per hectare, basal area, stand volume, quadratic mean diameter and mean height. Stand parameters determined by these models are average values; they do not provide predictions on the individual tree level (Pretzsch 2009). Size class models offer more detailed predictions than those offered by stand models about the structure of the stand and they use a class of trees as the basic unit for modelling (Vanclay 1994). The size class approach is a compromise between whole stand models and singletree models. Individual tree models use individual trees as the basic unit for modelling and they offer more detailed predictions than predictions offered by stand models and size class models (Vanclay 1994). The first yield studies regarding growth models were carried out 250 years ago, on evenaged pure stands. In the late 18^{th} and early 19^{th} century the earliest yield tables were established. They were a simple form of current examples. Pretzsch (2009) classifies development of growth models in four periods, namely, those in the form of a simple table with limited amount of data (from the last quarter of the 18^{th} century to the last quarter of the 19^{th} century), normal yield tables organized depending on the age and site index variables (from the last quarter of 19^{th} to the last quarter of 20^{th} century), yield models where growth relations are calculated via computeraided mathematical equations (from the first quarter to the last quarter of the 20^{th} century) and much more detailed stand simulation models that can be hardly operated even in computers with highspeed memory (from the last quarter of the 20^{th} century up to the present). Yield tables are divided into two groups: normal and variable density tables (Vanclay 1994). Normal yield tables can be established when the data are arranged according to age and site index; on the other hand, densityvariable yield tables can be formed when the data are arranged depending on stand density, in addition to age and site index. Normal yield tables provide estimates of expected yield tabulated by stand age and site index for ideal, fully stocked or "normal" forest stands (Weiskittel et al. 2011). Growth significantly varies as a function of stand density, as well as depending on stand age and site index. DensityVariable Yield Tables provide mean values regarding age, site index and density of stand elements; thus, they enable the determination of change in growth of stands with varying density. Fagus spp. is the most abundant and economically important hardwood genera in the temperate forests of the northern hemisphere (SaghebTalebi et al. 2011). There are 12 species under the beech genus; in Turkey, two of these grow naturally, namely, Oriental beech (Fagus orientalis Lipsky) and European Beech (Fagus sylvatica L.). Oriental beech is one of the principal forest tree species in economic terms in Turkey, covering a vast area. Oriental beech is naturally distributed in the Black Sea region, Aegean region, East Mediterranean and Marmara regions of Turkey, as well as in many other regions of the world (Davis 1982). According to the forest inventory conducted by the General Directorate of Forestry in 2012, the total area of forestlands in Turkey is equal to 22,342,935 ha; of which, 7,346,831 ha (33 %) consist of forests of broadleaved trees, and beech forests constitute 1,961,660 ha (9 %) (General Directorate of Forestry 2015). The yield tables developed in 1795 possibly represent the first serious approach to beech growth modeling (ÁlvarezGonzález et al. 2010). For the last 50 years, many yield tables and stand growth models have been developed for pure and natural beech stands in Europe (Dmitriev 1967, Carbonnier 1971, Marschall 1975, Mendlik 1983, Dittmar et al. 1986, Federal Agency of Forest Management of Russia 2008). The first study about growth relationships of Oriental beech in Turkey was presented by Kalıpsız (1962). In his study, he examined growth using data from singletree and stand level. Alemdağ (1963) studied certain growth relations regarding Oriental beech in the Tokat basin, and recommended a rotation age of 120 years pursuant to the technical maturation period. Asan (1987) carried out site index studies on managed, evenaged, one stored, normal and mediumstocked^{1} stands that constitute more than 50 % of Oriental beech in the Western and Central Black Sea regions. Carus (1998) made use of the site index table, as well as yield tables for normal and optimum formations, so as to put forth the development of Oriental beech in the Western and Central Black Sea regions, where these native species are the most widespread. In the study, normal yield tables revealed growth relations of normal density, evenaged and pure Oriental beech stands that are formed through natural ways without any intervention. Atıcı (1998), on the other hand, studied the growth of unevenaged Oriental beech forests. He established site index, pressure, bark coefficient, periodical diameter increment and fine doubleentry tree volume equation tables to be applied on beech forests of varying ages. Güvendi (2013) tried to determine absolute rotation ages on an ecological basis with regard to pure Oriental beech stands in the Sinop region. In the wake of his study, the absolute mean rotation age of pure Oriental beech forests in the lower zone was found to be 7298 years, while the same was determined to be 90102 years for those in the higher zone. The objectives of this study are (i) to establish a DensityVariable Yield Table by means of stand models applicable for pure and evenaged Oriental beech stands within the Sinop region, (ii) to compare these yield table parameters with previous ones on Oriental beech in Turkey and on European beech and (iii) to make a general evaluation of beech growth in Sinop. METHODS This study employs, as research material, data from 238 sample plots of evenaged and pure Oriental beech stands in Sinop (Sinop, Ayancık and Türkeli Forest Enterprises), the northernmost city in Turkey (Altun et al. 2012) (figure 1). All of the 238 sample plots are in natural grown stands. This study area ranges in latitude (North) from 41°27'00'' to 42°01'33'' and in longitude (East) from 34°16'25'' to 35°24'34''. During selection of the mentioned sampling area, we tried to include as many representative areas as possible, where stands of Oriental beeches are fairly balanced in terms of different ages, site index and density levels.
Figure 1. Map of Oriental beech distribution in Turkey (General Directorate of Forestry 2015)
The size of sampling areas in the hereby study varies between 400 m^{2} and 800 m^{2} depending on stand structure. In each sample area, diameter at breast height of each tree (d_{1.30}) was measured, as well as height and age of a certain number of trees that are equally distributed diameter classes. In addition to that, the tallest trees by the 100 tallest trees per hectare in each sample area were measured. The hereby study calculates mean diameter, mean height, top height and tree number, basal area and volume values, as well as stand age, site index and density level by using data from sampling areas, in order to establish stand growth models in form of densityvariable yield table. For sample areas, mean diameter weighted by basal area of stands is used as a basis for calculation of mean diameter. Stand mean height was calculated via mean height weighted by basal area. The number of trees per hectare was calculated by multiplying the number of trees thicker than 8 cm in sampling areas with hectare conversion coefficient. Stand basal area was found through addition of basal areas of single trees at the site and their multiplication for conversion to hectare. Stand top height was calculated by means of arithmetic mean of highest trees in the sampling area pursuant to 100tree per hectare calculation. The volume estimates in the sampling area made use of a volume equation developed by (Carus 1998) on Oriental Beech, given in equation number (1) below. Upon conversion of total single tree volumes in hectare via doubleentry tree volume equations, we attained stand volumes. In the model, "d" signifies diameter at breast height while "h" shows tree height. Site index of stands was calculated by site index models obtained through Generalized Algebraic Difference Approach (GADA). Hossfeld growth function was used to model the changes in the height values of trees at different ages according to age and other variables (base age and site index). These functions were rearranged by Cieszewski (2001) using GADA and were transformed into a model structure which can be used to make dominant height estimations for site and base ages. The site index model explains 97.6 % of the variance in height (R^{2}_{adj}). In the equation, the root mean square error (RMSE) = 1.07 m, bias was 0.49 m, while mean absolute error (MAE) is 0.78 m. The functional structure and parameter estimate value of the established site index model were given in equation number (2). In this model, h signifies top height estimated for a certain t age, t_{0} was standard age, whereas h_{0} shows top height at standard age of t_{0}, which in other words, is the site index. Stand site index can be calculated directly using a base age (t_{0}), stand age (t) and dominant height (h) determined via equality numbered 2 (Ercanlı et al. 2014). Relative Density Index (RD), developed by (Curtis 1982), is used for calculating stand density (Vanclay 1994). In the equation (equation 3), BA (m^{2} ha^{1}) shows total calculated basal area, and (cm) is the quadratic mean diameter, weighted by the basal area. Relative Density Index (RD): Descriptive statistics including mean, minimum, maximum and standard deviation of the plot characteristics, such as stand age, site index, stand density, quadratic mean diameter, mean height weighted by basal area, stand basal area, stand volume and number of trees per hectare were listed below (table 1).
Table 1. The minimum, maximum, mean, and standard deviation of main characteristics
Volume yield is calculated by adding the current standing volume (remaining stand volume) and the sum of the volume removed (removed stand volume) in all previous periods (Pretzcsh 2009). Stand volume and yield elements such as number of trees per hectare, basal area, volume, quadratic mean diameter and mean height of main stand are calculated as a function of stand age, site index and stand density. Regression equations using the Forward Selection Method, Backward Selection Method and Stepwise Selection Method techniques implemented using SPSS version 19.0 (SPSS Institute Inc. 2010). Removed stand volume and volume elements are necessary in order to determine total stand volume. The total stand volume is equal to the sum of the main stand volume at that age and the stand volume removed until the mentioned age (Kalıpsız 1962). Volume and volume elements of the removed stand can be directly fixed through continuous sampling areas, nevertheless, in the absence of continuous sampling areas, they can be estimated via temporary sampling areas (Vanclay 1994). In this study, removed stand elements are calculated using dying and dead trees data in the sampling area. This removed volume (the sum of the thinned volume) may amount to 2060 % of the total volume production depending on the selfthinning and the intensity of thinning operations (Vanclay 1994, Pretzcsh 2009). In addition to volume and volume elements of the main and removed stands on yield tables, in this study, we calculated other elements on yield table, including the current annual volume increment, the mean annual volume increment (periodic annual volume increment = PAI) and total stand volume. Mean annual increment (MAI) and current annual increment (CAI) are important for estimating the stand growth phase in the period among inventories, planning future stand treatment, analyzing thinning responses and calculating the quantity of removals (Pretzcsh 2009). Two statistical criteria obtained from the residuals were examined: the adjusted coefficient of determination (R^{2} adj) and root mean square error (RMSE). The model evaluation criteria were as follows:
Where: n is the number of observations; p is the number of parameters; y_{i}, ŷ_{i}, and ӯ_{i} are observed, predicted and mean values of the dependent variable, respectively. RESULTS The densityvariable yield tables were created for 20130 years, IIIIIIIVV site index, and 310 stand density. Not all table values were provided. In this study stand, age ranges from 24 to 131 years. Model results are valid for up to 130 years old. The other four stand variables in addition to the basal area value of old growth stands of 130yearold may not be accurately obtained with this model. The data evaluated in this study were obtained from the study conducted by Altun et al. (2013). The quadratic mean diameter (͞d_{q}), mean height (͞h_{q}), tree density (N), basal area (BA) and stand volume (V) of the main stand were calculated in densityvariable yield tables as functions of stand age (A), site index (SI) and stand density (RD) via regression equations (equations 610). Statistical information about the regression equations is provided in table 2. Subsequently, the change in the main and removed stand elements of the densityvariable yield tables as well as the other elements of the yield table were calculated on the basis of stands.
Table 2. Certain statistical information about growth relations of main stands.
*f (correction factors of logarithmic functions)
All coefficients in regression equations are found significant at significance level of P < 0.001. Among the developed regression models, stand basal area (R^{2}= 0.97) provided the highest relation, followed by models organized for mean height (R^{2}= 0.96), stand volume (R^{2}= 0.94), number of stand trees (R^{2}= 0.78) and mean diameter (R^{2}= 0.73), respectively. As the logarithmic values of dependent variables in equations 7, 8, 9 and 10 are estimated; it is necessary to extract the antilogarithm of these values after correction via f factor calculated under equation 11 (Baskerville 1972). The variance of mean diameter of main stands which are found using the aforesaid regression models, as well as of the number of trees per hectare and stand volume, depending on stand age, site index and stand density are given in figures 26. The comparison of yield tables were carried out using the volume and volume increment values at the stand density of 9 (9 RD) which is considered as normal in this study. Also, in this study, rotation age is taken as 100 years. The rotation age (100 years) is commonly used in these stands which are managed for production purposes. The change of some stand elements according to stand ages, site index and stand density are given in table 3.
Figure 2. Change in mean diameter according to (A) stand age and stand density for
Figure 3. Change of main stand volume according to (A) stand age and stand density for
Figure 4. Change of number of stems according to (A) stand age and stand density for
Figure 5. Change of current annual increment according to (A) stand age and stand density
Figure 6. Change of mean annual increment according to (A) stand age and stand density
Table 3. Change of some stand elements according to stand ages, site index and stand
SI site index, RD relative density, A age, quadratic median diameter, mean height weighted by
Mean diameter of pure Oriental beech stands increases as site quality is increased and stand grows older. The mean diameter falls as stand density rises (figure 2A and figure 2B). A hundred yearold stands of normal density (9 RD) have a mean diameter of 35.3, 32.6, 29.8, 27.0 and 24.2 cm, respectively, from high to low site index (IIIIIIIVV site classes). The mean diameters of 100yearold stands of the middle site index group (III site class) are 31.4, 30.7, 30.2 and 29.8 cm from low to high density (3579 RD) levels, respectively (table 3). As for evenaged Oriental beech stands, the main stand volume values increase depending on age, site index and density (figures 3A and 3B). Main stand volumes of 70 yearold stands at 9 stand density and site classes I, III and V are 502.4, 357.8 and 207.9 m^{3} respectively; for stands of 100 years volumes are 636.8, 453.6 and 263.6 m^{3}, and for 130 yearolds volumes are 723.6, 515.4 and 299.5 m^{3}, respectively. On the other hand, the main stand volumes of 70 yearolds at III, SI and density of 3579 RD are 119.2, 198.7, 278.3 and 357.8 m^{3}, while the volume for 100 yearolds are 151.1, 251.9, 352.7 and 453.6 m^{3}, and the volumes of 130 yearold stands are 171.1, 286.2, 400.8 and 515.4 m^{3}, respectively (table 3). According to the data model results main stand volume of 130 yearold at 9 stand density and 32 m site index is 723.6 m^{3} ha^{1} (table 3). This value is smaller than actual stand volumes of only 7 sample areas. Also, the maximum actual stand volume value is 994.6 m^{3} ha^{1} (table 1). Site index and stand density values of this sample area are 33.5 m and 10.26, respectively. It has a high volume value compared with the other stand volumes of sample areas. On the other hand, main stand volume of 130 yearold stands at 10.5 stand density and 32.0 m site index is 844 m^{3} ha^{1} for equation 10. As a result, we can say that the volume model is an accurate estimate for Oriental beech stands. For Oriental beech stands, the number of stems per hectare decreases as the stand grows older, and site index rises; the number, however, increases in parallel with stand density (figures 4A and 4B). Normal stand density and number of stems per hectare at site index groups of I, III and V are 3,116, 3,696 and 5,164 at 20 years; 734, 803 and 920 at 70 years; and 305, 369 and 456 at 130 years. Stems per hectare in middle site index class and at density levels of 3, 5, 7 and 9 are 1,141, 1,876, 2,694 and 3,696 at 20 years; 248, 407, 585 and 803 at 70 years; and 114, 187, 269 and 369, at 130 years (table 3). The current annual volume increment of Oriental beech stands reaches a peak at 3035 years at all site index groups and density levels, and begins to fall afterwards (figures 5A and 5B). As the site index gets better and stand density gets higher, current annual volume increment increases as well. For 35 yearold stands of normal density (9 RD), current annual increment changes as 11.72, 9.36, 7.87, 6.05 and 4.34 m^{3} ha^{1} from higher to lower site index (IIIIIIIVV site classes); while 35 yearold stands in middle site class (III site class) yielded values of 2.66, 4.48, 6.84 and 7.87 m^{3} ha^{1} from low to higher density (3579 RD), respectively (table 3). Mean annual volume increment peaks up at 6575 years at all site classes and density levels, before beginning to decline (figures 6A and 6B). Mean annual volume increments for site classes IIIIII peak up at 65, 70 and 75 years, respectively. Mean annual volume increment increases as site index gets better and stand density gets higher. For 70 yearold stands of normal density (9 RD), the mean annual increment is 8.42, 7.03, 5.66, 4.35 and 3.12 m^{3} ha^{1} from higher to lower site index (IIIIIIIVV site classes), while for 70 yearold stands in middle site index class (III site class), these figures change as 1.90, 3.16, 4.41 and 5.66 m^{3} ha^{1} from low to high density (3579 RD) (table 3). DISCUSSION The acquired results will enable the accurate determination of volume and volume elements of the oriental beech stands in addition to enabling the management of the forests at full capacity. Thus, maximum benefit will have been attained from forest soil. In addition, the optimal and actual tree stock that stands should have will be compared. Densityvariable yield tables, which are established under this study (SinopTurkey), were compared with normal yield tables for pure Oriental beech established by (Carus 1998) for Turkey, and with yield tables for European Beech in other countries (Ukraine, Dmitriev 1967; Sweden, Carbonnier 1971; Austria, Marschall 1975; Hungary, Mendlik 1983; Germany, Dittmar et al. 1986; Russia, Federal Agency of Forest Management of Russia 2008) (figures 79). For a more accurate result, the values in yield tables of stands of normal density in first site index category are compared. The normal density value of pure Oriental beech is taken as 9 (9 RD), taking into account the density in the studied sampling areas.
Figure 7. Comparison of Oriental beech stands for mean diameter (A) and basal area (B)
Figure 8. Comparison of Oriental beech stands for main stand volume (A) and number of
Figure 9. Comparison of Oriental beech stands for current (A) and mean annual increment (B)
At the end of our calculations, the quadratic mean diameters of Oriental beech stands of site index class I and normal density (9 stand density) found in our results are higher than those found in previous studies in stands aged between 35 and 55 (Dmitriev 1967, Carbonnier 1971, Marschall 1975, Mendlik 1983, Carus 1998), aged between 5565 (Dmitriev 1967, Carbonnier 1971, Marschall 1975, Carus 1998) and in stands aged between 65 and 120 (Marschall 1975, Carus 1998) (figure 7A). The basal area values in this study are lower than those from (Carus 1998, Federal Agency of Forest Management of Russia 2008) at every age, but higher than results obtained from all other studies (figure 7B). Although basal area values increased with increasing age, the periodic rate of basal area increment slows with older stand age. However, basal area growth continues. Main stand volumes of Oriental beech in site index class I and of normal density (9 stand density index) found in this study are lower than those found in previous studies for stands aged between 30 and 70 (Carus 1998, Federal Agency of Forest Management of Russia 2008) and for stands aged between 70 and 120 (Mendlik 1983, Carus 1998, Federal Agency of Forest Management of Russia 2008) (figure 8A). The stand volume values at which the mean annual volume increment values have reached maximum were determined in ascending order as Carbonnier (420.0 m^{3} at 110 years), Carus (428.0 m^{3} at 40 years), Dmitriev (429.0 m^{3} at 70 years), SinopTurkey (502.4 m^{3} at 70 years), Marschall (552.0 m^{3} at 90 years), Russia (552 m^{3} at 60 years), Dittmar (631.0 m^{3} at 110 years) and Mendlik (721.0 m^{3} at 100 years) (Figure 8A). The number of trees, calculated using stand models in this study, is higher than that found in previous studies done on stands aged between 3560 (Federal Agency of Forest Management of Russia 2008), aged between 60 and 75 (Dmitriev 1967, Carbonnier 1971, Mendlik 1983, Federal Agency of Forest Management of Russia 2008) and aged between 75 and 120 (Dmitriev 1967, Carbonnier 1971, Marschall 1975, Mendlik 1983, Federal Agency of Forest Management of Russia 2008) (figure 8B). The maximum point on the curve of CAI (the inflection point) is the point corresponding to the fastest change in yield. The peak ages of current annual volume increment values from low to high are as follows: Carus (before 20 years), SinopTurkey (3011.7 m^{3}), Russia (3017.6 m^{3}), Dmitriev (4012.3 m^{3}), Carbonnier (5011.6 m^{3}), Marschall (6014.0 m^{3}), Dittmar (6016.0 m^{3}) and Mendlik (7016.8 m^{3}), respectively (figure 9A). Increasing yield of stand is an important topic in forestry studies. For this, firstly, it is important to know where the intersection point for current and mean annual volume increment is situated on the curves. This intersection point corresponds to the age when the mean annual volume increment is at a maximum. This point of intersection denotes the age of most efficient volume production over time and also, the point in which they intersect is called the biological rotation age. Thus, it is necessary to know the ages at which the mean annual increments of stands are at a maximum. As it can be seen in Figure 9B, the mean annual volume increment values of the beech stands in Turkey reach the maximum value after the beech stands in Russia and before the beech stands in Austria, Hungary, Germany and Sweden. The peak age and values of general mean volume increment are, from low to high: Carus (4011.6 m^{3}), Russia (6015.6 m^{3}), Dmitriev (708.8 m^{3}), SinopTurkey (708.4 m^{3}), Marschall (909.1 m^{3}), Mendlik (10013.4 m^{3}) Carbonnier (1107.8 m^{3}) and Dittmar (11010.6 m^{3}) (figure 9B). If our objective is to attain maximum wood yield from the forest, the values at ages when the mean annual volume increment of stands is at a maximum should be considered as highly important. The average diameter values at ages when the mean annual volume increment values have reached their peak values are determined as Carus (13 cm at 40 years), SinopTurkey (26.0 cm at 70 years), Dmitriev (26.9 cm at 70 years), Russia (29.5 cm at 60 years), Marschall (32.5 cm at 90 years), Mendlik (42.5 cm at 100 years) and Carbonnier (47.6 cm at 110 years) with diameters ranging from thin to thick respectively (figure 7B). On the other hand, when 100 years is taken as basis, it is observed that the mean diameter values at this age vary from thin to thick diameter as 25.0 cm (Carus), 35.0 cm (Marschall), 35.3 cm (SinopTurkey), 38.0 cm (Dmitriev), 42.0 cm (Carbonnier), 42.5 cm (Mendlik) and 43.6 cm (Russia) (figure 7B). If our objective is to acquire quality logs with a certain diameter, the year when the trees in the stands reach these diameter values has a critical importance. For example, if we take 36 cm diameter value as the basis, the age when the trees reach this value can be listed in ascending order as 75 years (Russia), 86 years (Mendlik), 88 years (Carbonnier), 93 years (Dmitriev), 102 years (SinopTurkey) and 103 years (Marschall) (figure 7B). CONCLUSIONS Based on densityvariable yield tables derived in this study, it is now possible to make estimates according to a given age, site index class and stand density. The basal area, volume, mean diameter, and mean height of a main stand increases at the same age and density in parallel with site index; while number of trees fall under the same conditions. Number of trees, basal area and volume of main stands rise depending on stand density at the same age and site index class; whereas mean diameter and mean height decline. Current annual increment peaks at a certain age in every site index class and density levels before declining; besides, the increment rises as the site index grows better and stand density goes up. Total stand volume also increases with a better site index and higher stand density. The mean increment of general and main stand reaches the maximum level at the ages of 65 to 75 for every site index class and density, before falling. As site index prospers and stand density rises, the mean increment of the general and main stand goes up as well. In this study, the growth models of pure and evenaged Oriental beech stands are established based on data from sampling areas. This research enables estimation of stand growth elements for various ages, site indexes and density levels of Oriental beech stands in Sinop Province. In particular, the increment and growth amount of thinned interfered Oriental beech stands can be calculated based on the equation systems established in this study. Sinop Province, the northernmost point in Turkey, hosts the most important forests in the Black Sea region. Oriental beech stands in this province are of huge economic and biological importance. It is possible to protect and sustain the existence of Oriental beech stands only through forest management plans taking into account growth relations and stand structures in these forests. To that end, the recognition of growth relations regarding pure Oriental beech stands is crucial for preparation of forest management plans. The present study has determined the production potentials of pure and evenaged Oriental beech stands, and as such, has developed the required growth models for both management plans and silvicultural practices. The growth legality will bring about optimum utilization of ecological, economic and social functions of Oriental beech stands. NOTE ^{1} Medium stocked: crown closure or canopy closure = 4070 %; normal crown closure > 70 %; especially close = 100 %. ACKNOWLEDGEMENT The data for this study were obtained from Altun et al. (2012) project. The authors thank all who provided data. This work was supported by the TÜBÍTAK  TOVAG (The Scientific and Technological Research Council of Turkey  Research Committee of Agriculture, Forestry and Veterinary), Project Number: 107O752.
REFERENCES Alemdağ Ş. 1963. Some growth and yield in Oriental Beech (Fagus orientalis Lipsky.) Forests in Tokat and rotation age to be applied to the forest. Ankara, Turkey. Forestry Research Institute of Ankara. Technical Bulletins Number 12. Altun L, N Ay, B Serdar, N Köse. 2012. Determine the ecological based rotation length, anatomical, dendroclimatical, physical and mechanical characteristics of pure Oriental Beech (Fagus Orientalis Lipsky.) stands in Sinop Regional Directorate of Forestry. Trabzon, Turkey. TübitakTovag Project. Project Number: 107O752. ÁlvarezGonzález JG, A Zingg, KV Gadow. 2010. Estimating growth in beech forests: a study based on long term experiments in Switzerland. Annals Forest Science 67(3): 307307. Asan Ü. 1987. Site quality research of Oriental Beech (Fagus Orientalis Lipsky.) stands in West and Central Black Sea Region. Istanbul University Journal of Forestry Faculty A37(1): 106130. Atici E. 1998. Increment and growth in unevenaged beech (Fagus orientalis Lipsly.) forests. PhD Thesis. Graduate School of Natural and Applied Sciences. Istanbul University, Istanbul, Turkey. Baskerville G. 1972. Use of logarithmic regression in the estimation of plant biomass. Canadian Journal of Forest Research (2): 4953. Carbonnier C. 1971. Yield of beech in southern Sweden. Stockholm, Sweden. Royal College of Forestry. Studia Forestalia Number 91. 68 p. Carus S. 1998. Increment and growth in even aged beech (Fagus orientalis Lipsky.) forests. PhD Thesis. Graduate School of Natural and Applied Sciences. Istanbul, Turkey. Istanbul University. 360 p. Cieszewski CJ. 2001. Three methods of deriving advanced dynamic site equations demonstrated on inland DouglasFir Site Curves. Canadian Journal of Forest Research 31: 165173. Curtis RO. 1982. A simple index of stand density for Douglasfir. Forest Science 28(1): 9294. Davis PH. 1982. Flora of Turkey and the East Aegean Islands. Edinburgh, UK. Edinburgh University Press. 590 p. Dittmar O, E Knapp, G Lembcke. 1986. DDRBuchenertragstafel 1983. IFEBerichte aus Forschung und Entwicklung 4: 159. Dmitriev IP. 1967. Yield advancing of beech stands of western USSR  Moscow. Forest industry. p. 203205. (in Russian). Ercanlı İ, A Kahriman, H Yavuz. 2014. Dynamic baseage invariant site index models based on Generalized Algebraic Difference Approach for mixed Scotch Pine (Pinus sylvestris L.) and Oriental Beech (Fagus orientalis Lipsky) stands. Turkish Journal of Agriculture and Forestry 38 (1):134147. Federal Agency of Forest Management of Russia. 2008. Tables and models of growth and productivity of forests of major forest forming species of Northern Eurasia. Second edition, supplemented. Moscow, Russia. Federal Agency of Forest Management of Russia. 883 p. General Directorate of Forestry. 2015. Forets of Turkey. Ankara, Turkey. Turkey General Directorate of Forest Publications. 219 p. Güvendi E. 2013. Determination of ecology based rotation length of naturally distributed pure Oriental Beech (Fagus Orientalis Lipsky.) stands (Example of SinopAyancıkTürkeli). PhD Thesis. Trabzon, Turkey. Graduate School of Natural and Applied Sciences. Karadeniz Technical University. Kalıpsız A. 1962. Growth and yield of Oriental Beech (Fagus Orientalis Lipsky.) stands. Publications of General Directorate of Forestry, Number 339. Istanbul, Turkey. Marschall J. 1975. New yield tables for Austria. Allgemeine Forstzeitung 87(6):195200. Mendlik G. 1983. Beech yield table. ErdeszetiKutatasok 75:189198. Pretzcsh H. 2009. Forest dynamics, growth and yield. Berlin, Germany. Springer International. 664 p. SaghebTalebi K, Z Mirkazemi, R Akhavan, A Karimidoost, M Maghsoudloo, D Moghadasi. 2011. Some structural characteristics in the far east border of the distribution range of Oriental Beech (Fagus orientalis Lipsky) stands. In The 9th IUFRO International Beech Symposium. 1217 September, Göttingen, Germany. SPSS Institute Inc. 2010. IBM SPSS Statistics 19 Core System User's Guide. Chicago, UAS. SPSS. 426 p. Vanclay JK. 1994. Modeling forest growth and yield. Applications to mixed tropical forests. Wallingford, UK. CABI Publishing. 312 p. Weiskittel AR, DV Hann, JA Kershaw, JK Vanclay. 2011. Forest growth and yield modeling. Chicester, UK. Wiley  Blackweel. 430 p.
Recibido: 21.04.16

© 2017 •
Facultad de Ciencias Forestales, Universidad Austral de Chile
Teléfono/Fax: 56 63 221998 • Casilla 567 • Campus Isla Teja S/N • Valdivia • Chile
Email:
revistabosque@uach.cl