Next Article in Journal
Mechanical Characteristics of the Fine Roots of Two Broadleaved Tree Species from the Temperate Caspian Hyrcanian Ecoregion
Previous Article in Journal
Forwarder Productivity in Salvage Logging Operations in Difficult Terrain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Individual Distance-Independent Diameter Growth Models for Japanese Cedar (Cryptomeria japonica) Trees under Multiple Thinning Treatments

Forestry and Forest Products Research Institute, Forest Research and Management Organization, 1 Matsunosato, Tsukuba, Ibaraki 305-8687, Japan
*
Author to whom correspondence should be addressed.
Forests 2020, 11(3), 344; https://doi.org/10.3390/f11030344
Submission received: 24 January 2020 / Revised: 17 March 2020 / Accepted: 18 March 2020 / Published: 19 March 2020
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Understanding the tree growth process is essential for sustainable forest management. Future yields are affected by various forest management regimes such as thinning; therefore, accurate predictions of tree growth are needed under various thinning intensities. This study compared the accuracy of individual-level distance-independent diameter growth models constructed for different thinning intensities (thinning intensity-dependent multiple models: TDM model) against the model designed to include all thinning intensities (thinning intensity-independent single model: TIS model) to understand how model accuracy is affected by thinning intensity. We used long-term permanent plot data of Japanese cedar (Cryptomeria japonica) stands in Japan, which was gathered from four plots where thinning was conducted at different thinning intensities: (1) intensive (41% and 38% of trees removed at 25 and 37 years old, respectively), (2) moderate (38% and 34%), (3) light (32% and 34%), and (4) no thinning. First, we specified high interpretability distance-independent competition indices, and we compared the model accuracy both in TDM and TIS models. The results show that the relative spacing index was the best competition index both in TDM and TIS models across all thinning intensities, and the differences in the RMSE (Root mean square error) and rRMSE (relative RMSE) in both TDM and TIS models were 0.001–0.01 cm and 0.2–2%, respectively. In the TIS model, rRMSE varied with thinning intensity; the rRMSE was the lowest for moderate thinning intensity (45.8%) and the highest for no thinning (59.4%). In addition, bias values were negative for the TIS model for all thinning intensities. These results suggest that the TIS model could express diameter growth regardless of thinning intensities. However, the rRMSE had varied with thinning intensity and bias had negative values in the TIS model. Therefore, more model improvements are required for accurate predictions of long-term growth of actual Japanese cedar stands.

1. Introduction

Yield predictions are essential for sustainable forest management. Generally, future yields are affected by different forest management regimes, site conditions, other climate factors, and species [1,2]. Such growth processes are described by growth models [3]. Growth models can be classified into stand-level and individual-level models. The stand-level model, as the basis of growth prediction, has been studied for a long time [2]. This model can characterize yield and growth with a little information about the stand; therefore, many researchers have adapted this model for a variety of stands [4,5]. Individual-level models predict single-tree growth processes. These models incorporate flexible tree growth processes; thus, this model is robust in predicting tree growth under various conditions [3]. In addition, in contrast to the stand-level model, the individual-level model can predict the volume and other parameters of each diameter class within a stand by summarizing individual-tree estimates [6,7,8,9,10,11].
Individual-level models are further classified into distance-dependent and distance-independent models. Tree growth in a stand is restricted by stand age, stand condition, and competition with adjacent trees for light, water, and nutrients. This competition is measured by the competition index. Distance-dependent models account for distance between the subject tree and adjacent trees, so the model needs information regarding the trees’ positions in the stand. Generally, the distance-dependent competition index is defined as the weighted distance of competitor trees that are located within a certain distance from the subject tree. The distance-independent model uses a competition index that does not require the spatial information of trees. Rather, in this model, size ratios and stand density are used as a substitute for competition. In general, the accuracy of the distance-dependent model is higher than the distance-independent model because the distance-dependent model includes distance information regarding competitor trees [9,12,13,14]. On the other hand, some studies have indicated that the accuracy of the distance-dependent model is the same as the distance-independent model regardless of target species and region [15,16,17,18,19]. Therefore, the distance-independent model may sufficiently explain tree growth. The distance-dependent model is cost intensive and labor demanding because the competition index in the model needs detailed tree location information to generate calculations. If the distance-independent model, which does not require detailed tree location information, can sufficiently predict individual-tree growth, the yield prediction can be conducted in a cost-effective and labor-saving manner.
We should consider the effect of thinning on tree growth when predicting long-term tree growth. Several studies created diameter growth models under different thinning intensities [12,20,21]. For example, Uzoh and Oliver [20] designed a single diameter growth model using the site density index (SDI) and the basal area of larger trees than the subject tree (BAL). If a single model can explain the response of tree growth under different thinning intensities sufficiently well, we could make more efficient predictions. It is known that the dynamics of tree growth change under different thinning intensities [22,23,24,25]. For example, Mäkinen and Isomäki [22] reported that the tree height and diameter growth of Scots pine (Pinus sylvestris L.) were better in stands with higher thinning intensity. Thus, the effect of thinning intensity on tree growth may change under multiple thinning intensities. Model accuracy may differ under different thinning intensities; therefore, we need to construct a growth model for each thinning intensity. Previous studies have developed a single-tree growth model under different thinning intensities [22,23]; however, few studies have independently developed multiple tree growth models for each thinning intensity and then compared the prediction accuracies of single and multiple tree growth models. Even if a single model can predict tree growth under different thinning intensities, the model may not have the same precision and bias for stands with different thinning intensities. In such cases, the substantial precision and bias may differ between predicted and observed values when the model is applied to actual stands. Therefore, we need to examine the accuracy of both single and multiple models.
Our main objective is to compare the accuracy of the individual-level distance-independent diameter growth model constructed with each thinning intensity (thinning intensity-dependent multiple models: TDM model) to the accuracy of the model constructed with all thinning intensities (thinning intensity-independent single model: TIS model). To achieve this objective, we did the following: (1) selected the model effective among six distance-independent competition indices for predicting diameter growth of Japanese cedar (Cryptomeria japonica), (2) compared the accuracy of the best TDM and TIS models, and (3) examined how the accuracy of the best model changes under different thinning intensities.

