A Basal Area Increment-Based Approach of Site Productivity Evaluation for Multi-Aged and Mixed Forests

Accurate estimates of forest site productivity are essential for environmental planning and forest management. In this study, we developed a new productivity index, hereafter termed basal area potential productivity index (BAPP), to estimate site productivity for irregular and complex forests characterized by multi-aged, multi-species, and multi-layer stands. We presented the biological relevance of BAPP with its computational details. We also compared BAPP against basal area realized productivity (BARP) in order to verify the practicability and reliability of BAPP. Time-series data of the national forest inventory on 1912 permanent sample plots that were located in two main forest types and consisted of oak-dominated mixed forests and other broadleaf forests in northeast China were used to demonstrate the application of BAPP. The results showed that the value of BAPP for each sample plot was larger than or equal to the corresponding BARP value for each forest type. For appropriately managed stands with relatively better site conditions, the values of both BARP and BAPP were almost identical. The values of the difference between BAPP and BARP could therefore be used to effectively assess forest site productivity. Meanwhile, BAPP also provides much reliable and valuable information that can aid decision-making in forest management.


Introduction
Forest site productivity, which is a quantitative estimate of the potential of a site to produce wood or plant biomass, provides information that is essential for predicting the rate of change of wood or plant biomass production potential under different management regimes [1][2][3].Reliable information Forests 2017, 8, 119 2 of 18 on site productivity is an important prerequisite for forest management planning.The terms "site productivity" and "site quality" are often used interchangeably, and concepts of these terms were developed by forest managers to project the expected production potential of the trees of a certain species on a given site [4,5].
Various approaches have been used to estimate site productivity [1,[6][7][8], and among the methods, a dominant height-age model approach (i.e., the site index (SI)) is most commonly used, especially for planted stands [2,3,5,6].The SI is defined as the mean height of dominant and co-dominant trees in a stand at a specified reference age [5].The dominant height of a stand reflects the productivity of a fully stocked, even-aged stand, because the height growth of dominant trees is independent of stand density over a wide range of densities [1,9,10].Conceptually, the SI was developed for even-aged stands of a single species, but over the last two decades, it has been applied to multi-aged and mixed species stands as well [11][12][13].As a result, numerous drawbacks of using the SI have been identified, discussed, and reported [11].For instance, the SI is often not observable, because even-aged, well-stocked, free-growing, undisturbed, and undamaged dominant and co-dominant trees may not exist.This is a common situation in the degraded stands or stands managed using uneven-aged silviculture.Even if the trees are present, they may not be desirable species for estimating the SI.Small errors in the measurements of age or dominant height often lead to large uncertainties in SI prediction [9].
These limitations of estimating the SI have become much more noticeable in recent years as the scope of forest management has widened towards the creation of multi-purpose stands with multiple ages and species [1,2,10,12,14,15].Such a new paradigm of management has not only limited the application of the SI, but also led to the exploration of alternative approaches to evaluate site productivity of such stands [1,12,14], such as using either site and environmental factors [1,5,7,8,14] or stand-based measures: stand basal area, stand volume [10,14], and height-diameter models [12,16].Worldwide, there is an increasing interest in using different strategies to manage multi-aged and mixed species stands where other stand measures, rather than the SI, are often used to estimate site productivity [10,17].
In a natural forestry context, site productivity emphasizes the potential of timber or biomass production as a main indicator of a site.It is also used as a reliable indicator of sustainability and predicting the rate of change of the biomass or volume production [10,14].The maximum mean annual basal area increment (BAI) is often considered to be a more useful measure of site productivity than a stand height-based index [10,14].An index that represents stand volume or biomass production allows for more direct and meaningful comparison across the tree species, site types, and regions than a stand height-based index [10].Vanclay [6] proposed a growth index for complex and mixed tropical forests based on diameter increment, which was adjusted for tree size (diameter) and competition measures.This study also suggested that the annual BAI or volume increment might indicate the site productivity, especially for unmanaged stands.Berrill and O'Hara [10] proposed BAI-or volume increment-based site productivity indices, and demonstrated the applications of these indices in the individual-tree and stand growth models for a coastal redwood (Sequoia sempervirens Lamb.ex D. Don Endl.).
According to the Chinese National Forest Inventory program, China has approximately 122 million hectares of natural forests, which accounts for 64% of the total forested area; the standing volume of natural forests is approximately 12.3 billion m 3 , or 83 % of the total volume in the country [18].However, decreased forest productivity caused by inappropriate considerations of site quality during forest management is a very critical issue in China.Stand height-based or site and environmental factors-based approaches of evaluating site productivity of plantations were well documented [19][20][21], however, no general approaches that are applicable for multi-aged and mixed species stands are available in China.Thus, this study aims to: (i) develop a new productivity index (i.e., BAI-based index for estimating site productivity of multi-aged and mixed species stands); (ii) to provide computational details of this index; (iii) to verify and validate the index using the national forest inventory (NFI) data from two typical multi-aged and mixed species stands in northeast China; and (iv) to demonstrate the application of the new index by evaluating site productivity of the studied forests.

Study Area and Data
We used data acquired from a total of 1912 permanent sample plots (PSPs) of the Chinese NFI on two main multi-aged and mixed species stand types, including Mongolian oak (Quercusmongolica Fisch.ex Ledeb.)-dominatedstands (stand type I) and other broadleaf stands (stand type II) with low growth rates (the dominant species were Tilia tuan Szyszyl., Fraxinus mandshurica Rupr., Ulmuspumila L., and Betula costata Trautv.).Both forest types were distributed randomly across the middle and eastern parts of Jilin Province in northeast China (Figure 1).Each PSP was square with an area of 0.06 ha.Of the 1912 PSPs, 859 were allocated in stand type I and 1053 in stand type II.
Forests 2017, 8, 119 3 of 18 and (iv) to demonstrate the application of the new index by evaluating site productivity of the studied forests.

Study Area and Data
We used data acquired from a total of 1912 permanent sample plots (PSPs) of the Chinese NFI on two main multi-aged and mixed species stand types, including Mongolian oak (Quercusmongolica Fisch.ex Ledeb.)-dominatedstands (stand type I) and other broadleaf stands (stand type II) with low growth rates (the dominant species were Tilia tuan Szyszyl., Fraxinus mandshurica Rupr., Ulmuspumila L., and Betula costata Trautv.).Both forest types were distributed randomly across the middle and eastern parts of Jilin Province in northeast China (Figure 1).Each PSP was square with an area of 0.06 ha.Of the 1912 PSPs, 859 were allocated in stand type I and 1053 in stand type II.The PSPs represented a variety of different stand structures, stand densities, site qualities, tree sizes, ages, and growth conditions.The inventory was conducted partially or fully in 1994, 1999, 2004, and 2009, meaning that most sampled trees had three 5-year growth intervals.Within each PSP, we measured the diameter at breast height (dbh) of all standing living trees with dbh ≥ 5 cm.Six average trees were selected in each PSP, following the same sampling intensity as would be used for dominant height selection (i.e., proportion of 100 thickest trees per ha) [22].The total height of each selected tree was measured.The arithmetic mean of the trees was regarded as a stand mean height (H).In addition, the age of each selected tree was recorded by counting the growth rings on increment cores taken at 0.1 m above the ground [23], and their arithmetic mean value was also used as a stand age (T).The PSPs represented a variety of different stand structures, stand densities, site qualities, tree sizes, ages, and growth conditions.The inventory was conducted partially or fully in 1994, 1999, 2004, and 2009, meaning that most sampled trees had three 5-year growth intervals.Within each PSP, we measured the diameter at breast height (dbh) of all standing living trees with dbh ≥ 5 cm.Six average trees were selected in each PSP, following the same sampling intensity as would be used for dominant height selection (i.e., proportion of 100 thickest trees per ha) [22].The total height of each selected tree was measured.The arithmetic mean of the trees was regarded as a stand mean height (H).In addition, the age of each selected tree was recorded by counting the growth rings on increment cores taken at 0.1 m above the ground [23], and their arithmetic mean value was also used as a stand age (T).
A total of 5575 observations (2718 and 2857 from stand types I and II, respectively) were used to estimate site productivity.H, T, stand basal area (BA), species composition, and Reineke's stand density index (S) [24] for each PSP were calculated.The S was given by where β is a parameter which, as proposed by Tang et al. [25], we set to 1.605; D is the standard base diameter (recommended value 20 cm); and N D are the number of standing living trees and the arithmetic mean diameter at breast height of the stand, respectively.Summary statistics of measurements of the relevant stand characteristics are presented in Table 1.All PSPs were stratified into ten site groups using seven biogeoclimatic factors, including elevation, slope, aspect, slope position, soil depth, soil type, and humus layer, and we coded them by group = 1, 2, ..., 10 with smaller values indicating higher site productivities.Methodological details of dividing the PSPs into site groups are also found in Zhu et al. [26].The relationship between H and T for each of the ten site groups and two stand types is shown in Figure 2. A total of 5575 observations (2718 and 2857 from stand types I and II, respectively) were used to estimate site productivity.H, T, stand basal area ( BA ), species composition, and Reineke's stand density index ( S ) [24] for each PSP were calculated.The S was given by where  is a parameter which, as proposed by Tang et al. [25], we set to 1.605; D  is the standard base diameter (recommended value 20 cm); and N D are the number of standing living trees and the arithmetic mean diameter at breast height of the stand, respectively.Summary statistics of measurements of the relevant stand characteristics are presented in Table 1.All PSPs were stratified into ten site groups using seven biogeoclimatic factors, including elevation, slope, aspect, slope position, soil depth, soil type, and humus layer, and we coded them by group 1, 2, ...,10  with smaller values indicating higher site productivities.Methodological details of dividing the PSPs into site groups are also found in Zhu et al. [26].The relationship between H and T for each of the ten site groups and two stand types is shown in Figure 2.   Forests 2017, 8, 119 5 of 18

Basal Area Potential Productivity Index
In this study, a basal area potential productivity index (BAPP) was defined as the maximum annual BAI at a specific T for a stand type with a given site type.It can be obtained by seeking an optional Reineke's stand density index [24], ( S), in a given feasible region, that maximizes the following objective function: BAPP = arg max where S min and S max are the lower and upper limits of S, respectively; BAI T is the annual basal area increment at T; and f is the BAI-function with a predictor S, which depends on stand height and basal area growth models.
Based on Equation (1), BAPP has two primary characteristics: (i) Stability: BAPP is dependent on site type, stand type, and T. Thus, when these variables are quantified, BAPP can be determined.(ii) Maximization: BAPP reflects the maximum mean annual BAI of a particular site type, stand type, and T; therefore, the mean annual BAI of a stand is always less than or equal to BAPP.However, in theory, the mean annual BAI can be close, or even equal to BAPP for an appropriately managed stand.
Thus, BAPP can be effectively used to quantify site productivity for multi-aged and mixed stands.Therefore, a considerable amount of reliable information (e.g., site quality and potential for improving the BAI of an objective stand) could be provided by the values of difference between the basal area realized productivity (BARP) and BAPP, which is very important for decision-making in forestry.
Computational Conditions for BAPP Computational conditions for BAPP at a specific T based on Equation (1) as follows: (i) Determine the objective stand type; (ii) Develop a H-T model (Equation ( 2)) for the objective stand, and estimate the parameters using ordinary nonlinear least square (ONLS) regression: where T and H are stand age and stand mean height, respectively, and ΦH is the vector of the parameter estimates.It should be noted that for estimating BAPP based on NFI data, H instead of stand dominant height was applied in the whole calculation.This is because (a) stand dominant height was not measured in Chinese NFI data, especially for multi-aged and mixed forests; and (b) the correlation between H and stand dominant height is strong [27].
(iii) Develop a BA growth model (Equation ( 3)) for the objective stand, and estimate the parameters using ONLS regression: where BA and S are stand basal area and stand density index, respectively, and ΦG is the vector of the parameter estimates.
(iv) Determine the site group of a stand type (group); (v) Determine the specific age of a stand type (T 0 ); (vi) Determine a feasible region of the stand density index.
When the abovementioned conditions are given, Equation ( 1) can be a typical nonlinear programming problem.Applying these conditions, we used the combined method of a golden section and a dichotomy iteration algorithm (GS&DIA) [28,29] to compute the BAPP in Equation (1).

The GS&DIA Procedure
We implemented the GS&DIA for estimating BAPP as follows: Step 1: Set a feasible region of S [S min , S max ], and set a site group (group), a base age (T 0 ), a tolerance error (e), and iteration t = 1.The starting values of the stand density indices for four points of segmentation (S 3 , and S (t) 4 ), as determined by the golden section, are given by: Step 2: The BAIs of four points of segmentation (BAI 3 , and 3 , and S (t) 4 , respectively, at the t th iteration are obtained using subprogram A (see below), and a dichotomy iteration algorithm is applied. Step , and S (t+1) 4 : and one must go back to step 2 for the next iteration, otherwise the calculation stops.The final estimate of BAPP at T 0 and its corresponding stand density index are:

, and S (t)
4 values from the main program).Identical tree growth was assumed for an objective stand type.That is, tree numbers in the stand in the two succeeding stand ages are identical.A dichotomy iteration algorithm was applied to determine the BAI as below (Figure 3): Step 1: Compute the basal area of the stand, BA 0 , and number of trees, N 0 , at T 0 ; (i) Given T 0 , compute H 0 ; (ii) Given T 0 , S 0 , and H 0 , compute BA 0 using the following equation: (iii) Given BA 0 and S 0 , N 0 and the mean diameter at breast height, D 0 , of the stand at T 0 are obtained by solving the following equations: The estimates of D 0 and N 0 are D0 = 40000BA 0 πS 0 (20) 1.605   200/79 and N0 = G 0 /(πD 2 0 /40000).
Step 2: Compute the basal area of the stand (BA 1 ) at T 1 (T 0 + 1) years; We suppose that the tree numbers in the two succeeding stand ages (T 0 and T 1 ) are identical, which means that N 1 = N0 .
(i) Given T 1 and ΦH , compute H 1 using the following equation: (ii) Given H 1 , T 1 , N 1 , and ΦG , the basal area of the stand (BA 1 ), the stand density index (S 1 ), and the mean diameter at breast height (D 1 ) at T 1 are obtained by solving the following equations: If we solve D 1 directly, after simplifying the equations, the objective function is changed to the following form: while solving S 1 directly, the objective function is: The dichotomy iteration algorithm was employed to solve the above Equations ( 9) or ( 10) to obtain D 1 or S 1 , and then we computed BA 1 .
Step 3: Compute BAI using the following equation: We present a framework with all the computational steps through which BAPP was derived (Figure 3).The codes for the GS&DIA implemented in ForStat 2.2 (a statistical software with specific forestry analytical tools developed by the Chinese Academy of Forestry, Beijing, China, see [30]) are given in the Supplementary Code S1.All other calculations were also performed using ForStat 2.2 software.

BARP
The BARP is defined as the difference between the basal area of a stand for the two succeeding stand ages (T 0 and T 0 + 1) with a 1-year interval.However, in practice, it is impossible to conduct field surveys every year, especially for large-scale surveys.To address this problem, we proposed a reliable method for calculating the BARP.Two main functions of the BARP are: (i) To verify computational method of the BAPP.The values of BAPP must be greater than or equal to those of the BARP, otherwise the values of BAPP would be meaningless.
(ii) To evaluate the realized stand growth status using the following equation: where ∆ is the potential for improving the BAI of the objective stand.The computational procedures of BARP are similar to those in subprogram A for BAPP.It should be noted that, unlike the BAPP calculation, the stand characteristics, including T and S from field surveys, should be provided for the BARP calculation.20) ( ) Solve the following equations ( , , , ) / 40000 ( / 20) to obtain using dichotomy iteration algorithm.
and | e A framework showing all important computational steps to derive the basal area potential productivity index (BAPP).

