Spatial Forest Harvest Scheduling for Areas Involving Carbon and Timber Management Goals

Forest carbon sequestration has become an important ecological service for human society. Given the widespread attention paid to global climate change over the last few decades, a potential need has arisen to develop forest management plans that integrate carbon management and other spatial and non-spatial goals. The objective of this research was to develop a spatial forest planning process by which one could assess either a carbon stocks objective, a timber production objective, or a spatial objective related to the arrangement of forest management activities. This process was used to evaluate the maximization of (1) volume scheduled for harvest; (2) carbon stocks; and (3) spatial aggregation of the management activities through a utility function where all are equally weighted objectives. The process was employed for the development of 30-year plans for a forested landscape in northeast China that was approximately 120,000 ha in size. In addition, the sensitivity of the results with respect to four initial forest age structures was tested. Constraints mainly included those related to the need for an even flow of scheduled harvest volume and to the need to adhere to a maximum harvest opening size. The proposed scheduling process employed a simulated annealing algorithm to schedule harvests in an OPEN ACCESS Forests 2015, 6 1363 attempt to produce a high value of the utility function. Results showed that carbon stocks in the case study forests could significantly increase in the next 30 years under the proposed harvesting plans. Of the case study forest landscapes, the values of both the utility function and the computing time required were significantly different between different initial forest age structures (p < 0.05), i.e., the older forest landscape obtained the highest average solution value (0.6594 ± 0.0013) with the fastest processing speed (2.45 min per solution). For a fixed harvest level, the average carbon density (tons per hectare) at the end of planning horizon also increased by 4.48 ± 9.61 t/ha, 8.73 ± 10.85 t/ha, 2.99 ± 9.19 t/ha and 1.03 ± 9.77 t/ha when maximizing the total utility functions for the actual, young, normal and older forests, respectively, when compared those at their initial conditions. This heuristic spatial forest planning process can allow forest managers to examine a number of different management activities, for both timber production and carbon stocks, prior to selecting a preferred alternative.


Introduction
Forests play an important role in society's efforts to mitigate greenhouse gas emissions.Under the guidance of the United Nations Framework Convention on Climate Change Kyoto Protocol [1], many researchers have concentrated their attention on different scales of carbon inventories [2][3][4], models of carbon stocks by forest and carbon production [5][6][7], simulations of the effects of forest management [8][9][10], and analyses of the effects of climate change on carbon sequestration by forest types [11,12].However, more information is needed regarding greenhouse gas mitigation options, and in this context there is considerable interest in using forests for carbon sequestration efforts, and therefore in assessing how forest management can contribute to enhancing carbon sequestration efforts.
Carbon sequestration in forest ecosystems involves numerous components that include above-ground biomass carbon and below-ground soil carbon [13].Given the importance of forest carbon sequestration, it comes as no surprise that efforts have been put forth to value the potential additionality provided through changes in forest management practices [14].Several studies have examined the effects of varying rotations lengths [15,16], thinning intensities [8,10], and green-up and maximum opening size constraints [17][18][19] on the amount of timber produced and carbon sequestered (stocked) in forests.What all these studies have in common is that suitable forest management activities can enhance forest carbon sequestration efforts and timber production over the long-term.Therefore, it is valuable to determine which forest management activities (combinations of cutting regimes (if used), silvilcultural treatments, assumptions on the placement of activities, and management intensities) are more efficient and effective for these purposes for special management areas, and when to implement the selected activities in order to better reach timber production and carbon sequestration (or other) goals in the future.
Forests not only store carbon, they provide a wide range of economic, ecological and social value.Therefore, understanding ecosystem service interactions represents a challenge in land management decision-making [20].Generally speaking, the interaction can either be viewed as a trade-off or a synergy.As Cademus et al. [21] described, a trade-off occurs when production of one ecosystem service compromises the provision of other services.Conversely, a synergy occurs when the productions of multiple ecosystem services either increase or are provided at similar levels [21].However, a large-scale experiment that physically examines these types of forest management goals would be too costly and time-consuming to implement across a broad area [22].Therefore, a number of mathematical and heuristic methods have been used by forestry organizations over the past fifty years in an effort to support decision-making processes related to complex and interrelated forest management objectives [22], e.g., those that involve economic and commodity production, wildlife habitat, forest structure, biodiversity and other objectives.To our knowledge, some previous research efforts have examined the tradeoffs (or synergies) between commodity production goals and carbon stocks [23][24][25][26].These analyses generally indicated how changes in scheduled timber harvest levels could affect a land management organization's ability to address carbon-related goals (whether measured as a carbon stock or a carbon sequestration rate).Commodity production goals from a forest and the ability of a forest to store carbon are two areas of emphasis for land managers that are closely associated; therefore, depending on the management option employed, the outcomes of a forest plan can result in either positive or negative achievement of either area of emphasis.
However, most previous studies lack the spatial or temporal detail necessary to adequately evaluate the changes of forest landscapes and the spatial distribution of management activities.Nowadays, forest planning problems are becoming more complex because of the need to integrate spatial relationships in long-term planning process, to address certification requirements, aesthetic concerns, or issues related to fragmentation of older forests.One reason for this is that a spatial model for scheduling harvests requires characterization not only of the state of each stand, but also of the state of its neighbors, implying that nonlinear relationships between decision variables or integer variables will have to be used [27].While exact mathematical methods (i.e., mixed integer programming) can be used to address these types of problems, heuristic techniques are increasingly being employed in forest planning efforts as an alternative to produce near-optimal solutions with much shorter computational times.A variety of heuristic techniques can be used to address these problems, such as Monte Carlo integer programming (MCIP) [28], simulated annealing (SA) [27], tabu search (TS) [29], genetic algorithms (GA) [30], threshold accepting (TA) [31] and other hybrid algorithms [20,32].However, the performance of heuristics may vary depending on the complexity of the planning problem and rationality of the parameters.Previous work in this area has illustrated that a well-designed SA algorithm can produce forest plans as good (or better) than others heuristics for forest-level spatial planning problems with different scales [32,33].Bettinger et al. [32], for example, noted that SA was as good at developing forest plans as TA, great deluge (GD), TS with 1-opt and 2-opt moves and TS combined with GA, but was significantly better than TS with 1-opt moves (basic tabu search), GA and random search, when three increasingly difficult timber and wildlife planning problems were addressed.
The aim of this study was to develop and demonstrate a method which allows one to examine the synergies between timber production goals and carbon stocks goals.A utility function that includes timber production, carbon stocks and spatial considerations was developed as the objective.The process was subjected to an even flow of scheduled harvest volume and maximum size limit constraints.We described the results of the forest plans developed in terms of (a) scheduled timber harvest volume; (b) carbon stocks (standing, live trees) obtained across the forest; and (c) carbon density (amount of carbon per unit area) of each management unit.We hypothesized that for certain management options (combinations of cutting regimes (if used), silvilcultural treatments, assumptions on the placement of activities), land managers can achieve desired timber production goals as well as desired carbon management goals.To investigate how well the heuristic performs in developing a management plan in a realistic setting, it is applied to a national forest area of 123,293 ha in northeast China.We developed four sets of forest landscapes (one actual, and three hypothetical age class distributions with young, normal and older forests), and subsequently estimated the spatial dynamics of management activities, timber harvest levels and carbon stock levels associated with these age class distributions.

Materials
A number of structural and management variables can have significant influence in the forest management decision making process.However, when compared with other approaches, the use of age classes has a more explicit and credible application to forest management planning practices some parts of the world [34,35] (17,207 ha), 125 compartments were considered coniferous forests (44,927 ha), 59 compartments were considered broad-leaved forests (20,522 ha), 4 compartments were considered mixed broadleaf-conifer forests (1614 ha), and 8 compartments were represented by non-forest land (2737 ha).
In addition to the actual age class distribution (Figure 1), three other age class distributions were randomly assigned to the landscape database (young, normal and older forest), to assess effects on management strategies to maximize scheduled harvest volume or carbon stocks.The young forest represented an age class distribution of a typical forest in northeast China, where past natural disturbances and forest management practices have left many young, over-cut forests.The normal age class distribution was devised simply because many forest management organizations have this as a goal for their long-range forest plans.The older forest age class distribution was developed to represent an approximate, yet realistic inverse of the young forest age class distribution.All four models of the age class distributions of the forest landscapes used for this analysis have similar total areas, polygon sizes and others stand characteristics.

Methods
This research describes a forest planning process that integrates carbon management objectives into a wood production modeling environment that is subject to spatial and even-flow timber constraints.The methods consist of three distinct topical areas related to the planning process: (1) the growth and yield models of the five forest types; and (2) the forest planning problem formulation; and (3) the evaluation of the different age class distributions.

Forest Growth and Yield Models
The forest growth and yield models that the scheduling process used were mostly developed by the Department of Forest Management, Northeast Forestry University (NEFU-FM), and were based on data from the China National Forest Inventory (NFI) of the study area.These models include: average tree height (TH) growth models, average diameter at the breast height (DBH) growth models, site class index (SCI) models, stand density index (SDI) models, SDI dynamics models, stand basal area (BAS) growth models, stand volume (VOL) growth models, stand biomass (BIO) models and stand carbon (CAR) models (see details in [36,37]).Using these stand-and tree-level models, one can predict the forest growth and yield and forest response to different selective cutting intensities.

Forest Planning Formulation
The forest plans developed cover a 30-year time horizon that is divided into three 10-year planning periods.The objective function of the planning problem is one that seeks to maximize a utility function value that can include timber production, carbon stocks and the spatial pattern of management activities.

Objective Function
As we noted, the objective function for the planning problem involved simultaneously maximizing timber volume, carbon stocks and spatial aggregation of the management activities.The objective function was described as an additive utility model [38]: where U is utility function value, I is the number of objectives, ui is a scaled sub-priority function for management objective i, qi is the amount of objective variable i, and ai is the relative importance of objective variable i, and all objectives were assumed to be equally important, i.e., all ai in Equation (1) were 0.3333 throughout this paper.The three goals used in all utility functions were timber production, carbon stocks and spatial aggregation of management activities.Timber production and carbon stocks are predicted using the forest growth and yield models mentioned earlier.The clustering index of management activities of a planning period was calculated from [39].
where FSV is the forest spatial value, N is the number of compartments in the study area, NVi is the value of adjacency relationship for the i-th stand; Lik and Dik are the border length and the distance between the i-th stand and its k-th adjacency stand, Rik is a binary variable indicating whether (1) or not (0) the i-th stand and its k-th adjacent stand have an identical management activity (both planning period and selective intensity).
The sub-utility function for timber production was represented by two symmetric line-segments, one that increased linearly to a maximum value of 1.0 when a user-specific scheduled volume was achieved, and another that decreased linearly to a minimum value of 0.0 when the volume was equal to the user-specific scheduled volume multiplied by 2 (Figure 2A).For carbon stocks, the sub-utility function increased linearly to a maximum value of 1.0 which represented a user-specific carbon stock goal value, after which the sub-utility function was replaced by 1.0 (Figure 2B).For the forest spatial value, the sub-utility function increased linearly to a maximum value of 1.0, which represented the user-specific goal value (Figure 2C).In this research, the same set of sub-utility functions were applied to the four different forest landscapes.

Constraints
There were three main constraints within the planning formulation to guide the development of a forest plan.Two were designed as a proxy for some aspect of the decision-making process, where ecological and economic goals are recognized.Due to both the behavior of industrial landowners and the desires of state landowners, an even flow of harvest volume seems to play an important role for forest sustainable management.Rather than using a strict even-flow rule, the harvest volume constraints limited the ranges of the scheduled volume during each time period.
where VOLt and VOLt − 1 are the scheduled harvest volumes in time periods t and t − 1, bl and bu are the lower and upper limits harvest volume for each period, bl must be less than 1.0, and bu must be greater than 1.0.In this research, we assume the bl = 0.8, and bu = 1.2, respectively.The second constraint defined a maximum harvest size limit for each compartment n and its set of adjacent neighbors (Un) with a same management activity within each time period.Therefore, a size constraint not only can affect the scheduled volumes during each time period, but also affect the arrangement of management activities around each compartment scheduled during time period t.We defined the size constraint as: where again Un is the set of adjacent neighbors to unit n, k is the adjacent neighbors to unit n, Ak is the area of unit k, xnk is binary variable indicating whether unit k has the same management activity (including the harvest time and intensity) with the unit n, An is the area of unit n, Amax is the user-specific value.
The third constraint was a singularity constraint, which limited each management unit to one activity during the planning horizon: where T is the number of total planning periods, and xnt is a binary variable indicating whether (1) or not (0) the management unit n has been assigned a management activity during time period t.

Management Activities
The forests in our study area play an important role in the ecological function and economic development of the country.However, the quality and ecological function of these forests have been seriously degraded due to the adoption of an extensive farming strategy within the last 100 years.In order to restore and protect forest resources, the Heilongjiang Province Law prohibited final felling of forests on 1 April 2014.Based on the stipulations of the Laws, in this research we will assume that the management activities can only include four different selective cutting intensities (i.e., 0%, 10%, 20% and 30%) and each compartment can be managed (entered) only once during the planning time horizon (30 years).These four selective cutting intensities, in conjunction with the three planning periods, provide twelve potential management activities for each compartment.Given that there are 352 compartments, the total number of decision variables is 4224.

Evaluation of Different Age Class Distributions
When using heuristic techniques, one cannot be certain that the global optimum solution to a planning problem will be found, nor that the resulting solutions are even close to the global optimum [29,32].Therefore, the optimization process was repeated 20 times for each forest landscape (actual, young, normal, older).The maximum and minimum objective function values, and standard deviations and averages of these were compared as a self-validation of the results [40].As Barrett et al. [41] and Bettinger et al. [42] described, if a heuristic search was initiated from a random starting point, each resulting solution generated from a heuristic can be considered an independent sample from a larger population, for which one can assume independence and perform statistical tests.If normally distributed, a Students t-test can be used to quantitatively detect whether there are significant differences in the forest plans developed for the four heterogeneous forest landscapes.All tests are conducted at the α = 0.05 level.

Heuristic Search Algorithm
A heuristic search process, simulated annealing (SA), was developed and employed to schedule activities for each stand.SA is inspired by the process of annealing (cooling) of metals, and the motivation for it was first described in 1953 [43].As a search process, SA has shown to be useful in a number of important forestry problems, such as scheduling harvest [32,44], maintaining or developing wildlife habitat [31], maintaining core areas [45] and optimizing forest landscape values [46].Each of these efforts has shown that although SA cannot guarantee an optimal solution, it should be able to provide a number of good, feasible solutions within an acceptable period of time for large and complex planning problems.
SA is technically considered a s-metaheuristic neighborhood search technique [47] that examines a solution to a forest plan, then proposes a change to a single characteristic (harvest time or activity).As implemented here, we utilize 1-opt moves, whereby a forest plan is incrementally developed by changing the status (assignment of harvest time or activity) of a single management unit.SA always accepts every new solution which improves the value of the objective function.Non-improving moves are accepted based on a probability that is determined through a function where Temperature is the current temperature of the annealing process, Unew and Uold are the value of objective function for the current and previous best solutions, respectively.During the optimization process, the temperature cools according to a given cooling rate.At high temperatures the probability for accepting inferior solutions is high, and as the temperature decreases, the probability decreases.
In the forthcoming case study, the initial temperature, ending temperature, cooling rate, and the number of iterations per temperature level of SA were set to 1, 0.000001, 0.95 and 20, respectively.These parameters were based on several trial runs of the search process applied to the case study forest landscapes.

Results
The assumption regarding the initial age class distribution not only affected the forest plan values, but also the computing times necessary to develop forest plans (Table 1).The older forest produced the highest mean solution value to the utility function (0.6594) and the highest single maximum solution value (0.6623).The standard deviation for the older forest was also significantly lower than that for other three forest landscapes.A paired t-test indicated that there were significant differences (p < 0.01) among the four forest landscapes (age class distributions) in terms of the value of objective function.As for the average computing time needed for generating a single solution, plans for the older forest landscape were also developed faster (2.4483 minutes (min) per solution) and had the tightest standard deviations (0.7038 min) for the computing times than the other four forest landscapes.With a few exceptions, the young forest landscape needed approximately half an hour to produce a single solution, and the standard deviation was approximately 15 min.One reason for this is that most of the stands within the young forest landscape had much lower volumes available for harvesting, and thus the search process needed to explore more of the solution space to produce a good and feasible solution.Meanwhile, the higher standard deviations in the computing times also indicate that the quality of the initial solution can affect (significant shorten) the computing time, e.g., the fastest solution was generated in only 4.6667 min.However, the manner in which SA changed temperature also affected these results, as temperatures were only changed after 20 successful moves per temperature, rather than a pre-defined number of unsuccessful moves that have been used in other work to prevent long computational times.In all four of the case study forest landscapes, the total scheduled harvest volumes were very close to the target volume.However, the older forest volumes were significantly different than those provided by the young and actual forest landscapes (p < 0.01), as well as the normal forest landscape (p < 0.01) (Figure 3A).The scheduled harvest volume of three 10-year periods for each forest landscape all generally represented non-declining even flow, and gradually increased over time (Figure 3B).The forest spatial values were almost the same (approximately 180) for all four forest landscapes, which accounted for about 10% of the total FSV value (Figure 3C).The spatial distribution of harvest activities for the best solution of the 20 runs for each forest landscape is shown in Figure 4.For the harvest level desired for each of the four forest landscapes, the total carbon stocks increased due to the increases in the average stand age (p < 0.01, Figure 3D).In our synthesis, the different age class distributions seem to significantly affect the final carbon stock levels.illustrates the relative changes in carbon density for the four heterogeneous forest landscapes.The average increments of carbon density for the actual, young, normal and older forest landscapes were 4.48 ± 9.61 t/ha, 8.73 ± 10.85 t/ha, 2.99 ± 9.19 t/ha and 1.03 ± 9.77 t/ha, respectively, indicating appropriate forest management practices can effectively achieve the desired synergies between forest timber production and carbon stocks within different age class structures in the next 30 years.However, the number of stands with an increase in carbon sequestered for each forest landscape was significantly different (p < 0.05), i.e., which was reduced with the increases in average stand age.For the young forest landscape, approximately 80% of the stands increased their carbon levels.However, the percentage of the increase reduced to 65% for the normal forest landscape, and the value was much smaller (55%) for the older forest landscape.

Discussion
The major focus of this study was to develop and demonstrate a new method which can allow policy makers to predict future timber production and carbon stocks simultaneously.This was addressed by presenting and solving a forest planning problem with a desired high carbon stock level and a desired scheduled harvest volume, as well as an appropriate spatial distribution of forest activities.The planning problem was also subjected to even flow of harvest volume and maximum harvesting size constraints.In this case study, the three parts of the objective function were combined into a utility function using scaled sub-priority functions.The utility function was maximized, using simulated annealing, to obtain high-quality forest management plans for the next 30 years.Other methods capable of solving complex spatial planning problems could have been selected, and other methods for optimizing the search parameters may have been investigated [48].However, considering the effort of implementing the algorithm and running these problems, SA seemed to be an expedient choice.In addition, in other studies where heuristics have been compared, SA was shown to perform well [32,33].Further, we developed three typical hypothetical forest landscapes to estimate the effect on timber production and carbon stocks of different age class structures and therefore to assess the sensitivity of the results when the process was applied to other landscapes.
Technically, the approach presented in this study is immediately available for forest planning practice, because: (1) the harvest levels are in line with the actual harvest levels in the region and the average carbon level is 43.26 t/ha, which is close to the previously reported values (45.06 t/ha) [49].However, it should be noted that not all of the components of carbon stocks related to forests and forestry were included.For example, soil carbon sequestration, comprising as much as 50%-85% of current carbon stocks in the world's forests [15], plays an important role in the global carbon cycle.Unfortunately, incorporating below-ground carbon stocks into a forest planning process is challenging and not possible at this time, due to a lack of the necessary predictive models for the forests examined in this research.Further, the constraints we analyzed are required by law or are based on guidance from policy makers and evidence from recent management behavior, yet some technical issues may hinder their common use in the near future.Finally, the shape of the sub-utility functions can have a great influence on the level goals produced and their temporal development [38].In the case studies it was assumed that all goals use linear sub-utility functions, meaning the assumptions we employed may not completely match realistic situations.
Given the widespread attention paid to global climate change over the last two decades, a potential need has arisen to develop long-term forest management plans that consider climate change [50].For this purpose, either a set of process-based growth and yield models or adjustments to empirical growth and yield models (where growth is affected by climate) is needed.Along with changes in climate, the probability of occurrence for many unpredictable events (e.g., forest fires, disease and insect outbreaks) will likely change, and thus planning processes may need to incorporate natural disturbance regimes through the use of stochastic process methods [50].Another technical reason that may hinder the application of the process involves a subjective determination of the weights that are applied to the three portions of the objective function.The relative importance of the goals we examined may vary, and forest managers may want to develop forest plans that emphasize one or more of the three goals over the others.In weighting a goal, a planner may value preference (or avoidance) of the goal using the weighting term, rather than using a normalizing process to do so [51].Therefore, when one goal is weighted differently from the others, the search process may be biased [52].In our work, we assumed all three goals were of equal weight in the objective function, and that they were normalized appropriately to adequately reflect their equal importance.
The SA algorithm, like some other metaheuristic techniques, is a stochastic search technique and is particularly highly parameterized [31], so it is obvious that a single run of the model would be insufficient to assess the quality and the stability of the results [31][32][33].Therefore, with respect to the quality of results generated by the heuristic technique, we attempted several trial runs with varying combinations of annealing parameters to determine which would best suit the problem (as in [32,39,44]).The ones we presented in this work were deemed the most appropriate among several sets of others.In assessing the stability of results generated, we observed the objective function variability through 20 runs, and found that the results vary only from 0.20% to 1.34%, indicating adequate solution stability.However, each heuristic process also has its disadvantages.The most prominent is that these methods cannot ensure optimality [32,33,53].The quality of solutions generated by a forest planning model can be evaluated by comparing results with exact techniques for their "relaxed" problems, or by estimating a global optimum solution from a set of heuristic solutions, as proposed in Bettinger et al. [41].Due to the complexity of the objective (the spatial aspect in particular), we did not solve the problem exactly using non-linear mixed integer programming methods, Further, given problems that others have noted with the estimation of a global optimum (e.g., Bettinger et al. [29]) we did not pursue that course of action either, and thus our solution validation process can be considered a self-validation [39].
As may be evident in the methods of this research, the standard version of simulated annealing was used in this case study.Other forestry research efforts have shown that some enhanced heuristic methods can further explore the solution space through intensification and diversification and possibly improve the quality of solution values.For example, Heinonen and Pukkala [53] developed 2-opt heuristic search for random ascent, Hero, SA and TS, and found it can improve solution values by 0.61%, 1.35%, 0.25% and 0.38% for five different scales of spatial forest planning problems.Li et al. [33] combined two and three heuristic algorithms (among SA, TS, TA and the raindrop method), finding that more than 75% of the 3-algorithm and the best 2-algorithm (TA + TS) produced better solutions than the best standard heuristic in terms of mean and maximum solution values.In addition, one also can combine heuristic methods with exact techniques (e.g., linear programming) to produce a new method (e.g., Öhman and Eriksson [54]), where the exact techniques can be used for solving the non-spatial dimension, and the heuristic methods are used for solving the spatial dimension of the problems.Very recently, it was also observed that search reversion can help an s-metaheuristic solve a problem, such as the one described here, more effectively [55].

Conclusions
Opportunities for enhancing carbon sequestration planning research often center on forest policy makers who are seeking strategies to address global climate change and on forest landowners with multiple-use and less intensive timber production goals [56].In certain forested areas of the world, there may be opportunities to both reach the commodity production goals of a land management organization while also positively addressing carbon stock management goals.For example, of the three forest landscapes with hypothetical age class distributions used in the case study, all of them were shown to potentially obtain higher carbon stock levels in commercial forests while sufficiently satisfying a timber production objective.This work illustrated one planning process that could be employed to assess the trade-offs among goals and whether synergies might be achieved.However, the initial condition of a forest can influence the achievement of these goals, and therefore processes are necessary to assess whether these types of goals are compatible or incompatible.

Figure 1 .
Figure 1.Age class distribution of one actual and three hypothetical forests in a 120,293 ha landscape in the northeast China; (A) actual forest; (B) young forest; (C) normal forest; (D) older forest.

Figure 2 .
Figure 2. Shape of sub-utility functions for the planning formulation: (A) the timber production sub-utility increases linearly from 0.0 to 1.0 until a user-defined scheduled volume is achieved, after which it decreases linearly back to 0.0; (B) the carbon stock sub-utility function increases linearly from 0.0 to 1.0 until the user-defined carbon stock goal is achieved, after which it is 1.0; (C) the forest spatial value sub-utility function increases linearly from 0.0 to 1.0, which is achieved at the point of maximum forest spatial value for the landscape.

Figure 5
Figure5illustrates the relative changes in carbon density for the four heterogeneous forest landscapes.The average increments of carbon density for the actual, young, normal and older forest landscapes were 4.48 ± 9.61 t/ha, 8.73 ± 10.85 t/ha, 2.99 ± 9.19 t/ha and 1.03 ± 9.77 t/ha, respectively, indicating appropriate forest management practices can effectively achieve the desired synergies between forest timber production and carbon stocks within different age class structures in the next 30 years.However, the number of stands with an increase in carbon sequestered for each forest landscape was significantly different (p < 0.05), i.e., which was reduced with the increases in average stand age.For the young forest landscape, approximately 80% of the stands increased their carbon levels.However, the percentage of the increase reduced to 65% for the normal forest landscape, and the value was much smaller (55%) for the older forest landscape.

Figure 5 .
Figure 5. Changes in forest carbon levels for the four forest landscapes with different age class distributions; (a) actual forest; (b) young forest; (c) normal forest; (d) older forest.

by Larix gmelinii and mixed with Betula platyphylla and Populus davidiana. In this classification, 106 compartments were considered natural Larix gmelinii forests (36,286 ha), 50 compartments were considered natural Betula platyphylla forests
as in China, particularly with respect to partial harvesting activities.Therefore, the forest data used in this research involved one actual and three hypothetical forest landscapes with different age class distributions.The actual landscape database comprised 123,293 ha of national forest in northeast China.The Forest Inventory and Planning Center of Great Xing'an Mountain surveyed the forests using visual compartment inventory.The total area consisted of 352 stand compartments.

Table 1 .
Summary of solution quality and computing time of 20 runs for the four forest landscapes with different age class distributions.