Stocking Diagrams for Silvicultural Guideline in Korean Pines and Japanese Larch

: Appropriate management of stand density is necessary to avoid wasted growing space and overcrowding-induced self-thinning and therefore to optimize proﬁtability. We developed a Gingrich-style stocking chart for Korean red pine ( Pinus densiﬂora ), Korean white pine ( P. koraiensis ), and Japanese larch ( Larix kaempferi ) in Korea. Datasets for even-aged stands were categorized into two censored datasets via relative density based on species-speciﬁc stand density indices to assign adequate plots for tree allometry. Censored plot data for maximum density on full stocking were used to develop A-level stocking based on mean individual tree area. In censored plot data for minimum density on full stocking, individual trees represented open-grown trees, and a crown competition factor of 100 was proposed as B-level stocking. Based on parameters estimated from allometry, stocking diagrams comprising quadratic mean diameter, number of trees, and basal area were correctly expressed. A-level stocking at the same quadratic mean diameter revealed that Korean white pine had the most trees and largest basal area, while Japanese larch had the fewest trees and smallest basal area. In contrast, B-level stocking disclosed that Japanese larch had the most trees and basal area, whereas Korean white pine had the fewest trees and smallest basal area. The stocking diagrams suggest that silvicultural treatments for these species should be species-speciﬁc.


Introduction
Stand size-density relationships influence tree crowding and competition for resources such as light, water, and nutrients [1]. In forest management, one objective is to produce utility pole trees as rapidly as possible without wasting limited growing space. Hence, it is necessary to establish the growing space required for each target species. Poorly managed stand density may result in low profitability. The allocation of excessive growing space is a wasteful forest management practice and results in coarse stems and unwanted branches of poor quality. In contrast, competition among trees is extremely high in very dense stands, and the trees therein are prone to self-induced mortality [2].
Forest managers attempt to explain tree size-density relationships in order to control size, growth, and mortality [3,4]. Available tree-growing space is influenced by the degree of crowding within a stand and is usually defined by the number of trees, basal area, and volume per unit area. However, these factors are absolute and quantitative. Further, no single factor adequately defines the tree growing space. Two different stands with the same basal area or number of trees do not necessarily have the same amount of growing space for individual trees when their mean tree diameters differ [5].
The study area comprises pure stands of three major commercial conifers, namely Korean red pine, Korean white pine, and Japanese larch. It was located in Gangwon and North Gyeongsang Provinces of Korea. Permanent monitoring plots were established to develop growth and yield models for thinned stands with a squared plot design ranging in size from 0.04 ha (20 m × 20 m) to 0.09 ha (30 m × 30 m). Each experimental site consisted of three plots differing in size and thinning rate of plantation. The experimental sites varied in terms of tree number and age class. Detailed information is provided in Lee and Choi [27]. The plots have been measured at 3-year intervals since they were established in 2012-2014.
The plots were measured using the same inventory method for all three target species. The diameter at breast height (DBH) was measured at 1.2 m height using a diameter tape and was recorded in centimeters to the first decimal place. All trees ≥6 cm within each plot were measured except for young plantations, wherein the saplings were <6 cm in DBH and/or <1.2 m in height. Tree mortality was also determined and recorded. Crown width (CW) was measured in meters with a measuring tape to the first decimal place using two crown diameters from the second inventory (3 years after plot establishment) perpendicular to the slope direction. Crown data were collected from the longest branch under a normal crown shape in each of the four directions. The branches were approximately repositioned in pre-felling orientation using the overall leaf angles as a guide, according to the method of Choi et al. [28,29]. Unlike DBH, CW was determined for only 1/3 of the trees per plot. Trees were selected for CW measurements by sorting the DBH and by a systematic sampling method. When the number of trees available for CW measurement was <10 because of low density, then the minimum number of trees samples was set to 10. There were 252 plots for Korean red pine, 342 plots for Korean white pine, and 405 plots for Japanese larch. The numbers of Korean red pine, Korean white pine, and Japanese larch trees measured for DBH and CW were 1541, 2014, and 2921, respectively. After collecting and aggregating all these recurring inventory data, they were filtered for A-level and B-level development. The variables tested in this study were calculated on the basis of unit area or average. The number of trees per plot was converted to number of trees per hectare (TPH). The basal area per plot was converted to the basal area per hectare (BA). D q was calculated using the basal area as follows: where D q is the quadratic mean diameter (cm), BA is the mean basal area per tree (m 2 ), and π is a mathematical constant of~3.14159. D q simplified the interpretation of these relationships because yield tables and stand projections are relatively easy to follow [30]. Individual tree CWs were arithmetically calculated by averaging two perpendicular CWs per tree.

Maximum Density in Fully Stocked Stands (A-Level Stocking)
Overstocked stands undergo self-thinning to control their density. They revert to A-level stocking percentage wherein the trees are allocated just enough resources required for minimum growth [9]. A-level stocking is represented by an allometric relationship developed with regression models on normal stand data. Normal stands are undisturbed and even-aged. They have uniform tree area and no overstory canopy space. They nearly attain the maximum basal area and volume for each age and site productivity level [31]. However, the normal stand determinations made by the methods of Gingrich [6] and Chisman and Schumacher [32] for A-level stocking were dubious and difficult to Forests 2020, 11, 833 4 of 17 achieve because the stands evaluated for these determinations had exceptionally high or low mortality at each developmental stage [33]. The relative density or ratio of actual: maximum stand density is often used to evaluate the observed stand density and select normal stands. However, the selection criteria varied among studies [34].
Long [35] reported that the conifers in stands with >60% relative density showed the onset of self-thinning. This phenomenon was reported for longleaf pine by Shaw and Long [36]. Solomon and Zhang [37] established that mixed-softwood stands with relative density >0.7 were normal stands and manifested self-thinning. This criterion was recently applied to longleaf pine by Kara et al. [5]. Here, relative density >0.7 was adopted to ensure meticulous selection of normal stands [38]. The maximum stand density line used in the relative density calculations was adapted from Lee and Choi [27]. They used the same stand density indices for the same three species as those tested in the present study ( Figure 1a). Based on this criterion, 98, 118, and 128 plots were selected for Korean red pine, Korean white pine, and Japanese larch, respectively. The filtered data were used to develop an allometric equation for A-level stocking ( Table 1).
Forests 2020, 11, x FOR PEER REVIEW 4 of 17 stand density is often used to evaluate the observed stand density and select normal stands. However, the selection criteria varied among studies [34]. Long [35] reported that the conifers in stands with >60% relative density showed the onset of self-thinning. This phenomenon was reported for longleaf pine by Shaw and Long [36]. Solomon and Zhang [37] established that mixed-softwood stands with relative density >0.7 were normal stands and manifested self-thinning. This criterion was recently applied to longleaf pine by Kara et al. [5]. Here, relative density >0.7 was adopted to ensure meticulous selection of normal stands [38]. The maximum stand density line used in the relative density calculations was adapted from Lee and Choi [27]. They used the same stand density indices for the same three species as those tested in the present study ( Figure 1a). Based on this criterion, 98, 118, and 128 plots were selected for Korean red pine, Korean white pine, and Japanese larch, respectively. The filtered data were used to develop an allometric equation for A-level stocking ( Table 1). level stocking for Korean red pine, Korean white pine, and Japanese larch with maximum stand density lines by species. Sample data were constrained by plots with relative density >0.7 for A-level and relative density <0.5 for B-level. Calculations were made using maximum stand density lines derived from stand density indices in Lee and Choi [27].
is the quadratic mean diameter. white pine, and Japanese larch with maximum stand density lines by species. Sample data were constrained by plots with relative density > 0.7 for A-level and relative density < 0.5 for B-level. Calculations were made using maximum stand density lines derived from stand density indices in Lee and Choi [27]. D q is the quadratic mean diameter. B-level stocking is represented as the minimum density for full stocking. The entire mean individual tree area is occupied, and no available growing space is wasted. The crown growth of the trees under B-level stocking is not interrupted by competition among trees within the stand. Open-grown trees are commonly used to develop CW-DBH allometry for this level [39]. When these data cannot be measured, however, dominant forest-grown trees are used to evaluate CW-DBH allometry for B-level stocking [40,41]. In another study, allometric equations were also developed by using data for thinned stands [11]. The slope coefficient of the crown area-basal area relationship in forest-grown trees did not significantly differ from that for open-grown trees [2,42,43].
In this study, CW-DBH allometry by species was examined by using < 0.5 relative density. The data were filtered according to the number of tree samples, statistical summary and distribution, similarity of the allometry pattern with those of previous studies, and model fit statistics [38]. The numbers of plots selected for B-level stocking were 67 for Korean red pine, 98 for Korean white pine, and 129 for Japanese larch (Figure 1b). For the CW-DBH measurements within the filtered plots at low relative density, there were 543 trees in 50 plots, 755 trees in 68 plots, and 1128 trees in 105 plots, respectively. Summary statistics of the datasets for A-level and B-level stocking development are listed in Table 1.