Stand Mean Height and Stand Basal Area Growth Models
Based on the work of Zhu et al. [26], we used the Richards' height growth model (Equation ( 12)): where H and T are sample plot-level arithmetic mean height (m) and stand age (year), respectively; ˆ( , , ) is a parameter vector of the th i site group; and  is an error term.
The Richards' BA growth model from an integrated stand growth model, as proposed by Tang et al. [31], was applied to estimate the site productivity of stand types I and II: where BA is the stand basal area (m 2 ha −1 ); S is the stand density index; H and T are the same as in Equation ( 12); and ˆ( , , , ) ) is the parameter vector, in which 1i  is related to the site group.We estimated Equations ( 12) and ( 13) using ONLS and evaluated their fitting performance by coefficient of determination ( 2 R ):

Stand Mean Height and Stand Basal Area Growth Models
Based on the work of Zhu et al. [26], we used the Richards' height growth model (Equation ( 12)): where H and T are sample plot-level arithmetic mean height (m) and stand age (year), respectively; ΦHi = (a i , b i , c i ) T is a parameter vector of the i th site group; and ε is an error term.
The Richards' BA growth model from an integrated stand growth model, as proposed by Tang et al. [31], was applied to estimate the site productivity of stand types I and II: where BA is the stand basal area (m 2 ha −1 ); S is the stand density index; H and T are the same as in Equation (12); and ΦBAi = (β 1i , β 2 , β 3 , β 4 ) T (i = 1, ..., 10) is the parameter vector, in which β 1i is Forests 2017, 8, 119 9 of 18 related to the site group.We estimated Equations ( 12) and ( 13) using ONLS and evaluated their fitting performance by coefficient of determination (R 2 ): where y j and ŷj are the observed and predicted H or BA, respectively, for the j th observation, and y is the mean H or BA of the observations.