2. Materials and Methods

2.1. Data Collection

The data used in this study were obtained from Kudarukawa national forest, located in Kochi prefecture, western Japan (33°23’ N, 133°10’ E) (Figure 1). The study site is located in areas with an inclination of 5°10° at an elevation of 500 m. In 1972, four permanent plots were established in planted forests with different thinning intensities, which were: (1) intensive thinning, (2) moderate thinning, (3) light thinning, and (4) no thinning [26]. Each plot included approximately 200 Japanese cedar trees. The first thinning was conducted at 25 years old, and the second thinning was conducted at 37 years old (Table 1). Trees’ heights and diameters at breast height (DBH) in each plot were measured every 56 years from 14 years old to 61 years old. The DBH of each tree was measured at breast height position (1.2 m) from two directions using a caliper, whereas tree height was measured for approximately 30 selected trees with various DBH classes in each plot (Table 2). A tree number and breast height line were painted for each tree when tree measurements were conducted. Before our analysis, we compared the mean annual diameter growth for 47 years among four plots using a multiple-comparison test. The mean annual diameter growth of the unthinned plot was significantly smaller than that of the other thinned plots (Tukey–Kramer test, p < 0.05). Moreover, we estimated the heights of unmeasured trees using the diameter–height equation, created using the data of measured trees (R2 values in each thinning intensity were 0.88–0.99). To calculate the competition indices, we assumed that thinning was conducted in prior measurement periods in this study. In total, we used nine growing seasons of 1077 trees’ growth data.

2.2. Individual-Tree Diameter Growth Model

