Simulating Araucaria angustifolia (Bertol.) Kuntze Timber Stocks With Liocourt’s Law in a Natural Forest in Southern Brazil

: This paper presents a simulation of the regulation of Araucaria angustifolia (Bertol.) Kuntze timber stocks using Liocourt’s law. Although this species is currently protected by law, recent government initiatives are being considered to propose sustainable forest management practices by selecting small rural properties in Southern Brazil. Here, we simulate the applicability of Liocourt’s law in a typical rural property, the size of which is approximately 85 ha. Forest inventory measurements were conducted by estimating the diameter at the breast height ( d ), total height ( h ), and annual diameter increments of 308 trees that ﬁt the criteria of d ≥ 10 cm, distributed on 35 permanent plots of 400 m 2 each. As a result, a reverse J-shaped d distribution was found. On average, 303 trees can be found per hectare (ha). Local allometric equations showed their basal area (G) to be 21.9 m 2 · ha − 1 , and their commercial volume ( V ) to be 172 m 3 · ha − 1 , whereas Liocourt’s quotient (q) was 1.31. Based on these attributes, nine di ﬀ erent forest management scenarios were proposed by simulating a remaining basal area (Grem) of 10.0, 14.0, and 18.0 m 2 · ha − 1 , and Liocourt’s quotient was changed to 1.1, 1.3, and 1.5. All scenarios consider a d of 62.5 cm. In the less intensive scenario (i.e., q value = 1.5 and larger basal area of 18.0 m 2 · ha − 1 ) there is greater optimization of space, and higher economic return is ensured to rural producers due to the deﬁnition of shorter cutting cycles. This also allows a faster growth rate in both d and h for smaller trees, due to the higher incidence of light onto the lower canopy layer, increasing the natural regeneration implementation of other native species. Forest management should thus be considered a goal in addition to consumer market characteristics for deﬁning the ideal timber stock scenario.


Introduction
Mixed-species stands provide greater ecological and socio-economical services since they can be more resistant to natural disturbances [1,2], more stable with regard to climate change or disease control [3], and more productive than homogeneous stands [4,5]. Despite this, quantitative silvicultural guidelines facilitating the efficient management of mixed-species stands are still largely based on research information pertaining to pure even-aged stands [6].
Appropriate forest management practices serve to ensure that raw materials are exploited in a continuous and sustainable manner [7], i.e., by combining economic, social, and environmental benefits through efficient planning to allow the better management of the available resources, as exploiting more than the forest can produce at a given period may lead to the exhaustion of said resources, and removing less than what it can produce may lead to their underutilization and overpopulation, resulting in stagnation [8].
Growth and production models are important forest management planning tools, as they allow simulating the development of settlements and production under several management alternatives [9]. Cutting a single tree or a small group of trees distributed on a limited area creates spaces that can improve light conditions in the understory, promoting the growth of seedlings and young trees and improving the composition and structural diversity of plant species by forming habitats where they can coexist [10,11].
The regeneration and establishment of seedlings and their subsequent transition to larger classes must occupy the spaces left by the removal of bigger trees, with adequate intervals and intensity of cutting, which depend on the species' regeneration and growth rates [12]. Given this context, it is important to obtain information that allows the rational management of one of the Atlantic Forest Biome's most characteristic formations, found in the southern region of Brazil: Mixed Ombrophilous Forests (MOF) [13], also known as Araucaria Forests [14,15], with the implementation of plans and practices directed towards the species' characteristics.
MOF are characterized by the presence of Araucaria angustifolia (Bertol.) Kuntze (referred to here as A. angustifolia), as it forms the upper canopy layer and has a dominant character in this vegetation type. In spite of this, the middle and lower layers demonstrate large diversities of species, such as Ocotea porosa (Nees & Mart.) Barroso, Luehea divaricata Mart., Cedrela fissilis Vell., Dicksonia sellowiana, and Maytenus ilicifolia Mart. ex Reissek [16,17]. As A. angustifolia is currently considered an endangered species and protected by law [18], the management of forested areas is mostly forbidden. However, recent government initiatives are being considered to show that natural forests in Southern Brazil can be sustainable and provide environmental, social and economic benefits besides being an important part of the life of local farmers [19,20].
Some studies have already verified the potential of regulating the diameter at the breast height (d) of A. angustifolia using Liocourt' law [21] and quotient [22][23][24] to achieve sustainability by maximizing the growth of individual trees and, consequently, the production of timber [22][23][24]. However, all these studies fail to simulate changes of this quotient in association with the remaining basal area, which implies differences in cutting cycles and logging intensities. Such aspects must be clearly stated in sustainable forest management protocols by implementing reduced impact logging. More details and a better explanation of Liocourt's law and its practical applications are available in Kerr [25], Gove [26], and Govedar et al. [27]. Therefore, this study sought to simulate the regulation of A. angustifolia timber stocks using Liocourt's law in a typical rural property in the Southern region of Brazil. More specifically, it aimed to: (i) evaluate the forest's current diametric structure; (ii) define different forest management scenarios for the species based on Liocourt's law; (iii) evaluate the classification of wood assortment; (iv) provide insights for future forest management practices. We argue that sustainable forest management is the best way to promote the conservation of MOF and the renewal of this endangered species.
The remainder of this paper is organized as follows: Section 2 describes the rural property where this experiment was conducted and introduces the fundamentals of Liocourt's law. The protocol Forests 2020, 11, 339 3 of 14 followed in our experimental analysis is then presented. Section 3 shows the outcomes of the experiment and the simulation's results, whereas Section 4 discusses the experimental results. Finally, Section 5 summarizes the main conclusions of this study, and points to future directions.

Study Site
The study site is located in a typical rural property in Southern Brazil, of the size of which is 83.5 hectares. It is located in the municipality of Lages (Santa Catarina, Brazil), and its central geographic coordinates are 27 • 48' S and 50 • 19' W. More than 50 ha of the rural property's total size are covered by uneven-aged MOF, in which A. angustifolia is the dominant species. Although some authors consider A. angustifolia a pioneer species [28], there is also strong evidence that this is a climax species, which regenerates at different levels of light within the forest [29] and has seedlings that are extremely shade-tolerant [30]. The study site's climate characterization, according to the Köppen classification, is mesothermic humid (Cfb), without dry season and with temperate summers. The areas's altitude above sea level is approximately 987 m, with average annual temperature of 15.2 • C, and average annual precipitation of 1685.7 mm [31]. The predominant soil type in the region is Haplic Nitisol and Humic Cambisol, developed from basaltic rocks.

Data Collection and Preparation
A total of 35 permanent rectangular plots of 400 m 2 each were installed in the study site using the fixed area method, and distributed using a systematic sampling approach. The spacing between sample plots was 50 m by 100 m, resulting in a grid. In each sample plot, all trees with diameter at breast height (d) equal to or greater than 10 cm were measured. Sample intensity and sampling error were evaluated according to the methodology proposed by Sanqueta et al. [32].
In addition, two perpendicular increment cores were radially extracted at d from the 308 measured trees using a Pressler borer. The sampling included trees in all d classes, covering the entire diametric amplitude of the species. Such measurements allow quantifying the periodical increment in d. The increment core samples were fixed to supports, air dried and sanded in a laboratory in order to make it easier to visualize the tree rings. Following the delimitation of the tree rings, the samples were digitally scanned (1200 pixels resolution). After scanning, the last five tree rings (corresponding to the last five years of growth) were measured using Image Pro-Plus© (Media Cybernetics Inc., Silver Spring, MD, USA).

Balanced Frequency Distribution
After field data collection, the individuals were divided into diametric classes with regular 5-cm intervals to evaluate the forest's diametric structure. Meyer's equation [33,34] was adjusted to the sample units' data to represent the estimated frequency distribution per unit area as seen in Equation (1): where: Ni = frequency per hectare in diametric class i; d i = d class center, in cm; β 0 ; β 1 = estimated regression coefficients; e = base of natural logarithm; ε = residual error. The balanced frequency distribution was then applied according to Liocourt's law, demonstrating the proportionality between the number of trees of the current class and the subsequent one; this reason is called Liocourt's quotient, or "q" [35]. Based on the angular coefficient (β1) of Eqaution (1) adjusted to araucaria trees, it was possible to determine Liocourt's quotient: For the volume calculation, we followed the Kozak model [36] (Equation (3)), adjusted to the same area of the study by Costa et al. [37]. The model showed adjusted determination coefficient (R 2 adj.) = 0.980 and precision expressed by root mean square error in percentage (RMSE%) = 6.8.
where: d i = diameter in relative positions along the shaft in cm; h = total height in m; h i = relative height in m; d = diameter at breast height in cm; β 0 = 0.8944; β 1 = 1.0361; β 2 = 0.9994; β 3 = 0.0 (ns); β 4 = −0.1885; β 5 = 1.1267; β 6 = −0.2960; β 7 = −0.0420; (ns) = non-significant estimated regression coefficient at α = 5%, therefore, the model was adjusted without the respective coefficient; ε = residual error; ln = natural logarithm; p = inflection point considered in 1.3/h. However, because it is not analytically integrable, the Kozak model required the use of numerical integration techniques to obtain the partial volume (v i ) and total volume up to the crown (v tci ). In this way, with the reconstitution of the radial increment data at d of the last five years, the volume at the beginning and end of the selected period was obtained using the numerical integral approach. In order to accomplish this, a set of scripts were created in Visual Basic using Simpson's 1/3 rule method [38]. Then, the annual percentage of volume increments up to the crown (PAI% vtci ) for each tree was obtained by Equation (4), and later adjusted by Equation (5) to describe the variation of volumetric increments in all the diametric distribution of the forest.
where: PAI% vtci = annual percentage of volume increments up to the crown, with bark; vtci i+5 = total volume increments up to the crown at the end of the period in m 3 , with bark; vtci i = total volume increments up to the crown at the beginning of the period in m 3 , with bark; n = time period considered, in this case, five years.
where: PAI% vtci = annual percentage of volume increments up to the crown, with bark; d = diameter at breast height in cm; β 0 ; β 1 = estimated regression coefficients; ε = residual error. Following the adjustment of Equation (5), the average cutting rate was determined (PAIM% vtci ) by weighting the increment for the di class center, based on the volumes observed in the d i classes, according to Equation (6) [22]: where: PAIM% vtci = mean annual percentage of volume increments up to the crown, with bark; Vt i = total volume increments up to the crown in the d i class in m 3 ·ha −1 , with bark; PAI% vtcii = annual percentage of volume increments up to the crown in the d i class.

Cutting Cycle
The cutting cycle was estimated using Equation (7), which assumes that the volumetric annual growth of a tree or settlement accumulates according to the compound interest law [7]: where: Vt = total volume at the end of the period in m 3 ·ha −1 , with bark; Vr = remaining volume at the beginning of the period in m 3 ·ha −1 , with bark; i = annual growth rate; cc = cutting cycle in years.

Wood Assortments
Although the assortments may be redefined to any dimension with the use of Equation (3), for this simulation, the assortments with bark were divided into four scenarios, as follows: S1 = timber destined for sawmills with d ≥ 40 cm and length of 5.4 m; S2 = timber destined for sawmills with d ≥ 30 cm and ≤40 cm, and length of 2.7 m; S3 = timber destined for sawmills with d ≥ 20 cm and ≤30 cm, and length of 2.2 m; S4 = timber destined for the energy or wood industry with d ≤ 20 cm and the remaining available length.
While defining the assortments, the formation of timber stocks with larger diameter and length was prioritized, following the sequence S1, S2, S3, and S4.

Statistical Analysis
All the statistical analyses were processed using the Statistical Analysis System [39]. In the evaluation of the generated models, we considered as goodness of fit criteria the coefficient of determination (R 2 ) (Equation (8)) and the root mean square error (RMSE) (Equation (9)), in addition to the graphical analysis of the residuals.
where: y i = observed values;ŷ i = estimated values; n = number of observations; p = number of parameters.

Diametric Structure of the Forest
Our selected forest showed a reverse J-shaped diameter distribution, characterized by lower entry rate and by the presence of individuals in the smaller diametric classes (i.e., between 30 and 50 cm in diameter). The relative sampling error was 11.6% for the mean value of the volume, indicating that, in the sampling design adopted, the error was usually close to what was proposed in common forest management projects (E ar ≤ 10.0%) [40]. Thus, the characteristics of the population were inferred at this precision level.
The frequencies observed according to diametric classes ( Figure 1) indicate an unbalanced distribution. The observed basal area was 21.9 m 2 ·ha −1 , and the balanced frequency distribution for the entire rural property was predicted using Meyer's relation (Equation (1)), with the value of Liocourt's quotient (q) corresponding to 1.31. As a result, high adjustment and low error values can be noted.   (1.1), which is directly related to the number of individuals in the diametric classes. As a result, when the "q" value increases, the number of smaller trees in the forest also increases, whereas the number of larger trees decreases. Regarding the remaining basal area, the chosen values allowed the simulation of light (18.0 m 2 •ha −1 ), moderate (14.0 m 2 •ha −1 ), and strong (10.0 m 2 •ha −1 ) silvicultural intervention. The latter, when the current basal area (21.9 m 2 •ha −1 ) is reduced approximately by half, increases the number of candidate trees to be harvested in each diametric class. In summary, all these scenarios suggest possibilities of intervention in the forest while maintaining the balance of the current diametric structure.

Forest Management Scenarios Based on Liocourt's Law
The annual percentage of volume increments up to the crown, with bark (PAI%vtci), is presented in Figure 2a. It was determined according to the diametric classes, using (Equation (5)). The regression coefficients were quite significant, with R 2 = 0.4475 and RMSE = 0.5764. Figure 2b shows the graphical residual analysis as a function of the estimated values. In general, they showed a regular distribution, without the presence of bias in the predictions. Despite the low R 2 value, the model still explains almost 50% of the variation in the periodic increase in volume as a function of d. In this experiment, this value is not higher due to the lack of constancy in the increase in the growth rate of d in the last five years. However, this is a clear indication of the importance of forest management aimed at the conservation of this species.  The annual percentage of volume increments up to the crown, with bark (PAI% vtci ), is presented in Figure 2a. It was determined according to the diametric classes, using (Equation (5)). The regression coefficients were quite significant, with R 2 = 0.4475 and RMSE = 0.5764. Figure 2b shows the graphical residual analysis as a function of the estimated values. In general, they showed a regular distribution, without the presence of bias in the predictions. Despite the low R 2 value, the model still explains almost 50% of the variation in the periodic increase in volume as a function of d. In this experiment, this value is not higher due to the lack of constancy in the increase in the growth rate of d in the last five years. However, this is a clear indication of the importance of forest management aimed at the conservation of this species.

Forest Management Scenarios Based on Liocourt's Law
The mean volume per hectare found in the forest was 172.0 m 3 ·ha −1 , and the mean annual percentage of volume increments up to the crown (PAIM% vtci ) was 2.44%. Based on these values, it was possible to simulate different management scenarios (Table 1).  The mean volume per hectare found in the forest was 172.0 m 3 •ha −1 , and the mean annual percentage of volume increments up to the crown (PAIM%vtci) was 2.44%. Based on these values, it was possible to simulate different management scenarios (Table 1).   The values of Table 1 vary according to the values of "q" and Grem considered. When considering q = 1.1 and Grem = 10.0 m 2 ·ha −1 , the simulation shows that a volume of 82.7 m 3 ·ha −1 can be obtained every 27.2 years. By changing only the remaining Grem to 18.0 m 2 ·ha −1 , it is possible to obtain a volume of 11.3 m 3 ·ha −1 every 2.8 years. Thus, the optimal scenario for intervention may vary according to the forest management goals.
It is important to mention that this specific species is not explored commercially since the 90s due to the restrictive legislation [19]. Recently, there are efforts being made to approve an ongoing project for A. angustifolia in the state of Santa Catarina, to regulate the planting, preservation, sustainable management, and development of silviculture and its proven food resources [41].
It was observed that a larger volume per hectare (observed V·ha −1 ) is found between classes with d ranging from 22.5 to 57.5 cm (Figure 3a). Therefore, different scenarios keeping a q value of 1.5, and with varying Grem values tend to maintain a larger number of trees in smaller classes. Consequently, this leads to the cutting of a larger number of trees belonging to classes with higher d (Figure 3b). The opposite happens with a q value of 1.1.
The values of Table 1 vary according to the values of "q" and Grem considered. When considering q = 1.1 and Grem = 10.0 m 2 •ha −1 , the simulation shows that a volume of 82.7 m 3 •ha −1 can be obtained every 27.2 years. By changing only the remaining Grem to 18.0 m 2 •ha −1 , it is possible to obtain a volume of 11.3 m 3 •ha −1 every 2.8 years. Thus, the optimal scenario for intervention may vary according to the forest management goals.
It is important to mention that this specific species is not explored commercially since the 90s due to the restrictive legislation [19]. Recently, there are efforts being made to approve an ongoing project for A. angustifolia in the state of Santa Catarina, to regulate the planting, preservation, sustainable management, and development of silviculture and its proven food resources [41].
It was observed that a larger volume per hectare (observed V•ha −1 ) is found between classes with d ranging from 22.5 to 57.5 cm (Figure 3a). Therefore, different scenarios keeping a q value of 1.5, and with varying Grem values tend to maintain a larger number of trees in smaller classes. Consequently, this leads to the cutting of a larger number of trees belonging to classes with higher d (Figure 3b). The opposite happens with a q value of 1.1.

Classification of Wood Assortment
Considering the data shown in Table 2 for the diametric class of 27.5 cm, four trees per hectare can be harvested, with an estimated production of 1.6 m³•ha −1 . In this scenario, the trees had a total increment up to the crown of 15.3 m and 11.1 m. Therefore, 0.0% of the timber stocks' volume may be classified in assortments S1 and S2, 72.1% (corresponding to three units) in S3, and 27.9% or 0.4406 m³.cc in S4, referring to timber destined for the energy or wood industry. According to the same data, for a class center of 52.5 cm, the trees' height would be 19.1 m, with 51.1% (corresponding to one unit) of the timber stocks' volume being classified in assortment S1, 46.8% (corresponding to three units) in S2, 0% in S3, and 2.1% or 0.0437 m³•cc in S4. A total of 2.0796 m³ were thus categorized.

Classification of Wood Assortment
Considering the data shown in Table 2 for the diametric class of 27.5 cm, four trees per hectare can be harvested, with an estimated production of 1.6 m 3 ·ha −1 . In this scenario, the trees had a total increment up to the crown of 15.3 m and 11.1 m. Therefore, 0.0% of the timber stocks' volume may be classified in assortments S1 and S2, 72.1% (corresponding to three units) in S3, and 27.9% or 0.4406 m 3 ·cc in S4, referring to timber destined for the energy or wood industry. According to the same data, for a class center of 52.5 cm, the trees' height would be 19.1 m, with 51.1% (corresponding to one unit) of the timber stocks' volume being classified in assortment S1, 46.8% (corresponding to three units) in S2, 0% in S3, and 2.1% or 0.0437 m 3 ·cc in S4. A total of 2.0796 m 3 were thus categorized. Table 2. Simulation of activities in the observed/remaining stand, and of the cuts to be performed according to the wood assortment classification, with Liocourt's quotient (q) = 1.3, remaining basal area (Grem) = 14 m 2 ·ha −1 , and maximum desired diameter = 62.5 cm, for A. angustifolia (Bertol.) Kuntze in Lages, SC.

Observed Forest Stand Remaining Forest Stand Cut
Assortments

Discussion
This study simulates different forest management scenarios in a typical rural property within the MOF in which A. angustifolia is present. Several studies have shown this species' need for silvicultural intervention, since many diametric classes have already reached a growth stagnation point. Consequently, there is need to stimulate the growth of trees in lower diametric classes [42][43][44][45]. The complete lack of forest management contributes to the loss of diversity due to lower occurrence of sunflecks in the lower canopy layer and understory. This also contributes to the impoverishment of local landowners and farmers, and leads to antagonism in the conservation of this species [46]. Therefore, the choice of an appropriate remaining basal area is very important from the point of view of space optimization and vegetation recovery. It can be defined in relation to the observed frequencies per diametric class up to the maximum diameter desired [22]. In addition, the adoption of Liocourt's quotient allows identifying classes in which there is a deficit or surplus of trees [47], with the indication of interventions in the latter, and consequent increase in growth rates in the former. In this way, it is possible to regulate production without changing the forest's dynamics.
As for the variations in the values of Liocourt's quotient, these are associated with the rate of decrease in the number of trees in a class in relation to another. Thus, when the goal of management is to stimulate the development of trees in smaller classes, larger "q" values should be adopted (i.e., 1.5; Figure 1). However, if the goal is to obtain larger timber stocks in the future (e.g., S1 and 2), lower "q" values are recommended (e.g., 1.1), favoring the permanence of trees in larger d classes (Figure 1). Wu et al. [48] makes important considerations regarding the generation of large-diameter trees in the context of China's forests, and discusses the application of forest density management with appropriate thinning to meet the increasing demand for this type of tree.
Furthermore, other important implications of adopting either value of Liocourt's quotient and of Grem are associated with the logging cycle and logging intensity (Table 1). Maintaining a higher G results in shorter cutting cycles with less intensity. In contrast, a significant reduction in G leads to longer cutting cycles and more intensive interventions. Longhi et al. [20] assessed the experimental forest management of a secondary MOF in the state of Rio Grande do Sul. This site is just 250 km away from our study site. The referenced author suggested performing light selective cuts, reducing around 20% of the basal area per d class with short cutting cycles of about eight years, aiming to manage the forest with frequent cuts to induce a more productive structure.
As for the remaining volume, Figure 3 illustrates its distribution along the d classes as a function of the "q" value applied. In total, however, it is observed that the remaining volume tends to be greater with the application of a "q" value of 1.1 (Table 1). This may occur despite a larger number of trees being removed, as these are mostly thin trees.
In order to choose the ideal scenario, there are some factors to be considered. Interventions with shorter cutting cycles maintain a better basal area for every cutting (e.g., 18.0 m 2 ·ha −1 ), and offer a more regular income to landowners. In this type of forest management, because interventions are more frequent, there is also better space optimization, decreasing the competition's influence. Additionally, this type of intervention would ensure a more continuous forest coverage, promote a more conservative management, and reduce the antagonism that still exists in relation to this species in Southern Brazil [19]. Thus, the enrichment of the secondary forest of rural properties should be encouraged to increase the offer of edible seeds [49], while also taking the issue of higher interventions costs into account.
Favoring the species' natural regeneration should be especially considered when there is a Grem value that allows balancing the number of young plants [24]. In this case, it is also very important to set a higher q value by choosing to harvest a greater number of trees with larger d. This practice increases the incidence of light onto the lower canopy layer. Optimizing the value of d while recreating and maintaining aspects of the fauna's habitat [22] are also important factors in decision-making. In addition, with the possibility of producing different percentages of timber per assortment (Table 2), the regulation of production should also consider consumer market characteristics, as in some cases, opting for more intermediate (balanced) q values and maintaining Grem may be necessary as a safeguard against possible market/objective variances.
Finally, there are still many challenges to be faced when it comes to the forest management of A. angustifolia. Many of the methods applied around the world are still copies of forest management practices developed for monocultures. For this reason, advancements concerning the management of mixed forests bring great contributions to managers. The positive effects of mixed forest management have been diverse, such as greater diversity of natural regeneration [50], lower number of infestations, and reduction in the impact of insects on growth and yield [51,52], among others.
In spite of this, some studies also report the opposite, such as the increase in susceptibility to damaging biotic and abiotic agents due to the effects on the local microclimate, provision of fuel and resources to biotic and abiotic hazards, and changes in the individual physiology and development of trees [53]. These results, however, simply presents one more opportunity to develop silvicultural methods that are able to achieve forest management objectives while also minimizing the negative impacts, characterizing it as a permanent management alternative.
Many other factors still need to be considered in the context of mixed forest management, especially those related to the ecology of systems. These aspects, however, are beyond the scope of this research, and should be developed by future studies to ensure the proper implementation of the management of this species.

Conclusions
The forest fragment of the selected rural property has an inverted J tendency (observed diametric class frequency), which is very characteristic of this forest typology. The forest measurements show that it is necessary to regulate its diametric structure by cutting some trees with larger d to promote the opening of spaces and increase the incidence of light, consequently promoting the natural regeneration of this species. Silvicultural activities such as pruning and thinning ensure greater vitality to the remaining trees, improving the forest's structure and productivity, as well as the use and quality of wood assortments.
Considering the proposal of this study, it was possible to sustainably regulate the production of Araucaria wood in a natural forest located within a typical rural property in Southern Brazil. In the less intensive scenario (e.g., q value = 1.5 and larger basal area of 18.0 m 2 ·ha −1 ), there is greater space optimization, and higher economic return is ensured to the rural producer due to the definition of shorter cutting cycles (9.3 years). This also allows the faster growth of both d and h for smaller trees due to the higher incidence of light onto the lower canopy layer, thus increasing the regeneration rate of other native species.
Both the goal of forest management and consumer market characteristics need to be taken into account to define the ideal scenario. Additionally, it is necessary to consider the species' auto-ecology, its interaction with the wildlife, the site's peculiarities, among other factors not investigated in this research. Said factors should be addressed by future studies, so as to ensure the proper implementation of the forest management of this species in the Southern region of Brazil.