Parameter Estimation
The parameter estimates of Equations ( 12) and ( 13) fitted by OLS regression are listed in Table 2.All parameter estimates for both models were significant (p < 0.05).Based on the fit statistics, the R 2 of Equation ( 12) was 0.987 and 0.978 for stand types I and II, respectively.For Equation ( 13), R 2 for the same stands was 0.956 and 0.945, respectively.This shows that both models fitted the data adequately well.12) and ( 13) for stand types I and II; a, b, c are the parameters of Equation ( 12), and β 1 , β 2 , β 3 , β 4 are the parameters of Equation (13).

BAPP Calculation
BAPPs under ten site groups for stand types I and II were obtained using GS&DIA based on Equations ( 12) and ( 13) with the known parameters estimates (Table 2).Using stand type I as an example, we presented the computational details of BAPP.The application conditions for computing BAPP of a stand type I are listed below: (i) Equation ( 12) was used as a H growth model, and its parameter estimates are given in Table 2; (ii) Equation ( 13) was used as a BA growth model, and its parameter estimates are listed in Table 2; (iii) Ten site groups denoted by group = 1, ..., 10 were applied.(iv) The feasible range of T ranged from 10 to 150 years.(v) A 50-year base age (T 0 ) was used.(vi) The feasible region of S ranged from 30 to 3000.
The BAPP at T was obtained by seeking an optimal stand density index, S, from the given feasible range (30, 3000) that maximized the objective function (Equation (1)).For example, Figure 4 shows the relationship of annual BAI and S for the feasible region 30 ≤ S ≤ 3000 when T = 20 years (young), 50 years (middle-aged), 100 years (sub-mature), and 150 years (mature) for the first site group (group = 1).For each specific T (T = 20, 50, 100, and 150 years), the relationship of BAI and S exhibited a single curve; thus, the maximum annual BAI (culmination point) on this curve is a BAPP at T.
(iv) The feasible range of T ranged from 10 to 150 years.(v) A 50-year base age ( 0 T ) was used.
(vi) The feasible region of S ranged from 30 to 3000.
The BAPP at T was obtained by seeking an optimal stand density index, S  , from the given feasible range (30, 3000) that maximized the objective function (Equation (1)).For example, Figure 4 shows the relationship of annual BAI and S for the feasible region 30 3000 S   when T = 20 years (young), 50 years (middle-aged), 100 years (sub-mature), and 150 years (mature) for the first site group ( group = 1 ).For each specific T ( T = 20, 50, 100, and 150 years), the relationship of BAI and S exhibited a single curve; thus, the maximum annual BAI (culmination point) on this curve is a BAPP at T .The values of S for the culmination points decreased as T increased (Figure 4), meaning that S  was negatively correlated with T .In addition, BAPP and T were also negatively correlated.The connection of all culmination points at different values of T within its feasible region (e.g., 10 150 T   ) formed a three-dimensional curve of BAPP (the green line in Figure 4).This curve depicts the correlation between BAPP and S  at different values of T for the multi-aged oak- dominated mixed stands (stand type I).Three-dimensional curves of BAPP could be projected as twodimensional ones based on T , which would be simpler.This simplified curve is defined as a BAPP curve for stand type I (the blue line in Figure 4).For different site groups in each stand type, there were negative correlations between BAPP and T .For a specific T , BAPP and site group ( group=10 ) were also negatively correlated (i.e., BAPP decreased with increasing site group (Figure 5)).These results are consistent with the definition of The values of S for the culmination points decreased as T increased (Figure 4), meaning that S was negatively correlated with T. In addition, BAPP and T were also negatively correlated.The connection of all culmination points at different values of T within its feasible region (e.g., 10 ≤ T ≤ 150) formed a three-dimensional curve of BAPP (the green line in Figure 4).This curve depicts the correlation between BAPP and S at different values of T for the multi-aged oak-dominated mixed stands (stand type I).Three-dimensional curves of BAPP could be projected as two-dimensional ones based on T, which would be simpler.This simplified curve is defined as a BAPP curve for stand type I (the blue line in Figure 4).
For different site groups in each stand type, there were negative correlations between BAPP and T. For a specific T, BAPP and site group (group = 10) were also negatively correlated (i.e., BAPP decreased with increasing site group (Figure 5)).These results are consistent with the definition of the site groups, where the first site group denotes the highest site quality and the last site group denotes the lowest site quality of the stands.the site groups, where the first site group denotes the highest site quality and the last site group denotes the lowest site quality of the stands.

