The most widely accepted method of evaluating site productivity is site index. In spite of some important restrictions it is still a useful concept in both forest research and management. One of the most important challenges when using site index is an age trend manifested by a negative correlation between site index and stand age. Age trend may result from the inappropriateness of site index models. In this paper we develop a new approach for assessing age bias in site index models. Field data collected from 311 sample plots established in Norway spruce stands in the Polish region of the Carpathians formed the basis of this study. In the proposed approach the appropriateness of site index models is assessed by analyzing the existence of age trends in residuals of geocentric site index prediction models. Using the developed approach we demonstrated that when significant correlations exist between residuals of site index prediction models and stand age, it likely indicates the existence of an age trend and thus the inappropriateness of site index model. To remedy this situation, we demonstrated that the observed age trend can be quantified and utilized in new, non-biased, site index models.
Site ProductivityAge TrendClimate ChangeHeight GrowthSite Index ModelForest InventoryIntroduction
Site productivity is a quantitative estimate of the potential of a site to produce plant biomass (Skovsgaard & Vanclay 2008). Forest site productivity remains a fundamental variable in forestry (Bontemps & Bouriaud 2013). Reliable estimates of potential site productivity are essential for sustainable forest management, as it is often a key criteria when considering site- and species-specific decisions concerning species composition, silvicultural treatments, determining allowable cut, rotation period and forecast of timber yield (Splechtna 2001, Pretzsch et al. 2008). To be most useful for modeling and prediction, a measure of site quality must be quantitative, that is expressed by a number (Burkhart & Tomé 2012).
The most commonly used and widely accepted method of evaluating site productivity is site index (Skovsgaard & Vanclay 2008, Sharma et al. 2011). In spite of the fact that in many cases site index is not always sufficient to characterize site productivity (Bontemps & Bouriaud 2013), it is still one of the most widely used indicators of wood production and is useful in both forest research and management (Skovsgaard & Vanclay 2013). Therefore, the construction of site index models remains a fundamental task for site productivity differentiation (Skovsgaard & Vanclay 2008, Skovsgaard & Vanclay 2013, Sharma et al. 2011).
One of the most important challenges in the site index concept is an age trend manifested by a negative correlation between site index and stand age, which has been frequently reported (Elfving & Tegnhammar 1996, Socha 2008, Albert & Schmidt 2010, Nothdurft et al. 2012, Sharma et al. 2012, Yue et al. 2014). In most studies problems associated with age trend in site index resulted from changes in site conditions caused by non-static factors: nitrogen deposition, increasing CO_{2} concentration and climate change (Elfving & Tegnhammar 1996, Hasenauer et al. 1999, Schadauer 1999, Fries et al. 2000, Solberg et al. 2009, Nothdurft et al. 2012, Sharma et al. 2012, Bontemps & Bouriaud 2013, Yue et al. 2014). Highly productive forests are usually managed under short rotation, and consequently older stands may be over-represented on infertile sites, since they require a longer time to reach a certain goal of wood production which also can cause an age trend in site index estimates (Tegnhammar 1992, Nothdurft et al. 2012, Yue et al. 2014).
In a study conducted in Sweden, Johansson (1995) found significant variation in site index curves of forest stands due to mineral soil type and geographical location. Effect of genetic variability on the height development of trees and resulting age trends were also reported by Adams et al. (2006). Genetic variability has been identified as the factor modifying the trajectory of site index curves by Buford & Burkhart (1987), who found that the site index curves developed for different provenances differed in their parameters describing the asymptote for height growth. Climate, soil and genotype-by-environmental interactions can favor the adaptation of species to local ecological conditions, causing variable height patterns for the same species among ecoregions (Monserud & Rehfeldt 1990, Alvarez-González et al. 2005). Significant differences in the pattern of tree heights can also occur for the same species, even within a given site index class (Monserud 1984). García (2010) reported that three sources of height growth development variability can be identified: (a) between sites, (b) within sites, and (c) observation error. The differences between site index model trajectories with age may also result both from the method and data used in deriving site equations (Cieszewski 2001, Cieszewski & Strub 2007).
In order to avoid age bias in site index prediction, Sharma et al. (2012) included age as one of the predictor variables in site index prediction models for Norway spruce in Norway and obtained an increase in the proportion of explained variance. Including age as one of the variables reduces bias of site index prediction model; however, in such cases the problem with age trend in the site index model still remains unsolved. A different approach was used by Bøhler & Bernt-Håvar (2011), which proposed a linear correction of site index models in order to take into account the inappropriacy of site index model for Norway spruce resulting from age trend in high altitude semi-natural old stands in Norway.
In summary, an accurate site index model should take into consideration age trends irrespectively of its source, and ignoring these potential biases will cause erroneous assessment of site productivity. Therefore, there is a need to have a practical procedure which allows managers to assess if particular site index models are sufficiently reliable for describing the site productivity of given tree species in areas of interest. In the case of a lack of adequate models, it follows there is a further need to develop more appropriate models by fitting base-age invariant dominant height equation. The decision as to whether existing site index models are sufficiently accurate or whether new site-specific site index models are required for local site conditions should be based on consistent and well developed criteria. However, in the forestry literature the aspect of adequacy of site index models which is very important both from the practical and scientific view point has rarely been analysed. In most cases authors suggested inappropriateness of site index models based on simple comparisons of site index curves obtained from stem analysis or permanent sample plots with site index models or tables. These comparisons assume that the collected tree or plot observations data effectively describes changes in stand height with age. However, the appropriateness of such assumptions is not valid in many cases and just because the measured and modeled height does not match, it does not mean the site index equation is bad.
In this paper we develop a new approach for the assessment of site index models. In the proposed approach the appropriateness of site index models is assessed by analyzing the existence of age trends in residuals developed using forest inventory data and environmental characteristics. According to the developed procedure, a significant correlation between residuals of site index prediction models and stand age indicates the existence of a bias resulting from an age trend and thus the inappropriateness of site index model. We used residual values of site index prediction models instead of a direct correlation of site index models with age in order to remove frequently observed effects of over-representation of older stands on infertile sites or other types of unequal distribution of sample plots in given ages on different environmental gradients. In cases where site index models have a demonstrated age trend we proposed a correction procedure which takes into consideration the age trend found during the validation process by incorporating it into the new site index model.
Material and methodsA case study area and sample design
Norway spruce stands in the Polish part of Western Carpathians (18^{o }48′ 50″ E and 19^{o }58′ 58″ E longitude, 49^{o }23′ 52″ N and 49^{o }41′ 3″ N latitude) formed the basis of this study. Field data were collected from 311 sample plots located in even-aged spruce stands aged 40-157 years. The stands were located between 540-1360 m a.s.l. Sample plots were located in a regular (1250×1250 m) grid and were circular ranging in size from 0.02-0.10 ha. The lowest number of trees per plot was 8 and the maximum 62, with an average of 27. Tree height and diameter at breast height were measured for all trees in the sample plots.
Top height of individual plots was calculated from the mean height of the 100 largest trees per hectare. The age of each plot was determined from at least six dominant trees by counting the number of rings on increment cores extracted using an increment borer at a height of between 20-30 cm above the ground surface, adding 3 years to reach the point of extraction. In general, the difference in age among individual trees in the sample plots was less than 3-5 years.
At each plot the characteristics of topography and the geological substratum were also assessed. Topography was described by the elevation above the sea level (E) measured with the use of an altimeter, the aspect (A) measured in degrees with a compass, the slope (S) and the landform category (Lc). Landform category was derived from a 10×10 m digital terrain model using a combination of topographic position indexes for the scale of 300 and 2000 meters. Geological substratum (Gs) was estimated for individual plots using the geological map in the scale of 1:50 000. For the modeling purpose, the aspect values were transformed as follows (eqn. 1):
K_{B} = sin \left [ A + (90 - A_{max}) \right ] + 1
where K_{B} is the transformed value of aspect, varying from 0 to 2, A is the aspect expressed by azimuth (degrees), A_{max} is the assigned direction (azimuth) of maximum importance, the value A_{max} = 45° was assumed.
In the analysis the relationships between the site index and predictor variables were investigated using generalized additive models (GAM, Hastie & Tibshirani 1990 - eqn. 2).
where Y is the dependent variable, X_{i} are the predictor variables, G is the link function, f_{1}, f_{2}, …, f_{n} are smoothing spline functions estimated from the data, β is a free term, and the errors ε are assumed to have constant variance and mean of 0.
GAM’s were selected as the tool for analysis because they are ranked as one of the most sufficient modeling techniques for the development of site index prediction models (Aertsen et al. 2010). The relationship between residuals of GAM models and age of the stands was analyzed using standard linear models.
The procedure of validation of site index models
The proposed method of validation and correction of site index models was applied to three site index models used in practice for Norway spruce in Poland.
(1) An anamorphic site index model for Norway spruce in Poland (SI_{B}) developed by Bruchwald et al. (1999) using the stem analysis data collected from 548 Norway spruce stands located at the lowland area of North-eastern and mountain area of southern Poland (eqn. 3):
SI_{B} = \frac{H}{A_{T}}
where SI_{B} denotes the site index with a base age of 100 years, H is the top height of the stand, and A_{T} is an age-dependent function defined as (eqn. 4):
(2) The site index model describing Lorey mean height growth trajectories from Schwappach (1943) yield tables (SI_{S}) developed with the use of base age specified method (eqn. 5):
where H_{L} is the Lorey mean height at the age T.
(3) The site index model for Norway spruce in Polish Mountains (SI_{M}) developed using the GADA method based on height growth trajectories obtained from stem analysis of 162 individual trees measured from 53 sample plots located in the Polish Carpathians (Socha 2011 - eqn. 6):
where H is the top height at the age T and R is defined as (eqn. 7):
R = -15.03036 + H + \left [(-15.03036+H)^2 + {\frac{2 H \cdot 66568.71}{ T^{1.95817}}} \right ]^{0.5}
The proposed validation procedure of the site index models described above is outlined as follows. Initially, each site index model was used for the calculation of site index values for all of the sample plots. Two approaches were then developed to assess the adequacy and validity of these site index estimates.
First, the relationship between the site index calculated on the basis of eqn. 3, eqn. 5 or eqn. 6 and the available environmental variables is modeled using the developed GAM. The result of this analysis is a site index prediction model describing the site index as a function of a number n of site characteristics X_{n} (eqn. 2). The adequacy of site index models to the local conditions is validated based on the relationship between the residuals from the model fitted in step 2, and stand age. The existence of a significant correlation between the residuals from the site index prediction model and stand age implies that there are systematic errors in the site index determination.
Secondly, we demonstrated an alternative approach using the age of the stands as one of the independent variables in the site index prediction model (SIPM). In this case the significance of “age” as an independent variable in SIPM demonstrates the inadequacy of site index model. As a result, the height of stands in base ages will be under- or overestimated. The “+” sign of the parameter connected with the variable “age” in SIPM model indicates that site index for young stands is underestimated and overestimated for older stands. In the opposite case a “-” implies overestimation of site index for young and overestimation for old stands.
Procedure of temporary adjusting of site index model to local site conditions by incorporating age trend
Assuming that a site index model is proven invalid, we proposed a simple procedure to adjust the model to the local growth conditions. Using the regression model describing the relationship between residuals of the site index prediction model and the stand age, it is possible to determine the error in the predicted site index as a function of stand age. The relationship between the residuals of the site index prediction model and the stand age can be described both using a simple linear regression model (eqn. 8):
where δ_{SIPM} is the error of site index prediction model, β_{0} is the model parameter, T denotes the stand age and ε the random residual error of the equation.
The site index of the stand at a given age is then corrected by subtracting the correction factor calculated as the difference between the calculated SI and the error described by eqn. 9, which subsequently removes the age trend from site index model (eqn. 9):
SI_{cor} = SI - \delta _{SIPM}
where SI_{cor} denotes the corrected value of site index, SI denotes the value of the site index calculated on the basis of the validated model, and δ_{SIPM} denotes the error of the site index prediction model calculated from eqn. 8.
The above procedure allows systematic errors from inadequate models due to age trend to be minimized. It should be noted, however, that the eliminated systematic errors are also free from the influence of unbalanced distribution of sample plots across different site conditions. Using the proposed approach, it is also possible to adjust invalid site index models using a correction factor I_{cor}, defined as (eqn. 10):
I_{cor} = \left (I \cdot {\frac{SI_{E}}{SI_{cor}}} \right )
where SI_{E} denotes the observed (existing) value of site index, and (eqn. 11):
I =\frac{H_{T}}{H_{BA}}
where H_{T} is the height of the stand at specified age T, and H_{BA} is the height of the stand in a base age.
Corrected value of the stand height trajectory (H_{cor}) in a certain age is calculated by its multiplication by the correction factor (I_{cor}) according to the following equation (eqn. 12):
H_{cor} = I_{cor} \cdot SI_{E}ResultsValidation of Bruchwald’s site index model (SI<sub>B</sub>)
First, a local site index GAM model describing SI_{B} as a function of site characteristics was developed. It was found that the site index calculated using the Bruchwald et al. (1999) model was significantly correlated with site characteristics including altitude, product of slope and elevation (S × K_{B}), geological layer and landform category (Tab. 1), and explains 71.8% of the variability of site index (SI_{B}). Residuals from this site index prediction model were significantly correlated with the stand age (T), which explain over 34% of the variability (Fig. 1, Tab. 2). This relationship was described by linear regression model (eqn. 13):
\delta _{SI_{B}} = -0.06561 \cdot (T-100)
where δ_{SIB} denotes the error of site index estimation using the site index prediction model, and T the age of the stand.
Using the second approach that included “age” into site index prediction model, the proportion of explained variance significantly increased to 81.6% (Tab. 2).
Applying our approach, the observed significant correlation between residuals of the site index prediction model and the stand age, as well as the significance of the age variable in the SIPM age trend, implied an age trend bias and indicated the inadequacy of the SI_{B} model under local site conditions. This inadequacy results in systematic errors in site index estimation.
Validation of site index model by Schwappach (SI<sub>S</sub>)
Site indexes SI_{S} were calculated using the site index model by Schwappach (1943) and used as dependent variables in a GAM site index prediction model. In this case, the variance explained by the model (80.4% - Tab. 1) is substantially higher than the variance of the site index prediction model developed using SI_{B} as the dependent variable. By including the variable “age” in this site index prediction model, the variance only slightly increases to 81.6%. Also, stand age was a significant variable when added to the SIMP equation. However, residuals of the site index prediction model are weakly correlated with the stand age (T) explaining only about 17% the total variation (Fig. 2, Tab. 2).
Validation of the site index model for Norway spruce Polish ountains (SI<sub>M</sub>)
Lastly, the site index prediction GAM model for SI_{M} (eqn. 5) as function of environmental factors was fitted (Tab. 1). Again, the validation of the SI_{M} model was based on the analysis of the relationship between residuals of the site index prediction GAM model (Tab. 1) and stand age. In this case, no significant correlation was detected between the residuals of site index prediction model developed with the use of SI_{M} and stand age (Fig. 3, Tab. 2).
Considering stand age in site index prediction model with SI_{M} as the independent variable also showed no increase in the explained variance. As a result, in this case, there was no systematic error in site index estimation related to age of the stands (Fig. 3, Tab. 2), thus the site index is stable along the stand age.
The example of temporary correction of SI<sub>B</sub> model to local conditions
The results of the analysis indicated a strong age trend in site index resulting from the use of the SI_{B} model in local conditions. Therefore, the proposed correction procedure equation describing the dependence of expressed residuals of SIPM on stand age was used to adjust the height development patterns of the SI_{B} model. Examples of the correction are graphically presented for individual SI curves 20m, 25m, 30m, 35m in Fig. 4. In general, in the case of medium and good sites showing a site index over 25m, there is close similarity between the corrected SI_{B} and the local SI_{M} site index models (Fig. 4). Only the least productive stands had large differences (Fig. 4).
Height increment analysis indicates that, in the case of SI model developed for local conditions of Polish mountains for stands up to 45 years old, the height growth rate is higher than that assumed in the model SI_{B}, whereas for stands over 50 years the height growth rate is lower than expected by the validated model. Especially, large differences in site index curves are observed for the oldest stands (> 100 years - Fig. 4). The above mentioned differences are the reason of the overestimation of site index for younger stands and underestimation of site index for stands in ages older than base age (100 years).
Discussion
Inappropriateness of site index models manifested by age trend has been widely observed (Spiecker et al. 1996, Subedi et al. 2009, Bøhler & Bernt-Håvar 2011, Socha & Orzel 2011). The existence of correlations between site index and stand age have been frequently reported and most often an age trend has been interpreted as driven by changes in site conditions (Elfving & Tegnhammar 1996, Kiviste 1999, Bravo-Oviedo et al. 2010, Nothdurft et al. 2012, Yue et al. 2014). However, other important sources of site index model inappropriateness are reported, such as local differences in top height trajectories resulting from specific site conditions (Johansson 1995), genetic variability (Buford & Burkhart 1987, Adams et al. 2006) and inappropriate data, model or methods used in development of site index model (Cieszewski 2001, Cieszewski & Strub 2007, García 2010).
Site index estimated on the basis of stand height and age is only an approximation of an unobservable site variable, which express the site conditions for given tree species. As a consequence of the aforementioned definition, it can be assumed that:
Environmental drivers and the associated plant physiological response remain constant over the life of the stand, without changes in genetic structure, forest management practices or any other modifying growth factors. Given the same site conditions, site index values calculated for young stands should be equal to site index of older stands. Therefore, in such cases the site index curves describing change in height with age for young and old stands should be similar.
In the case of deteriorated or improved site conditions (such as those caused by improved climatic conditions, CO_{2} fertilization or nitrogen deposition, genetic improvement, change in forest management practices or other factors influencing growth rate), younger stands may grow faster and therefore trajectories describing change in stand height with age for younger stands are different from those for older stands (Nothdurft et al. 2012, Yue et al. 2014). As the consequence of the change in site productivity over time, models describing change in stand height with age are not explicitly linked to site index. In such cases site index models should take into account time trend in site productivity in order to obtain comparable site indexes for sites characterized with comparable site productivity, irrespectively of the age of stand growing in a given location. Otherwise, the use of site index is restricted only to stands in a specific age range.
The proposed procedure assumes that the residuals from the SIPM describing the site index as the function of site properties should not be correlated with stand age. Site index prediction models were used in this case because direct correlation between the site index and stand age may also result from the irregularity in localization of stands of varying ages on different site conditions, which is commonly observed (Nothdurft et al. 2012). It is particularly likely in mountainous areas, where stands at the highest altitudes under poorer site conditions are frequently older than stands at lower altitudes in richer sites. This irregularity results from the unequal accessibility to stands located at different altitudes and possibly from differences in the rotation age.
Most authors agree that observed age trend in site index estimates are the result of changing site conditions (Nothdurft et al. 2012, Pretzsch et al. 2012, Pretzsch et al. 2014, Waring et al. 2014, Yue et al. 2014). However, in this paper more localized site index models developed for the Carpathians did not revealed correlations with the age of stands. Likely, this is caused by the local model built using data from trees growing under environmental changes. Moreover, a large component of the stem data used to develop the site index model comes from relatively young (60-70 years old) stands (Socha 2011). Thus, most trees were grown under the ameliorating effect of nitrogen deposition, increased temperature and CO_{2} concentration, and this may result in marked changes of site conditions as compared with the early 1950s (Yue et al. 2014). Likewise, in lowland Scots pine stands in southern Poland (Socha & Orzel 2011), site indexes calculated using site index models developed for Poland (Bruchwald et al. 2000) and Schwappach’s (Schwappach 1943) yield tables indicated strong correlation between site index and the age of stands. However, after the use of site index model developed for local site conditions, substantial overestimation of site indexes for young stands disappeared (Socha & Orzel 2011).
Strong correlation between site index and age of stands was apparent only in the case of site index model SI_{B}. We assume that the most probable cause of this is that this model was developed for the whole population of Norway spruce in Poland, including both for mountain and lowland sites. It could be expected that in mountain areas trees in young stands initially grow relatively faster than in older stands. Therefore, the use of site index models appropriate for lower elevations applied to high mountain forests may result in an age trend. Comparably different growth trajectories for high- and low-elevation stands were observed between subalpine fir, Engelmann spruce, lodgepole pine and black spruce in British Columbia (Klinkal et al. 1996, Chen & Klinka 2000, Nigh et al. 2002). Evidence of different growth in mountain Norway spruce stands was also observed in Norway by Bøhler & Bernt-Håvar (2011).
The most problematic issue with the proposed validation procedure is the difficult separation of the effect of change in site conditions from the effect of using inappropriate site index models on observed age trends in site indexes. However, regardless of the source of the site index age trends, the correlation between residuals of site index prediction models and the age of stands indicates inappropriate site index predictions, which makes the proposed approach useful in both cases. Moreover, both in the case of changing site productivity and inappropriate site index model, the usefulness of site index predictions is questionable.
Climate-based site index models may also represent an improved alternative, which take into account changes in site parameters such as temperature, precipitation, drought length, etc. (Bravo-Oviedo et al. 2008, Monserud et al. 2008). However, both traditional site index models describing site index as the function of stand age and height, and site index models including climatic parameters as additional independent variables in the case of changing site conditions, gives the assessment of the site index dependent on the age of the stand.
Most of the recently published papers on forest productivity demonstrate the high probability of further changes in site conditions around the world (Coops et al. 2007, Coops et al. 2010, Albert & Schmidt 2010, Nothdurft et al. 2012). Therefore, one of the most important aspect for forest management is to attain a reliable evaluation of site index models appropriateness and its correction to local growth conditions.
Conclusions
Age trends may result from inappropriateness of site index models, often driven by: (a) local differences in top height trajectories resulting from specific site conditions, forest management practice or genetic variability; (b) inappropriate model or methods used in the development of the site index model; and (c) changes in site conditions. Relationship between the residuals of site index prediction models describing site index as a function of environmental factors and age of the stands may be used for the validation and correction of site index models in forest inventory. The existence of a correlation between residuals of site index prediction models and the age of stands implies the occurrence of systematic errors in site index determination. Such model inadequacy results in the under-/overestimation of the height of the stand in base ages. The regression model describing the relationship between the residuals of site index prediction model and the age of the stand may be used to detect the age trend in site index model, and to make it more suitable to local site conditions.
Acknowledgments
This study was financially supported by the Polish National Science Center, Grant NN-309-066939. JS acknowledges the support by the University of British Columbia (Vancouver, BC, Canada) and the University of Agriculture in Krakow (Poland).
ReferencesAdams JP, Matney TG, Land SB, Belli KL, Duzan HWIncorporating genetic parameters into a loblolly pine growth-and-yield model. Canadian Journal of Forest Research 36: 1959-1967.2006Aertsen W, Kint V, Van Orshoven J, Ozkan K, Muys BComparison and ranking of different modelling techniques for prediction of site index in Mediterranean mountain forests. Ecological Modelling 221: 1119-1130.2010Albert M, Schmidt MClimate-sensitive modelling of site-productivity relationships for Norway spruce (Picea abies (L.) Karst.) and common beech (Fagus sylvatica L.). Forest Ecology and Management 259: 739-749.2010Alvarez-González JG, Ruiz-González AD, Rodríguez-Soalleiro R, Barrio-Anta MEcorregional site index models for Pinus pinaster in Galicia (northwestern Spain). Annals of Forest Science 62: 115-127.2005Bøhler F, Bernt-Håvar ØEstimation of site-index in old, semi-natural stands of Norway spruce at high altitude. In: Proceedings of the Conference “Forest Management and Silviculture in the North - Balancing Future Needs”. Stjørdal (Norway) 6-8 Sep 2011, Book of Abstracts. Rapport Fra Skog Og Landskap, Norway, pp. 53.2011Bontemps JD, Bouriaud OPredictive approaches to forest site productivity: recent trends, challenges and future perspectives. Forestry 87: 109-128.2013Bravo-Oviedo A, Tomé M, Bravo F, Montero G, Del Río MDominant height growth equations including site attributes in the generalized algebraic difference approach. Canadian Journal of Forest Research 38: 2348-2358.2008Bravo-Oviedo A, Gallardo-Andrés C, Del Río M, Montero GRegional changes of Pinus pinaster site index in Spain using a climate-based dominant height model. Canadian Journal of Forest Research 40: 2036-2048.2010Bruchwald A, Dudek A, Michalak K, Rymer-Dudzinska T, Wróblewski L, Zasada MModel wzrostu dla drzewostanów swierkowych [A growth model for Norway spruce stands] . Sylwan 143 (1): 19-31. [in Polish]1999Bruchwald A, Michalak K, Wróblewski L, Zasada MWzrost wysokosci sosny w róznych regionach Polski [Height growth of Scots pine in various regions of Poland] . In: “Przestrzenne zróznicowanie wzrostu sosny” [Spatial variation in Scots pine growth] (Bruchwald A ed). Fundacja Rozwój SGGW, Warszawa, Poland, pp. 77-83. [in Polish]2000Buford MA, Burkhart HEGenetic improvement effects on growth and yield of Loblolly pine plantations. Forest Science 33: 707-724.1987Burkhart HE, Tomé MModeling forest trees and stands. Springer-Verlag, Dordrecht, The Netherlands, pp. 471.2012Chen HYH, Klinka KHeight growth models for high-elevation subalpine fir, Engelmann spruce, and Lodgepole pine in British Columbia. Western Journal of Applied Forestry 15: 62-69.2000Cieszewski CJThree methods of deriving advanced dynamic site equations demonstrated on inland Douglas-fir site curves. Canadian Journal of Forest Research 31: 165-173.2001Cieszewski CJ, Strub MParameter estimation of base-age invariant site index models: which data structure to use? A discussion. Forest Science 53: 552-555.2007Coops NC, Coggins SB, Kurz WMapping the environmental limitations to growth of coastal Douglas-fir stands on Vancouver Island, British Columbia. Tree Physiology 27: 805-815.2007Coops NC, Hember R, Waring RHAssessing the impact of current and projected climates on Douglas-fir productivity in British Columbia, Canada, using a process-based model (3-PG). Canadian Journal of Forest Research 40: 511-524.2010Elfving B, Tegnhammar LTrends of tree growth in Swedish forests 1953-1992: an analysis based on sample trees from the national forest inventory. Scandinavian Journal of Forest Research 11: 26-37.1996Fries A, Lindgren D, Ying CC, Ruotsalainen S, Lindgren K, Elfving B, Karlmats UThe effect of temperature on site index in western Canada and Scandinavia estimated from IUFRO Pinus contorta provenance experiments. Canadian Journal of Forest Research 30: 921-929.2000García ODynamical implications of the variability representation in site-index modelling. European Journal of Forest Research 130: 671-675.2010Hasenauer H, Nemani RR, Schadauer K, Running SWForest growth response to changing climate between 1961 and 1990 in Austria. Forest Ecology and Management 122: 209-219.1999Hastie TJ, Tibshirani RJGeneralized additive models. CRC Press, Boca Raton, FL, USA, pp. 352.1990Johansson TSite index curves for Norway spruce plantations on farmland with different soil types. Studia Forestalia Suecica 198: 1-19.1995Kiviste ASite index change in the 1950s-1990s according to Estonian forest inventory data. In: Proceedings of the International Seminar “Causes and Consequences of Accelerating Tree Growth in Europe” (Karjalainen T, Spiecker H, Laroussinie O ed). Nancy (France) 14-16 May 1998. EFI Proceedings no. 27, Joensuu, Finland, pp. 87-100.1999Klinkal K, Carter RE, Chert HYHHeight growth-elevation relationships in subalpine forests of interior British Columbia. The Forestry Chronicle 72: 193-198.1996Monserud RA, Yang Y, Huang S, Tchebakova NPotential change in lodgepole pine site index and distribution under climatic change in Alberta. Canadian Journal of Forest Research 38: 343-352.2008Monserud RAHeight growth and site index curves for inland Douglas-fir based on stem analysis data and forest habitat type. Forest Science 30: 943-965.1984Monserud RA, Rehfeldt GEGenetic and environmental components of variation of site index in Inland Douglas-fir. Forest Science 36: 1-9.1990Nigh GD, Krestov PV, Klinka KHeight growth of black spruce in British Columbia. The Forestry Chronicle 78: 306-313.2002Nothdurft A, Wolf T, Ringeler A, Böhner J, Saborowski JSpatio-temporal prediction of site index based on forest inventories and climate change scenarios. Forest Ecology and Management 279: 97-111.2012Pretzsch H, Biber P, Schütze G, Uhl E, Rötzer TForest stand growth dynamics in Central Europe have accelerated since 1870. Nature Communications 5: 4967.2014Pretzsch H, Dieler J, Seifert T, Rötzer TClimate effects on productivity and resource-use efficiency of Norway spruce (Picea abies [L.] Karst.) and European beech (Fagus sylvatica [L.]) in stands with different spatial mixing patterns. Trees 26: 1343-1360.2012Pretzsch H, Grote R, Reineking B, Rötzer T, Seifert SModels for forest ecosystem management: a European perspective. Annals of Botany 101: 1065-1087.2008Schadauer KSite index and site characteristics of Norway spruce according to data of the Austrian Forest Inventory. Report 171, Mitteilungen der Forstlichen Bundesversuchsantalt, Wien, Austria, pp. 38.1999Schwappach AErtragstafeln der wichtigeren Holzarten [Yield tables for the most important forest tree species] . Druckerei Merkur, Prague, Czech Republic, pp. 49-55.1943Sharma RP, Brunner A, Eid TSite index prediction from site and climate variables for Norway spruce and Scots pine in Norway. Scandinavian Journal of Forest Research 27: 619-636.2012Sharma RP, Brunner A, Eid T, Oyen BHModelling dominant height growth from national forest inventory individual tree data with short time series and large age errors. Forest Ecology and Management 262: 2162-2175.2011Skovsgaard JP, Vanclay JKForest site productivity: a review of the evolution of dendrometric concepts for even-aged stands. Forestry 81: 13-31.2008Skovsgaard JP, Vanclay JKForest site productivity: a review of spatial and temporal variability in natural site conditions. Forestry 86: 305-315.2013Socha JEffect of topography and geology on the site index of Picea abies in the West Carpathian, Poland. Scandinavian Journal of Forest Research 23: 203-213.2008Socha JSite index curves for Norway spruce on mountain habitats. Sylwan 155: 816-826.2011Socha J, Orzel SDynamiczne krzywe bonitacyjne dla drzewostanów sosnowych Puszczy Niepolomickiej [Dynamic site index curves for Scots pine stands in NiepoÅomice primeval forest] . Sylwan 155: 301-312. [in Polish]2011Solberg S, Dobbertin M, Reinds GJ, Lange H, Andreassen K, Fernandez PG, Hildingsson A, De Vries WAnalyses of the impact of changes in atmospheric deposition and climate on forest growth in European monitoring plots: a stand growth approach. Forest Ecology and Management 258: 1735-1750.2009Spiecker H, Mielikäinen K, Kohl M, Skovsgaard JGrowth trends in European forests. Springer-Verlag, Berlin, Heidelberg, New York, pp. 41-88.1996Splechtna BEHeight growth and site index models for Pacific silver fir in southwestern British Columbia. Journal of Ecosystem Management 1: 1-14.2001Subedi N, Sharma M, Parton JAn evaluation of site index models for young Black spruce and Jack pine plantations in a changing climate. Climate Change Research Report, Ministry of Natural Resources, Ontario, Canada, pp. 15.2009Tegnhammar LOn the estimation of site index for Norway spruce. Report 53, Department of Forest Survey, Swedish University of Agricultural Sciences, Uppsala, Sweden, pp. 259.1992Waring R, Coops N, Mathys A, Hilker T, Latta GProcess-based modeling to assess the effects of recent climatic variation on site productivity and forest function across Western North America. Forests 5: 518-534.2014Yue C, Mäkinen H, Klädtke J, Kohnle UAn approach to assessing site index changes of Norway spruce based on spatially and temporally disjunct measurement series. Forest Ecology and Management 323: 10-19.2014
Relationship between residual values (obs-pred) of the site index prediction model based on validated site index model SI_{B} (Bruchwald et al. 1999) and the age of stands.
Relationship between residual values (obs-pred) of site index prediction model developed for site indexes calculated on the basis of Schwappach (1943) site index model SI_{S} and age of the stands.
Relationship between residual values (obs-pred) of site index prediction model based on local site index model SI_{M} for Polish mountains (Socha 2011) and stand ages.
Actual (SI_{B} - dotted line) and corrected (SI_{Bcor} - solid line) site index curves drawn on the basis of the SI_{B} model by Bruchwald et al. (1999) against a background of site index curves drawn with the use of local site index model SI_{M} (grey line) for Polish mountains (Socha 2011).
Site index prediction models describing the relationship between site indexes calculated using the assessed site index models used for Norway spruce in Poland (SI_{B}, SI_{S} and SI_{M}) and selected predictor variables developed using generalized additive models techniques.
Dependent variable - site index model(model source)
Predictor variables
AdjustedR^{2}
SI_{B} (Bruchwald et al. 1999)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs)
71.8
SI_{B} (Bruchwald et al. 1999)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs), stand age (T)
81.6
SI_{S} (Schwappach 1943)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs)
80.4
SI_{S} (Schwappach 1943)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs), stand age (T)
81.6
SI_{M} (Socha 2011)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs)
81.2
SI_{M} (Socha 2011)
elevation above sea level (E), product of slope and aspect (S×K_{B}), landform category (Lc), geological substratum (Gs), stand age (T)
81.5
Characteristics of the models describing the relationship between residuals of site index prediction models and the age of stands (eqn. 7).