We used the linear-mixed effects model to characterize Japanese cedar diameter growth. We defined individual diameter growth of the TDM model as follows:
l n ( G I ,   i , j + 1 ) = α 0 + α 1 D B H I , i , j + α 2 D B H I , i , j 2 + α 3 A g e I ,   i , j + φ i + ε I , i , j ,
where G I , i , j is the growth of the subject tree i in the I th plot in the period between the jth measurement and the subsequent measurement. D B H I , i , j and D B H I , i , j 2 are DBH and the DBH squared, respectively, and A g e I , i , j is the stand’s age. Generally, tree diameter growth is affected by individual diameter size [27]. Moreover, tree diameter growth follows a sigmoid curve, indicating the existence of restraint by aging and increasing DBH [3]. To express this restraint, we used a simple quadratic function of DBH. α 0 , …, α 3 are each parameters, and φ i is a random parameter for subject tree i. ε I , i , j is the error term produced by the normal distribution with mean 0 and variance σ 2 ( ε I , i , j ~ Normal [ 0 , σ 2 ] ). The annual diameter growth in the period between two successive measurements was expressed using DBH, D B H 2 , and the age at the beginning of the period. A full TDM model derived by adding each distance-independent competition index C I I , i , j is given by
l n ( G I ,   i , j ) = α 0 + α 1 D B H I , i , j + α 2 D B H I , i , j 2 + α 3 A g e I ,   i , j + α 4 C I I ,   i , j + φ i + ε I , i , j ,
where α 0 ,   , α 4 are fixed effects parameters, and C I I , i , j is the distance-independent competition index of the Ith plot in subject tree ith at jth measurements. We estimated the parameters of Equation (2) for each of the four plots, with varying thinning regimes. Thus, we obtained a total of 20 fixed effects parameters for the TDM model. Corresponding to Equations (1) and (2), the TIS model was defined as follows:
l n ( G I ,   i , j ) = α 0 + α 1 D B H I , i , j + α 2 D B H I , i , j 2 + α 3 A g e I ,   i , j + φ I + φ I , i + ε I , i , j ,
l n ( G I ,   i , j ) = α 0 + α 1 D B H I , i , j + α 2 D B H I , i , j 2 + α 3 A g e I ,   i , j + α 4 C I I ,   i , j + φ I + φ I , i + ε I , i , j ,
where φ I , i is the random parameter for the Ith plot in subject tree i. The TIS model is a single model that expresses diameter growth for all plots because the model includes the plot-level random effect φ I . In this study, we used six distance-independent competition indices as described in the following section. The generalized liner mixed effects model was used in the lme4 package [28] in R version 3.6.0 [29]. Also, we used a Correction Factor (CF = exp((MSE)/2, where MSE is the mean square error) to predict expected diameter growth. We applied the CF to correct bias occurring in the log transformation [30].

2.3. Distance-Independent Competition Index

We used six distance-independent competition indices: basal area of trees larger than the subject tree (BAL), diameter ratio (DR), basal area ratio (BR), cumulative distribution function (CDF), relative spacing index (Sr), and SDI according to the definition of Sun et al [11].

2.3.1. Basal Area of Trees Larger than the Subject Tree (BAL)

BAL is defined as
BAL I , i , j = π 4 D B H c 2 ,
where BAL I , i , j is the sum of the basal area of competitor trees for subject tree i in Ith plot at the jth measurement, and D B H c is competitor trees that have larger diameter than subject tree i in the Ith plot at the jth measurement [14,31].

2.3.2. Diameter Ratio (DR)

DR is defined as
DR I , i , j = D B H I , i , j Q M D I , j ,
where DR I , i , j is the DR of Ith plot in subject tree i at the jth measurement. D B H I , i , j is the DBH of subject tree i at the jth measurement, and Q M D I , j is the quadratic mean diameter of the Ith plot at the jth measurement [32].

2.3.3. Basal Area Ratio (BR)

BR is defined as
BR I , i , j = B A I , i , j B A ¯ I , j = D B H I , i , j 2 Q M D I , j 2 ,
where BR I , i , j and B A I , i , j are the BR and BA of the Ith plot in subject tree i at the jth measurement, respectively. B A ¯ I , j is the mean BA of the Ith plot at the jth measurement. In this study, we used the square of Equation (5) [11].

2.3.4. Cumulative Distribution Function (CDF)

CDF is usually used to express the diameter distribution in the stand. Also, CDF expresses the relative position ( F I , i , j ) of the Ith plot in subject tree i at the jth measurements as follows:
ln ( F I , i , j ) = k 0.5 N I , j ,
where k is the rank of the size of subject tree i ordered from smallest to largest in the Ith plot. N I , j is the number of trees per hectare in the Ith plot at the jth measurement.

2.3.5. Relative Spacing Index (Sr)

Sr is a stand density index that is calculated as the ratio of distance between stand trees and mean height as follows:
S r I , j = 100 2 H ¯ I , j N I , j ,
where S r I , j and H ¯ I , j are Sr and mean height in the Ith plot at the jth measurement, respectively [33].

2.3.6. Site Density Index (SDI)

Similarly, SDI is expressed by the number of trees in the stand and quadratic mean diameter as follows [34]:
SDI I , j = N I , j ( Q M D I , j 25 ) c ,
where c is a constant defined by a self-thinning model. In this study, we used c = 1.486, which was evaluated in a Japanese cedar stand [35,36].

2.4. Model Fitting and Evaluation

We evaluated the goodness of the fit of the model for each competition index based on the marginal R2, conditional R2, Akaike’s information criteria (AIC), and ΔAIC. Marginal R2 and conditional R2 were calculated using observed and predicted diameter growth. Marginal R2 was calculated with only the fixed effect, and conditional R2 was calculated with both fixed and random effects [37]. ΔAIC is the difference between the AIC of the best model and each other model. Larger values of marginal R2 and conditional R2 indicate better model accuracy; conversely, a lower value of AIC indicates better model accuracy. Additionally, to evaluate the effectiveness of the competition indices under different thinning intensities, we calculated the root mean square error (RMSE), relative RMSE (rRMSE), bias (mean error), and variance error using the best model. rRMSE is calculated as RMSE divided by mean diameter of observed diameter growth, and a smaller rRMSE indicates better prediction accuracy [38]. Accuracy is higher when bias is closer to 0. Negative bias values indicate underestimation, and vice versa.

3. Results

3.1. The Goodness of Fit for Distance-Independent Competition Index

In the TDM model, the best model included Sr for each thinning intensity based on marginal R2, conditional R2, AIC and ΔAIC (Table 3). Then, the marginal R2 and conditional R2 of each thinning intensity were 0.36–0.49 and 0.62–0.70, respectively (Table 3). The second-best model included BAL in intensive and light thinning, and SDI in moderate thinning and no thinning (Table 3). Coefficients of the TDM model in all thinning intensities showed that DBH, DBH2, age, and Sr were 0.174 to 0.213, −0.002 to −0.001, −0.056 to −0.030, and 0.161 to 0.209, respectively (Table 4). The TIS model showed that the best model included Sr (Table 5). Then, for the best model, marginal R2 and conditional R2 were 0.42 and 0.67, respectively (Table 5). Also, the coefficients of the best model showed that DBH, DBH2, age, and Sr were 0.177, −0.002, −0.040, and 0.163, respectively (Table 6).

3.2. Comparison of Model Accuracy between TDM Model and TIS Model

Visual assessment of the relationships between actual and predicted DBH growth showed little difference in prediction performance between TDM and TIS models of each thinning intensity (Figure 2). RMSE and rRMSE also differed little for both TDM and TIS models of all thinning intensities (Table 7). The differences in RMSE and rRMSE among TDM and TIS models were 0.001–0.01 and 0.2–2%, respectively (Table 7). The bias of the TDM and TIS models specifically differed in no thinning, with a difference of 0.032. Variance decreased with decreasing thinning intensity, and the lowest variances for no thinning in both TDM and TIS models were 0.039 and 0.033, respectively.

3.3. Model Accuracy under Different Thinning Intensities

For the TDM model, the lowest RMSE value was 0.203 for no thinning, and the highest RMSE value was 0.249 for intensive thinning (Table 7). rRMSE was the lowest for moderate thinning (47.8%), and the highest for no thinning (60.6%). Bias values were negative for all thinning intensities, indicating that the predicted values were underestimated compared with the observed values. For the TIS model, RMSE values showed that the lowest value was 0.199 for no thinning, and the highest value was 0.250 for intensive thinning (Table 7). rRMSE values were lowest for moderate thinning (45.8%), and highest for no thinning (59.4%). Bias values were negative for all thinning intensities, indicating that the predicted values were underestimated compared to the observed values.

4. Discussion

4.1. Model Selection for TDM and TIS Models

Several studies that used distance-independent competition indices showed that the best competition indices were BAL, SDI, and CDF [6,10,20]. Goodness of fit of competition indices may be species-specific. In this study, Sr was the best competition index both in the TDM and TIS models (Table 3 and Table 5). This result indicates that individual level competition is expressed as other explanatory variables such as DBH and DBH squared because Sr is expressed as stand level competition. Moreover, a previous study showed that Sr was related to crown length [39]. Tree crown is widely used for predicting tree growth, and there is high correlation between crown and tree growth [9,40]. Sr may indirectly reflect information of the tree crown; therefore, the best models might perform well. Moreover, tree height is suited to express site productivity. Sr could reflect differences in productivity among each study plot, therefore increasing the accuracy of the model. In the TDM model, the second-best models were different for each thinning intensity (Table 3). BAL was selected in intensive and light thinning, and SDI was selected in moderate and no thinning. This result suggests that BAL and SDI may also be possible to express diameter growth of Japanese cedar under different thinning intensities.

4.2. Diameter Growth Model Accuracy for TDM and TIS Models

Previous studies evaluated model accuracy using indexes such as RMSE, rRMSE, and AIC [9,10,41]. For example, Palahí et al. [41] constructed an individual-tree distance-independent diameter growth model for Scots pine (Pinus sylvestris L.) in northeast Spain using long-term growth data. They calculated model accuracy using rRMSE, and the rRMSE in their model was 64.1%. On the other hand, in our study, the rRMSE values of TDM and TIS models were 45.8–60.6% (Table 7). Thus, our model accuracy was not inferior in comparison with previous studies.
There was little difference in RMSE and rRMSE among TDM and TIS models (Table 7). Also, the best competition index was the same for both TDM and TIS models (Table 3, Table 4, Table 5 and Table 6). These results indicate that the TIS model can explain diameter growth under different thinning intensities. If the TIS model can explain the response of diameter growth under different thinning intensities, even a single model could estimate tree growth in stands with various thinning intensities. rRMSE showed that there were differences in the TIS model under each thinning intensity (Table 7). Specifically, the rRMSE values for the model of intensive thinning and no thinning were higher than those of other thinning intensities, indicating low model accuracy. The model may not suit specific thinning intensities, because the TDM model had a similar tendency. On the other hand, both rRMSE and bias were smallest with moderate thinning; therefore, the model can roughly predict DBH growth for forests with moderate thinning intensities.
There were not significant differences in diameter growth among three thinned plots in our dataset, as shown in “2.1. Data Collection”. Therefore, our analysis could not examine in detail the behavior of TDM and TIS models under conditions in which there are large differences in diameter growth. On the other hand, there were significant differences in diameter growth among thinned plots and no thinning plot in our dataset. In addition, the TIS model accurately predicted diameter growth in thinned and no thinning plots, suggesting that the TIS model is effective for the multiple plots with largely different diameter growth. However, this suggestion is not based on actual data from the multiple plots with largely different diameter growth. Thus, the discussion here requires further analysis.

4.3. Application for Forest Management

Our results suggest that the TIS model can express diameter growth under various thinning intensities. In this study, the model predicts individual-tree diameter growth. When the model is applied to actual stands, the diameter size is predicted based on tree diameter sizes in the prior period. The long-term growth is repeatedly predicted. In addition, stand-level growth is predicted by accumulating individual tree-level growth. If the error is large, such as in intensive and no thinning, error will propagate. The more we predict long-term growth, the more the error will increase between actual stands and predicted values. Also, bias values were negative for all thinning intensities in the TIS model, indicating the predicted diameter growth was underestimated (Table 7). Underestimation of tree diameter is problematic for optimal operation planning, such as determining thinning intensity. Stands with low tree diameter usually require intensive thinning to allow residual trees to grow well; thus, underestimated tree diameter may lead to unnecessary thinning intensity. For those reasons, we need to improve models in which the error is constant with thinning intensity and bias values are not negative.

5. Conclusions

This study compared the model accuracy of TDM and TIS models with different thinning intensities. For both the TDM and TIS models, the selected competition index was the same and there was little difference in model accuracy. These results suggest that the TIS model could express diameter growth regardless of thinning intensity. This result contributes to yield prediction techniques for various thinning intensities. There are little differences in diameter growth among thinned plots in our study site. We need to consider model construction under the variety of thinning intensities. On the other hand, the TIS model had varying error under each thinning intensity, and the model accuracy was low in intensive thinning and no thinning. Also, bias values were negative; thus, predicted diameter size may be underestimated when the model is adapted to actual stands. These results suggest that our model has limitation regarding adaptation in unusual thinning intensities or other sites. Therefore, we need to improve the model to generate accurate predictions of the long-term growth of actual stands. We focused on individual tree growth in terms of tree competition in this study. Recently, some researchers have mentioned that site productivity is affected by climate change [42,43]. Thus, we need to consider factors such as temperature, precipitation, and carbon concentration when predicting long-term individual or forest stand growth.

Author Contributions

K.F., T.N. and F.K. conceived and designed the study. F.K. provided field measurements data. All authors contributed to data analysis. K.F and T.N. wrote the paper. All authors have read and agreed to the published of the manuscript.

Funding

This research was funded by research grant #20184 of the Forestry and Forest Products Research Institute.

Acknowledgments

We appreciate the academic editor and anonymous reviewers for their constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fujimori, T. Ecological and Silvicultural Strategies for Sustainable Forest Management; Elsevier: Amsterdam, The Netherlands, 2001. [Google Scholar]
  2. Davis, L.S.; Johnson, K.N.; Bettinger, P.; Howard, T. Forest Management: To Sustain Ecological, Economic, and Social Values; McGraw Hill: New York, NY, USA, 2001. [Google Scholar]
  3. Weiskittel, A.R.; Hann, D.W.; Kershaw, J.A.; Vanclay, J.K. Forest Growth and Yield Modeling; John Wiley and Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  4. Bose, A.K.; Weiskittel, A.; Kuehne, C.; Wagner, R.G.; Turnblom, E.; Burkhart, H.E. Does commercial thinning improve stand-level growth of the three most commercially important softwood forest types in North America? For. Ecol. Manag. 2018, 409, 683–693. [Google Scholar] [CrossRef]
  5. Lussetti, D.; Kuljus, K.; Ranneby, B.; Ilstedt, U.; Falck, J.; Karlsson, A. Using linear mixed models to evaluate stand level growth rates for dipterocarps and Macaranga species following two selective logging methods in Sabah, Borneo. For. Ecol. Manag. 2019, 437, 372–379. [Google Scholar] [CrossRef]
  6. Wykoff, W.R. A basal area increment model for individual conifers in the northern Rocky Mountains. For. Sci. 1990, 36, 1077–1104. [Google Scholar]
  7. Monserud, R.A.; Sterba, H. A basal area increment model for individual trees growing in even- and uneven-aged forest stands in Austria. For. Ecol. Manag. 1996, 80, 57–80. [Google Scholar] [CrossRef]
  8. Hökkä, H.; Alenius, V.; Penttilä, T. Individual-tree basal area growth models for Scots Pine, Pibescent Birch and Norway Spruce on Drained peatlands in Finland. Silva Fenn. 1997, 31, 161–178. [Google Scholar] [CrossRef]
  9. Contreras, M.A.; Affleck, D.; Chung, W. Evaluating tree competition indices as predictors of basal area increment in western Montana forests. For. Ecol. Manag. 2011, 262, 1939–1949. [Google Scholar] [CrossRef]
  10. Quiñonez-Barraza, G.; Zhao, D.; De Los Santos Posadas, H.M.; Corral-Rivas, J.J. Considering neighborhood effects improves individual dbh growth models for natural mixed-species forests in Mexico. Ann. For. Sci. 2018, 75, 78. [Google Scholar] [CrossRef] [Green Version]
  11. Sun, S.; Cao, Q.V.; Cao, T. Evaluation of distance-independent competition indices in predicting tree survival and diameter growth. Can. J. For. Res. 2019, 49, 440–446. [Google Scholar] [CrossRef]
  12. Wimberly, M.C.; Bare, B.B. Distance-dependent and distance-independent models of Douglas-fir and western hemlock basal area growth following silvicultural treatment. For. Ecol. Manag. 1996, 89, 1–11. [Google Scholar] [CrossRef]
  13. Miyamoto, A.; Amano, M. Difference of growth by the spatial distribution of individual trees and growth factor. Bull. FFPRI 2002, 1, 163–178. [Google Scholar]
  14. Tenzin, J.; Tenzin, K.; Hasenauer, H. Individual tree basal area increment models for broadleaved forests in Bhutan. Forestry 2017, 90, 367–380. [Google Scholar] [CrossRef]
  15. Takata, K.; Kobayashi, S. Relationships between diameter or basal-area growth and competition factors for Cryptomeria japonica D. DON stands after the last thinning. J. Jpn. For. Soc. 1983, 65, 113–118. [Google Scholar]
  16. Daniels, R.F.; Burkhart, H.E.; Clason, T.R. A comparison of competition measures for predicting growth of loblolly pine trees. Can. J. For. Res. 1986, 16, 1230–1237. [Google Scholar] [CrossRef]
  17. Biging, G.S.; Dobbertin, M. Evaluation of Competition Indices in Individual Tree Growth Models. For. Sci. 1995, 41, 360–377. [Google Scholar]
  18. Kahriman, A.; Şahin, A.; Sönmez, T.; Yavuz, M. A novel approach to selecting a competition index: The effect of competition on individual-tree diameter growth of Calabrian pine. Can. J. For. Res. 2018, 48, 1217–1226. [Google Scholar] [CrossRef]
  19. Kuehne, C.; Weiskittel, A.R.; Waskiewicz, J. Comparing performance of contrasting distance-independent and distance-dependent competition metrics in predicting individual tree diameter increment and survival within structurally-heterogeneous, mixed-species forests of Northeastern United States. For. Ecol. Manag. 2019, 433, 205–216. [Google Scholar] [CrossRef]
  20. Uzoh, F.C.C.; Oliver, W.W. Individual tree diameter increment model for managed even-aged stands of ponderosa pine throughout the western United States using a multilevel linear mixed effects model. For. Ecol. Manag. 2008, 256, 438–445. [Google Scholar] [CrossRef]
  21. Roberts, S.D.; Harrington, C.A. Individual tree growth response to variable-density thinning in coastal Pacific Northwest forests. For. Ecol. Manag. 2008, 255, 2771–2781. [Google Scholar] [CrossRef]
  22. Mäkinen, H.; Isomäki, A. Thinning intensity and growth of Scots pine stands in Finland. For. Ecol. Manag. 2004, 201, 311–325. [Google Scholar] [CrossRef]
  23. Mäkinen, H.; Isomäki, A. Thinning intensity and growth of Norway spruce stands in Finland. Forestry 2004, 77, 349–364. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, S.Y.; Chauret, G.; Swift, D.E.; Duchesne, I. Effects of precommercial thinning on tree growth and lumber quality in a jack pine stand in New Brunswick, Canada. Can. J. For. Res. 2006, 36, 945–952. [Google Scholar] [CrossRef]
  25. Del Río, M.; Calama, R.; Cañellas, I.; Roig, S.; Montero, G. Thinning intensity and growth response in SW-European Scots pine stands. Ann. For. Sci. 2008, 65, 308. [Google Scholar]
  26. Kitahara, F. Summary of Field Measurements Data of Japanese Sugi in Kudarukawa National Forest; Annual Report 56; Shikoku Research Center Foreatry and Forest Products Research Institute: Kochi, Japan, 2015; pp. 33–34. [Google Scholar]
  27. Tanaka, K. A Stochastic Model of Diameter Growth in an Even-Aged Pure Forest Stand. J. Jpn. For. Soc. 1986, 68, 226–236. [Google Scholar]
  28. Bates, D.; Mächler, M.; Bolker, B.M.; Walker, S.C. Fitting linear mixed-effects models using lme4. arXiv 2014, arXiv:1406.5823. [Google Scholar]
  29. R Core Team R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 8 November 2019).
  30. Sprugel, D.G. Correcting for Bias in Log-Transformed Allometric Equations. Ecology 1983, 64, 209–210. [Google Scholar] [CrossRef]
  31. Wykoff, W.R.; Crookston, N.L.; Stage, A.R. User’s Guide to the Stand Prognosis Model; General Technical Report INT-133; United States Department of Agriculture: Washington, DC, USA, 1982. [Google Scholar]
  32. Curtis, R.O.; Marshall, D.D. Technical Note: Why Quadratic Mean Diameter? West. J. Appl. For. 2000, 15, 137–139. [Google Scholar] [CrossRef] [Green Version]
  33. Wilson, F.G. Numerical expression of stocking in terms of height. J. For. 1946, 44, 758–761. [Google Scholar]
  34. Reineke, L.H. Perfecting a stand-density index for even aged forests. J. Agric. Res. 1933, 46, 627–638. [Google Scholar]
  35. Nishizono, T.; Tanaka, K. Does the relationship between quadratic mean diameter and stem density in old thinned and unthinned Cryptomeria japonica forests deviate from a power function? J. For. Plan. 2012, 18, 63–76. [Google Scholar]
  36. Kunisaki, T. Selecting a convenient index of stand density for the management of dense plantations of young Sugi (Cryptomeria japonica). Bull. Iwate Univ. 2013, 44, 1–18. [Google Scholar]
  37. Nakagawa, S.; Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 2013, 4, 133–142. [Google Scholar] [CrossRef]
  38. Tanaka, S.; Takahashi, T.; Nishizono, T.; Kitahara, F.; Saito, H.; Iehara, T.; Kodani, E.; Awaya, Y. Stand Volume Estimation Using the k-NN Technique Combined with Forest Inventory Data, Satellite Image Data and Additional Feature Variables. Remote Sens. 2014, 7, 378–394. [Google Scholar] [CrossRef] [Green Version]
  39. Kanazawa, Y.; Kiyono, Y.; Fujimori, T. Crown Development and Stem Growth in Relation to Stand Density in Even-Aged Pure Stands (II) Clear-length model of Cryptomeria japonica stands as a function of stand density and tree height. J. Jpn. For. Soc. 1985, 67, 391–397. [Google Scholar]
  40. Mäkinen, H.; Isomäki, A. Thinning intensity and long-term changes in increment and stem form of Norway spruce trees. For. Ecol. Manag. 2004, 201, 295–309. [Google Scholar] [CrossRef]
  41. Palahí, M.; Pukkala, T.; Miina, J.; Montero, G. Individual-tree growth and mortality models for Scots pine (Pinus sylvestris L.) in north-east Spain. Ann. For. Sci. 2003, 60, 1–10. [Google Scholar]
  42. Aussenac, G. Interactions between forest stands and microclimate: Ecophysiological aspects and consequences for silviculture. Ann. For. Sci. 2000, 57, 287–301. [Google Scholar] [CrossRef]
  43. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Seidl, R.; Delzon, S.; Corona, P.; Kolström, M.; et al. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