BAPP Verification
The BARPs of all the PSPs for both stand types were obtained using Equations ( 12) and ( 13) according to the computational method of BARP presented in Section 2.3.Figures 6 and 7 show the BAPPs and BARPs of all PSPs for stand types I and II, respectively.The BAPP value of each PSP was larger than or equal to the value of the corresponding BARP for each site group.
For young and middle-aged stands, the differences of BARPs from BAPPs for stand types I and II were large.However, such difference decreased with increasing T .For each site group, the BARP values were much smaller than those of the BAPP, and BARPs tended to be closer to the corresponding BAPPs when a stand was getting more and more matured stages.In addition, for a specific T , BARP was always smaller than BAPP.However, when stands had higher site productivities and were managed more appropriately, both BARPs and BAPPs were almost identical (Figure 6 and 7).These results suggest that BAPP could be used to assess site productivity effectively.

BAPP Verification
The BARPs of all the PSPs for both stand types were obtained using Equations ( 12) and ( 13) according to the computational method of BARP presented in Section 2.3.Figures 6 and 7 show the BAPPs and BARPs of all PSPs for stand types I and II, respectively.The BAPP value of each PSP was larger than or equal to the value of the corresponding BARP for each site group.

Application of BAPP in Forest Management
The BAPP is used to evaluate site productivity effectively.Based on the values of difference of the BARP from BAPP, three important questions for forest management could be answered during decision-making in forestry.As an example, we illustrated its application in the forest stands of Wangqing Forestry Bureau.
(i) How much could BAI of a stand be improved?
The potential improvement of the BAI of a stand at time T for each site group is given by For young and middle-aged stands, the differences of BARPs from BAPPs for stand types I and II were large.However, such difference decreased with increasing T. For each site group, the BARP values were much smaller than those of the BAPP, and BARPs tended to be closer to the corresponding BAPPs when a stand was getting more and more matured stages.In addition, for a specific T, BARP was always smaller than BAPP.However, when stands had higher site productivities and were managed more appropriately, both BARPs and BAPPs were almost identical (Figures 6 and 7).These results suggest that BAPP could be used to assess site productivity effectively.

Application of BAPP in Forest Management
The BAPP is used to evaluate site productivity effectively.Based on the values of difference of the BARP from BAPP, three important questions for forest management could be answered during decision-making in forestry.As an example, we illustrated its application in the forest stands of Wangqing Forestry Bureau.
(i) How much could BAI of a stand be improved?
The potential improvement of the BAI of a stand at time T for each site group is given by where BAI T is the potential basal area increment at time T, and BAPP T and BARP T are the basal area potential productivity and realized productivity, respectively, at time T. Using stand type I as an example, ∆BAI T at different values of T (T = 10, ..., 150) for the ten site groups (group = 1, ..., 10) are displayed in Figure 8.The ∆BAI T (m 2 ha −1 ) values at the base age (T = 50 years) for the ten site groups (group = 1, ..., 10) were 0.31, 0.22, 0.34, 0.26, 0.29, 0.32, 0.28, 0.24, 0.18, and 0.19, respectively.Similar results could be attained at forest management level-Wangqing Forestry Bureau (Figure 9).area potential productivity and realized productivity, respectively, at time T .Using stand type I as an example, Similar results could be attained at forest management level-Wangqing Forestry Bureau (Figure 9).(ii) When should the site quality of a stand be improved?
The optimal time for improving the site quality of a forest stand (e.g., via new or improved forest management treatments) could be determined with the help of the curves of the relationship between T BAI  and T (Figure 8).This figure shows that .This clearly indicates that any tending operations should be applied to a stand before it becomes 50 years old.The stands could then be substantially improved when appropriate measures were applied in an appropriate range of stand ages.
(iii) What stand is given priority for improvement?
In forestry, stands with the best or worst site conditions are usually selected for improving the productivity.However, potential improvements of these stands are commonly ignored.In this study, the values of BAI  could be used to determine which stands should have priority for quality improvement.For example, stands with large values of BAI  might have small BARP and, thus, they would represent high-priority stands for quality improvement.In contrast, stands with smaller values of BAI  would have lower priority under proper management.Using stand type I as an example, according to the BAI  (m 2 ha −1 ), the priorities of stand improvements at the base age ( (ii) When should the site quality of a stand be improved?
The optimal time for improving the site quality of a forest stand (e.g., via new or improved forest management treatments) could be determined with the help of the curves of the relationship between ∆BAI T and T (Figure 8).This figure shows that ∆BAI T decreased with increasing T for each site group.However, ∆BAI T first decreased rapidly before T = 50 and then relatively slowly after T = 50.This clearly indicates that any tending operations should be applied to a stand before it becomes 50 years old.The stands could then be substantially improved when appropriate measures were applied in an appropriate range of stand ages.
(iii) What stand is given priority for improvement?
In forestry, stands with the best or worst site conditions are usually selected for improving the productivity.However, potential improvements of these stands are commonly ignored.In this study, the values of ∆BAI could be used to determine which stands should have priority for quality improvement.For example, stands with large values of ∆BAI might have small BARP and, thus, they would represent high-priority stands for quality improvement.In contrast, stands with smaller values of ∆BAI would have lower priority under proper management.Using stand type I as an example, according to the ∆BAI (m 2 ha −1 ), the priorities of stand improvements at the base age (T = 50 years) for the ten site groups are: group = 3 (∆BAI = 0.34) > group = 6 (∆BAI = 0.32) > group = 1 (∆BAI = 0.31) > group = 5 (∆BAI = 0.29) > group = 7 (∆BAI = 0.28) > group = 4 (∆BAI = 0.26) > group = 1 (∆BAI = 0.23) > group = 2 (∆BAI = 0.22) > group = 10 (∆BAI = 0.19) > group = 9 (∆BAI = 0.18).The values of ∆BAI can also be used to determine the priorities of the stands with different values of T.
Based on the abovementioned functions and provisions, BAPP could be applied for site productivity evaluation of stands at a large scale.Figure 9 shows the distribution of the potential improvement of BAI for stand type I at current stand age of each sub-compartment in the Wangqing Forestry Bureau, which is a main distribution region of this stand type in Jilin province based on the NFI data acquired in 2009.We could directly obtain the site quality class of any locations in the Wangqing Forestry Bureau in BA growth for stand type I in Figure 9.In addition, we could also know the potential improvement of the BAI for each sub-compartment with the management priority being their quality improvement.Therefore, the combination of the information given in Figure 9 and the abovementioned three questions could be useful for decision-making in forest management.group = 9 ( BAI  = 0.18).The values BAI  can also be used to determine the priorities of the stands with different values of T .
Based on the abovementioned functions and provisions, BAPP could be applied for site productivity evaluation of stands at a large scale.Figure 9 shows the distribution of the potential improvement of BAI for stand type I at current stand age of each sub-compartment in the Wangqing Forestry Bureau, which is a main distribution region of this stand type in Jilin province based on the NFI data acquired in 2009.We could directly obtain the site quality class of any locations in the Wangqing Forestry Bureau in BA growth for stand type I in Figure 9.In addition, we could also know the potential improvement of the BAI for each sub-compartment with the management priority being their quality improvement.Therefore, the combination of the information given in Figure 9 and the abovementioned three questions could be useful for decision-making in forest management.

Discussion
Estimating forest site productivity is necessary for effective forest management and useful for evaluating basic site conditions for ecological field studies [32].The high variability in silvicultural and ecological conditions for multi-aged and mixed species forests makes it difficult to use the dominant height-age index or the SI.For such conditions, as mentioned in the introduction section of this article, it is necessary to use alternative indices, such as basal area growth or diameter growth as productivity indices, and to integrate them into the stand structure or typology in the analyses [32].In this study, we first proposed a new BAI-based index (i.e., BAPP) for measuring the site productivity of multi-aged and mixed species forests.This index is the maximum mean annual BAI at a given T and site condition for a specific stand, and the computational details of BAPP are presented.The BARP of mixed species stands was then used to verify this index.In addition, the application of BAPP was also demonstrated using the NFI time-series data from two types of multiaged and mixed species forests in northeast China.The results showed that BAPP could be used to

Discussion
Estimating forest site productivity is necessary for effective forest management and useful for evaluating basic site conditions for ecological field studies [32].The high variability in silvicultural and ecological conditions for multi-aged and mixed species forests makes it difficult to use the dominant height-age index or the SI.For such conditions, as mentioned in the introduction section of this article, it is necessary to use alternative indices, such as basal area growth or diameter growth as productivity indices, and to integrate them into the stand structure or typology in the analyses [32].In this study, we first proposed a new BAI-based index (i.e., BAPP) for measuring the site productivity of multi-aged and mixed species forests.This index is the maximum mean annual BAI at a given T and site condition for a specific stand, and the computational details of BAPP are presented.The BARP of mixed species stands was then used to verify this index.In addition, the application of BAPP was also demonstrated using the NFI time-series data from two types of multi-aged and mixed species forests in northeast China.The results showed that BAPP could be used to assess forest site productivity effectively.Meanwhile, BAPP also provides much reliable information for decision-making in forestry.
The SI derived from stand dominant height-dominant age model is empirical and based on two assumptions: (i) height growth is highly dependent on site quality and independent of competition effects or stand density effects; and (ii) height growth over time is asymptotic, reaching a maximum that is defined by site quality [33].In practice, when dominant or co-dominant trees were selected as site trees to estimate the SI, an assumption was made that these trees have always been dominant or co-dominant and would continue to be so [34,35].However, this assumption might be difficult to meet because dominant trees observed at the time of sampling might have been suppressed before.The dominant height-age relationship approach for site productivity estimation of multi-aged and mixed species forests was based upon following three assumptions [12]: (i) decreasing tree taper (D divided by H) is associated with increasing site productivity [12]; (ii) stand density does not affect the height-diameter relationship of the dominant and co-dominant trees in the uneven-aged stands; and (iii) height growth over time is asymptotic, whereas diameter is not.The first and the third assumptions generally hold [12].However, the second assumption is very difficult to meet because stand density is a very important factor related to competition effects that is well correlated with stand growth.As per BAPP, there are not any restrictions or assumptions to be made on height growth.This is because BAPP is a local maximum of BAI in a stand feasible region of S at a given T and it is obtained from empirical stand growth models by optimizing an iterative algorithm in theory.In addition, the effects of stand density on BAI also were included in the BAPP calculation.Therefore, relative to SI, BAPP is a more reliable and effective index for site quality evaluation of multi-aged and mixed species forests.Furthermore, BAI is a reliable indicator of sustainability of biomass or volume productions in mixed species stands [6,10,14].
In this study, both H and BA growth models were needed to determine BAPP.For a H growth model, T is the only predictor.This is useful, because adding other tree and stand predictors would increase forest inventory costs and the computational complexity of BAPP.Even though the Richards' height growth model (Equation ( 12)) showed high fitting accuracy for stand types I and II, other model forms, such as exponential, logistics, Schumacher, and power functions, could also be used to develop height growth models [9,12,[32][33][34][35][36][37].The predictors in the BA growth model were restricted to T, S, and H when computing BAPP.In addition, the functional relationship of the BAI obtained from the BA growth model and S should be non-monotonic.The BA growth model (Equation ( 13)) in this study had the above-referenced characteristics.Similar to the results of the integrated stand growth model [31], this model (Equation ( 13)) also shows good performance for BA prediction.
Stand density index, S is negatively correlated with T (Figure 4), and this is consistent with the laws of stand growth.In forestry, the mean diameter at breast height of a stand increases and stand density decreases during self-thinning as T increases, which also causes S to decrease.Self-thinning and tree growth are the two common phenomena of the stands, and self-thinning is a natural process in which the numbers of trees per unit area decrease as the mean size of the trees increases over time [25,38].It is a process intrinsic not only to the forest types used in this study, but also to all forest types and plant communities whose compositions and structures are influenced by competitive interactions among the individuals [25,38].Stand attributes, such as mean height and mean diameter at breast height, are extremely difficult to measure when self-thinning occurs, which is more obvious in mixed species stands [25,38].Therefore, in this study, the annual BAI of a stand was obtained based on the assumption that the density of the stand would remain stationary.In reality, the stands selected for estimating BAPP are relatively stable at a chosen base age.The self-thinning phenomenon is also not obvious for these stands, and therefore, the assumption of stand growth for the same trees is feasible for estimating BAPP.An additional self-thinning algorithm of stands needs to be developed if a modeler wants to include self-thinning in the BAPP computation.We are in the process of developing such algorithms.
The differences of BAPPs among the ten site groups were quantified by parametrizing the H and BA growth models (Figures 4-6).The site quality of a stand decreased with increasing site groups (from group = 1 to group = 10).For both stand types, the values of BAPP at a specific T decreased with increasing site groups, which further suggests that the BAPP proposed in this study could be suitable for effectively assessing site productivity (Figure 5).In addition, at present, because of the state of abandonment and decay of the stands, it is necessary to establish a priority area where the resources should be allocated to improve the stands [39], and BAPP may be a key tool for this purpose.The utility of any approach for estimating the site productivity of a multi-aged and mixed species stand is largely dependent on its application conditions.In reality, the application conditions for BAPP are easily met for different stands.The BAPP could provide a quick and simple approach for quantifying site productivity for multi-aged and mixed stands.The main drawback of this approach is that the accuracy of the site productivity depends largely on the prediction accuracies of the H and BA growth models, as these models cause an uncertainty regarding the computation of BAPP.However, this should not be a major issue because of the rapid development of software and techniques that can be easily applied for modeling forest growth and site productivity.In addition, since BAPP was defined as the maximum annual BAI with the optimal stand density, it could guide silvicultural practices to adjust the stand density close to the optimal one.This needs to be further studied in the future.
It is worth noting that, because of the uncertainty in climate change (e.g., extreme climate), the predictions using empirical H and BA growth models (e.g., Equations ( 12) and ( 13)) should not be made very far ahead of the current conditions [40].For making a short-term site quality evaluation (such as ten years or a shorter time interval), Equations ( 12) and ( 13) with the estimates of parameters could be directly applied because the structure of multi-aged and mixed forests and the site-related variables would not be obviously changed [41].However, for making a long-term site quality evaluation, besides the growth of dominant species in the multi-aged and mixed forests spatially affected by climate change, the distribution of habitats and the structure of the multi-aged and mixed forests could also substantially change [42,43].Therefore, when Equations ( 12) and ( 13) are used, sampling surveys for predicting H and BA of each multi-aged and mixed forest in a study area should be conducted spatially and periodically, and used to recalibrate and validate the existing models.

Figure 1 .
Figure 1.Location of the permanent sample plots (PSPs) in two main multi-aged and mixed species stands (859 PSPs in oak-dominated stands (stand type I) and 1053 PSPs in other broadleaf stands (stand type II)) distributed across the middle and eastern parts of Jilin Province in northeast China.

Figure 1 .
Figure 1.Location of the permanent sample plots (PSPs) in two main multi-aged and mixed species stands (859 PSPs in oak-dominated stands (stand type I) and 1053 PSPs in other broadleaf stands (stand type II)) distributed across the middle and eastern parts of Jilin Province in northeast China.

Figure 2 .
Figure 2. The relationships of stand mean height and stand age for ten site groups (group = 1, 2, …,10) for oak-dominated stands (stand type I) and other broadleaf stands (stand type II).

Figure 2 .
Figure 2. The relationships of stand mean height and stand age for ten site groups (group = 1, 2, . . ., 10) for oak-dominated stands (stand type I) and other broadleaf stands (stand type II).

Figure 3 .
Figure 3.A framework showing all important computational steps to derive the basal area potential productivity index (BAPP).

Figure 4 .
Figure 4.The distribution of the annual basal area increment ( BAI ) and stand density index ( S ) for the first site group ( group = 1 ) at the stand ages of T = 20, 50, 100, and 150 years for oak-dominated stands.

Figure 4 .
Figure 4.The distribution of the annual basal area increment (BAI) and stand density index (S) for the first site group (group = 1) at the stand ages of T = 20, 50, 100, and 150 years for oak-dominated stands.

Figure 5 .
Figure 5.The clusters of stand basal area potential productivity (BAPP) curves with group = 1,...,10 for oak-dominated stands (stand type I) and other broadleaf stands (stand type II).

Figure 5 .
Figure 5.The clusters of stand basal area potential productivity (BAPP) curves with group = 1, ..., 10 for oak-dominated stands (stand type I) and other broadleaf stands (stand type II).

Forests 2017, 8 , 119 12 of 18 Figure 6 .
Figure 6.The basal area potential productivity (BAPP) and basal area realized productivity (BARP) of each permanent sample plot against stand age T for group = 1,...,10 for stand type I. BAI is the annual basal area increment.

Figure 6 .
Figure 6.The basal area potential productivity (BAPP) and basal area realized productivity (BARP) of each permanent sample plot against stand age T for group = 1, ..., 10 for stand type I. BAI is the annual basal area increment.

Figure 6 .
Figure 6.The basal area potential productivity (BAPP) and basal area realized productivity (BARP) of each permanent sample plot against stand age T for group = 1,...,10 for stand type I. BAI is the annual basal area increment.

Figure 7 .
Figure 7.The basal area potential productivity (BAPP) and basal area realized productivity (BARP) of each permanent sample plot against stand age T for group = 1,...,10 for stand type II.BAI is the annual basal area increment.

Figure 7 .
Figure 7.The basal area potential productivity (BAPP) and basal area realized productivity (BARP) of each permanent sample plot against stand age T for group 1, ..., 10 for stand type II.BAI is the annual basal area increment.

Figure 8 .
Figure 8.The potential improvement of the annual basal area increment ( BAI ) for each permanent sample plot of oak-dominated stands for the ten site groups ( group = 1,...,10 ), where T is the stand age.

Figure 8 .
Figure 8.The potential improvement of the annual basal area increment (∆BAI) for each permanent sample plot of oak-dominated stands for the ten site groups (group = 1, ..., 10), where T is the stand age.

Figure 9 .
Figure 9.The distribution of the potential improvement of the basal area increment of stand type I at the current stand age of each sub-compartment in the Wangqing Forestry Bureau, which is a main distribution region of this stand type in Jilin province based on the national forest inventory data acquired in 2009.

Figure 9 .
Figure 9.The distribution of the potential improvement of the basal area increment of stand type I at the current stand age of each sub-compartment in the Wangqing Forestry Bureau, which is a main distribution region of this stand type in Jilin province based on the national forest inventory data acquired in 2009.

Table 1 .
Summary statistics of the stand variables for two main multi-aged and mixed species stand

Table 1 .
Summary statistics of the stand variables for two main multi-aged and mixed species stand