Minimum Tree Area in Fully Stocked Stands (for A-Level Stocking)
Tree area was traditionally used to determine A-level stocking and stocking percentage. This technique was originally developed by Chisman and Schumacher [32]. A certain area is allocated to the trees growing inside the area. The sum of all individual tree areas is converted to a certain area, and the following allometric equation is proposed for numeric interpretation: where TA is the tree area (such as 10, 000 m 2 ); n is the number of trees per unit area; d i is i th individual tree DBH; and β 0 , β 1 , and β 2 are coefficients. This equation was applied by Gingrich (1967) to determine A-level stocking via regression modeling. However, the response variable was a uniformized constant set to a unit area such as 10,000 m 2 [7]. Hence, both terms in Equation (2) were divided by the number of trees (n), and both the response and explanatory variables were modified as follows: where all variables and parameters indicated were previously defined. Using Equation (3), a stocking diagram was developed for an Eastern cottonwood-silver maple-American sycamore forest [7].
To express A-level stocking on the stocking diagram, the allometry was not limited to the type of tree area equation. A log-log relation between TPH and D q was recently applied for longleaf pine [5].
In the present study, the tree area method was selected for A-level stocking as it directly describes the spatial concept within an area according to the crown area for B-level stocking [38]. The former equation was modified to solve for multicollinearity between independent variables and improve model fitting: where TA is the mean individual tree area (m 2 ) equal to TA divided by n. All other variables and coefficients are the same as those previously defined. The allometric equation indicates that TA is proportional to the square of D q . In the present study, this equation was applied without any statistical processing or validation issues. Using the estimated coefficients, TPH was explained as follows via D q determined from TA: where % stocking is the stocking level expressed as a percentage (%). All other variables and coefficients are the same as those previously defined. Equation (5) was applied to the stocking diagram to describe the relationship between the average tree diameter and the TPH as stocking percentage.

Crown Area on Open-Grown Sites (for B-Level Stocking)
B-level stocking is the minimum TPH estimated from a fully occupied growing space and defined as crown competition factor = 100 [39]. In the plots, the individual trees with these characteristics were used in the CW-DBH allometry model (Table 1). CW-DBH allometry is described by linear, power, logarithmic, and other equations [3,44]. Here, a linear equation was used because it effectively describes open-grown trees [2,5,7]. The following allometric equation was applied to individual open-grown trees for B-level stocking: where all variables and coefficients are the same as those previously defined. The equation was successfully fitted without examining any other factors. It was based on the premise that the relationship between DBH and CW for open-grown trees is not influenced by site productivity or stand age [32,44]. Based on the estimated parameters, the allometry of DBH and CW resulted in a relationship between average tree diameter and number of trees at the onset of full stocking: where TPH O is the number of trees per hectare at the onset of full stocking that reaches the limit of the open-growing condition (trees·ha −1 ). All other variables and coefficients are the same as those previously defined. TPH O was computed from the crown area and represents B-level stocking consisting of open-grown trees. It can also be defined as crown competition factor = 100. Equation (7) was applied to determine the B-level stocking percentage in the stocking diagram.

Statistical Analysis
The model parameters were estimated by linear regression for Equation (4) (square root-transformed), and for Equation (6) (original arithmetic form). Each experimental site had a nested design because there were three plots per site, and they were measured thrice after establishment. Linear or nonlinear mixed-effects models were developed for the forest trees and stands to account for multivariate multilevel random plot effects and covariances among repeated measurements [45,46]. By regarding these factors as random effects, a linear mixed-effect model was formulated as follows: where TA is the square root of the mean individual tree area, b 0 and b 1 are fixed-effect parameters, u ei , u ei.pj , u ei.pj.tk are random-effect parameters, and ε ei.pj.ml and ε ei.pj.tk.ml are the measurement errors for each equation that were not explained by fixed-effect or random-effect parameters. In Equation (8) for the plot-level analysis, u ei was incorporated because of the nested plot design and u ei.pj was incorporated because of the longitudinal data structure. In Equation (9) for the tree-level analysis, u ei and u ei.pj accounted for the spatial correlations among the individual trees within the same area, while u ei.pj.tk explained the longitudinal data structure [47]. All parameters were estimated with the PROC MIXED module in SAS v. 9.4 [48]. Parameter signs and values were checked, and the significance level was determined from the p-value to validate the model. The coefficient of determination (R 2 ) was calculated to evaluate the percentage variation explained by the linear mixed-effects model; R 2 = 1 − (sum of squared errors/Total sum of squares). The root mean squared error (RMSE) and standard deviation (SD) of the residuals were calculated to evaluate their spread. Scatterplots of the residuals (differences between the measured and predicted values) were made to identify and discriminate undesirable and abnormal distributions such as biased trends and heteroscedasticity.

Parameter Estimates and Tree Allometric Descriptions
The parameters of Equations (8) and (9) fit well with the range of suitable allometric relationships. All fixed-effect slope parameters were statistically significant (p < 0.01). In Equation (8), the slope and intercept parameters were highest for Japanese larch. As the slope parameter was larger for Japanese larch, TA increased more with D q , and TPH decreased more than other species in A-level stocking in Equation (5). The relatively large TA for Japanese larch caused its TPH to be lower than that of Korean red pine or Korean white pine at the same stocking percent. The variance of the residuals by species not explained by the fixed and random effects were in the range of 0.013-0.022, and there were no marked differences among them. R 2 > 0.968 and RMSE < 0.110 for all three species.
The slope parameter in the CW-DBH model was largest for Korean white pine and smallest for Japanese larch (Equation (9)). At the same DBH, then, the crown area was larger for Korean white pine and smaller for Japanese larch. For B-level stocking, there were relatively fewer TPH O for Korean white pine and relatively more TPH O for Japanese larch, as calculated using Equation (7) with the fixed-effect parameters in Table 2. The relationship between DBH and CW influenced both the crown area and TPH O for B-level stocking. The observed random effects in the CW-DBH model varied with species. The residual variances were 0.840 for Korean red pine, 0.603 for Korean white pine, and 0.789 for Japanese larch. R 2 was highest and RMSE were lowest for Korean white pine. All RMSE values were lower than the SD for each species with a meaningful model fitting (Table 1). In the residual plots of both models (Equations (8) and (9)), the residuals were symmetrically distributed, and undesirable patterns such as heteroscedasticity were not detected within the random distribution ( Figure 2). varied with species. The residual variances were 0.840 for Korean red pine, 0.603 for Korean white pine, and 0.789 for Japanese larch. was highest and RMSE were lowest for Korean white pine. All RMSE values were lower than the SD for each species with a meaningful model fitting (Table 1). In the residual plots of both models (Equations (8) and (9)), the residuals were symmetrically distributed, and undesirable patterns such as heteroscedasticity were not detected within the random distribution ( Figure 2).  (8) and (9) in Table 2.
The scatterplots of and used for modeling are displayed with regression lines developed for A-level stocking (Figure 3a). The original distribution for the data used in A-level stocking corresponded well to the squared values of the predicted values by regression (Equations (4) and (8) (8) and (9) in Table 2.
The scatterplots of D q and TA used for modeling are displayed with regression lines developed for A-level stocking (Figure 3a). The original distribution for the data used in A-level stocking corresponded well to the squared values of the predicted values by regression (Equations (4) and (8)). TA were 3.5 m 2 for Korean red pine, 3.5 m 2 for Korean white pine, and 4.1 m 2 for Japanese larch at 10 cm D q . At 40 cm D q , they were 22.7 m 2 , 22.3 m 2 , and 25.3 m 2 , respectively. The calculated TA resulted in the maximum TPH for the fully stocked stand (Equation (5)). The order of TPH by species for the A-level stocking was the opposite of the order by species for TA. This finding was consistent with the carrying capacity sequence by species for the stand density index [27]. This TPH discrepancy amongst species is explained by the differences in slope parameter estimated from TA-D q for each species (Table 2).    Table 2.
The TA used in B-level stocking was larger than that used for the A-level stocking at the same D q because the former represented open-grown trees (Figure 3b). The TA were 7.6-8.7 m 2 at 10 cm D q and 35.1-44.2 m 2 at 40 cm D q . These TA were larger than those used in A-level stocking. TA for B-level stocking were exponentially higher than that for A-level stocking. A similar pattern was reported for a previous study, but the absolute sizes differed [7].
We used the individual tree data with CW-DBH measurements for the plots shown in Figure 3b to display species-specific tree allometry between DBH and CW with scatterplots and regression lines for B-level stocking (Figure 4). There was a positive linear correlation between DBH and CW ( Table 2). The linearity was strongest for Korean white pine and weakest for Japanese larch. Neither the slopes nor the intercept parameters noticeably differed among species. The predicted CW were 2.8 for Korean red pine, 2.9 m for Korean white pine, and 2.8 m for Japanese larch in 10 cm DBH and 7.3 m, 7.6 m, and 7.0 m, respectively, in 40 cm DBH (Figure 4d). CW increased with DBH and influenced the TPH O developed for B-level stocking on the stocking diagram structure.    Table 2.

Illustration and Comparison of Stocking Diagram Applications
The stocking diagram was plotted from the allometric relationship and its estimated parameters in the fitting dataset. The average tree diameter and TPH for the A-level and B-level stockings were calculated using Equations (5) and (7) and the fixed-effect parameters of Table 2. The relationships between BA and the other variables were calculated using Equation (1). The final stocking diagram by species consisted of D q , TPH, and BA with A-level and B-level stocking ( Figure 5). The stocking percentage was calculated from the ratio of the A-level stocking percentage to the relative percentage range of 20-120% (Equation (5)). The average tree diameter was in the range of 12-32 cm so that it could encompass the Korean small and large tree diameter criteria of 16 cm and 30 cm, respectively [49].  The stocking diagram was categorized as overstocked and fully stocked stand conditions based on the A-level stocking percentage and as fully stocked and understocked stand conditions based on the B-level stocking percentage and the allometric equations ( Figure 5). The stocking diagram for Japanese larch differed from that for Korean red pine and Korean white pine in terms of general shape. The two pine species had different absolute densities (TPH and BA) at the same D q . In terms of 100% stocking of 32 cm D q , the quantitative absolute densities for Korean white pine were 636 trees·ha −1 TPH and 51.1 m 2 ·ha −1 BA, and for Japanese larch they were 558 trees·ha −1 TPH and 44.9 m 2 ·ha −1 BA. Japanese larch had a lower quantitative absolute density than either Korean red pine or Korean white pine. In 12 cm D q , the 100% stockings of Japanese larch were 2013 trees·ha −1 TPH and 22.8 m 2 ·ha −1 BA. This carrying capacity was lower than that for the other two species. This discrepancy originated from the differences among the species in terms of their slope parameters for the TA allometry model ( Table 2).
The A-level stocking and stocking percentage were affected by the TA and, by extension, the relationships among the variables. Similarly, the B-level stocking percentage was influenced by the CW-DBH allometry model. Korean red pine and Korean white pine had large CW and, therefore, large crown areas (Figure 4). The B-level stockings (crown competition factor = 100) in 32 cm D q were 343 trees·ha −1 and 313 trees·ha −1 TPH O , respectively (Figure 5a,b). In contrast, the Japanese larch had a relatively smaller CW and, by extension, a smaller crown area. Its B-level stocking in 32 cm D q was 370 trees·ha −1 TPH O (Figure 5c). This characteristic slightly differed across the D q range in the stocking diagram because each species had a different intercept and slope parameters ( Table 2). In 12 cm D q , the TPH O of B-level stocking were 1321 trees·ha −1 for Korean red pine, 1219 trees·ha −1 for Korean white pine, and 1311 trees·ha −1 for Japanese larch. The relationship between TPH O and D q affected the differences in BA among species in the B-level stocking. The ranges of BA in 12-32 cm D q were 14.9-27.5 m 2 ·ha −1 for Korean red pine, 13.8-25.2 m 2 ·ha −1 for Korean white pine, and 14.8-29.7 m 2 ·ha −1 for Japanese larch. The B-level stocking % across D q were in the ranges of 55-57%, 49-52%, and 65-66% for Korean red pine, Korean white pine, and Japanese larch, respectively ( Figure 5). The stocking diagrams for the three species were generally well represented and had similar shapes. Nevertheless, the stocking levels based on D q , TPH, and BA clearly differed among the three species.

Suitability of Gingrich-Style Stocking Chart
We compared several Gingrich-style stocking charts developed in North America against that which was prepared in the present study. We validated the form and suitability of the stocking diagram designed here because this type of chart has not yet been implemented for the same species in Korea. Here, the overall shapes of D q , TPH, and BA resembled those reported in Gingrich [6]. Our absolute TPH for A-level stocking was 200 trees·ha −1 lower than that of the previous study. Hence, a relatively lower BA was shown in the present study. However, this discrepancy was attributed to the presence of hardwoods in the upland regions of the United States. Although these regions might substantially differ from those of the present study in terms of species, geography, climate, and stand origin and structure, our stocking diagrams generally presented a similar stocking level to that of other studies with reference to general plant ecology or silvicultural practices. The A-level stocking reported in Larsen et al. [7] for midwestern U.S. bottomland forests was similar to that reported in the current study, especially for Japanese larch and especially in the quadratic mean diameter range 8-12 inches D q (~20-30 cm D q ). For the A-level stocking in this D q range for Japanese larch, there were 1089 trees·ha −1 TPH and 34.2 m 2 ·ha −1 BA at 20 cm D q , and 615 trees·ha −1 TPH and 43.5 m 2 ·ha −1 BA at 30 cm D q (Figure 5c).
The highest comparability was determined for a recently plotted stocking diagram for longleaf pine (Pinus palustris Mill.) forest [5]. The stocking charts for the Korean pines were analogous to the stocking chart for longleaf pine after direct comparison of the metric unit charts of the previous study (Figure 5a,b). The stocking diagram of Korean red pine was highly similar to that of longleaf pine (693 trees·ha −1 TPH and 49.0 m 2 ·ha −1 BA at 30 cm D q on A-level stocking in the current study). As there was high comparability even at the lower D q (≤20 cm) despite the differences in species and region, the stocking diagrams for the pine species in the present study can provide feasible stocking ranges for practical application.
On the other hand, the B-level stocking of Larsen et al. [7] was as low as 30%, while the B-level stocking calculated in the present study was 55-57% for Korean red pine, 49-52% for Korean white pine, and 65-66% for Japanese larch. The B-level stocking of longleaf pine was <40% and markedly differed from that of the present study [5]. However, it was analogous to the result of Gingrich [6], whose B-level stocking was~60%. In this way, B-level stocking either resembled or differed from the stocking of the current work depending on the species and stand condition. The allometric equation for CW-DBH of young Japanese larch stands provided similar estimates for relatively smaller DBH [50]. Hence, B-level stocking was deemed appropriate for the present study considering that it serves as the marginal density level after thinning. Thus, absolute quantitative comparisons are undesirable in situations where the regions, species, data characteristics, and management objectives differ [7].
The stand density management diagram enables accurate density management via allometry and is suitable for practical field applications [12]. It allows forest managers to compute the remaining variable based on the other two variables and to plan silvicultural operations such as thinning by interpolating the amount of thinned BA from the y-axis. Moreover, combinations of any two of the three variables (D q , TPH, and BA) facilitates designation of the relative stocking as a major feature.

Practicability and Applicability of the Stand Density Management Diagram
In the attempt to plot a stand density management diagram for Korea, an allometric line of Korean red pine in the Gangwon area was proposed by using the maximum stand density line [18]. Relative to the present study, similar stocking and maximum TPH were suggested at 30 cm D q ; however, compared to our study, the maximum TPH was lower as the D q decreased at <30 cm D q . This approach was also tested on Japanese larch stand density management [19]. The density chart for that study resembled that of the present study. It indicated 100% stocking of 22 cm D q in 1000 trees·ha −1 (TPH). In general, the density charts developed in the past [18,19] were similar to those plotted here. Nevertheless, the stand structures of the previous studies markedly differed from those of the present time in terms of age class, mean diameter, stand volume, and other variables. Stand density management diagrams calculated for our study reflect current stand condition characteristics. The diagrams used in the present study were evaluated to be feasible and practicable for field application. In contrast, those used in the previous study were plotted on a log scale and lacked B-level stocking; therefore, they can be difficult to interpret.
The maximum stand density lines used in previous stand density management diagrams [22][23][24] developed using the methodology proposed by Ando (1968) were slightly lower than the A-level stocking for the species in the present study. This method was recently applied to propose a Korean red pine density guideline based on National Forest Inventory data [25,26]. TPH corresponded to a relative yield index of unity (1) and referred to maximum stand density. The TPH in the previous study resembled that of the A-level stocking in the present study. The adequate relative yield index of 0.84 recommended for density management by Park et al. [26] corresponded to 550 trees·ha −1 (TPH) at 30 cm DBH and was contained within the fully stocked density area of the present study. Stocking diagrams for density management are nonetheless required because the previous study differed from the current work in terms of National Forest Inventory data sources and the stocking diagram approach used. Moreover, considering that the seed tree method is used for the natural regeneration of Korean red pine in Korea, both the B-level stocking percentage and the stocking level may be implemented under this guideline for silvicultural practices or recruitment [5].
Empirical yield tables for each species were provided by the National Institute of Forest Science (NIFoS) and were consulted to validate stocking diagrams developed in our study. The DBH, TPH, and BA in the empirical yield table provided by NIFoS all fell within the range of the fully stocked densities by species in the stocking diagram of the present study. Hence, the stocking diagrams effectively represented the stocking level characteristics of the empirical stands [51,52]. The stocking diagrams for all three species were effective and easily applicable as they corresponded to those of the previous studies. In earlier studies on longleaf pine and upland hardwoods, it was reported that 70-80% stocking caused density-induced self-thinning [36,53]. The stocking percentage was~70% for all three species after the empirical yield table data were plotted [51,52]. In the present study, <70% stocking based on the consideration of the stocking level of density-induced mortality and post-thinning stand growth was recommended following the silvicultural treatments such as pre-commercial thinning or commercial thinning.
The difference in TPH by D q on the x-axis increased with decreasing D q for all three species. Thus, variation in maximum TPH on the stocking level increased with decreasing D q . In other words, the carrying capacity at higher density increases with decreasing tree diameter [5]. In cases of equal stand density such as BA, large trees are less prone to crowding than small trees [2]. Thus, when stocking is not considered, silvicultural treatments such as thinning may waste growing space in a stand. When stocking meets the stand management objectives, it may effectively control density in the field because it accounts for the ratio of relative crowding (based on the observed growing space) to the maximum density criterion line. In terms of percentage stocking, although the stocking diagrams of Korean red pine and Korean white pine seemed to be similar to each other, they must be provided separately for each species. This is not only because of the obvious difference in B-level stocking, but also based on the characteristics of each species such as tree physiology, shade tolerance, and afforestation method [22,24,49,51,52]. The models developed in the present study were validated against other models reported for species in Korea. Cross-checking of the datasets from the other reports and the trajectories of stand density should be verified in future work. However, based on our models and comparisons against models for species inside and outside of Korea, the recommended stocking percentage is to be higher than B-level stocking by species and lower than 70% stocking in general after thinning. For this reason, the type of stocking chart developed in the present study was evaluated for use as a practical stand density management diagram.

Conclusions
To develop A-level and B-level stocking with a Gingrich-style stocking chart, we used data from permanent plots for three species and applied the relative density based on the maximum stand density line. A-level stocking was calculated from the allometric relationship of TA-D q . B-level stocking was calculated from the allometric relationship of CW-DBH. According to the estimated parameters, TA was largest for Japanese larch and smallest for Korean white pine. Hence, the descending order of maximum allowable fully stocked density (TPH and BA) was Korean white pine, Korean red pine, and Japanese larch. In contrast, the crown area determined from CW-DBH was largest for Korean white pine and smallest for Japanese larch. Therefore, the B-level stocking, generated by a crown competition factor of 100, resulted in the lowest TPH O for Korean white pine and the highest TPH O for Japanese larch.
In our stocking diagrams, the ranges of overstocked, fully stocked, and understocked densities differed among the species. TPH and BA for the A-level stocking were Korean white pine, Korean red pine, and Japanese larch in the order across all D q . In view of the references and other guidelines for Korea and the results of the present study, the recommended post-thinning stocking percentage is <70%. B-level stocking resulted in the highest TPH O and BA for Japanese larch followed by Korean red pine and Korean white pine. The B-level stocking was also assessed as a reference line for the creation of open-grown canopy space because Korean red pine stands may be established in Korea by the seed tree method. Based on the stocking diagrams developed in our study, forest managers are encouraged to manage the stand density by generally keeping the stocking percentage less than 70% after silvicultural treatments. The graphic illustrations provided by the stocking diagrams are useful and convenient diagnostic tools for the practical management of stand density.