Figure 1. Location of the study site and permanent plots.
Figure 1. Location of the study site and permanent plots.
Forests 11 00344 g001
Figure 2. The relationship between predicted and observed diameter growth under different thinning intensities for the (1) TDM model and (2) TIS model. Note: Predicted diameter growth was calculated by the best model of TDM model according to Table 4, and the best model of TIS model according to Table 6.
Figure 2. The relationship between predicted and observed diameter growth under different thinning intensities for the (1) TDM model and (2) TIS model. Note: Predicted diameter growth was calculated by the best model of TDM model according to Table 4, and the best model of TIS model according to Table 6.
Forests 11 00344 g002
Table 1. Summary of thinning intensities in study plots.
Table 1. Summary of thinning intensities in study plots.
Thinning Rate at 25 Years Old (%) Thinning Rate at 37 Years Old (%)
VolumeNumber of treesVolumeNumber of trees
(1) Intensive thinning26412938
(2) Moderate thinning22382334
(3) Light thinning20322334
(4) No thinning----
Table 2. Summary of measurement data in study site.
Table 2. Summary of measurement data in study site.
Age (Year)
14192329344045515661
(1) Intensive thinning (0.116 ha)
Mean height (m)11.313.916.720.723.226.028.830.432.033.3
Mean DBH (cm)14.517.218.825.527.632.334.437.338.940.8
Stand density (trees/ha)25092509245712411241879724707707707
(2) Moderate thinning (0.123 ha)
Mean height (m)10.412.915.420.922.725.828.530.331.633.6
Mean DBH (cm)13.917.219.125.127.332.134.736.939.041.0
Stand density (trees/ha)22362236215412281228919805805780780
(3) Light thinning (0.106 ha)
Mean height (m)11.314.216.921.923.126.329.130.431.833.4
Mean DBH (cm)15.618.920.726.228.033.135.637.538.640.6
Stand density (trees/ha)21892189213213491349887887887877868
(4) No thinning (0.113 ha)
Mean height (m)10.213.115.618.920.824.425.726.928.329.5
Mean DBH (cm)14.116.918.821.023.425.826.728.430.431.8
Stand density (trees/ha)2469246924072248200917701770168114871416
Note: DBH, Diameter at breast height.
Table 3. Evaluation of thinning intensity-dependent multiple (TDM) models using each distance-independent competition index.
Table 3. Evaluation of thinning intensity-dependent multiple (TDM) models using each distance-independent competition index.
Marginal R2Conditional R2AICΔAIC
(1) Intensive thinning
Null0.290.532611126
BAL0.350.63252338
DR0.330.56256983
BR0.310.562602117
CDF0.290.532618133
Sr0.360.6224850
SDI0.320.59254661
(2) Moderate thinning
Null0.300.592596247
BAL0.370.672475126
DR0.340.612553204
BR0.330.612585236
CDF0.300.592602253
Sr0.420.7023490
SDI0.360.652463114
(3) Light thinning
Null0.280.572481230
BAL0.380.70231766
DR0.360.622396144
BR0.330.622441189
CDF0.280.572488237
Sr0.410.6922510
SDI0.340.632371119
(4) No thinning
Null0.350.583957363
BAL0.430.683749155
DR0.420.653824230
BR0.380.623934340
CDF0.350.583964371
Sr0.490.6835940
SDI0.480.68361420
Note: BAL, Basal area of trees larger than the subject tree; DR, Diameter ratio; BR, Basal area ratio; CDF, Cumulative distribution function; Sr, Relative spacing index; SDI, Site density index.
Table 4. Summary of parameter estimation using the relative spacing index (Sr) for the TDM model.
Table 4. Summary of parameter estimation using the relative spacing index (Sr) for the TDM model.
ParametersDescriptionEstimateSEp-Value
(1) Intensive thinning
α0intercept−5.5710.340<0.001
α1DBH0.2000.014<0.001
α2DBH2−0.0020.000<0.001
α3Age−0.0560.003<0.001
α4Sr0.1610.013<0.001
(2) Moderate thinning
α0intercept−6.2600.346<0.001
α1DBH0.1740.012<0.001
α2DBH2−0.0010.000<0.001
α3Age−0.0430.003<0.001
α4Sr0.2000.012<0.001
(3) Light thinning
α0intercept−6.6500.368<0.001
α1DBH0.1990.014<0.001
α2DBH2−0.0020.000<0.001
α3Age−0.0440.003<0.001
α4Sr0.2090.013<0.001
(4) No thinning
α0intercept−6.8010.307<0.001
α1DBH0.2130.014<0.001
α2DBH2−0.0020.000<0.001
α3Age−0.0300.003<0.001
α4Sr0.1950.010<0.001
Table 5. Evaluation of the thinning intensity-independent single (TIS) model using each distance-independent competition index.
Table 5. Evaluation of the thinning intensity-independent single (TIS) model using each distance-independent competition index.
Marginal R2Conditional R2AICΔAIC
Null0.330.5811578803
BAL0.400.6711080305
DR0.360.6211393618
BR0.350.6111504729
CDF0.330.5811587811
Sr0.420.67107750
SDI0.370.6311205430
Table 6. Summary of parameter estimation using Sr for TIS model.
Table 6. Summary of parameter estimation using Sr for TIS model.
ParametersDescriptionEstimateSEp-Value
α0intercept−5.5870.156<0.001
α1DBH0.1770.007<0.001
α2DBH2−0.0020.000<0.001
α3Age−0.0400.002<0.001
α4Sr0.1630.005<0.001
Table 7. Evaluation of TDM and TIS models using Sr under different thinning intensities.
Table 7. Evaluation of TDM and TIS models using Sr under different thinning intensities.
RMSE (cm)rRMSE (%)Bias (cm)Variance (cm)
TDMTISTDMTISTDMTISTDMTIS
(1) Intensive thinning0.2490.25054.254.4−0.092−0.1040.0540.052
(2) Moderate thinning0.2300.22047.845.8−0.074−0.0680.0470.044
(3) Light thinning0.2370.22851.149.1−0.078−0.0710.0500.047
(4) No thinning0.2030.19960.659.4−0.049−0.0810.0390.033

