Forest Management of Pinus pinaster Ait. in Unbalanced Forest Structures Arising from Disturbances—A Framework Proposal of Decision Support Systems (DSS)

: Forests assume a great socioeconomic and environmental importance, requiring good management decisions to value and care for these natural resources. In Portugal, forest land use accounts for 34.5% of the continental area. The softwood species with the highest representation is maritime pine ( Pinus pinaster Ait.). Traditionally, the species is managed as pure and even-aged stands for timber production, with a rotation age of 45 to 50 years. Depending on the initial stand density, the stands are thinned 2 to 4 times during the rotation period. Disturbances associated with forest ﬁres have a negative impact on the age structure of stands over time, as they result in a narrow range of stand ages. This age homogenization over large forest areas increases with the recurrence and size of forest ﬁres, bringing new challenges to forest management, namely the difﬁculty in ensuring the long-term sustainability of the wood supply. The problem aggravates with the increasing demand pressure on pine wood. This article aims to suggest a framework of DSS for Pinus pinaster that can effectively support the management of forest areas under these circumstances, i.e., narrow age ranges and high demand of harvested timber volume. A communal woodland area in the Northern region of Portugal affected by forest ﬁres was selected as a study case. The Modispinaster model was used as the basis of the DSS, to simulate growth scenarios and interventions along the optional rotation period. Two clear-cut ages were considered: 25 and 40 years. The results obtained were the input data for an integer linear programming (ILP) model to obtain the plan that maximizes the volume of timber harvested in the study area, during the planning horizon. The ILP model has constraints bounding the area of clearings, and sustainability, operational and forestry restrictions. The computational results are a powerful tool for guidance in the decision-making of scheduling and forecasting the execution of interventions determining the set of stands that are exploited according to the different scenarios and the period in which the clear-cut is made throughout the planning horizon. Considering all constraints, the solution allows a balanced extraction of a total of 685 m 3 · ha − 1 , over the 50-year horizon, as well as the representation of all age classes at the end of the planning period.


Introduction
Forests play important roles and functions, including watershed protection, prevention of soil erosion and mitigation of climate change. Forests provide habitats for animals and livelihoods for humans and, after oceans, they are the world's largest storehouses of carbon. Inside the four major categories of ecosystem services of provisioning, regulating, cultural and supporting services, one major benefit of the forest, which is the main focus of this study, is wood production and the supply of timber. This activity is directly linked to the production chain (saw milling, carpentry for construction, wood packaging, furniture and other wood products) and can be of particular economic relevance in many countries. While the broad role of the functions offered by forests justifies adequate forest management and planning, the growing demand for wood on a global and national scale makes this an essential requirement, as it allows for verification of the possibility of obtaining a wood supply over a time horizon while ensuring the sustainability of the wood supply. This planning procedure goes beyond traditional "forest regulation," leading to the need for intelligent management of unbalanced forest structures [1] and of specific and advanced tools.
The application of computerized tools designated as decision support systems (DSS) emerged at the beginning of the 1980s. DSS are computerized technological resources capable of supporting business and organizational decision-making activities by compiling useful information through data analysis, model study, and support tools designed to assist decision makers in identifying and solving problems [2][3][4][5][6][7][8]. They aim to provide the best decision-making process by improving the efficiency and productivity of managers and professionals [9]. The application of DSS to decision-making in forestry has been increasing over the years, embracing different problems and goals (e.g., [10,11] and references herein).
DSS presents particular interest when optimized management solutions for forestlevel planning are the aim. Forest DSS can help forest managers in management planning by building knowledge and capacity through the use of advanced techniques and tools (e.g., growth and yield simulators, optimization techniques, tools of communication and graphical representation, among others), that helps to provide support in the decisionmaking process.
This research paper exposes the advantages of using simulation models and optimization models, under a framework of a DSS, for testing alternative management in imbalanced forest structures using as a case study, the management of maritime pine (Pinus pinaster Ait.) forests in Portugal.
Maritime pine is the most represented softwood species in Portugal [12], presenting a high economic, social and ecological interest. Nevertheless, the availability of the resources is below demand levels [13]. The incidence and recurrence of forest fires have strongly impacted the availability of the resources and its future sustainability. According to Centro Pinus [13], pine wood consumption in 2020 was 4.1 million m 3 without bark. In 2018, this consumption was over 4.2 million m 3 , where the availability of this wood was only 2.25 million m 3 , which supports the fact that wood consumption is higher than its availability and there is a deficit of wood for consumption. In accordance with the possibility of the existing cutting down in the forest, the structural deficit of maritime pine wood amounted to 56.6% of wood consumption in 2020. This deficit tends to worsen as Portugal's wildfire risk continues to increase [14], so proper management planning of the forest production function takes on a high importance. Ideal management options considering to climate change expected impacts are not among the ones that can be currently implemented in the common lands as there are two major identified difficulties: presence of large number of naturally regenerated dense areas and lack of policies to favor (costly) adaptive silvicultural measures (e.g., reforestation with improved material; see PEN-CAFoRRCOST Action, http://www.pen-caforr.org/for overall options accessed on 8 March 2021).
It is agreed that the use of DSS helps to efficiently support the management of forest resources. Examples of DSS applied in Portugal are AGFLOR [15] and web-based SADfLOR [16], the former being designed to help assess the impacts of policy changes on regional land use patterns whereas the latter is a complete DSS that integrates an information system, a supply-driven stand level forest simulator and several optimization tools to support forest management planning. In terms of DSS for maritime pine forest, Garcia-Gonzalo et al. [17] applied SADfFLOR to generate strategic management plans for the Leiria National Forest in Portugal aimed at the production of wood. Almost 80% of the state managed area of the Leiria National Forest are even-aged homogeneous stands of maritime pine spanning a large range of ages. Other relevant references are [18,19] The suggested framework of a DSS applied to forest management of maritime pine stands here presented is supported by real data of forest inventory. It concerns a forest area located in the Perimeter Florestal of Barroso, in Ribeira de Pena, in the Northern region of Portugal. This area is made up of common lands, called baldios, which are co-managed by the local communities and the Institute for Nature and Forest Conservation (ICNF). Although these forests provide several products, services, and other benefits, primary incomes are directly related to timber removal following thinning and clear-cutting. With time, the forest disturbances have been altering the forest structure, resulting in a narrow range in terms of stand age. Moreover, the stands that are naturally regenerated after forest fires usually present a reduced or absence of management.
To help with the decision, growth and yield models are available for use in the region [20][21][22][23], one of them, the ModisPinaster [22] available as growth simulator of freely use and used in the DSS framework. During the last two decades, efforts have been made to implement the use of these tools, through workshops and dissemination events. Further, optimization models have been developed for more efficient management of these forest resources [24][25][26][27].
With this research, the authors want to highlight the importance of DSS and optimization models as a powerful tool to provide decision support in forest management of areas that are structurally unbalanced, showing limited age differentiation, due to fire disturbances and being under pressure from wood utilization. Major differences in forest structure from the case studied by Garcia-Gonzalo et al. (2013) [17], for maritime pine, in Portugal, namely with regard to the range of ages that in [17] is large (nine 10-year age classes are used, between 0-10 and >80 years) and on growth development of the stands, justified the design of a new approach for the forest area considered in the present research.
The DSS framework (Section 2) considers two scenarios with distinct rotation length with the aim of maximizing the volume of timber harvested to meet the interests of the local community. One of the scenarios considers a rotation length of 40 years, in accordance with the traditional silvicultural models. The alternative scenario considers a reduced rotation length of 25 years, approaching the age of harvest of new silvicultural models for the species, suggested in [28]. Managing the stands with shorter rotation lengths can be considered both as an adaptive measure to account for the forest fires risk and a means of providing viable yields from shorter term harvesting.

Study Area Characteristics
The study area is located in North of Portugal, in the Tâmega valley, in the county of Vila Real, Figure 1. It focuses on the communal area of Santo Aleixo, which belongs to the national forest area of Perímetro Florestal do Barroso, occupying 439 ha [29,30].
The available information for each stand was collected on forest inventories and refers to stand area (A, ha); stand age (t, years); number of trees (N, trees/ha); stand basal area (G, m 2 ·ha −1 ); average diameter (dg, cm); dominant diameter (d dom , m) and dominant height (h dom , m). The variable dg is calculated as the quadratic mean diameter measured at the height level of 1.30 m of the living trees. The variable d dom corresponds to the average diameter of the largest trees in a basis of the largest 100 trees per hectare. The average height of those dominant trees is the dominant height. A summary of the data is shown in Table 1.  Legend: A-stand area (ha); t-stand age (yr); N-number of trees per hectare (trees·ha −1 ); G-basal area per hectare (m 2 ·ha −1 ); dg-quadratic mean diameter measured at the height level of 1.30 m (cm); d dom -dominant diameter (cm); h dom -dominant height (m); Min-data minimum; Mean-data average; Max-data maximum; and sd-data standard deviation.

Simulation Scenarios Using ModisPinaster
Growth simulators provide useful information about stand development and production under traditional or user specified silvicultural scenarios, which can be used to support decision making. In this work, the Model with Diameter Distribution for Pinus pinaster (ModisPinaster) was selected as the appropriate simulator as it was developed to be used in the overall region of Tâmega's Valley, where the case study is located. Mod-isPinaster [21,22] is a model developed for maritime pine, that simulates stand growth and yield, predicts the occurrence of mortality associated with competition or damages from extreme weather events, such as windrow and uprooting, and presents flexibility to select alternative silvicultural procedures to support forest management. It is a model that is available for free use, in a friendly interface, on the Computer Aided Projection of Strategies in Silviculture (CAPSIS) platform (http://www.inra.fr/capsis accessed on 8 March 2021) [31]. A comprehensive explanation of ModisPinaster's components and examples of its use can be seen in a reference by Fonseca et al. [22].
The minimum input data required to the initialization of ModisPinaster is information about the variables stand age (t, yr); number of trees per hectare (N, trees·ha −1 ); basal area per hectare (G, m 2 ·ha −1 ); dominant diameter (d dom , cm) and dominant height (h dom , m), as defined in Section 2.1.
After initializing ModisPinaster with input data, simulations are run on an annual basis from the initialization point up to a maximum value of 65 years. The model allows testing alternative approaches to density regulation (manual or automatic prescriptions based on spacing factors and density indices, in addition to the usual density measures of number of trees or basal area per hectare) and provides diverse information as output, such as stem volume, biomass (aboveground, per component, and root), carbon, and energy, among other characteristics. The output data is accessible at any age within the simulation period. The authors selected as output information the standing volume and the volume of material that is removed in thinning. These volume values serve as input data in the proposed DSS. We considered a stand density regulation model based on the stand density index (SDI), as defined for maritime pine by Luis and Fonseca (2004) [32], ranging from 35% (the lower limit of complete occupation) to 55% (a value close to the lower limit of self-thinning, or natural mortality by competition). That is, when the stand reaches a maximum stipulated value of SDI equal to 55%, the model simulates thinning from below, so that the remaining stand has an SDI equal to 35%.
The simulations were performed for two different scenarios, presented in Figure 2. Scenario 1 consists of maintaining the stand until 25 years of age and, at that age, performing a clear-cut. For stands older than 25 years at the first period, Scenario 1 assumes that a clear-cut is performed in the first year of simulation, and this simulates the growth of a new regenerated stand for another 25 years and a clear-cut is performed again. Scenario 2 consists of letting the stand develop to 40 years of age and then performing a clear-cut at that age. Thinning is only carried out on stands that were managed according to Scenario 2 because, in the case of Scenario 1, as the age of rotation age is 25 years, the inventory data showed that the density of the stands was below the stipulated threshold for thinning.

Proposed DSS Structure and Components
It is intended to create a growth simulator with an interactive interface, which facilitates the user's understanding and is simple to apply and manipulate the results obtained. Through the insertion of the output data obtained with the ModisPinaster model, it is intended that this DSS framework provides management strategies in different time periods and allows the identification of interventions to be performed in the study area. The main objective of this system is to obtain optimal management plans, to optimize the resources provided by this type of forests, considering environmental and sustainability features.
The proposed DSS is composed of an external component, the ModisPinaster model, which performs the initial data processing (simulation of growth over time and of specific silvicultural interventions, the thinnings) and of internal components which are the data from ModisPinaster, the optimization model and the user interface. Figure 3 presents the structure and components of the proposed DSS. The objective of the DSS model is to maximize the removed timber. Therefore, it is essential to survey management alternatives and to choose the one that best fits the problem, providing the user with a solution corresponding to the maximum achievable volume, while assuring a specified set of constraints.

Optimization Model
This section presents the optimization model to obtain the management plan for a given study area aiming to maximize the harvested timber volume, considering constraints that limit the areas of clearings, and also sustainability, operational and silvicultural constraints. The practices that are considered are thinning and clear-cuts following one of the two scenarios, Scenario 1 and Scenario 2, presented in Section 2.2. The problem consists of determining which set of stands are harvested under Scenario 1 and which are harvested under Scenario 2. In addition, the year/period in which the first clear-cut is performed is also identified.
Let T = {1, 2, · · · , n} be the planning horizon with n one-year periods, P be the stands set, and K be the set of age classes.
We assign the following parameters for each stand j ∈ P: A j -area of stand j (ha); t jage (years) of stand j in the first period; C l j and D l j -timber volume obtained by, respectively, clearcutting and thinning stand j in period l (m 3 ); S j ⊂ T-set of periods to perform a thinning at stand j.
It is also considered parameter A max = 10 ha which is the prescribed value in Portugal for the maximum area that can be cut in the same year and the following two years. To formulate constraints on the maximum clear-cut area and the green-up constraints, consider set l of all possible clusters that cannot be harvested as a whole at period l, with l ∈ T, and which are minimal [27,33]. Each of these clusters is a contiguous group of stands, whose total area is greater than A max and that does not contain any similar cluster in the sense that if a stand is excluded from a cluster then the remaining could be cut at once. If r is the number of stands in such a cluster then at most r-1 of them can be harvested during the green-up period (3 years). It should be noted that two stands are adjacent when both share a boundary that is not a discrete set of points. The minimal infeasible clusters are obtained recursively, starting from clusters with only one stand, and growing them one stand at each time, until the maximum area is met.
We define the following variables: x 1 l j (clearcutting Scenario 1 binary variable) taking value 1 if a clearcutting is performed, adopting Scenario 1, in stand j in period l and 0 otherwise; x 2 l j (clearcutting Scenario 2 binary variable) taking value 1 if a clearcutting is performed, adopting Scenario 2, in stand j in period l and 0 otherwise; z 1 j (Scenario 1 binary variable) taking value 1 if stand j is harvested according Scenario 1 and 0 otherwise; z 2 j (Scenario 2 binary variable) taking value 1 if stand j is harvested according Scenario 2 and 0 otherwise.
In the following is presented the base forest management optimization model, FMO, which aims to maximize the volume of wood harvested during the planning horizon considering only two groups of constraints: the ones that impose that each stand is managed by one of the two scenarios with the suitable ages, constraints (2)-(8), and the ones that limit the area of clearings, constraints (9).
Subject to: The objective function (1) corresponds to maximization of the volume of timber removed over the planning horizon.
Constraints (2) guarantee that each stand is harvested under one of the two scenarios, 1 or 2. Constraints (3) identify, for each stand, the harvested scenario adopted. Constraints (3) link variables x and z, ensuring that if variable z 1 j takes the value one then there is at least one clear cut of stand j under Scenario 1, otherwise, no clear-cutting takes place over the planning horizon, under this scenario. Moreover, if variable z 2 j takes the value one then there is one clear-cut of stand j under Scenario 2, otherwise, no clear-cutting takes place over the planning horizon under this scenario. Note that a stand can be subjected to a maximum of two clear-cuts whether Scenario 1 or Scenario 2 is adopted.
Constraints (4)- (7) impose conditions on the date of the first clear-cut, considering a tolerance tol years. For a stand with initial age greater or equal to 25 and less than 40, constraints (4), impose that if Scenario 1 is adopted the first clear-cut must take place in the first period or in the following tol periods. Constraints (5) imposes that if a stand reaches 25 years for an intermediate period and is managed under Scenario 1, then it must be cut in that period or in the following tol periods, as long as it is within the planning horizon. Constraints (6) impose that if a stand reaches 40 years in period l of the planning horizon and is managed according to scenario 2, then it must be cut in that period or the following tol periods. Constraints (7) impose that if the initial age is greater or equal to 40 years then the first clear-cut must take place in the first period or in the following periods within the tolerance tol, independently of being managed according to scenarios 1 or 2.
Constraints (8) impose an interval of at least four years between interventions, that is, it has to be at least four years between a thinning and a clear-cut.
Constraints (9) are known as path constraints. These constraints prevent the formation of clear-cuts whose areas are greater than 10 ha and ensure the green-up period requirement [25,26]. Constraints (10) establish the variables domain.
Besides constraints (2)-(10), after some computational tests, it was concluded that, although the green-up constraints contribute to ensuring sustainability, they are not sufficient, it would be necessary considering a group of sustainability constraints, which will be described next.
Additional constraints on the final average age and age structure will be included. Let I l j be an integer age variable representing the age, in years, of stand j in period l. It holds: The following constraint imposes that the average age of the forest stands at the end of the planning horizon should be at least Ft years: where I n j is the integer variable corresponding to the age of stand j in the last year of the planning horizon. The obtained model, including this final age constraint, is denoted by FMO_FA.
To ensure sustainability, constraints on the age structure of the forest will be also included. Constraints (15) ensure that, on a period l, the area in each age class is within (1 − ∆ t ) and (1 + ∆ t ) times the target area for each class, Ta: where w l jk are age class binary variable taking value 1 if stand j in period l belongs to age class k and 0 otherwise. The obtained model, including the class age constraint, is denoted by FMO_CA.
The following constraints, (16)- (19), should also be included to state the age class that stand j belongs at each period l. Considering the age class k ∈ K be given by k = [L k , U k ], we have: ∑ k∈K w l jk = 1, ∀j ∈ P , l ∈ T w l ji ∈ {0, 1} ∀j ∈ P, l ∈ T, i ∈ K As operational constraints, a balanced income over the planning horizon is guaranteed. To do so, the planning horizon is divided into several consecutive periods, T 1 , · · · , T d , with equal length and it is ensured that the relative difference in the removed volume between two consecutive periods does not exceed ∆: where VolumeP k is the volume of timber removed in the k-th period, T k , given by: The obtained model, including balanced revenue of volume constraints, is denoted by FMO_RV.
To demonstrate the dynamics of a DSS, the original model was adjusted according to the achieved results. In a first approach, model FMO was considered, where neither the age constraints, nor the revenue regularity were admitted. These sets of constraints were gradually included, resulting in the five models presented in Table 2. In order to solve the models, we used the optimization software Xpress 7.2 (available at http://www.fico.com/Xpress accessed on 11 January 2021). Computations were performed on a computer equipped with an i5-2500 3.3 GHZ/3.7 Turbo CPU and 8 GB of RAM. The branch-and-bound algorithm was allowed to run for two hours at the most.

Results
In this section the computational results for the Santo Aleixo baldio are presented.

Remarks about the Model Parameters
The set P has 65 stands, characterized on Table 1. n = 50 was considered as the number of years in the planning horizon, and five age classes, K = {[0, 10], [11,20], [21,30], [31,40], [41, 65]} were chosen. The upper bound on the last class is in accordance with the maximum age that the stands can achieve during the planning horizon.
Regarding the tol parameter, constraints (4)- (7), it should be noted that there was a need to make the date of the first clear-cut more flexible due to the clear-cut area constraints because there are several adjacent stands with the same characteristics. After several computational tests, to gets a feasible model, a margin of tol = 14 years to perform the first clear-cut is considered, that is, a stand that is in condition to be cut in period l (reached 25 years being harvested according to Scenario 1 or reached 40 years being harvested according to Scenario 2) should be clear-cut in this period or the fourteen following periods. After the first clear-cut, the stand starts a new development phase of natural regeneration and is cut again when it reaches 25 years (for Scenario 1) or 40 years (for Scenario 2).
Furthermore, based on computer experiments, it was considered F t = 20 for the average age of the stand in the last year, constraints (14).
To get feasible models, due to the age structure of the study area in the beginning of the planning horizon, constraints (15) are only included for the last period, l = n, with ∆ t = 0.4. Moreover, given that the latest cutting age is 40 years, we will consider at the end of the planning horizon only the first four area classes in K and, as the total area is 439 hectares, it is considered that T a = 100 ha.
On constraints (20), ten 5-year periods, T 1 , . . . , T 10 , and ∆ = 0.4 were considered, values that allow one to find the optimal solution with acceptable variations in the volume of removed timber.

Obtained Results with the Framework DSS
In Figure 4 is presented, as a bar chart, the volume of timber removed per ha, in each five-year period during the planning horizon for formulations FMO and FMO_all. In Table 3, the optimal solutions for the base model, FMO, and the model with all constraints incorporated, FMO_all, are characterized in terms of the number of stands and area managed under each one of the scenarios, as well as the removed timber volume by clearcutting and by thinning, the total removed volume, and the final average age of the maritime pine stands in baldio of Santo Aleixo. The results for the intermediate formulations, FMO_FA, FMO_CA and FMO_RV, are presented on Table A1 in the Appendix A. Tables 4 and 5 characterize by each 5-years of the planning horizon the number of stands and the area managed under scenario 1 and under scenario 2, specifying the practice undertaken, clear cut or thinning, for models FMO and FMO_all, respectively. Table 4. Characterization of FMO results according to the optimal solution during the planning horizon, discriminated by scenario.     Figures 6 and 7 illustrate the spatial representation of the optimal solution obtained with the FMO-all model. The cartography presented concerns about the interventions made during the first ten periods, with the indication of the period in which it takes place. Figure 6 shows the clear cuts and Figure 7 shows the thinning.

Discussion
Analysing the results presented in Section 3, in particular Tables 3-5 and Figures 4 and 5, it can be concluded that although the total volume of harvested timber is smaller with model FMO-all, it is the one that provides a result more coherent with a sustainable management. Similar results have been achieved also in [34], even with the connection with other spatial environmental requirements. In fact, in the optimal solution of model FMO_all the area by age class varies more evenly over the planning horizon, Figure 5, and there are less differences between the removed volume on each 5-year period.
It should be noted that, as the initial age of the stands' study area has a range of values with reduced amplitude it was only possible to obtain feasible problems by imposing constraints (15) only at the end of the planning horizon, which will ensure sustainability in later periods. During the planning horizon, the volume regularity constraints already guarantee revenue over the entire planning period. Regarding the value of the average age of the stand in the last year, F t , computational tests were carried out and, in the case of including revenue regularity constraints, (20), the problem becomes infeasible for values of F t greater than 20. With no constraints (20) one can consider F t = 30 but, in this case, in the 5-year period from 36 to 40 there is no revenue and in the period from 41 to 45 there is only 743 m 3 , which does not meet the interests of the populations that are interested in receiving some regular dividend.
Also, if we include in the case F t = 30 and the constraints (15) we obtain an unfeasible problem. It should be noted that the fact that the area under study has a very unbalanced distribution in terms of areas by age class at the beginning, Table 1 and Figure 5, force the choice of ∆ and ∆ t parameters greater than what could be anticipated, as well as a value of F t lower than expected. However, the model enables to obtain a balanced area class structure in the last year of the planning horizon, which ensures even more balanced planning afterwards. This is the great added value of the work presented.
In the optimal solution of FMO, Table 4, all 47 stands managed under Scenario 1 are subject to two clear cuts during the planning horizon and nine of them undergo thinning in the first decade. The 18 stands managed under Scenario 2 are subject to one clear cut during the planning horizon and all of them are thinned at least once and at most three times. In the optimal solution of FMO_all, Table 5, only two of the stands managed by Scenario 1 are thinned, and all 40 stands managed under Scenario 1 are subject to two clear cuts during the planning horizon. The 25 stands managed under Scenario 2 are subject to one clear cut during the planning horizon and only three of them are not thinned while the remaining stands suffer one to three thinnings.
The inclusion of sustainability and operational constraints leads to a reduction in the amount of timber volume harvested by 0.3% in the case of model FMO_FA, 0.7% in the case of FMO_CA, 3.8% in the case of FMO_RV, and 4.0% in the case of model FMO_all. A reduction on the timber volume was also reported by [17] with the DSS that the authors proposed for maritime pine in Coastal Region of Portugal adding sustainability constraints.
Analysing the solution obtained from FMO_FA and FMO_CA, although the stand age concerns are considered, the revenues are not regular. During the 5-year period from 41 to 45 the amount of extracted volume is greatly reduced, 1133 m 3 , while in period from 11 to 15 year it is very high, 80,114 and 89,339 m 3 . This situation is not the most favorable for the local communities who aim to obtain a regular yield over time. Hence the need to include the constraints (20) to prevent large differences in the volume harvested each five-year period, leading to the formulation of model FMO_RV. This model provides wood removals that are more balanced over time.
Analysing Tables 3 and A1, it can be concluded that using any of the models, Scenario 1 is the most adopted, both in number of stands and in the intervened area. This difference is more evident with model FMO, where 70.9% of the study area is managed under Scenario 1 and the remaining 29.1% of the area is managed under Scenario 2. For the other models around 65% of the area is managed under Scenario 1.
Regarding to the differences between the removed timber volume in thinning and clearcutting, in any of the models, the volume of material removed in clearcutting is much higher, as it would be expected (Tables 3 and A1). In the case of model FMO, the total volume of timber harvested is 313,214.2 m 3 , of which 83.7% is from clearcutting and the remaining is originated from thinnings. For the model FMO_all, the total volume of timber harvested is 300,717.8 m 3 , of which 80.8% is removed in clear-cut. For the other models, around 83% of the volume of wood harvested was due to clearcutting. Figures 6 and 7 show the proposed interventions with FMO-all between periods 1 to 10. There are no interventions in periods 3 and 9. There are two stands, 142b1 and 142a2, which are thinned on period 1 and will be clear cut during period 7. In Figure 6 the effect of application of the constraints (9) can be seen, which prevent the creation of clearings with an area larger than 10 ha and guarantees the green-up requirements with an exclusion period of 3 years. A clear-cut is performed at period 1 in stand 141b1, and 3 years later an adjacent stand, 141b2 is clear cut. Stands 144a2 and 146a1 are clear cut on period 2 and only after at least three years will be adjacent stands clear cut: 146a2 clear cut on period 5; 142a3 on period 7; 144a1 and 144a3 on period 10; 147a1 on period 8.
The proposed model is consistent with the management in short rotation and with the specification of the problem that was intended to be studied: it guarantees a balanced age class structure at the end of the planning horizon; it contemplates restrictions on the maximum area of clearings that are in accordance with what is stipulated in Portugal, which is a maximum of 10 hectares; the operational constraints ensure a sustainable management plan that meets the interests of the communities.
Next steps for future development of the proposed model and decision support framework should be incorporating the constraints expressing the different ecosystem services which can change the optimal solution significantly. Even that, the demand for multiple ecosystem services is increasing [35].
In terms of generalization to other cases of study, the proposed approach of a DSS framework can be easily extended or adjusted to application to other regions or forest species, provided that the adequacy of the input information is ensured. This practical implementation encompassed the use of the ModisPinaster model, from which information was acquired for insertion into the proposed optimization model. ModisPinaster was developed for application to maritime pine forests of the northern region of Portugal but can be used to support the simulation of pine forests in other regions and countries (Spain, France). As it works to obtain information to support to the DSS, the suitability for other areas (preferably validated before its use) does not limit the interest of the presented DSS. For other species, appropriate forest growth model specific to those species should be selected. The optimization models developed can be used and/or adapted to other regions, regardless of the species, where scenarios that contemplate sustainability constraints such as those indicated in this model are to be studied.

Conclusions
The research presented refers to forest management planning of stand areas that have grown after rural fire disturbance showing an unbalanced age distribution observed at the beginning of the horizon planning with a narrow range of stand age range.
The selected case study represents not only a real situation but also the expected pattern in the age structure of stands to occur naturally following forest disturbances, such as those associated with recurrent forest fires in Portugal, and generalizable to other causes.
The optimization model created allowed to essay different management alternatives, providing the best one and making a good selection of the available information, in a relatively short period of time. Some difficulties arose during the development of the model because, being in the presence of a common land where the majority of the stands in the study area were more than 25-30 years old, the implementation of interventions became quite restrictive in the sense that if no rules were imposed on the implementation of interventions, the planning horizon would end with a very low average age. As output information, the users have a table indicating what to do in each stand and throughout the planning horizon (we chose not to include it in this article so as not to make it too long), and maps that provide a spatial and temporal perception of the common land area and the interventions to be performed.
The application performed makes the forest management process more efficient, ensuring the satisfaction of legal requirements, the sustainability of resources and the balance of revenues. Further developments should enlarge the optimization problem in order to take into account the broad role and multiple functions of forests. For pine forests, two of the relevant functions that might be of importance for the DSS are the production of resin and the carbon sequestration.

Appendix A. Complementary Results
This section presents results for the intermediate models, FMO_FA (model with constraints bounding the final average age), FMO_CA (model with area by classes of age regulation constraints), and FMO_RV (model with balanced revenue of volume constraints).
In Table A1, the optimal solutions for models FMO_FA, FMO_CA and FMO_RV are characterized in terms of the number of stands and area managed under each one of the scenarios, as well as the removed timber volume by clearcutting and by thinning, the total removed volume, and the final average.