Next Article in Journal
Impact of Environmental Conditions on Wood Anatomical Traits of Green Alder (Alnus alnobetula) at the Alpine Treeline
Previous Article in Journal
Insular Cycas micronesica Habitats Respond Similarly to Aulacaspis yasumatsui Invasion, Regardless of Co-Occurring Consumers
Previous Article in Special Issue
Individual Tree Height Increment Model for Quercus mongolica Secondary Forest in the Northeastern China Using Generalized Nonlinear Two-Level Mixed-Effects Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating Forest Site Quality Using the Biomass Potential Productivity Approach

1
College of Forestry, Shanxi Agricultural University, Jinzhong 030800, China
2
Research Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
3
Institute of Forestry, Tribhuwan University, Kathmandu 46000, Nepal
4
College of Mathematics and Statistics, Xinyang Normal University, Xinyang 464000, China
*
Author to whom correspondence should be addressed.
Forests 2024, 15(1), 23; https://doi.org/10.3390/f15010023
Submission received: 14 October 2023 / Revised: 30 October 2023 / Accepted: 19 December 2023 / Published: 21 December 2023
(This article belongs to the Special Issue New Theories and Methods in Tree and Stand Measurement and Modeling)

Abstract

:
Biomass productivity is of great significance for the evaluation of forest quality, which is important for the improvement of forest management. We propose the computational methods of biomass potential productivity (BPP) and biomass realistic productivity (BRP), both of which provide reliable practical guides for predicting forest growth under multi-aged, multi-species, and multi-layered canopy conditions. We used 2222 national forest inventory plots that were measured in four consecutive periods in the Jilin Province for this purpose. We analyzed and verified the computational methods of BPP based on the BRP and evaluated its practical significance. The results showed that growth models of the stand height, stand basal area, and stand biomass of four forest types (pure larch forest, larch broadleaf mixed forest, Mongolian oak pure forest, and Mongolian oak broadleaf mixed forest) fit adequately, BPP was greater than BRP, and this difference decreased with an increasing stand age, suggesting that the potential productivity of the middle-aged and young forest was higher than that of the mature forest, although the difference is minimal. In addition, the realistic productivity of stands with better site quality was close to the potential productivity, which is consistent with the biological significance of the potential productivity of the biomass. The degree of difference between the potential productivity of the biomass and the realistic productivity of biomass also decreases with the decline in site quality, and it can be termed as the potentially improved stand biomass. The BPP model was able to perform well in both the pure and mixed forests. The BRP not only verifies the rationality of the BPP but can be also used to quantify the forest site quality, which is helpful for evaluating forest growth and informed decision making in forestry.

1. Introduction

Site quality evaluation involves an assessment of the stand productivity potential of a forest stand with a certain method adopted, which is a necessary theoretical basis for the full utilization of the productive potential of the land and achieving scientific management of forest resources [1]. Since the 1920s, site type and site quality assessments have been receiving widespread attention [2]. So far, significant theoretical and technical progress has been made within research on forest site quality evaluation, and many evaluation methods have been proposed, such as the index plant method [3,4], the site index method [5,6], the growth index method [7,8,9], and the advanced exponential method [10]. On the one hand, traditional forest site quality evaluation is mainly aimed at plantation forests. For natural forests with a complex stand structure, an irregular matching of tree height class, age structure, and other factors, the method of forest site quality evaluation could be difficult to meaningfully unity [11]. On the other hand, the current research methodology assesses only the realistic growth of the forest, which essentially reflects the real productivity of the stand. In contrast, it is difficult for forest sites to reach maximum productivity under natural conditions, as there could be disturbances caused by human and natural factors. Researchers have attempted to accurately evaluate forest quality by constructing ecological models, but data acquisition (climate, soil, etc.) problems increase the difficulty of subsequent large-scale applications [12,13]. The remote sensing approach, which helps identify and distinguish the stand features (species composition and structure, etc.) is also commonly applied for site quality assessment [14]. However, the results can be substantially affected by the environment and accuracy. The existing methods of site quality evaluation are still limited in practice, and thus, we need a comprehensive method that is suitable for evaluation of site quality on a large scale.
For assessing forest site quality on a large-scale, data availability is the basis for developing forest quality indicators. Since forest biomass is a major driver of forest structure and diversity, its accumulation or decrease determines the forest site quality, and hence, forest biomass can be used as a reliable indicator of site productivity [15,16,17,18]. Accurate computation of forest biomass productivity and its dynamics can have a reference value, which can be used for forest development, management, and utilization, and such an estimate may also solve the ecological problems [19,20,21]. The forest site quality depends on the forest and growth potential [22]. Biomass potential productivity (BPP) is the maximum annual biomass growth that can be achieved at a given stand age for the given site type and stand type [23].
To overcome the difficulties mentioned above, this study assumes that under the same site conditions and stand type (similar species composition) with a similar stand structure and stand density, there would be an approximately similar growth in the stand height, stand basal area, and stand biomass. Based on the concept of basal area potential productivity and volume potential productivity, this study proposes a biomass potential productivity that is applicable to multi-aged and multi-layered mixed forests [24,25]. The methods of forest site quality evaluation presented in this study can provide practical guidance for scientific forest management [26,27]. The ninth national forest resources inventory in the Jilin Province, which is one of the key provinces for forest development in China [28], shows that the forest area is 7.85 million hectares, and the standing volume is 101.296 million cubic meters. Thus, scientific evaluation of forest site quality in the Jilin Province is conducive to raising social awareness of improvements to the current state of forests. Forest site quality evaluation in the Jilin Province is considered very important, as this helps improve the prediction performance of future stand productivity and can also lay the methodological foundation for effective forest management strategies, the sustainable development of China’s future forest ecosystems, and achieving carbon neutrality. This study aims to (1) propose a new forest site quality evaluation index: biomass potential productivity; (2) provide detailed computational methods of biomass potential productivity; (3) demonstrate the application of the computational methods using real data of the larch and Mongolian oak in the Jilin Province as an example; and (4) verify and validate the proposed methods and models. The methods and models proposed in this paper will provide both theoretical and practical knowledge that will be useful for forest site quality evaluation—one of the necessary informational inputs for scientific forest management.