Share and Cite

MDPI and ACS Style

Fukumoto, K.; Nishizono, T.; Kitahara, F.; Hosoda, K. Evaluation of Individual Distance-Independent Diameter Growth Models for Japanese Cedar (Cryptomeria japonica) Trees under Multiple Thinning Treatments. Forests 2020, 11, 344. https://doi.org/10.3390/f11030344

AMA Style

Fukumoto K, Nishizono T, Kitahara F, Hosoda K. Evaluation of Individual Distance-Independent Diameter Growth Models for Japanese Cedar (Cryptomeria japonica) Trees under Multiple Thinning Treatments. Forests. 2020; 11(3):344. https://doi.org/10.3390/f11030344

Chicago/Turabian Style

Fukumoto, Keiko, Tomohiro Nishizono, Fumiaki Kitahara, and Kazuo Hosoda. 2020. "Evaluation of Individual Distance-Independent Diameter Growth Models for Japanese Cedar (Cryptomeria japonica) Trees under Multiple Thinning Treatments" Forests 11, no. 3: 344. https://doi.org/10.3390/f11030344

APA Style

Fukumoto, K., Nishizono, T., Kitahara, F., & Hosoda, K. (2020). Evaluation of Individual Distance-Independent Diameter Growth Models for Japanese Cedar (Cryptomeria japonica) Trees under Multiple Thinning Treatments. Forests, 11(3), 344. https://doi.org/10.3390/f11030344

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop