Abstract
Tropical forest plantations play an important role in meeting global wood demand. While research has highlighted the ecological potential of mixed-species plantations, studies on the economic viability and management of such plantations are largely missing in the context of tropical plantation forestry. In this study, we estimated the economic potential and optimized the management of commercial mixed-species plantations of four tree species native to Central America (Dalbergia retusa, Dipteryx oleifera, Hieronyma alchorneoides, and Vochysia guatemalensis) and Teak (Tectona grandis). We combined the forest growth model 3-PGmix and detailed economic data for two plantation sites in Costa Rica to optimize the management of 11 different mixtures using a genetic optimization algorithm. We found that several of the modeled mixed-species stands can be highly profitable with net present values (NPV) up to 4821.2 USD/ha at an 8% discount rate, and internal rates of return up to 17% (under excellent site conditions). This indicates that the most profitable mixtures (e.g. of V. guatemalensis-D. oleifera or T. grandis-D. oleifera on excellent sites) could compete economically with conventional monoculture plantations such as Teak monocultures. Further, mixed stands can be managed based on the same simple even-aged management approaches currently applied in monoculture plantations if the specific management parameters are adapted. The optimized management parameters also lead to improved NPV of the modeled stands under alternative valuation assumptions but are site-specific. In the present study, we only considered financial benefits from timber production. However, establishing mixed-species plantations in the tropics could provide a wide range of ecosystem services including climate change mitigation and biodiversity protection.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Tropical forest plantations play a disproportionately important role in supplying the world with wood products. In the tropics, forest plantations cover only 3% (131 million ha) of the total forested area but recently produced 46.3% of the total global and 63.7% of the tropical supply of industrial roundwood (FAO 2020; Jürgensen et al. 2014; Payn et al. 2015).
Tropical forest plantations are dominated by a handful of tree species grown in monocultures, including species of Eucalyptus spp., Acacia spp., Pinus spp., and Tectona grandis, while mixed-species plantations are restricted almost exclusively to forest plantations planted for ecological restoration or protective purposes (Gunter et al. 2013; Pancel and Köhl 2016). This is despite the fact, that researchers have repeatedly highlighted the possible silvicultural benefits of mixed-species plantations, such as increased productivity and reduced overall risk (Kelty 2006; Knoke et al. 2008; Liu et al. 2018; Messier et al. 2021). Several studies have shown that some tropical mixed-species plantations can potentially achieve higher stand productivity than their monoculture peers, if species are carefully selected (Ewel et al. 2015; Le et al. 2020; Mayoral et al. 2017; Piotto 2008; Schnabel et al. 2019). Some rare economic studies for tropical mixed-species plantations concluded that mixed-species plantations could achieve IRR between 7.37 and 15.64% in Costa Rica (Piotto et al. 2010; Streed et al. 2006). Most recently, Sinacore et al. (2022) estimated mixed-species reforestation with (partly) native tree species to yield IRR between 4 and 8% in Panama, with mixed-species stands partly outperforming their monoculture peers. Meanwhile, Trujillo-Miranda et al. (2021) found that mixed-species reforestations in Mexico would not be economically viable at any interest rate without receiving subsidies, while monocultures were viable at low interest rates (< 4%).
Despite the apparent silvicultural and economic potential of mixed-species forest plantations, the uptake by forest plantation practitioners has been limited. The present preference towards monoculture forest plantations might be driven by the simplicity of managing monoculture systems, a lack of management recommendations for mixed-species plantations, perceived higher costs of mixed-species plantations, the hesitance of giving up monoculture plantation systems with a long-proven track record of attractive financial returns, and the use of simplified economic models that are favorably biased towards monoculture plantations (Knoke et al. 2008; Nichols et al. 2006). So far, the aforementioned economic studies offer little insight and guidance into how tropical mixed-species forest plantations should be managed.
In this study, we aim to address several of these barriers to mixed-species plantations by (1) determining the economic potential of tropical mixed-species plantations, (2) optimizing plantation management regimes of mixed-species plantations, (3) determining whether forest management practices currently applied in monoculture plantations can be applied to mixed-species plantations, and (4) exploring the sensitivity of optimized management regimes and their associated economic potential to valuation assumptions, and other management constraints. To meet these research objectives, we will pursue a simulation–optimization approach. Simulation–optimization approaches are widely used to solve forest modeling problems in plantation forestry related to, for example, stand establishment, thinning, and harvesting (e.g. Jayaraman and Rugmini (2008); Jin et al. (2017); Mabvurira and Pukkala (2002); Seppänen and Mäkinen (2020)). A central part of this approach is the prediction of the effects of different silvicultural treatments through simulation, and evaluation according to the management objectives.
Just recently, Nölte et al. (2022) published a forest growth model for 5 tropical tree species in Central America. The model allows for full-rotation growth predictions of single- and mixed-species stands under varying site and climatic conditions and hence allows for the management optimization of complex mixed-species plantation systems.
Methods
To apply the simulation–optimization approach we constructed a forest economic model that can simulate the economic performance of tropical mixed-species plantations under different silvicultural treatments. We then used this model to optimize stand management. Our study was based on two plantation sites in Costa Rica and five tree species grown in mixed-species stands. This study was performed using R (R Core Team 2019).
Case study: Mixed-species plantations in Costa Rica
In the twentieth century, Costa Rica experienced a sharp decrease in forest cover from above 60% before 1960 to below 30% by the 1980s (Kleinn et al. 2002). Since then the forest cover has rebounded (FAO 2020). While only 1–2% of the land is covered by forest plantations, 75% of the harvested wood (407,046 m3 in 2014) was derived from forest plantations of mainly T. grandis and Gmelina arborea monocultures (Barrantes Rodríguez and Ugalde Alfaro 2019; REDD/CCAD-GIZ-SINAC 2015; SIREFOR 2016).
Study sites and tree species
We conducted our analysis for two sites in northern Costa Rica. The first site, “Las Delicias” is characterized by tropical moist forest conditions according to the Holdridge Life Zone System (Kohlmann et al. 2007) with an average of 2221 mm/m2 precipitation yearly, mean temperatures ranging between 21 and 35 °C and a pronounced dry period (Fig. 1). The second site, “La Virgen”, is characterized by tropical wet forest conditions with a yearly precipitation of 3549 mm/m2, mean temperatures between 20 and 32 °C, and no pronounced dry period. Both sites have previously been used as cattle pastures and partly reforested with mixed-species stands of two or three species. In the mixtures, fast-growing pioneer species are mixed tree-by-tree with intermediate and/or slow-growing species, or intermediate species are mixed with slow-growing species. We selected V. guatemalensis, H. alchorneoides, D. oleifera, Dalbergia retusa, and T. grandis as our study species. V. guatemalensis, H. alchorneoides, D. oleifera, and D. retusa are native species that are currently important timber species in the region or have been so in the past, and have been considered promising species for commercial native species plantation forestry (Cordero 2003a, 2003b, 2003c; Murillo 2018; González J 2018; SIREFOR 2021). We included T. grandis as it is the most common plantation-grown species in the region (SIREFOR 2021).
Forest economic model
The forest economic model roughly consists of four parts: forest growth modeling, timber valuation, cost and subsidies estimation, and economic performance evaluation. Figure 2 shows an overview of the model components and data flows. A detailed description of our forest economic model can be found in Supplementary Information (SI) 1.
Forest growth model
For our forest growth model, we used 3-PGmix. 3-PGmix is an extension of the 3-PG model (Physiological Principles Predicting Growth) developed by Waring and Landsberg (1997). While the original 3-PG model allows for growth predictions of single-species stands, the 3-PGmix model by Forrester and Tang (2016) was modified to allow for growth predictions of mixed-species stands through vertical and horizontal differentiation in the canopy. 3-PG calculates the total carbon fixed from absorbed photosynthetically active radiation, corrected for mediating effects from soil and climatic conditions, and stand age. To predict the growth of our species of interest we used the 3-PGmix parameters published by Nölte et al. (2022).
By default, the 3-PGmix model can make predictions for one rotation (i.e. planting to harvesting), and management options are limited to pre-determined age-based thinning and harvesting. We extended the management options to include basal area-based thinning, and target mean dbh (diameter at breast height) harvesting. Furthermore, we allowed for the replanting of trees after harvesting to model multiple rotations.
Basal area-based thinning is initiated when the modeled stand basal area (BA) exceeds a specified BA thinning threshold. The modeled stand is then thinned to a target BA. In modeled mixed-species stands, the species with the largest species-level BA is targeted for thinning. Target mean dbh harvesting is initiated when the species-level modeled mean dbh exceeds a set target dbh.
We included the option to replant harvested trees either after the modeled stand is clear-cut or underplanting trees under the remaining stand. Underplanting of harvested species was integrated by applying a stand density threshold for underplanting based on the relative spacing of the remaining stand. When the remaining modeled stand’s relative spacing after the harvest of a tree species is higher than the applied underplanting threshold (i.e. the stand density is lower) the species is replanted. Otherwise, the replanting is postponed and reevaluated at the next management intervention. We assume that the trees are pruned to a stem height of 2.5 m, 5 m, and 8 m when the tree heights exceed 5 m, 10 m, and 12 m (Pérez et al. 2003).
Timber valuation
For the timber valuation, we first calculated the merchantable volume based on the modeled volume of the harvested trees. In Costa Rica, the commercial stem height of plantation trees is usually 7 m and stems are usually sold in two sections of 3.5 m, each with a minimum diameter of 10 cm. We assumed a bark thickness of 1 cm (Berrocal et al. 2018). We applied a Weibull distribution to convert the mean dbh from the 3-PGmix output into a diameter distribution (see SI 2). For each diameter class, we calculated merchantable volume up to a stem height of 7 m. If the lower diameter of the stem pieces was below 10 cm, or if the stem pieces were not pruned we did not consider them to be merchantable.
We used the most recent timber prices reported by the Costa Rican Oficina Nacional Forestal (ONF) and local price reports (Banco Central de Costa Rica 2021; Barrantes Rodríguez and Ugalde Alfaro 2020; Dario Extra 2014; Fordaq.com 2016; Robbins 2012; Vega and Chaves 2011; World Bank 2021a).
A majority of wood harvested Costa Rica is used for pallets (47%), construction (21%), or exported as roundwood or pallets (22%) (Barrantes Rodríguez and Ugalde Alfaro 2019). In the official price reports, timber is graded by stem diameter (for roundwood), and – in some instances – by origin (plantation or natural forest) (Barrantes Rodríguez and Ugalde Alfaro 2020). Timber prices are mainly influenced by two important quality characteristics: (1) stem diameter, and (2) stem heartwood content (Barrantes Rodríguez and Ugalde Alfaro 2020; Pérez and Kanninen 2005, 2003). For each modeled stem piece, we incorporated price corrections to address these quality characteristics. A detailed description of the timber valuation process and price corrections can be found in SI 1.
Costs and subsidies
We retrieved establishment and management costs from Vallejos (2019). The detailed cost data originate from 12 mid-size to large Costa Rican plantation forest enterprises and were collected from 2008 until today (Guevara-Bonilla 2011; Guevara-Bonilla and Murillo 2021). We applied fertilization and weeding costs in the first 5 years after planting. Costs for tilling and fertilization during stand establishment were only applied during the first rotation. As the modeled stems are sold standing, we only applied costs for tree selection and marking during harvesting. Finally, we assumed that the value of the modeled wood harvested must exceed 1000 USD/ha at stumpage for the operation to be profitable for the harvesting company. If this minimum threshold was not met, the difference was added as a harvesting cost.
Finally, we applied species-dependent reforestation subsidies offered by the state during the first 5 years of the forest establishment between a total sum of 1239.2–2096 USD/ha (MINAE 2020). A detailed overview of all costs and subsidies can be found in SI 1.
Economic performance
Based on the cashflows of the modeled stands, we calculated the net present value (NPV) for one planning period of 80 years (i.e. 2020–2100). We used a discount rate of 8% (Cubbage et al. 2020, 2007). NPV is the sum of the discounted cash flows of an asset with a definite planning horizon. We also calculated the IRR. IRR is the discount rate at which the NPV is 0 and is commonly used to communicate the economic performance of forest assets. We calculated NPV using Eq. 1. IRR was determined using the R package “jrvFinance” (Varma, 2021).
where T is the total number of simulation periods (in years), Rt and Ct are all revenues and costs USD/ha) occurring in each simulation period t, and r is the discount rate.
Baseline scenarios for optimization
Planting schemes
We simulated the study species grown in mixed-species stands. The modeled stands were planted at 800 stems/ha, where up to three species of different growth characteristics are mixed within a single modeled stand. An overview of the planting schemes modeled on each site can be found in Table 1. T. grandis and D. retusa need a seasonal dry period to thrive and were hence not planted in La Virgen.
Site fertility
In 3-PG, site quality is represented in terms of the Fertility Rating (FR) parameter. FR is a relative measure of a site's species-specific productivity adjusted for soil water and climatic growth influences. For our selected species, site preferences are still partly poorly understood, and sensible generalized site fertility scenarios challenging to formulate. Therefore, we used the species-specific FR model-fitted for our study sites by Nölte et al. (2022) (Table 2). Nölte et al. (2022) included FR as a parameter to be determined during the model calibration, i.e. FR was determined such, that the modeled growth data best fitted the observed data on the study sites. Similar approaches were applied in e.g. Forrester and Tang (2016) and Gonzalez-Benecke et al. (2014). It should be noted, that 3-PG growth predictions are highly sensitive to the chosen FR. We will address this issue further in the discussion.
Baseline management scenarios
To assess whether our optimization leads to improved performance we formulated 4 baseline management scenarios (Table 3) by reviewing recent literature on the management of tropical timber plantations (Griess and Knoke 2011; Ladrach 2004; Nölte et al. 2018; Onyekwelu et al. 2011; Pérez and Kanninen 2005). In two of the management scenarios, we applied age-based thinning and harvesting, while we applied basal area-based thinning in the other two.
Management optimization
We aimed to determine the thinning and harvesting regimes that offer the highest NPV. We optimized the modeled stand NPV based on two different thinning and harvesting approaches: (1) an age-based approach, and (2) an approach based on stand growth performance. In the age-based approach, we optimized thinning ages and stem numbers, and optimal harvesting ages. We allowed for 3 thinnings in mixed-species stands of two species, and 4 thinnings in mixed-species stands of three species (Eq. 2). In the second approach, we optimized the upper BA threshold for thinning, the target BA after thinning, and the target mean dbh for harvesting (Eq. 3). We further included the stand density threshold for underplanting as a parameter for optimization.
where TA is the applied thinning ages, TS is the applied number of stems remaining after thinning, HA is the applied harvesting ages, and RS is the applied relative spacing threshold for underplanting.
where BA1 is the applied upper basal area threshold to initiate thinning, BA2 is the applied target basal area after thinning, TD is the applied target mean dbh.
We relied on a genetic optimization algorithm using Differential Evolution (Mullen et al. 2011). A detailed description of the optimization and the applied parameter bounds can be found in SI 1. We ran the optimization for each scenario with 100 iterations.
Sensitivity analysis
To assess whether the optimized management parameters lead to improvements compared to the baseline management scenarios under alternative valuation and site assumptions we applied 7 alternative sensitivity scenarios (Table 4).
Further, to maintain the mixed-species character of the modeled stands we introduced a revenue share constraint. The constraint prescribes that each species must contribute a certain minimum percentage of the modeled stand-level discounted revenues from timber sales. We applied the three different constraint scenarios so that each species in a mixture contributes at least 5%, 10%, and 20% of the discounted revenue from timber sales. The constraint ensures that at least some trees of all species reach merchantable sizes. We reran the optimization under the application of the revenue share constraints.
Results
We ran optimizations to find management parameters for thinning and harvesting that would maximize the economic performance of the studied planting schemes. We applied both an age-based management approach and a management approach based on growth performance (i.e. BA and dbh).
Economic performance and productivity
First, we found that all mixtures except H. alchorneoides-D. retusa were economically viable under optimized management on our two study sites at a discount rate of 8%. The most profitable mixtures could support required rates of return above 15%.
Figure 3 shows (A) the discounted cash flows and NPV for the studied mixed-species plantations when applying the optimized growth performance-based management parameters (see also Tables 5 and 6), (B) the mean annual increments (MAI) of stem biomass, and (C) merchantable stem volume. The performance outcomes under optimized age-based management were largely similar (but slightly lower in pairwise comparison) and can be found in SI 3.
On average an NPV of 2640.1 USD/ha was achieved across all mixtures. The mixture of V. guatemalensis-D. oleifera (4821.2 USD/ha and 3433.4 USD/ha) shows the highest NPV of all mixtures in La Virgen and Las Delicias respectively. IRR varied between 4.7 and 16.7%.
Most mixtures accumulated stem biomass and merchantable volume at a rate of 10–15 tons/ha/year and 4–22 m3/ha/year respectively. In Figs. 2B and C it is evident that under (unconstrained) optimized management a single species could dominate the mixed stands in terms of biomass and volume, effectively transforming the modeled stands into single-species stands.
Optimized management parameters
Secondly, we found that the optimized management parameters varied strongly between the different mixtures. Generally, under optimized management, the mixed stands were managed in conventional rotations of even-aged trees (i.e. no underplanting takes place). The modeled stands were kept below a BA of 16–31 m2/ha, with higher stand densities in stands where V. guatemalensis is dominant, and lower stand densities in mixtures with D. oleifera and T. grandis (Tables 5 and 6). Between 150–300 trees/ha of each species are maintained until the end harvest is reached. Due to space constraints, we report the optimized age-based parameters for age-based thinning and harvesting in SI 3.
Thinning
When examining the BA thinning thresholds roughly 5 distinguishable thinning schemes emerged: (1) in mixtures where D. oleifera was dominant (i.e. when mixed with D. retusa or H. alchorneoides) thinning was initiated when the modeled stand-level BA exceeded 16.1–18.8 m2/ha and the modeled stand was thinned to 14–16.5 m2/ha (i.e. BA of the dominant species was lowered by 1.8–2.3 m2/ha during thinning); (2) in mixtures where V. guatemalensis was dominant (i.e. when mixed with D. retusa and/or H. alchorneoides) thinning was initiated when stand-level BA exceeded 28.1–30.6 m2/ha and the stand was thinned to 23.3–27.1 m2/ha (i.e. BA of the dominant species was lowered by 3.5–5.3 m2/ha); (3) when V. guatemalensis occured in a mixture with D. oleifera the stand was thinned from 19.9–23.9 to 13.4–19 m2/ha (4.8–6.5 m2/ha of the dominant species were removed); (4) when T. grandis was mixed with D. oleifera the stand was thinned from 18.6 m2/ha to 13.2 m2/ha, and (5) when T. grandis was mixed with D. retusa the stand was thinned from 22.9 to 17.1 m2/ha.
Harvesting
Target mean dbh varied between 13.1 and 54.8 cm for D. retusa, 25.7–35.1 cm for D. oleifera, 10.1–23 cm for H. alchorneoides, 34.6–36.9 cm for T. grandis, and 36.2–43.8 cm for V. guatemalensis. However, it should be noted, that the large target mean dbh around 50 cm or higher for D. retusa was the result of high competition-related mortality (in 3-PG small trees die first). As competition-related mortality removed individuals from the lower part of the diameter distribution, the mean dbh was shifted upwards.
Underplanting
Underplanting occurred in 7 of the 15 optimized mixtures. The minimum relative spacing threshold for underplanting varies between 0.15 and 3.96 For reference, with a minimum relative spacing threshold of 1, the mean tree height of the remaining stand cannot exceed 7 m if the stocking is above 200 tree/ha, or 14 m if the stocking is above 50 tree/ha. Accordingly, underplanting was likely to occur when only D. retusa remains in a stand (D. retusa generally showed a slow height growth and low maximum height), or in mixtures of V. guatemalensis-H. alchorneoides. Otherwise, replanting only occurred after clear-cut events.
Improvements and sensitivity analysis
Finally, we found that applying the optimized management parameters lead to improvements in the economic performance of all mixtures compared to the best-performing baseline management scenarios (under sensitivity scenario “MAIN”) (Fig. 4). The optimized management parameters also lead to improvements in NPV in all valuation scenarios except H. alchorneoides-D. retusa under “V.LIST”. However, when alternative site fertility ratings were applied the optimized management parameters did not consistently improve the economic performance.
Maintaining the mixed-species character of stands
When enforcing a minimum revenue share of each admixed species, the economic performance decreased consistently for all mixtures (Fig. 5). On average, increasing a species' minimum revenue share by 1 percentage point decreased the stand-level NPV 119 USD/ha.
Discussion
In this study, we aimed to determine the economic potential and optimized management of mixed-species plantations in the tropics based on a case study in Costa Rica. To our knowledge, we are the first to combine detailed economic and operational data with a process-based forest growth model that can integrate species-interaction effects in the context of tropical mixed-species plantation management.
Economic potential and management of tropical mixed-species plantations
We found, that most of the modeled planting schemes were profitable at a discount rate of 8%, with IRR ranging between 9.9 and 16.7%. Only the mixture H. alchorneoides-D. retusa was not economically viable, presumably due to poor site suitability. Earlier studies had found that mixed-species plantations in Central America could yield internal rates of return 4–15.6% (Piotto et al. 2010; Sinacore et al. 2022; Streed et al. 2006). For monocultures of selected native species and T. grandis in Panama, Griess and Knoke (2011) and Paul et al. (2015) estimated IRR of 10–15% and 2.8–12.6% respectively. Cubbage et al. (2020) in reviewing the economic performance of forest plantations across the globe IRR of tropical species in Central and South America between 8.1 and 21.8%, with plantations of Gm. arborea and T. grandis (the two most important plantations species in Costa Rica) reaching IRR up to 24.5% and 18.6% respectively. Our findings are within the range of previous economic performance estimates for tropical mixed-species plantations and indicate that mixed-species plantations can be economically competitive with conventional monocultures. However, our estimated IRR and NPV might be considered as upper boundaries for the economic performance of mixed-species plantations with the selected species, as applied fertility ratings for D. oleifera and T. grandis are exceptionally high and we assumed the mixtures are managed to achieve their optimal performance. On less fertile sites or under suboptimal management the economic performance could be considerably lower.
The second objective of our study was to determine how mixed-species plantations should be managed to achieve optimal economic outcomes. Our results indicate that both age-based and growth performance-based management approaches can achieve similar economic outcomes. Our optimized thinning parameters were highly variable from mixture to mixture with upper BA thinning thresholds ranging from 16.1 to 30.6 m2/ha, target BA after thinning between 8.6 and 27.1 m2/ha, and BA being lowered by 1.8–7.5 m2/ha per thinning. The optimized BA thinning parameters reflect the typical dichotomy between maximizing stand-level timber volume and maximizing individual tree stem volume. In our study, this dichotomy is partly driven by the integration of heartwood content into the valuation process. D. oleifera, D. retusa, H. alchorneoides, and T. grandis all produce distinguishable (and valuable) heartwood. As heartwood content increases with increasing stem diameter, optimized management for mixtures with only these four species will tend to favor individual tree growth by maintaining relatively low stand densities (upper BA threshold between 16.1 and 22.3 m2/ha with 1.8–7.5 m2/ha removed per thinning). V. guatemalensis does not produce distinguishable heartwood and optimized management favored maintaining higher stand densities in mixtures where V. guatemalensis was dominant (upper BA threshold between 28.4 and 30.6 m2/ha with 3.5–5.3 m2/ha removed per thinning). In mixtures where V. guatemalensis occurred with D. oleifera (both species have a high economic potential) a midway was found (upper BA threshold between 19.9 and 23.9 m2/ha with 4.8–6.5 m2/ha removed per thinning). Pérez and Kanninen (2005) in a comprehensive study for T. grandis applied BA thresholds for thinning 20–23 m2/ha to 12–17 m2/ha (the upper threshold for thinning and target BA after thinning respectively) when high individual tree growth is the management objective, and 22–25 m2/ha to 13–18 m2/ha when high stand growth is the main management objective.
Target mean dbh for harvesting generally varied between 26 and 44 cm for D. oleifera, T. grandis, and V. guatemalensis. These values lay within the range of diameters at harvest for plantation-grown species reported in other studies (Griess and Knoke 2011; Pérez and Kanninen 2005). D. retusa had a target mean dbh between 15 and 20 cm (when correcting for mortality-induced shifts of the mean dbh). The species is mainly used for artisanal work or fine carpentry and even wood pieces of small dimensions might fetch considerable prices. Due to its slow growth, achieving a target dbh of 20 cm can take more than 30 years in a mixed-species plantation.
Underplanting of harvested tree species generally did not occur under optimized management. Instead, replanting occurred after clear-cutting.
Overall, it should be noted, that the optimization approach we applied in this study approximates the global optimum of the objective function (i.e. NPV). However, whether the global optimum is reached at the end of the optimization is unknown. Often, (forest) economic optimization studies prefer to use models that can be optimized mathematically. This ensures both that the global optimum of the optimization problem is reached and allows for detailed static analyses but requires a substantial simplification of complex real-world biological and economic processes (Knoke et al. 2008).
Sensitivity and uncertainties
The results from our sensitivity analysis showed that using our optimized management parameters improved economic performance consistently also when alternative valuation assumptions are applied. However, under alternative fertility scenarios, the optimized parameters did not consistently lead to improvements in the modeled stands' economic performance, compared to the baseline management scenarios. This highlights that the optimized management parameters are site-specific, i.e. specific to the applied FR. On sites with diverging FR, the optimized management parameters from this study might not applicable, and the economic potential of such sites might diverge substantially from the ones estimated in this study. This is demonstrated in the sensitivity analysis, where only 4 of 15 modeled stands are economically viable under generic median site fertility (all FR = 0.5). This is especially important for the species D. oleifera and T. grandis which showed exceptionally high FR on the two study sites. On poorer sites, reaching the target mean dbh for harvesting might take much longer than suggested in this study (when applying a growth performance-based management strategy), or might not be reached within the planning period at all. Generally, poorer sites can be expected to have longer optimal rotation times (due to the lower opportunity costs for delaying cashflows from future rotations) and smaller target mean dbh for harvesting. For broader applications, the optimized management parameters need to be generalized in future research. The data underlying the calibration of the growth model used in this study included examples of commercial plantations with very poor FR (between 0 and 0.3) (Nölte et al. 2022), where stand growth stagnated with increasing stand age (e.g. T. grandis growth stagnated at a mean dbh around 20 cm). Such low FR can either indicate nutrient-poor soils or unsuitable species selection. While our model allows for the optimization of such stands, they are unlikely to be economically viable regardless of the applied management parameters. With improved knowledge of species-specific site requirements, more site-suitable species should be selected for the reforestation of such sites in the future.
Secondly, the applied timber prices in our study remain an important source of uncertainty. For D. oleifera and D. retusa plantation-grown wood is not yet available on the market and prices hence refer to wood from natural forests. While we corrected timber prices for two important quality characteristics of tropical high-value timber (heartwood and stem diameter), it remains to be seen whether these price estimates materialize once plantation-grown wood from D. oleifera and D. retusa enters the market.
Finally, changing interest rates pose a crucial uncertainty in our study due to the long planning horizon of 80 years. As indicated by the IRR, 4 of 15 mixtures would remain economically viable even if interest rates increased to 15%.
Potential beyond timber production
In the present study, we assessed mixed-species plantations exclusively based on their ability to generate revenues through timber sales. However, other ecosystem services are becoming increasingly important objectives of plantation establishment and management.
For example, reforestation in the tropics has been suggested as a powerful strategy for combating climate change through carbon sequestration (Bastin et al. 2019; Griscom et al. 2017). However, past reforestation efforts have been criticized for the excessive establishment of exotic species monocultures with a low long-term carbon sequestration potential (Lewis et al. 2019). Establishing mixed-species reforestations has been suggested as an alternative to such monocultures in carbon plantings (Hulvey et al. 2013). Further, mixing species in commercial plantations might also be used as a risk mitigation strategy against climate change uncertainties, such as shifts in tree species distributions and site suitability (Allen et al. 2010; Anadón et al. 2014; Baumbach et al. 2021; Messier et al. 2021).
In general, mixed-species plantations are expected to perform better than monocultures when multiple ecosystem services and biodiversity protection are taken into account as management objectives (Carnus et al. 2006; Di Sacco et al. 2021; Hulvey et al. 2013; Liang et al. 2016; Lindenmayer et al. 2003; Stephens and Wagner 2007).
Conclusion
Overall, this study found that mixed-species plantations in Costa Rica could be economically viable and even highly profitable if managed appropriately and planted on adequate sites. Even-aged management approaches currently applied in mono-specific plantations could be applied successfully to mixed-species plantations if management thresholds and target values are adapted to individual species combinations and site conditions. Interestingly, also tree species currently grown almost exclusively in monocultures, such as T. grandis, showed a high economic potential when grown in mixed stands. Taking a broader selection of management objectives into account—such as climate change and risk mitigation, and biodiversity protection—might highlight further advantages of establishing mixed-species plantations as an alternative to current monoculture practices in tropical plantation forestry.
References
Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, Kitzberger T, Rigling A, Breshears DD, Hogg EH, Ted Gonzalez P, Fensham R, Zhang Z, Castro J, Demidova N, Lim JH, Allard G, Running SW, Semerci A, Cobb N (2010) A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Ecol. Manage. 259:660–684. https://doi.org/10.1016/j.foreco.2009.09.001
Anadón JD, Sala OE, Maestre FT (2014) Climate change will increase savannas at the expense of forests and treeless vegetation in tropical and subtropical Americas. J Ecol 102:1363–1373. https://doi.org/10.1111/1365-2745.12325
Banco Central de Costa Rica, 2021. US Dollar purchase and sale rates [WWW Document]. URL https://gee.bccr.fi.cr/indicadoreseconomicos/cuadros/frmvercatcuadro.aspx?CodCuadro=400&Idioma=2&FecInicial=2020/01/01&FecFinal=2020/12/31&Filtro=0 (accessed 9.9.21).
Barrantes Rodríguez A, Ugalde Alfaro S (2019) Usos y aportes de la madera en Costa Rica. San José, Costa Rica.
Barrantes Rodríguez A, Ugalde Alfaro S (2020) Precios de la madera para las especies mas comercializadas -Primer semestra del 2020. San José, Costa Rica.
Bastin J-F, Finegold Y, Garcia C, Mollicone D, Rezende M, Routh D, Zohner CM, Crowther TW (2019) The global tree restoration potential. Science. 365(80):76–79. https://doi.org/10.1126/science.aax0848
Baumbach L, Warren DL, Yousefpour R, Hanewinkel M (2021) Climate change may induce connectivity loss and mountaintop extinction in Central American forests. Commun Biol 4:1–13. https://doi.org/10.1038/s42003-021-02359-9
Berrocal A, Gaitan-Alvarez J, Moya R, Fernández-Sólis D, Ortiz-Malavassi E (2018) Development of heartwood, sapwood, bark, pith and specific gravity of teak (Tectona grandis) in fast-growing plantations in Costa Rica. J for Res. https://doi.org/10.1007/s11676-018-0849-5
Carnus JM, Parrotta J, Brockerhoff E, Arbez M, Jactel H, Kremer A, Lamb D, O’Hara K, Walters B (2006) Planted forests and biodiversity. J for 104:65–77. https://doi.org/10.1093/jof/104.2.65
J Cordero (2003a). Dipteryx oleifera. In: J Cordero, D Boshier (Eds.), Arboles de Centroamerica. OFI-CATIE.
JCordero (2003b). Dalbergia retusa. In: J Cordero, D Boshier (Eds.), Arboles de Centroamerica. OFI-CATIE.
JCordero (2003c) Hyeronima alchorneoides. In: J Cordero, D Boshier (Eds.), Arboles de Centroamerica. OFI-CATIE.
Cubbage F, Mac Donagh P, Sawinski J, Rubilar R, Donoso P, Ferreira A, Hoeflich V, Olmos VM, Ferreira G, Balmelli G, Siry J, Báez MN, Alvarez J (2007) Timber investment returns for selected plantations and native forests in South America and the southern United States. New for 33:237–255. https://doi.org/10.1007/s11056-006-9025-4
Cubbage F, Kanieski B, Rubilar R, Bussoni A, Olmos VM, Balmelli G, Donagh PM, Lord R, Hernández C, Zhang P, Huang J, Korhonen J, Yao R, Hall P, Del La Torre R, Diaz-Balteiro L, Carrero O, Monges E, Thu HTT, Frey G, Howard M, Chavet M, Mochan S, Hoeflich VA, Chudy R, Maass D, Chizmar S, Abt R (2020) Global timber investments, 2005 to 2017. For Policy Econ 112:102082. https://doi.org/10.1016/j.forpol.2019.102082
Di Sacco A, Hardwick KA, Blakesley D, Brancalion PHS, Breman E, Cecilio Rebola L, Chomba S, Dixon K, Elliott S, Ruyonga G, Shaw K, Smith P, Smith RJ, Antonelli A (2021) Ten golden rules for reforestation to optimize carbon sequestration, biodiversity recovery and livelihood benefits. Glob Chang Biol. https://doi.org/10.1111/gcb.15498
Ewel JJ, Celis G, Schreeg L (2015) Steeply Increasing Growth Differential Between Mixture and Monocultures of Tropical Trees. Biotropica 47:162–171. https://doi.org/10.1111/btp.12190
Dario Extra, 2014. Policía confirma existencia de mafia del cocobolo [WWW Document]. URL https://www.diarioextra.com/Noticia/detalle/243440/policia-confirma-exi (accessed 9.9.21).
FAO (2020) Global Forest Resources Assessment 2020. FAO. https://doi.org/10.4060/ca9825en
Fordaq.com, 2016. VENTA DE MADERA COCOBOLO DALBERGIA RETUSA [WWW Document]. URL https://madera.fordaq.com/srvAuctionView.html?AucTIid=18056636 (accessed 9.9.21).
Forrester DI, Tang X (2016) Analysing the spatial and temporal dynamics of species interactions in mixed-species forests and the effects of stand density using the 3-PG model. Ecol Modell 319:233–254. https://doi.org/10.1016/j.ecolmodel.2015.07.010
Forrester DI (2020) 3-PG User Manual. Birmensdorf, Switzerland
Murillo O (2018) ¿Cuáles especies forestales nativas debemos priorizar en el país? Ambientico 267:4–9
González J, E., (2018) Especies forestales nativas para la reforestación comercial en las principales biorregiones de Costa Rica. Ambientico 267:17–22
Gonzalez-Benecke CA, Jokela EJ, Cropper WP, Bracho R, Leduc DJ (2014) Parameterization of the 3-PG model for Pinus elliottii stands using alternative methods to estimate fertility rating, biomass partitioning and canopy closure. For Ecol Manage 327:55–75. https://doi.org/10.1016/j.foreco.2014.04.030
Griess VC, Knoke T (2011) Can native tree species plantations in Panama compete with Teak plantations? An Economic Estimation New for 41:13–39. https://doi.org/10.1007/s11056-010-9207-y
Griscom BW, Adams J, Ellis PW, Houghton RA, Lomax G, Miteva DA, Schlesinger WH, Shoch D, Siikamäki JV, Smith P, Woodbury P, Zganjar C, Blackman A, Campari J, Conant RT, Delgado C, Elias P, Gopalakrishna T, Hamsik MR, Herrero M, Kiesecker J, Landis E, Laestadius L, Leavitt SM, Minnemeyer S, Polasky S, Potapov P, Putz FE, Sanderman J, Silvius M, Wollenberg E, Fargione J (2017) Natural climate solutions. Proc Natl Acad Sci U S A 114:11645–11650. https://doi.org/10.1073/pnas.1710465114
Guevara-Bonilla M (2011) Productividad y costos operacionales de la chapea mecanizada en plantaciones jóvenes de Acacia mangium en la zona norte de Costa Rica. Rev for Mesoam Kurú 8:27–30
Guevara-Bonilla M, Murillo O (2021) Productividad y costos del primer raleo en plantaciones de Acacia mangium Willd en la zona norte de Costa Rica. Rev Mesoam Kurú 16:55–61. https://doi.org/10.18845/rfmk.v16i42.5539
Gunter S, Weber M, Stimm B, and MR (2013) Silviculture in the Tropics.
Hulvey KB, Hobbs RJ, Standish RJ, Lindenmayer DB, Lach L, Perring MP (2013) Benefits of tree mixes in carbon plantings. Nat Clim Chang 3:869–874. https://doi.org/10.1038/nclimate1862
Hung TT, Almeida AC, Eyles A, Ratkowsky D, Lam VT, Mohammed C (2019) Maximising growth and sawlog production from Acacia hybrid plantations in Vietnam. New for 50:785–804. https://doi.org/10.1007/s11056-018-09699-1
Jayaraman K, Rugmini P (2008) Optimizing management of even-aged teak stands using a growth simulation model: A case study in Kerala. J Trop for Sci 20:19–28
Jin X, Pukkala T, Li F, Dong L (2017) Optimal management of Korean pine plantations in multifunctional forestry. J for Res 28:1027–1037. https://doi.org/10.1007/s11676-017-0397-4
Jürgensen C, Kollert W, Lebedys A, FAO Jim Carle, Peter Holmgren (2014) Planted Forests of the Future: Sustainable Management Pathways to achieve ecosystem restoration, 1–14.
Kelty MJ (2006) The role of species mixtures in plantation forestry. For Ecol Manage 233:195–204. https://doi.org/10.1016/j.foreco.2006.05.011
Kleinn C, Corrales L, Morales D (2002) Forest area in costa rica: a comparative study of tropical forest cover estimates over time. Environ Monit Assess 73:17–40
Knoke T, Ammer C, Stimm B, Mosandl R (2008) Admixing broadleaved to coniferous tree species: a review on yield, ecological stability and economics. Eur J for Res 127:89–101. https://doi.org/10.1007/s10342-007-0186-2
Kohlmann B, Solís Á, Elle O, Soto X, Russo R (2007) Biodiversity, conservation, and hotspot atlas of costa rica: a dung beetle perspective (coleoptera: scarabaeidae: scarabaeinae). Zootaxa 1457(1):1–34
Ladrach WE (2004) Harvesting and comparative thinning alternatives in Gmelina arborea plantations. New for 28:255–268. https://doi.org/10.1023/B:NEFO.0000040952.31476.03
Le HD, Smith C, Herbohn J, Nguyen H (2020) A comparison of growth, structure and diversity of mixed species and monoculture reforestation systems in the philippines a comparison of growth, structure and diversity of mixed. J Sustain for 00:1–30. https://doi.org/10.1080/10549811.2020.1767145
Lewis S, Wheeler CE, Mitchard ETA, Koch A (2019) Regenerate natural forests to store carbon. Nature 568:25–28
Liang J, Crowther TW, Picard N, Wiser S, Zhou M, Alberti G, Schulze ED, McGuire AD, Bozzato F, Pretzsch H, De-Miguel S, Paquette A, Hérault B, Scherer-Lorenzen M, Barrett CB, Glick HB, Hengeveld GM, Nabuurs GJ, Pfautsch S, Viana H, Vibrans AC, Ammer C, Schall P, Verbyla D, Tchebakova N, Fischer M, Watson JV, Chen HYH, Lei X, Schelhaas MJ, Lu H, Gianelle D, Parfenova EI, Salas C, Lee E, Lee B, Kim HS, Bruelheide H, Coomes DA, Piotto D, Sunderland T, Schmid B, Gourlet-Fleury S, Sonké B, Tavani R, Zhu J, Brandl S, Vayreda J, Kitahara F, Searle EB, Neldner VJ, Ngugi MR, Baraloto C, Frizzera L, Bałazy R, Oleksyn J, Zawiła-Niedźwiecki T, Bouriaud O, Bussotti F, Finér L, Jaroszewicz B, Jucker T, Valladares F, Jagodzinski AM, Peri PL, Gonmadje C, Marthy W, O’Brien T, Martin EH, Marshall AR, Rovero F, Bitariho R, Niklaus PA, Alvarez-Loayza P, Chamuya N, Valencia R, Mortier F, Wortel V, Engone-Obiang NL, Ferreira LV, Odeke DE, Vasquez RM, Lewis SL, Reich PB (2016) Positive biodiversity-productivity relationship predominant in global forests. Science. https://doi.org/10.1126/science.aaf8957
Lindenmayer DB, Hobbs RJ, Salt D (2003) Plantation forests and biodiversity conservation. Aust for 66:62–66
Liu CLC, Kuchma O, Krutovsky KV (2018) Mixed-species versus monocultures in plantation forestry: Development, benefits, ecosystem services and perspectives for the future. Glob Ecol Conserv 15:e00419. https://doi.org/10.1016/j.gecco.2018.e00419
Mabvurira D, Pukkala T (2002) Optimising the management of Eucalyptus grandis (Hill) Maiden plantations in Zimbabwe. For Ecol Manage 166:149–157. https://doi.org/10.1016/S0378-1127(01)00666-1
Mayoral C, van Breugel M, Cerezo A, Hall JS (2017) Survival and growth of five Neotropical timber species in monocultures and mixtures. For Ecol Manage 403:1–11. https://doi.org/10.1016/j.foreco.2017.08.002
Messier C, Bauhus J, Sousa-Silva R, Auge H, Baeten L, Barsoum N, Bruelheide H, Caldwell B, Cavender-Bares J, Dhiedt E, Eisenhauer N, Ganade G, Gravel D, Guillemot J, Hall JS, Hector A, Hérault B, Jactel H, Koricheva J, Kreft H, Mereu S, Muys B, Nock CA, Paquette A, Parker JD, Perring MP, Ponette Q, Potvin C, Reich PB, Scherer-Lorenzen M, Schnabel F, Verheyen K, Weih M, Wollni M, Zemp DC (2021) For the sake of resilience and multifunctionality, let’s diversify planted forests! conserv. Letter. https://doi.org/10.1111/conl.12829
MINAE (2020) Resolución Ministerial: R-055–2020-MINAE. MINISTERIO DE AMBIENTE Y ENERGÍA, DIRECCIÓN DE ASESORÍA JURÍDICA.
Moya R, Bond B, Quesada H (2014) A review of heartwood properties of Tectona grandis trees from fast-growth plantations. Wood Sci Technol 48:411–433. https://doi.org/10.1007/s00226-014-0618-3
Mullen KM, Ardia D, Gil DL, Windover D, Cline J (2011) DEoptim: An R package for global optimization by differential evolution. J Stat Softw 40:1–26. https://doi.org/10.18637/jss.v040.i06
Nichols JD, Bristow M, Vanclay JK (2006) Mixed-species plantations: Prospects and challenges. For Ecol Manage 233:383–390. https://doi.org/10.1016/j.foreco.2006.07.018
Nölte A, Meilby H, Yousefpour R (2018) Multi-purpose forest management in the tropics: Incorporating values of carbon, biodiversity and timber in managing Tectona grandis (teak) plantations in Costa Rica. For Ecol Manage 422:345–357. https://doi.org/10.1016/j.foreco.2018.04.036
Nölte A, Yousefpour R, Cifuentes-Jara M, Piotto D, Murillo O, Zúñiga P, Hanewinkel M (2022) Broad-scale and long-term forest growth predictions and management for native, mixed species plantations and teak in Costa Rica and Panama. For Ecol Manage 520:120386. https://doi.org/10.1016/j.foreco.2022.120386
Onyekwelu JC, Stimm B, Evans J (2011) Plantation Forestry (Review). In: Günter S, Weber M, Stimm B, Mosandl R (eds) Silviculture in the Tropics. Springer, London, pp 399–454
Pancel L, Köhl M (2016) Tropical forestry handbook, 2nd edn. Springer, Berlin Heidelberg, pp 204–206
Paul C, Griess VC, Havardi-Burger N, Weber M (2015) Timber-based agrisilviculture improves financial viability of hardwood plantations: a case study from Panama. Agrofor Syst 89:217–235. https://doi.org/10.1007/s10457-014-9755-9
Payn T, Carnus JM, Freer-Smith P, Kimberley M, Kollert W, Liu S, Orazio C, Rodriguez L, Silva LN, Wingfield MJ (2015) Changes in planted forests and future global implications. For Ecol Manage 352:57–67. https://doi.org/10.1016/j.foreco.2015.06.021
Pérez D, Kanninen M (2005) Stand growth scenarios for Tectona grandis plantations in Costa Rica. For Ecol Manage 210:425–441. https://doi.org/10.1016/j.foreco.2005.02.037
Pérez LD, Viquez E, Kanninen M (2003) Preliminary pruning programme for Tectona grandis plantations in Costa Rica. J Trop for Sci 15:557–569
Pérez LD, Kanninen M (2003) Heartwood, sapwood and bark content, and wood dry density of young and mature teak (Tectona grandis) trees grown in Costa Rica. Silva Fenn 37:45–54. https://doi.org/10.14214/sf.511
Piotto D (2008) A meta-analysis comparing tree growth in monocultures and mixed plantations. For Ecol Manage 255:781–786. https://doi.org/10.1016/j.foreco.2007.09.065
Piotto D, Craven D, Montagnini F, Alice F (2010) Silvicultural and economic aspects of pure and mixed native tree species plantations on degraded pasturelands in humid costa rica. New for 39:369–385. https://doi.org/10.1007/s11056-009-9177-0
Pretzsch H (2019) Grundlagen der Waldwachstumsforschung, Grundlagen der Waldwachstumsforschung. Springer, Berlin Heidelberg
R Core Team, 2019. R: A language and environment for statistical computing.
REDD/CCAD-GIZ-SINAC, 2015. Inventario Nacional Forestal de Costa Rica 2014–2015. San José, Costa Rica.
Rizanti DE, Darmawan W, George B, Merlin A, Dumarcay S, Chapuis H, Gérardin C, Gelhaye E, Raharivelomanana P, Kartika Sari R, Syafii W, Mohamed R, Gerardin P (2018) Comparison of teak wood properties according to forest management: short versus long rotation. Ann for Sci. https://doi.org/10.1007/s13595-018-0716-8
Robbins, S., 2012. Altos precios de exótica madera impulsan tala ilegal en Panamá [WWW Document]. InSight Crime. URL https://es.insightcrime.org/noticias/noticias-del-dia/altos-precios-de-exotica-madera-impulsan-a-la-tala-ilegal-en-panama/ (accessed 9.9.21).
Schnabel F, Fichtner A, Schwarz JA, Nock CA, Dănescu A, Bauhus J, Potvin C (2019) Drivers of productivity and its temporal stability in a tropical tree diversity experiment. Glob Chang Biol. https://doi.org/10.1111/gcb.14792
Seppänen, P., Mäkinen, A., 2020. Comprehensive yield model for plantation teak in Panama. Silva Fenn. 54, 1–25. https://doi.org/10.14214/sf.10309
Sinacore K, García EH, Howard T, van Breugel M, Lopez OR, Finkral AJ, Hall JS (2022) Towards effective reforestation: growth and commercial value of four commonly planted tropical timber species on infertile soils in Panama. New for. https://doi.org/10.1007/s11056-022-09906-0
SIREFOR (2016) REPORTE ESTADÍSTICO FORESTAL 2014–2015. Costa Rica.
SIREFOR, 2021. Estadísticas de aprovechamiento forestal en Costa Rica [WWW Document]. URL https://www.sirefor.go.cr/Sirefor/estadisticas (accessed 8 Apr. 21)
Stephens SS, Wagner MR (2007) Forest plantations and biodiversity: A fresh perspective. J for 105:307–313. https://doi.org/10.1093/jof/105.6.307
Streed E, Nichols JD, Gallatin K (2006) A Financial analysis of small-scale tropical reforestation with native species in Costa Rica. J for 104:276–282
Tenorio C, Moya R, Salas C, Berrocal A (2016) Evaluation of wood properties from six native species of forest plantations in Costa Rica. Bosque 37:71–84. https://doi.org/10.4067/S0717-92002016000100008
Trujillo-Miranda AL, Toledo-Aceves T, López-Barrera F, Günter S (2021) Tree diversity and timber productivity in planted forests: Pinus patula versus mixed cloud forest species. New for 52:177–195. https://doi.org/10.1007/s11056-020-09787-1
Vallejos J (2019) El cultivo del pilón (Hyeronima alchorneoides Allemão) en Costa Rica. Instituto Tecnológico de Costa Rica
Varma J (2021) Package ‘jrvFinance.’
Vega RP, Chaves GS (2011) Precios de la Madera en Costa Rica 2011 y tendencias de las principales especies comercializadas.
Víquez E, Pérez D (2005) Effect of Pruning on Tree Growth. Yield, and Wood Properties of Tectona Grandis Plantations in Costa Rica 39:381–390
Waring RH, Landsberg JJ (1997) A generalised model of forest productivity using simplified concepts of radiation-use efficiency, carbon balance and partitioning. For Ecol Manage 95:209–228
World Bank, 2021a. Consumer Price Index [WWW Document]. URL https://data.worldbank.org/indicator/FP.CPI.TOTL (accessed 9.9.21).
World Bank, 2021b. GDP per capita, PPP (current international dollars) [WWW Document]. URL https://data.worldbank.org/indicator/NY.GDP.PCAP.PP.CD (accessed 9.27.21).
Funding
Open Access funding enabled and organized by Projekt DEAL. The research was funded by Deutsche Bundesstiftung Umwelt with grant number AZ 670 34502/01–33/0.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors have no competing interests to declare that are relevant to the content of this article.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pinnschmidt, A., Yousefpour, R., Nölte, A. et al. Economic potential and management of tropical mixed-species plantations in Central America. New Forests 54, 565–586 (2023). https://doi.org/10.1007/s11056-022-09937-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11056-022-09937-7