2. Materials and Methods

2.1. Study Area

The forest area selected for this study lies in the northeast of China, which is a key forestry province (121°38′–131°19′ E, 40°50′–46°19′ N). The province has 45.2% forest cover, which includes 8,792,800 hectares of forest land and 1086 million cubic meters of accumulation. The reported mean annual temperature of the Jilin Province is 2–6 °C, and mean annual precipitation is 765.9 mm. The forest types are mainly medium-aged temperate coniferous and broad mixed forests, with a wide variety of tree species and forest site productivities involved.

2.2. Data Materials

We used Chinese national forest inventory (NFI) data from the Jilin Province, collected in 1999, 2004, 2009, and 2014, which consist of 2222 plots that are distributed across the four main multi-aged and mixed-species stand types. Of the 2222 plots, 384 were allocated in pure larch forest (stand type I), 258 in larch broadleaf mixed forest (stand type II), 821 in Mongolian oak pure forest (stand type III), and 759 in Mongolian oak broadleaf mixed forest (stand type IV). The distribution of the sample plots is shown in Figure 1.

2.3. Methods

2.3.1. Computing Sample Plot-Level Variables

We determined the number of trees, mean breast height, basal area, stand density index, biomass per plot, and stand density based on the basic NFI data. The lowest breast height diameter of a tree is >5 cm. The individual tree biomass was estimated using the method suggested by the State Forestry Administration (LY/T2258-2014) [29]. That is, individual tree biomass is the sum of trunk, bark, branches, leaves, and below-ground biomass, and the total biomass of the sample plot is the sum of all the individual tree biomasses (Table 1).
We summarized the relevant stand characteristics and various indicators (Table 2).

2.3.2. Site Grouping Based on Stand Height Growth

Due to the existence of different stand heights of the same age, we classified each stand type into the following five site groups based on the age range of the stand without considering the effects of environmental factors and management practices [24,30] as follows:
(1)
Determination of age group for each plot:
We divided the age of the sample plots in 5-year intervals.
A c = int ( a g e / 5 ) + 1
A c = { k i f min a g e + ( k 1 ) D a g e a g e < min a g e + k D a g e m i f min a g e + ( k 1 ) D a g e a g e max a g e k = 1 , 2 , m 1 k = m
where A c is the age grouping of the plot; a g e is the stand age of the plot; int ( a g e / 5 ) is an integrating function with parameter a g e / 5 , it is a vectorial mode of operation with an interval of 5, so the initial value is set to 5; m is the number of groupings across the age span of the sample ( m = 5 ) ; min a g e is the minimum age; max a g e is the maximum age; D a g e is the differences between age groups ( max a g e min a g e ) / m . After computation of (Equations (1) and (2)), a sample plot can only belong to one age group.
(2)
Determination of the initial height group of each sample plot:
For each plot, the stand height hierarchy at the beginning of the period was determined in each age group to which it belongs.
H i l = { k i f min H i A + ( k 1 ) D H i A H i < min H i A + k H i A H i l = n i f min H i A + ( k 1 ) D H i A H i < max H i A k = 1 , 2 , n 1 k = n
where H i is the stand height; H i l is the stand height classification of the plot, which was assumed to be l ; min H i A is the minimum stand height for the age group to which the plot belongs; max i A is the maximum stand height for the age group to which the plot belongs; n is the stand height group ( n = 5 ) ; D H i A is the difference in stand height within the age group to which the sample site belongs ( max H i A min H i A ) / n .
(3)
Determination of height-age model form:
H = f h ( A ge | F , L )
where f h ( ) is the height-age model; H is the stand height; A ge is the stand age; F is the stand type (I, II, III, IV); L is the site group (1, 2, 3, 4, 5).

2.3.3. Biomass Potential Productivity and Realistic Productivity

Stand productivity includes both the realistic productivity and potential productivity. The biomass potential productivity ( B P P ) refers to the biomass increment hat seeks a maximum objective function from the given stand density index ( S ) interval under the conditions of growth of stand basal area and stand biomass corresponding to a given stand type ( F ) and the corresponding parameter estimates, age (expressed by A g e 0 ) and site group ( L ):
B P P = f ( A g e 0 , S , Φ ^ G , Φ ^ B | F , L ) S [ S min , S max ]
where A g e 0 is the specific stand age, Φ ^ G and Φ ^ B parameters for stand basal area and stand biomass growth models, respectively. S min and S max is the set a feasible region of S ; F is the stand type (I, II, III, IV); L is the site group (1, 2, 3, 4, 5). The specific calculation steps are as follows:
Step 1, we need to have 5 known conditions:
(1)
Given the site group ( L ) of the stand type ( F ). Site group ( L ) was known from the Section 2.3.2.
(2)
Given the specific stand age to be computed ( A g e 0 )
(3)
Given the density index S search interval [ S min , S max ] , the feasible region of S is assumed to vary from 30 to 3000 [24].
(4)
Form of basal area growth model corresponding to the specific site group ( L ) of the stand type ( F ) is
G = f G ( A g e , S , Φ ^ G | F , L )
where G is the basal area growth model of the stand type ( F ); A g e is the specific stand age; S is the stand density index of the specific stand age; Φ G is the model parameter; F is the stand type (I, II, III, IV); L is the site group (1, 2, 3, 4, 5).
(5)
Form of biomass growth model corresponding to the specific site group ( L ) of the stand type ( F ) is
B = f B ( A g e , S , Φ ^ B | F , L )
where B is the biomass growth model of the stand type ( F ); A g e is the specific stand age; S is the stand density index of the specific stand age; Φ ^ B is the parameter vector to be estimated; F is the stand type (I, II, III, IV); L is the site group (1, 2, 3, 4, 5).
Step 2, finding a reasonable stand density index (S) that maximizes biomass annual growth (BI).
(1)
Given a feasible region of S ( S min , S max ) ,a stand age ( A g e 0 ), a stand type ( F ), a site group ( L ), the error term ( ε ), and iteration t = 1. Computing of initial values of stand density index ( S ) at 4 points ( S 1 A g e 0 , S 2 A g e 0 , S 3 A g e 0 , S 4 A g e 0 ) using the golden section method [24].
(2)
By computing B 0 at A g e 0 and B 1 at A g e 1 = A g e 0 + 1 , we can get biomass annual growth ( B I 1 A g e 0 , B I 2 A g e 0 , B I 3 A g e 0 , B I 4 A g e 0 ).
B I = B 1 B 0 = f B ( A g e 1 , S 1 , Φ ^ B ) f B ( A g e 0 , S 0 , Φ ^ B )
where f B is the biomass growth model; S 1 is the assumed stand density index at A g e 1 ; S 0 is the assumed stand density index at A g e 0 .
(3)
Finding the S that made the B I maximized.
Step 3, Computing biomass potential productivity (BPP).
B P P = ( B I 2 A g e 0 + B I 3 A g e 0 ) / 2
where B I 2 A g e 0 and B I 3 A g e 0 are the S 2 A g e 0 and S 3 A g e 0 corresponds to the maximizes biomass growth.
The detailed procedure of the above B P P algorithm is shown in the Supplementary Information S1.
The biomass realistic productivity ( B R P ) is the continuous annual growth of the stand biomass computed from the stand age ( A g e 0 ), stand density index ( S ), stand type ( F ) and site groups ( L ), combined with the basal area growth model, and the biomass growth model. The computational procedures of B R P are similar to B P P . It should be noted that, unlike the B P P calculation, the stand characteristics, S from field surveys, should be provided for the B R P calculation. The stand density index ( S 0 ) for biomass realistic productivity in A g e 0 is the measured data from the sample plots. Theoretically, the realistic productivity is lower or closer to the potential productivity. The formula is as follows:
B R P = B I R 1 B I R 0 = f B ( A g e 1 , S 1 , Φ ^ B | F , L ) f B ( A g e 0 , S 0 , Φ ^ B | F , L )
where B R P is the biomass realistic productivity; B I R 1 is the realistic biomass annual growth at A g e 1 ; B I R 0 is the realistic biomass annual growth at A g e 0 ; S 1 is the stand density index at A g e 1 ; S 0 is the realistic stand density index at A g e 0 . F is the stand type (I, II, III, IV); L is the site group (1, 2, 3, 4, 5).
The detailed procedure of the above B R P algorithm is shown in the Supplementary Information S2.

2.3.4. Computing the Available Potential Improved Stand Biomass

The potential for improving the biomass productivity is obtained from the difference between the potential biomass productivity and the realistic biomass productivity.
B P I = B P P B R P
where B P I is the available improved potential stand biomass; B P P and B R P are the biomass potential productivity and biomass realistic productivity corresponding to a certain stand age.
All computations in this study were implemented in R 4.3.0.

2.3.5. Modeling and Parameter Estimation

The Chapman-Richards function was selected as a basic model for fitting stand height data [25], and the model form is:
H ( F , L ) = 1.3 + a i ( F , L ) [ ( 1 e b i ( F , L ) A g e ) ] c i ( F , L ) + ε ( F , L )
where H ( F , L ) is stand height to the F t h stand type for the L th site group; A g e is the stand age for the plots; a i ( F , L ) , b i ( F , L ) , c i ( F , L ) are the parameters corresponding to the model; ε ( F , L ) is the error term.
The stand basal area growth model adopts the functional form by Chapman-Ricards [25]:
G ( F , L ) = a j ( F , L ) [ 1 exp ( b j ( F , L ) ( S / 1500 ) c j ( F , L ) A g e ) ] d j ( F , L ) + ε ( F , L )
where G ( F , L ) is the stand basal area to the F t h stand type for the L th site group (m2·ha−2); S is the stand density index; A g e is the stand age; a j ( F , L ) , b j ( F , L ) , c j ( F , L ) and d j ( F , L ) are the parameters corresponding to the model; ε ( F , L ) is the error term.
The biomass growth model we used here is exactly the same model form as the basal area growth model:
B ( F , L ) = a k ( F , L ) [ 1 exp ( b k ( F , L ) ( S / 1500 ) c k ( F , L ) A g e ) ] d k ( F , L ) + ε ( F , L )
where B ( F , L ) is the stand biomass to the F t h stand type for the L th site group (t2·ha−2); S is the stand density index; A g e is the stand age; a k ( F , L ) , b k ( F , L ) , c k ( F , L ) and d k ( F , L ) are the parameters corresponding to the model; ε ( F , L ) is the error term.

2.3.6. Model Evaluation

We estimated models (12), (13) and (14) using the R package and evaluated their fitting performance using the coefficient of determination ( R 2 ) and root mean square error ( RMSE ).
R 2 = 1 i = 1 k j = 1 n i ( y i j y ^ i j ) 2 i = 1 k j = 1 n i ( y i j y ¯ i j ) 2
R M S E = i = 1 k j = 1 n i ( y i j y ^ i j ) 2 / n
where: y i j is the j th actual value of the i th sample plot; the j th estimated value of the i th sample plot; y ¯ i j is the mean H or B A or B of the observations; k is the number of classifications; n is the number of sample plots in the k th category. n is the number of all sample plots.

2.3.7. Reference Age

Reference age refers to an age after which the growth of the stand height tends to be stable, and the age at which the mean growth of the stand height or timber volume is the largest. The reference age of mixed forests is determined by dominant tree species or tree species more frequently occurring in a stand. According to the LY/T 2415-2015 “Technical Regulations for the Preparation of Status Index Tables”, the reference age for both larch and Mongolian oak are 30 years [31].

3. Results

3.1. Parameter Estimates for Growth Models of Stand Height, Basal Area, and Biomass

Based on the NFI data of four forest types in the Jilin Province, the height–age model, basal area growth model, and biomass growth model of forest stands models (12), (13), and (14) were fitted, and parameter estimates were obtained (Table 3). All the parameter estimates were significant (p < 0.05), and each model described more than 70% variations in the response variables of the interest (Table 4). The fitting effect of each growth model was good, and the results were feasible.
Using models (12), (13) and (14) and their parameter estimates, we plotted stand height, basal area and biomass as a function of stand age ( Y e a r ) for different stand types at different site groups. Basal area and biomass are determined by stand age and stand density. For showing the relationship between basal area, biomass, and stand age more clearly in a 2-dimensional plan view, we assumed stand density ( S ) to be 1000. As we can be seen in Figure 2, models (12), (13) and (14) all exhibited a clear division by site groups, and the models are reasonable.

3.2. Computation of BPP by Stand Type

We plotted the BPP by stand type in the age range between 5 and 150 years (Figure 3), and there appeared to be a unimodal BPP curve for each site group and stand age in different age groups, which is inversely correlated. That is, the potential productivity decreased with stand age. The BPP of each stand type decreased with an increasing site group, and BPP was the highest when the site group was 1 (site group 1 stands for the highest site quality and 5 stands for the lowest site quality). There was a negative relationship between the potential productivity and the stand class for a given age (i.e., the BPP decreased with an increasing site group).

3.3. BRP and BPI for Each Stand Type

The BRP reflects the realistic biomass accumulated in a stand during a year [27]. The realistic productivity of biomass was computed based on the stand age and stand density index. The computed biomass realistic productivity (BRP) and the available improved potential stand biomass (BPI) are shown in Figure 4. The trends of BRP and BPI seem to be similar for each stand type, with BRP showing an overall increasing trend with the stand age followed by decreasing trend, while the BPI decreases with stand age. The biomass potential productivity (BPP) is far from being brought into play at the middle-aged and young stage, which suggests that the earlier the forest quality is improved, the higher the stand productivity would be.

3.4. Verification of BPP and Evaluation of Realistic Stands

The BPP and BRP for four forest types are shown in Figure 5, Figure 6, Figure 7 and Figure 8. The difference between the BPP and BRP corresponding to each stand is large, the degree of difference is gradually decreasing as the stand age increases, and they are very close when the stand reaches the near-mature or mature stage. An identical pattern was observed for all site groups.

3.5. Site Quality Evaluation at Reference Age

The potential productivity at a reference age can reflect the growth potential when the stand height tends to be stable and the mean growth is at its maximum (Table 5). In terms of the potential productivity of the biomass at reference age, the potential productivity levels of the four stand types were different at different site groups, and all showed that the potential productivity was the highest at site group 1 and the lowest at site group 5. Assuming a reference age of 30 for larch pure forest as an example, the potential productivity of the biomass corresponding to group 1 is 3.642 t·ha−2·year−1, which is 23% higher than that of site group 5 (2.966 t ha−2 year−1). Thus, site group 1 has a higher potential productivity and better stand quality.

4. Discussion

An accurate estimate of site productivity is an important part of formulating effective forest management strategies, and it is also the core basis for evaluating the basic status of sites in the ecological environment [32,33]. Biomass potential productivity and realistic productivity refer, respectively, to the maximum biomass accumulated in a stand and the realistic biomass accumulated in a stand during a year. We used the biomass potential productivity as a site quality evaluation index and ignored the traditional practice of selecting stands with good or poor site quality for site quality assessment. As mentioned in Section 1, we have combined the stand classification and site quality evaluation, introduced the detailed BPP computational processes, and have verified the rationality of this index by BRP, confirming a precise method that is applicable to forest site evaluation. The potential productivity of the biomass can be used as an indicator of site quality, providing reliable quantitative information for forest management decisions.
In the traditional methods of site quality evaluation, the dominant height at any reference age is chosen as the reference index. This method can be inaccurate when tree height growth is constrained by stand density or competition [34,35]. The potential productivity of the biomass as an evaluation index takes into account the influence of stand density and sets the maximum value of stand density, which is more reasonable for reflecting the actual situation of forest complexity in multi-aged and multi-layered mixed forests. In our study, the biomass potential productivity decreased with stand age in each forest type, and this decreasing pattern is consistent with the characteristics of the potential productivity (Figure 3). Because young stands would have a higher number of trees or a high stand density, consequently, there would be a higher potential productivity. That is, at a point in time, the potential productivity is a local extreme of the successive annual growth of the biomass. If the site conditions are favorable and the management practices employed are appropriate, the realistic productivity of the stand can approach or reach the potential productivity (Figure 4). The BPI is the available improved potential stand biomass. This indicator can also rationalize the priority of different stands for upgrading. For stands of the same age, the higher the BPI is, the higher the priority of promotion is; conversely, the lower the priority of promotion is (Figure 5, Figure 6, Figure 7 and Figure 8). This is consistent with the trend of site quality in terms of the potential productivity evaluated in terms of the stand basal area approach [24]. This indicates that the management of young and middle-aged forests with good site quality should be emphasized in the process of forest nurturing. Our proposed concept of biomass potential productivity is reasonable, and our method for quantifying the site quality is feasible. In fact, the selection of evaluation indices is highly related to the purpose of forest site evaluation. The method of the basal area potential productivity emphasizes the stand diameter growth, which is more suitable for the evaluation of site quality aimed at producing stands with large-sized timber. However, for the special-purpose forests and public ecological welfare forests, where emphasis is placed on the forest biomass growth, the biomass productivity index can be intuitive and comprehensive to reflect the overall forest ecosystem quality. Therefore, it is necessary to choose the appropriate evaluation method according to the objectives.
The models we propose in the article are reliably applicable to other stand types in terms of species of interest and forested regions. The computation of the biomass potential productivity requires, in addition to the given biomass growth model and its parameter estimates, a stand height growth model and a basal area growth model and their parameter estimates. All the growth models included in the biomass potential productivity model should have high accuracies in order to accurately estimate the forest site productivity. In our study, the whole-forest holistic model form was chosen, and both the stand basal area and stand biomass growth models were fitted with high accuracy for the actual stand conditions of the multi-aged and multi-layered mixed forests [36]. In the practical application, the stand basal area growth model or biomass growth model may be chosen based on the available stand type data. In conclusion, the biomass potential productivity approach is flexible enough to be applicable to a large-scale forestry, as data used to construct the biomass growth model are based on a large number of NFI plots, which cover the extensive forests (Figure 1). Our methods and models can provide both theoretical and practical knowledge, and their applicability will be useful for a larger scale site quality evaluation and even for tree species suitability evaluation.
We conducted this study under the assumption of equal annual growth of biomass regardless of the tree species and their growth behaviors. This is an ideal situation that is difficult to actually achieve. We were searching for the maximum value of the objective function in a given search interval for any stand density, which enables the computation of the potential productivity. The result obtained was that young stands would have a higher stand density, and thus have a higher potential productivity of the biomass. In contrast, near-mature and mature stands would have a lower stand density and lower biomass potential productivity. We may construct the segmentation function to solve this problem by assuming that a given forest type has the highest production potential at a given reference age. The maximum value of the objective function can be found for a given interval of the stand density, so that the potential productivity in the segmentation function can be obtained. In the actual forest growth, the natural process of self-thinning occurs when the age increases, and an increased age is accompanied by a decrease in the stand density of young and middle-aged stands. At the reference age, the stand density would be relatively stable [31,37]. Theoretically, the self-thinning model might improve the accuracy of the biomass potential productivity model if included; however, this requires the construction of the self-thinning model, which demands long-term time series data. The steps of computing biomass potential productivity may become more complicated when a self-thinning model is included. In our future studies, we are going to consider calculating the biomass potential productivity using the self-thinning model. In addition, as ecological problems and climate issues have attracted significant global attention, climate factors need to be included in the forest site quality evaluation methods [38,39].

5. Conclusions

In this study, we computed the potential productivity of the stand basal area and biomass for four stand types in the Jilin Province using four consecutively measured NFI data. We presented the growth model forms and their parameter estimates for stand height, basal area, and biomass, each of which had a reasonable fitting accuracy and model stability. The potential productivity is the local extreme value of the continuous annual growth of the forest biomass, and the realistic productivity of a stand can only be less than or equal to the potential productivity. If the site condition is good, the realistic productivity can approach or reach the potential productivity if the stand is in a good condition and the appropriate management measures are applied to the stands. In addition, the potential productivity of the middle-aged and young stands was significantly higher than that of the mature stands, and stands with better growth had a higher potential productivity than those with poorer stand growth. When evaluating the forest site quality, the biomass potential productivity at the reference age or close to the mature stand stage can often be used, which could provide a relatively high accuracy. Our potential biomass productivity concept is reasonable, and the method for its computation is feasible and applicable to multi-aged and multi-layered mixed forests. Forest managers should pay attention to the development of young and middle-aged forests with better site quality for the improvement of forest production. Our study provides technical support for scientific forest management, which will have significant implications for the improvement of forest management and promotion of the modern forestry.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/f15010023/s1, Supplementary Information S1: Method of biomass potential productivity (BPP); Supplementary Information S2: Method of biomass realistic productivity (BRP). All computations in this study were implemented in R4.3.0. The sample site factor computations such as biomass were implemented in the foreign package, the openxlsx package, the dplyr package, and the tidyr package. Stand height, basal area, and biomass growth models were used in the nlme package. The plotting was carried out with the ggplot package in R4.3.0 (www.r-project.org) and orgin 2022 software (www.origin.com/store/odc/store-plugin/views/checkoutLoading.html, accessed on 18 December 2023).

Author Contributions

X.Y., R.P.S., G.D., J.G. and L.F. (Linyan Feng). conceived the study. X.Y., L.F. (Liyong Fu), and L.P. performed the analysis and wrote the initial draft of the paper. All authors contributed to interpreting the results and the improvement of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (Grant Number: 2022YFD2201005) and Accurate estimation of carbon sequestration potential of Plantation (Grant Number: DL2022052001L).

Data Availability Statement

The data used for the analysis and graphing in this paper are the 4-period national forest inventory data provided by the State Forestry and Grassland Administration (SFA) (https://www.forestdata.cn/index.html, accessed on 18 December 2023).

Acknowledgments

We also appreciate the valuable comments and constructive suggestions from three anonymous referees and the Associate Editor that helped improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, Y.; Wu, B.; Liu, Y.; Ju, H.; Wang, C. Research Progress of Site Quality Evaluation. World For. Res. 2012, 25, 47–52. [Google Scholar]
  2. Carmean, W.H. Forest Site Quality Evaluation in The United States. Adv. Agron. 1975, 27, 209–269. [Google Scholar]
  3. Zas, R.; Alonso, M. Understory vegetation as indicators of soil characteristics in northwest Spain. For. Ecol. Manag. 2002, 171, 101–111. [Google Scholar] [CrossRef]
  4. Wilson, S.M.; Pyatt, D.G.; Malcolm, D.C.; Connolly, T. The use of ground vegetation and humus type as indicators of soil nutrient regime for an ecological site classification of British forests. For. Ecol. Manag. 2001, 140, 101–116. [Google Scholar] [CrossRef]
  5. Kimberley, M.; West, G.G.; Dean, M.G.; Knowles, L.R. The 300 index—A volume productivity index for radiata pine. N. Z. J. For. 2005, 50, 13–18. [Google Scholar]
  6. Elfving, B.; Kiviste, A. Construction of site index equations for Pinus sylvestris L. using permanent plot data in Sweden. For. Ecol. Manag. 1997, 98, 125–134. [Google Scholar] [CrossRef]
  7. Mailly, D.; Gaudreault, M. Growth intercept models for black spruce, jack pine and balsam fir in Quebec. For. Chron. 2005, 81, 104–113. [Google Scholar] [CrossRef]
  8. Pretzsch, H. Forest Dynamics, Growth and Yield; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  9. Berrill, J.P.; O’Hara, K. Estimating site productivity in irregular stand structures by indexing the basal area or volume increment of the dominant species. Can. J. For. Res. 2014, 44, 92–100. [Google Scholar] [CrossRef]
  10. Guo, J.P.; Wang, J.R. Comparison of height growth and growth intercept models of jack pine plantations and natural stands in northern Ontario. Can. J For. Res. 2006, 36, 2179–2188. [Google Scholar] [CrossRef]
  11. Skovsgaard, J.A.; Vanclay, J.K. Forest site productivity: A review of the evolution of dendrometric concepts for even-aged stands. Forestry 2008, 81, 13–31. [Google Scholar] [CrossRef]
  12. Roxburgh, S.H.; Karunaratne, S.B.; Paul, K.I.; Lucas, R.M.; Armston, J.D.; Sun, J. A revised above-ground maximum biomass layer for the Australian continent. For. Ecol. Manag. 2019, 432, 264–275. [Google Scholar] [CrossRef]
  13. Czucz, B.; Molnar, Z.; Horvath, F.; Nagy, G.G.; Botta-Dukat, Z.; Toeroek, K. Using the natural capital index framework as a scalable aggregation methodology for regional biodiversity indicators. J. Nat. Conserv. 2012, 20, 144–152. [Google Scholar] [CrossRef]
  14. Xiao, Y.; Ouyang, Z.Y.; Wang, L.Y.; Rao, E.; Jiang, L.; Zhang, L. Spatial patterns of ecosystem quality in Inner Mongolia and its driving forces analysis. Acta Ecol. Sin. 2016, 36, 6019–6030. [Google Scholar]
  15. Liu, W.; Lu, F.; Luo, Y.; Bo, W.; Kong, L.; Zhang, L.; Liu, B.; Ouyang, Z.; Wang, X. Human influence on the temporal dynamics and spatial distribution of forest biomass carbon in China. Ecol. Evol. 2017, 7, 6220–6230. [Google Scholar] [CrossRef]
  16. Ouyang, Z.; Zheng, H.; Xiao, Y.; Polasky, S.; Liu, J.; Xu, W.; Wang, Q.; Zhang, L.; Xiao, Y.; Rao, E.; et al. Improvements in ecosystem services from investments in natural capital. Science 2016, 352, 1455–1459. [Google Scholar] [CrossRef]
  17. Ji, B.Y.; Tao, J.X.; Wang, W.W. Forest site quality evaluation based on biomass data from Continuous Forest Inventories. J. Southwest For. Univ. 2012, 32, 45–50. [Google Scholar]
  18. Kukrety, S.; Wilson, D.C.; D’Amato, A.W.; Becker, D.R. Assessing sustainable forest biomass potential and bioenergy implications for the northern Lake States region, USA. Biomass Bioenergy 2015, 81, 167–176. [Google Scholar] [CrossRef]
  19. Becker, D.R.; Moseley, L. A supply chain analysis framework for assessing state-level forest biomass utilization policies in the United States. Biomass Bioenergy 2011, 35, 1429–1439. [Google Scholar] [CrossRef]
  20. Chen, Q.; Chen, Y.H.; Wang, M.J.; Jiang, W.G.; Hou, P.; Li, Y. Ecosystem quality comprehensive evaluation and change analysis of Dongting Lake in 2001–2010 based on remote sensing. Acta Ecol. Sin. 2015, 35, 4347–4356. [Google Scholar]
  21. Sprugel, D.G. Density, biomass, productivity, and nutrient-cycling changes during stand development in wave-regenerated balsam fir forests. Ecol. Monogr. 1984, 54, 165–186. [Google Scholar] [CrossRef]
  22. Ding, Z.; Li, R.; O’Connor, P.; Zheng, H.; Ouyang, Z. An improved quality assessment framework to better inform large-scale forest restoration management. Ecol. Indic. 2021, 123, 107370. [Google Scholar] [CrossRef]
  23. LY/T3315-2022; Technical Regulation for Forest Site Quality Evaluation. China Forestry Industry Standards. National Forestry and Grassland Administration: Beijing, China, 2022. (In Chinese)
  24. Fu, L.; Sharma, R.P.; Zhu, G.; Li, H.; Hong, L.; Guo, H.; Duan, G.; Shen, C.; Lei, Y.; Li, Y.; et al. A Basal Area Increment-Based Approach of Site Productivity Evaluation for Multi-Aged and Mixed Forests. Forests 2017, 8, 119. [Google Scholar] [CrossRef]
  25. Lei, X.; Fu, L.; Li, H.; Li, Y.; Tang, S. Methodology and applications of site quality assessment based on potential mean annual Increment. Sci. Silvae Sin. 2018, 54, 116–126. (In Chinese) [Google Scholar]
  26. Mensah, A.A.; Holmström, E.; Nyström, K.; Nilsson, U. Modelling potential yield capacity in conifers using Swedish long-term experiments. For. Ecol. Manag. 2022, 512, 120162. [Google Scholar] [CrossRef]
  27. Duan, G.X.; Zheng, Y.L.; Hong, L.; Song, X.; Fu, L. A Potential Productivity-Based Approach of Site Quality Evaluation for Larch Pure Forest and Birch-Aspen Mixed Forest. Sci. Silvae Sin. 2022, 58, 9. (In Chinese) [Google Scholar]
  28. National Forestry and Grassland Administration. China Forest Resources Report: The Ninth National Forest Resources Inventory; China Forestry Publishing House: Beijing, China, 2019. (In Chinese) [Google Scholar]
  29. LY/T2258-2014; Technical Regulation on Methodology for Tree Biomass Modeling. China Forestry Industry Standards. National Forestry and Grassland Administration: Beijing, China, 2014. (In Chinese)
  30. Tang, S.; Meng, C.H.; Meng, F.R.; Wang, Y.H. A growth and self-thinning model for pure even-age stands: Theory and applications. For. Ecol. Manag. 1994, 70, 67–73. [Google Scholar] [CrossRef]
  31. LY/T2415-2015; Technical Regulations for Drafting of Site Index Table. China Forestry Industry Standards. National Forestry and Grassland Administration: Beijing, China, 2015. (In Chinese)
  32. Heiderman, R.R.; Kimsey, M.J., Jr. A species-specific, site-sensitive maximum stand density index model for Pacific Northwest conifer forests. Can. J. For. Res. 2021, 51, 1166–1177. [Google Scholar] [CrossRef]
  33. Aertsen, W.; Kint, V.; Muys, B.; Van Orshoven, J. Effects of scale and scaling in predictive modelling of forest site productivity. Environ. Model. Softw. 2012, 31, 19–27. [Google Scholar] [CrossRef]
  34. Chen, Y.F. Research on standard age affecting analysis of site quality evaluation. For. Res. Beijing 2010, 23, 283–287. [Google Scholar]
  35. Adame, P.; Ca Ellas, I.; Roig, S.; Río, M.D. Modelling dominant height growth and site index curves for rebollo oak (Quercus pyrenaica Willd.). Ann. For. Sci. 2006, 63, 929–940. [Google Scholar] [CrossRef]
  36. Tang, S.Z.; Li, Y.; Fu, L.Y. Statistical Foundation for Bio-Mathematical Models, 2nd ed.; Higher Education Press: Beijing, China, 2015; 435p. (In Chinese) [Google Scholar]
  37. Tang, S.; Meng, F.R.; Meng, C.H. The impact of initial stand density and site index on maximum stand density index and self-thinning index in a stand self-thinning model. For. Ecol. Manag. 1995, 75, 61–68. [Google Scholar] [CrossRef]
  38. Popkin, G. How much can forests fight climate change? Nature 2019, 565, 280–283. [Google Scholar] [CrossRef]
  39. Bontemps, J.-D.; Bouriaud, O. Predictive approaches to forest site productivity: Recent trends, challenges and future perspectives. Forestry 2014, 87, 109–128. [Google Scholar] [CrossRef]
Figure 1. Overview of the study area and NFI plot locations.
Figure 1. Overview of the study area and NFI plot locations.
Forests 15 00023 g001
Figure 2. Growth curves for different site groups. For stand type, see definition in Table 2. For Site group, see definition in Section 2.3.2.
Figure 2. Growth curves for different site groups. For stand type, see definition in Table 2. For Site group, see definition in Section 2.3.2.
Forests 15 00023 g002
Figure 3. The biomass productivity potential (BPP) of different forest types on different site groups. For stand type, see definition in Table 2. For site group, see definition in Table 3.
Figure 3. The biomass productivity potential (BPP) of different forest types on different site groups. For stand type, see definition in Table 2. For site group, see definition in Table 3.
Forests 15 00023 g003
Figure 4. Biomass realistic productivity (BRP) and the available improved potential stand biomass (BPI) of each stand type. For stand type, see definition in Table 2.
Figure 4. Biomass realistic productivity (BRP) and the available improved potential stand biomass (BPI) of each stand type. For stand type, see definition in Table 2.
Forests 15 00023 g004
Figure 5. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type I (pure larch forest). For site group, see definition in Table 3.
Figure 5. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type I (pure larch forest). For site group, see definition in Table 3.
Forests 15 00023 g005
Figure 6. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type II (larch broadleaf mixed forest). For site group, see definition in Table 3.
Figure 6. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type II (larch broadleaf mixed forest). For site group, see definition in Table 3.
Forests 15 00023 g006
Figure 7. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type III (Mongolian oak pure forest). For site group, see definition in Table 3.
Figure 7. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type III (Mongolian oak pure forest). For site group, see definition in Table 3.
Forests 15 00023 g007
Figure 8. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type IV (Mongolian oak broadleaf mixed forest). For site group, see definition in Table 3.
Figure 8. The biomass potential productivity (BPP) and biomass realistic productivity (BRP) of each sample plot against stand age (year) for group = 1, …, 5 for stand type IV (Mongolian oak broadleaf mixed forest). For site group, see definition in Table 3.
Forests 15 00023 g008
Table 1. Computation of stand-level variables.
Table 1. Computation of stand-level variables.
VariableFormula
Mean of the diameter at breast height (cm) D = 1 n D i 2
Number of trees per hectare N = n p a
Stand basal area (m2 ha−2) B A = ( π D i 2 / 40 , 000 ) / p a
Stand density index S = N ( D / 20 ) 1.605
Stand biomass (t ha−2) B = B P p a
Note: n is the number of the trees per sample plot; D i is the diameter at breast height; P a is the plot area (0.067 hectares); and B p is the total biomass per sample plot.
Table 2. Description statistics of stand variables by stand type.
Table 2. Description statistics of stand variables by stand type.
Stand TypeVariableMin.Max.MeanSD
IStand height (m)4.1023.0012.193.99
Stand age (year)8.0067.0028.6011.88
Stand density index61.851057.62456.77221.52
Stand basal area (m2 ha−2)10.43259.4658.8935.79
Stand biomass (t ha−2)5.12126.0829.2117.64
IIStand height (m)4.0023.0012.724.216
Stand age (year)7.0068.0030.6213.70
Stand density index5.971006.8509.83237.48
Stand basal area (m2 ha−2)0.11150.0327.3127.21
Stand biomass (t ha−2)0.04219.3149.0445.36
IIIStand height (m)1.5022.0013.053.74
Stand age (year)5.00163.0062.5826.69
Stand density index137.051670.96694.60265.11
Stand basal area (m2 ha−2)3.4355.0520.098.43
Stand biomass (t ha−2)20.37485.65136.0567.17
IVStand height (m)4.0024.0014.273.58
Stand age (year)8.00169.0063.1227.20
Stand density index149.981614.57682.70261.69
Stand basal area (m2 ha−2)3.6553.4820.028.79
Stand biomass (t ha−2)20.38421.27135.8770.54
Note: stand type I is the pure larch forest; stand type II is the larch broadleaf mixed forest; stand type III is the Mongolian oak pure forest; stand type IV is the Mongolian oak broadleaf mixed forest.
Table 3. Parameter estimates of growth models by stand type.
Table 3. Parameter estimates of growth models by stand type.
Stand TypeSite Group Model (12)Model (13)Model (14)
a i b i c i a j b j c j d j a k b k c k d k
I123.190.0360.9343.180.034.710.21528.030.0000110.090.11
219.8142.11505.12
316.7940.86493.13
413.5338.89454.20
59.7837.69430.16
II122.580.051.34105.140.000015.180.11685.030.0000210.020.13
219.04104.09673.58
316.48100.27645.42
414.0998.47576.05
59.8897.28549.29
III120.220.031.0896.760.000029.330.11539.570.0000713.990.08
217.4194.01516.80
315.1192.26498.15
412.5687.47442.82
59.3384.89419.56
IV122.740.0190.6864.700.0025.810.17506.840.00582.540.42
219.5363.09485.45
317.2161.88460.61
414.6361.27458.08
510.9759.71447.23
Note: stand type I is the pure larch forest; stand type II is the larch broadleaf mixed forest; stand type III is the Mongolian oak pure forest; stand type IV is the Mongolian oak broadleaf mixed forest. “—” means same parameters.
Table 4. Evaluation indices of the growth models by stand type.
Table 4. Evaluation indices of the growth models by stand type.
Stand TypeModel (12)Model (13)Model (14)
R 2 R M S E R 2 R M S E R 2 R M S E
I0.97260.700.97941.030.832414.63
II0.97250.700.97941.040.748234.24
III0.96360.700.97231.400.857125.37
IV0.96250.670.97871.270.926419.12
Note: R2 is the coefficient of determination; RMSE is the root mean square error.
Table 5. Biomass potential productivity by stand type for each site group.
Table 5. Biomass potential productivity by stand type for each site group.
Stand TypePotential ProductivitySite Group
12345
IBiomass (t ha−2 year−1)3.6423.4833.4003.1322.966
IIBiomass (t ha−2 year−1)5.5525.4595.2304.6684.452
IIIBiomass (t ha−2 year−1)7.7637.4367.1676.3716.037
IVBiomass (t ha−2 year−1)7.6567.3336.9586.9206.756
Note: For stand type, see definition in Table 2. For site group, see definition in Table 3.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, X.; Feng, L.; Sharma, R.P.; Duan, G.; Pang, L.; Fu, L.; Guo, J. Evaluating Forest Site Quality Using the Biomass Potential Productivity Approach. Forests 2024, 15, 23. https://doi.org/10.3390/f15010023

AMA Style

Yan X, Feng L, Sharma RP, Duan G, Pang L, Fu L, Guo J. Evaluating Forest Site Quality Using the Biomass Potential Productivity Approach. Forests. 2024; 15(1):23. https://doi.org/10.3390/f15010023

Chicago/Turabian Style

Yan, Xingrong, Linyan Feng, Ram P. Sharma, Guangshuang Duan, Lifeng Pang, Liyong Fu, and Jinping Guo. 2024. "Evaluating Forest Site Quality Using the Biomass Potential Productivity Approach" Forests 15, no. 1: 23. https://doi.org/10.3390/f15010023

APA Style

Yan, X., Feng, L., Sharma, R. P., Duan, G., Pang, L., Fu, L., & Guo, J. (2024). Evaluating Forest Site Quality Using the Biomass Potential Productivity Approach. Forests, 15(1), 23. https://doi.org/10.3390/f15010023

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