Skip to main content
  • Research article
  • Open access
  • Published:

Adaptive economic and ecological forest management under risk

Abstract

Background

Forest managers must deal with inherently stochastic ecological and economic processes. The future growth of trees is uncertain, and so is their value. The randomness of low-impact, high frequency or rare catastrophic shocks in forest growth has significant implications in shaping the mix of tree species and the forest landscape. In addition, the fluctuations of wood prices influence greatly forest revenues.

Methods

Markov decision process models (MDPs) offer a rigorous and practical way of developing optimum management strategies, given these multiple sources of risk.

Results

Examples illustrate how such management guidelines are obtained with MDPs for combined ecological and economic objectives, including diversity of tree species and size, landscape diversity, old growth preservation, and carbon sequestration.

Conclusions

The findings illustrate the power of the MDP approach to deal with risk in forest resource management. They recognize that the future is best viewed in terms of probabilities. Given these probabilities, MDPs tie optimum adaptive actions strictly to the state of the forest and timber prices at decision time. The methods are theoretically rigorous, numerically efficient, and practical for field implementation.

Background

Forest planning and decision making models are often deterministic. They assume that the future is known, or they reduce uncertain variables to their expected values, with the hope that the results will work, at least “on average”. However, the future is highly uncertaina, and it can only be described in terms of probabilities, be they “objective probabilities” derived from factual data, or “subjective probabilities” reflecting a personal belief (de Finetti 1937).

There is a large literature concerning the role of risk in forest ecosystems (see e.g. Kant and Alavalapati 2014, p. 307–369). But there are few applications of Markov decision process models (MDPs) in forestry, a paucity also observed in the general operations research literature (White 1993). Nevertheless, MDPs are powerful and general methods that deserve more attention in forest management. This was recognized early on in Lembersky and Johnson (1975), and Lembersky (1976) who dealt with applications to the management of forest plantations. Their pioneering work was followed by several applications to the uneven-aged management of various forest types for financial objectives (Kaya and Buongiorno 1987), ecological objectives (Zhou and Buongiorno 2006), and combined financial and ecological objectives (Lin and Buongiorno 1998; Rollin et al. 2005; Zhou et al. 2008a, b).

The purpose of this paper is to describe how stochastic simulation and MDPs can be used to deal with decision making under risk in forestry. The data and the specific applications described below deal with mixed loblolly pine (Pinus taeda L.)-hardwood stands in the southern United States. However, the methods are general and can be applied to a variety of forest ecosystems, ranging from even-aged monospecific plantations to mixed-species uneven-aged forests, and from single stands to large forest areas. The objective throughout is to develop simple adaptive management guidelines that can be readily applied in the field.

The paper is organized as follows. The next section presents the methods: stochastic forest growth models and their reduction to Markov chains which are then used for prediction and optimization. This is followed by the results which show first some of the consequences of (wrongly) predicting the evolution of forest ecosystems without taking into account the effects of random ecological, climatic, or economic shocks. Next are shown predictions of the long-term evolution of forest ecosystems with and without management and the implications for economic and ecological criteria. And, MDP results optimizing decisions adaptively for discounted or undiscounted economic and ecological objectives, and combinations thereof. The discussion and conclusion deal with limitations of the methods and potential improvements.

Methods

The guiding principle of the methods is to represent the evolution of a forest ecosystem subject to various impact variables, such as prices and catastrophic disturbances, with Markov chains: probability matrices describing the frequencies of transition of each variable between discrete states. After the transition probabilities between the system states have been determined, the Markov model is used to predict future system states subject to specific management policies. Furthermore, decision process models (MDPs), optimization techniques based on Markov chains, are used to determine the best management policy for a particular objective. A policy is a set of rules that prescribe a decision for each observed system state at decision time. Various adaptive strategies are obtained in this way depending on the objective function which may deal with discounted or undiscounted criteria and constraints that may limit the decision domain.

Markov forest model

The different possible states of a forest stand are described in a practical manner by the level of basal area in trees of different sizes and species. Basal area per hectare, the sum of the cross sectional area of the trees, typically measured at breast height, is one of the simplest measures of forest stand density and composition. To be practical, the number of tree size and species categories should be kept lowb. For example Table 1 shows the definition of forest stand states used here, with the basal area in three tree size classes (small, medium, large), and two species categories (softwoods, hardwoods). With two levels of basal area, high or low, indicating below average or above average basal area in each species and size category, there are 26 = 64 possible stand states. The example stand state #13 in Table 1, denoted (001, 100) has high basal area in large pine trees and small hardwood trees, and low basal area in the other categories (Zhou and Buongiorno 2006).

Table 1 Definition of forest stand states according to basal area level by species and tree size

In MDPs, the evolution of stands over time is described by T = [p(s’|s)], a matrix of yearly transition probabilities between stand states s and s’. These probabilities cannot be established directly by direct observation due to lack of sufficient data in each stand. Simulation methods are used instead. A typical stochastic simulation model of stand growth suitable for this purpose has the form:

$$ {\mathbf{y}}_{t+\mathbf{1}}={\mathbf{G}}_t{\mathbf{y}}_t+{\mathbf{i}}_t+{\boldsymbol{\upvarepsilon}}_t $$
(1)

where y t  = [y ijt ] is the vector of the number of trees per unit area of species i and size j at time t. G t is a matrix of parameters describing the probability that a tree of a particular species stays alive and grows into a higher size class from t to t + 1. The vector i t refers to recruitment, the number of trees of a particular species that enters the smallest diameter class from t to t + 1. Both G t and i t vary depending on the stand condition, y t , through the stand density that affects tree growth and recruitment. ε t is a vector of random disturbances.

In the following applications the parameters of model (1) are estimated from observations on permanent sample plots in the mixed loblolly pine (Pinus taeda L.)-hardwood forests of the Southern United States (Schulte et al. 1998)c. The differences between the deterministic predictions of model (1), and the observations on the plots give observations on the random shocks ε t due to ice storms, wind, insect outbreaks, abnormal weather, etc.… that have affected forest growth during the observation period (Zhou and Buongiorno 2004).

To get the transition probabilities between stand states p(s’|s) model (1) is used to predict the future stand state, s’, of a random initial stand in state s, by bootstrapping a random shock ε t from the set of the observed shocks. This is repeated a sufficient number of times to obtain stable estimates of the transition probabilities p(s’|s) (Zhou 2005, p. 51).

Due to the short time typically covered by plot observations, the transition probability matrix T = [p(s’|s)] obtained in this manner reflects mostly small scale, high frequency disturbances. In addition, catastrophic, low frequency events, such as hurricanes must be added where needed, based on weather and geographic data. Table 2 shows the non-zero elements of the 64 × 64 transition probability matrix computed in this way and used in the results of this paper. The entry in row and column 000,000 means that if a stand has now low basal area in all species and size categories, the probability is 0.797 that it will be in the same state next year. The entry in row 000,000 and column 000,001 means that the same stand has a 0.027 probability of growing into a stand with high basal area in the largest size class. Table 2 also indicates that regardless of its current state, there is a 0.025 probability that the stand will be next year in state 000,000 due to a catastrophic event, a hurricane in this example (Zhou and Buongiorno 2006).

Table 2 Transition probabilities between stand states

The resulting T matrix is then used to predict future states:

$$ {\boldsymbol{\Pi}}_{t+1}={\boldsymbol{\Pi}}_t\mathbf{T}\kern0.75em t = 0,\ 1,\ 2\dots $$
(2)

where Π t  = [π st ] is the vector of probabilities of stand state s at time t, and the steady-state probabilities, Π = Π t  = Π t + 1 are obtained by solving the system of equations:

$$ \boldsymbol{\Pi} \left(\mathbf{I}-\mathbf{T}\right)=\mathbf{0} $$
(3)

where Π = [π s ] is the vector of long-term steady state probabilities for each stand state s, and I is the identity matrix. As equation (3) indicates, the long-term probabilities are independent of the initial state. Each probability is also interpreted as the fraction of a large forest landscape that is in each state s in the long run when only natural disturbances occur. As indicated below, an equation equivalent to (3) can also be used to predict the steady state probabilities when a particular management regime is followed indefinitely.

Markov price model

The financial returns from forest management depend critically on the price of wood. Perfect prediction of future prices is not possible, but probable price levels and changes can be drawn from the behavior of past prices. In parallel with the procedure for forest stands, the method describes market fluctuations by transition probabilities between price levels. For example, based on past data prices may be classified as low, medium, or high. And the transition probabilities between price levels are obtained from a price model such asd:

$$ {Q}_t=\alpha +\beta {Q}_{t-1}+{\nu}_t $$
(4)

where Q t is the price in year t and α, β are parameters estimated by regression analysis of past data. The matrix of transition probabilities between price levels is M = [p(m’|m)] where p(m’|m) is the probability that next year’s price level is m’ given the current price level m. These probabilities are obtained by simulations and repeated sampling from the distribution of the residual v t . An example of price transition probabilities used later in this paper is in Table 3 (Zhou 2005). Given a currently low price level (less than $84 · m−3), the probability that it will be low following year is 0.82. The probability that it will be medium (between $84 · m−3 and $94 · m−3) is 0.10, and the probability that it will be high (above $94 · m−3) is 0.08. Table 3 illustrates the high price autocorrelation (price stickiness): there is a high probability that the price stays at the same level from year to year.

Table 3 Annual transition probability matrix between price levels

System states and transition probabilities

Given transition probability matrices between stand states T = [p(s’|s)] and price levels M = [p(m’|m)], the transition probabilities of the forest-price system are:

$$ \mathbf{S}=\left[\mathrm{p}\left(\mathrm{j}\Big|\mathrm{i}\right)\right]=\left[\mathrm{p}\left(\mathrm{s}^{\prime}\Big|\mathrm{s}\right)\mathrm{p}\left(\mathrm{m}^{\prime}\Big|\mathrm{m}\right)\right] $$
(5)

where i is the current system state (stand state s and price level m), and j is next year’s system state (stand state s’ and price level m’). For example, with the 64 stand states in Table 2 and the three price levels in Table 3, there are 192 system states and 192 × 192 transition probabilities.

Decisions and immediate rewards

Decisions consist in moving instantly from one stand state to another by harvesting some of the trees. It is assumed that market prices are exogenous, i.e. the decisions have no effect on prices. A decision produces a stand state with less basal area, less carbon sequestered in the stand, and a different composition of tree species and size. It also generates revenues that depend on the volume of timber harvested and on the price at decision timee. A policy consists of a set of decisions applied systematically to each stand-price state. If the policy is applied to an entire forest, decisions alter the distribution of stand states and thus the landscape diversity, and possibly the part of the forest that stays in an old-growth (late-seral) state.

Ecological criteria include a variety of indices concerning the state of the stand and of the landscape induced by a policy. Although several diversity indices are possible (Magurran 1988), Shannon’s index (Shannon 1948) was used here throughout. For example, applied to tree species diversity Shannon’s index is:

$$ H=-{\displaystyle \sum_{a=1}^n{f}_a} \ln \left({f}_a\right) $$
(6)

where H is the tree-species diversity of a stand state, f a is the fraction of trees of species a, and n is the number of species. The fraction f a is preferably based on the basal area of trees rather than their number, to give more importance to the large trees. With this index, species diversity is between 0 and ln(n), and it is at the maximum when f a is the same for all species. The tree size diversity of a stand state is calculated in a similar way.

The forest landscape diversity is also expressed similarly with f s , the fraction of the forest landscape in state s. Another useful criterion is the fraction of the forested landscape in “old growth” state. For some forest types, old growth or late-seral stands have been defined in previous studies (Hummel and Calkin 2005). Here, old-growth is defined as the most frequent states that develop in the steady state predicted with equation (3) in the absence of harvest (Zhou and Buongiorno 2006). With the data in Table 2, the five old-growth states are 001,011; 001,111; 011,111; 101,111; and 111,111, which all have high basal area in the largest softwoods and hardwood trees.

The data on the expected timber volume and ecological criteria by stand state are obtained during the simulations with model (1) that give the transition probabilities. The immediate ecological reward of a decision is the characteristic of the state induced by the decision, for example, its tree diversity. The immediate financial reward is equal to the change in volume obtained by harvesting the stand from the current state to another, multiplied by the price level at decision time. The carbon sequestered was estimated from the amount of growing stock left after harvest, assuming 1.24 t · m−3 of CO2e for pulpwood and 1.57 t · m−3 for sawtimber (AFC 2014).

In sum, each immediate reward, monetary or ecological, is summarized by a vector V = [v ik ] where v ik is the immediate reward when the stand-price system is in state i and the decision is k (which changes the stand state or leaves it intact).

Consequences of management policies

A policy, c, is described by a particular decision matrix:

$$ {\mathbf{D}}_{\mathrm{c}}=\left[{\mathrm{p}}_{\mathrm{c}}\left(\mathrm{k}\Big|\mathrm{i}\right)\right] $$
(7)

where p c (k|i) is the probability of harvesting a stand to state k, given the stand-price state i under policy c. If the harvest policy is deterministic (the same action is taken all the time in each state) the probability of a particular decision is 0 or 1. To investigate the consequences of continuing the current management, independently of the price level, D c is based on the frequency at which harvests have changed the state of the sample plots during the observation period.

The transition probabilities between stand-price states given a particular policy are then obtained from:

$$ {\mathbf{E}}_{\mathrm{c}} = \left[{\mathrm{p}}_{\mathrm{c}}\left(\mathrm{j}\Big|\mathrm{i}\right)\right]=\left[{\mathrm{p}}_{\mathrm{c}}\left(\mathrm{k}\Big|\mathrm{i}\right)\mathrm{p}\left(\mathrm{j}\Big|\mathrm{k}\right)\right] $$
(8)

where p c (j|i) is the probability of ending in stand-price state j after one period from stand-price state i under policy c, and p(j|k) is the probability of ending in stand-price state j after one period from stand state k. Substituting E c for T in equation (3) gives the probability vector \( {\boldsymbol{\Pi}}_c=\left[{\pi}_i^c\right] \) where \( {\pi}_{\mathrm{i}}^{\mathrm{c}} \) is the steady state probability of stand-price state i under policy c.

Predicting expected undiscounted rewards

The expected undiscounted reward, R, such as the expected tree species diversity, or the expected annual timber revenue, that results in the long run by following a specific policy is:

$$ \mathbf{R}={\mathbf{V}}^{\hbox{'}}{\boldsymbol{\Pi}}_c $$
(9)

where \( \mathbf{V}\hbox{'}=\left[{r}_i^c\right] \) is a row vector of the immediate reward in stand-price state i under policy c. For financial returns \( {\mathrm{r}}_{\mathrm{i}}^{\mathrm{c}} \) is the timber income when the stand-price state is i and the policy is c. The expected long-run diversity of tree species and size, the expected basal area and carbon sequestered, and the expected annual production under a particular policy are obtained in similar fashion. The expected long-run forest landscape diversity is based on equation (6) where p a is replaced by the steady state probability of each stand state, \( {\pi}_{\mathrm{s}}^{\mathrm{c}} \), under the policy c.

The expected interval between harvests (or cutting cycle), in years, is the inverse of the yearly probability of a harvest. Let h i c = 1 if the decision called for by policy c is a harvest when the stand-price state is i. Then, the annual probability of a harvest is:

$$ {\pi}_h={\displaystyle \sum_i{\pi}_i^c}{h}_i^c $$
(10)

And the expected cutting cycle is:

$$ \mathrm{C}=\frac{1}{\pi_{\mathrm{h}}} $$
(11)

Predicting expected discounted rewards

Some criteria are typically discounted, future outcomes receiving less weight than current ones. This is especially the case for timber revenues, but ecological criteria may also be discounted (Boscolo et al. 1997; Howarth 2009). For a particular management policy, c, the present value of a criterion is obtained by solving the following system of equations (Hillier and Lieberman 2005, p. 785):

$$ {V}_i^c={r}_i^c+\frac{1}{1+\rho }{\displaystyle \sum_j{p}_c\left(j\Big|i\right){V}_j^c}\kern1.44em \forall i $$
(12)

where \( {\mathrm{V}}_{\mathrm{i}}^{\mathrm{c}} \) is the present value of the criterion of interest over an infinite horizon for an initial stand-price state i, when the harvest policy c is followed, with an immediate reward \( {\mathrm{r}}_{\mathrm{i}}^{\mathrm{c}} \) and a yearly interest rate ρ. In the results presented below, the interest rate was set at 3.08% per year, the real yield on AAA corporate bonds from 2001 to 2013, using the consumer price index as deflator (U.S. Government 2014). This choice implies a long-term, low-risk investment by a conservative investor.

Then, the average expected present value per ha over an entire forest with a particular initial distribution of stands and probability of price level is:

$$ PV={\displaystyle \sum_{i=1}^{192}{V}_i^c{\pi}_i^0} $$
(13)

where \( {\uppi}_{\mathrm{i}}^0 \) is the initial probability of stand-price state i.

Optimizing management policies

In parallel with the predictions of the effects of specific management policies, the optimization of policies differs depending on whether the rewards are discounted or undiscounted.

Optimizing expected discounted rewards

The policy that maximizes an expected discounted reward over an infinite horizon, such as the net present value of the timber income, without any other constraint, is found by solving the following linear programming program, originally due to d’Epenoux (1963).

$$ \begin{array}{l}\underset{y_{ij}}{ \max }{\displaystyle \sum_i{\displaystyle \sum_k{r}_{ik}{y}_{ik}}}\\ {}\mathrm{subject}\ \mathrm{t}\mathrm{o}:\\ {}{\displaystyle \sum_k{y}_{jk}-\frac{1}{1+\rho }}{\displaystyle \sum_i{\displaystyle \sum_kp\left(j\Big|i,k\right){y}_{ik}={\beta}_j\kern1.75em \forall j}}\\ {}{y}_{ik}\ge 0\forall i,k\end{array} $$
(14)

where y ik is the total discounted probability of being in state i and making decision k over an infinite horizonf. β j is the probability that the system starts in stand-price state j. This probability must be strictly positive, i.e. β j  > 0 and \( {\displaystyle \sum_{\mathrm{j}}{\beta}_{\mathrm{j}}}=1 \). After finding the best solution, \( {y}_{\mathrm{ik}}^{*} \), the best decisions are obtained from:

$$ {D}_{ik}=\frac{y_{ik}^{*}}{{\displaystyle \sum_k{y}_{ik}^{*}}}\kern0.75em \forall i,k $$
(15)

where D ik is the probability of making decision k when the stand-price state is i. The optimal policy is deterministic (D ik  = 0 or 1), and it is independent of the initial condition, i.e. of the β j ’s (Hillier and Lieberman 2005, p. 921), although the value of the objective function, the maximum discounted value of the reward, depends very much on the initial conditiong.

Optimizing expected undiscounted rewards

Another linear program gives the policy that maximizes an unconstrained undiscounted reward (Manne 1960). Examples of undiscounted criteria are the expected species diversity observed yearly in steady state, or the expected annual wood harvest.

$$ \begin{array}{l}\underset{y_{ij}}{ \max }{\displaystyle \sum_i{\displaystyle \sum_k{r}_{ik}{z}_{ik}}}\\ {}\mathrm{subject}\ \mathrm{t}\mathrm{o}:\\ {}{\displaystyle \sum_k{z}_{jk}-}{\displaystyle \sum_i{\displaystyle \sum_kp\left(j\Big|i,k\right){z}_{ik}=0\kern2.5em \forall j}}\end{array} $$
(16)

where z ik is the steady-state probability of stand-price state i and decision k.

Given the best solution, \( {z}_{\mathrm{ik}}^{*} \), the corresponding best decisions are:

$$ {D}_{{}_{ik}}^{\hbox{'}}=\frac{z_{ik}^{*}}{{\displaystyle \sum_k{z}_{ik}^{*}}}\kern0.75em \forall i,k $$
(17)

The best decisions obtained with equation (17) are still deterministic (D’ ik  = 0 or 1) and independent of the initial stand-price state. And, the maximum value of the objective function, such as the maximum expected undiscounted species diversity over an infinite time horizon, is also independent of the initial condition.

Multiple objectives

Policies that best meet multiple objectives simultaneously are obtained by modifying the objective function of models (14) and (16), or/and by adding constraints. In the results shown below the models were kept linear by expressing both the objective function and the constraints in undiscounted or discounted terms.

For example, in a discounted framework, a policy that maximizes the net present value of timber income, while keeping a specific fraction of the forest area in old growth state, is obtained with a model consisting of equations (14)–(15) and the additional constraint:

$$ {\displaystyle \sum_i{\displaystyle \sum_k{o}_{ik}{y}_{ik}}}\ge \alpha {O}^{*} $$
(18)

where o ik = 1 if state k is an old-growth state, o ik = 0 otherwise. O* is the maximum, unconstrained, discounted value of the fraction of old growth, obtained with model (16) to (17), and α is the desired fraction of this maximum that must be maintained by the policy.

As an example of an undiscounted framework, the maximum expected annual harvest that keeps the sequestered carbon dioxide equivalent (CO2e) at a specific level is obtained by solving model (16)–(17) with r ik equal to the amount of timber harvest resulting from decision k in state i, and the additional constraint:

$$ {\displaystyle \sum_i{\displaystyle \sum_k{c}_{ik}{y}_{ik}}}\ge {C}^{*} $$
(19)

where c ik is the amount of CO2e stored in the stand state that results from decision k in state i, and C* is the desired expected value of the stored CO2e.

Results

Effects of disturbances on predicted stand growth

Figure 1 illustrates the importance of recognizing random shocks in predicting the evolution of forest ecosystems. The data result from long-term simulations of a forest stand of mixed hardwoods-loblolly pines in the south of the United States, growing without human intervention (Zhou and Buongiorno 2004). The figure shows long-term predictions of total stand basal area, and basal area of pine trees. The continuous line refers to predictions obtained with the deterministic part only of model (1), ignoring the random shocks, ε t . The jagged shows instead the predictions obtained with the random shocks reflecting the high-frequency, low impact disturbances derived from the plot data, without catastrophic events.

Figure 1
figure 1

Simulated total basal area (a) and basal area of pine trees (b) in an uneven-aged loblolly pine stand with a deterministic model (smooth line) or a stochastic model with low-impact high-frequency shocks but no catastrophic event (jagged line). Source: Zhou and Buongiorno (2004).

While the total basal area is qualitatively similar for both models, the predictions of the basal area of pine trees are very different. Ignoring the random shocks, pines totally disappear in about 300 years. Instead, when the random disturbances are taken into account the basal area of the pine trees is never less than 20 m2 · ha−1, about a third of the total basal area. Forest scientists agree that stochastic disturbances “can be a major determinant of forest structure” (Oliver and Larson 1990). In the case of loblolly pine, recurrent fires are typical in the south of the United States. Lobllolly pines are more fire resistant than hardwoods. Fires also cause openings that favor pine regeneration. The simulations with stochastic shocks are thus more likely to predict correctly the long-term evolution of loblolly pine-hardwood stands than the deterministic version.

Effects of catastrophes on forest landscape

In addition to the high-frequency, low-impact disturbances, the forests are also subject to low-frequency catastrophes, such as hurricanes in the region investigated here. The long-term effects of catastrophic events on the forest landscape were predicted with the Markov model in Table 2, with and without the 0.025 probability of a catastrophe appearing in the first column. The fraction of the landscape that would be present in the long run (steady state) in each of the 64 possible states was predicted with equation (3).

Figure 2a shows the prediction without catastrophic events, while Figure 2b shows the prediction with catastrophes. Both figures show the fraction of the landscape area occupied in the steady state by each of the five old-growth states, and the aggregate fraction occupied by the 59 other states. Without catastrophic disturbances, 81% of the landscape was in an old growth state, and 38% consisted of stands in state #16 (001,111) with high basal area in hardwoods in all three size categories, and high basal area in large pine trees. However, when catastrophic events were taken into account, only 29% of the landscape was in an old-growth state and the dominant old-growth state #16 occupied only 12% of the landscape area. In sum, the long-run expected diversity of the forest landscape measured with Shannon’s index (6) was 80% higher due to catastrophic disturbances than it would have been without them.

Figure 2
figure 2

Expected distribution of forest area in steady state, without catastrophic disturbances (a) and with catastrophic disturbances (b).

Economic and ecological consequences of current management

The economic consequences of the current management, were expressed by the net present value (NPV) of timber income over an infinite horizon. The NPV was derived with equation (12) with the initial distribution of stands reflected in the plot data and a low initial price levelh. Other consequences, such as the annual wood production and the diversity indices were undiscounted expected values in steady state, obtained with equation (9), while equation (11) gave the expected length of the cutting cycle. The results in Table 4 compare the consequences of the current management with those arising from forest growth alone without management.

Table 4 Long-term effects of the current management on economic return and ecological criteria, compared with the outcome without management

Continuing the current management indefinitely generated an NPV of $1,339 · ha−1 in timber revenues and an average expected production of 1.9 m3 · ha−1 · y−1 over the entire forest. Due to this harvest the expected average basal area was 2 m2 · ha−1 lower with the current management than it would be without human intervention.

Using the diversity indices without management as a basis of comparison, the current management increased tree species diversity slightly (1%). Like catastrophic events, management also increased landscape diversity, though to a lesser extent (3%). However, by harvesting large trees it decreased the expected tree size diversity by 2%.

The fraction of old-growth in the forest maintained by the current management was 11% lower than it would be in a natural forest, and the value of the average carbon sequestration was 42 t · ha−1 lower. All these outcomes stemmed from harvests in individual stands that occurred at average intervals of 7 years under the current management.

Optimizing financial objectives

The policy that maximized the expected NPV was obtained by solving model (14)–(15), with a uniform distribution of initial stand-price states (i.e. β j = 1/192  j). This best policy was then used to obtain with equations (12)–(13) the expected NPV given the stand states distribution in the region of interest, as reflected by the plot data, and the current price level (“low” in 2014 according to the definitions in Table 3).

The results in Table 5 show how the best decisions adapted to the stand state and price level observed at decision time. For example, for a stand in state #64 (111,111), at low price, the best decision was to cut the stand to state #52 (110,011) by reducing the basal area of the large pine trees and of the small hardwood trees. If instead the current price was average, the best decision was to cut the stand to state #3 (000,010) which meant reducing further the basal area of the small and medium pine trees. Last, if the price was high, the best decision was to reach stand state #1 (000,000) by reducing the basal area to a low level in all species and size categories. In contrast, under current management, the plot data revealed that a stand in this same state #64 (111,111) was left intact 85% of the time and cut to state #44 (101,011) only 15% of the time.

Table 5 Management policy that maximizes the net present value of harvested wood, depending on stand state and price level

The expected consequence of this best financial policy for timber production and the various ecological indicators is in Table 6. The expected NPV and annual timber production were more than four times higher than with the current policyi. However, the expected basal area, and the diversity of tree species and size were all lower. The landscape diversity was less than a third of its level with the current policy, and no old growth was left in the entire landscape. The amount of CO2e sequestered was less than half, while the frequency of harvests was about double that with the current policy.

Table 6 Long-term effects of maximizing economic returns (NPV) from harvested wood on NPV and ecological criteria, compared with the outcomes due to the current management

Optimizing ecological objectives

Two criteria were used as applications of policies that best meet ecological objectives: maximizing tree species diversity, or carbon sequestration. The results in Table 7 show that at most, tree species diversity could be 7% higher than with the current management. Maximizing tree species diversity would not hurt financial returns, as in fact the expected NPV was $713 · ha−1 higher than with the current management as it led to substantially higher harvest of sawtimber and pulpwood. This more intensive harvest, however, reduced the fraction of the forest in old-growth state by 5 percentage points. Accordingly, the landscape diversity was reduced by approximately 30%, and the amount of CO2e sequestered was 52 t · ha−1 less. The tree size diversity index was also reduced slightly.

Table 7 Long-term effects of maximizing expected tree species diversity or CO 2 e sequestered, compared with the outcome due to the current management

Figure 3a shows the expected fraction of the total forest area in four stand states, starting in the current condition and evolving over a century with the policy that maximized the expected tree species diversity. The states in Figure 3a dominated the steady-state landscape when state #1 (000,000), with low basal area in all species and size categories occupied 20% of the forest landscape, state #7 (000,110) with high basal area in small and medium size hardwoods occupied 14.5%, state #8 (000,111) with in addition high basal area in large hardwoods occupied 14%, and state #16 (001,111) with in addition high basal area in large pine trees occupied 13%. With the tree species maximization policy, state #1 quickly over reached its steady-state level, while the other states approached it more slowly.

Figure 3
figure 3

Predicted changes of the expected fraction of forest area in selected stand states (a), and of the expected tree species diversity index (b), given the initial distribution of stand states revealed by the forest inventory and the policy that maximized tree species diversity. State 1 = (000,000), state 7 = (000,110), state 8 = (000,111), state 16 = (001,111).

Figure 3b shows how the tree species diversity index of the entire forest evolved from its initial level to the steady state with the policy that maximized expected tree species diversity. The diversity index is relative to the steady-state diversity of the unmanaged forest. With this policy it took approximately 30 years for the species diversity index to reach 90% of its maximum steady-state value of 108%.

Maximizing carbon sequestration

The amount of CO2e sequestered in the trees living biomass in steady state could be increased by 42 t · ha−1 compared to the current management (last column of Table 7), but in doing so, wood production would be very low (0.2 m3 · ha−1 · y−1). In fact, the level of CO2e sequestered in this way would be almost the same as that obtained by letting the forest grow naturally without any harvest (Table 4, last column). Indeed, the hands-off option maybe superior as it kept 29% of the landscape in old-growth state against only 15% when maximizing CO2e sequestration, at the small opportunity cost of $119 · ha−1 in foregone NPV from timber revenuesj. The other criteria were similar to those obtained by natural forest growth.

Constrained financial and ecological objectives

Figure 4 shows the results of a series of optimizations in which the objective was to maximize the NPV of timber revenues given its initial stands distribution and current low price, subject to increasing constraints on the fraction of old growth maintained on the forestk. Figure 4 documents the opportunity cost, in terms of the NPV of foregone revenues, of higher levels of old growth. While approximately $5,300 · ha−1 of NPV could be obtained without any old growth, this was reduced to about $1500 · ha−1 when 22% of the landscape was kept in old growth. The marginal cost of old growth increased at increasing levels of old growth, from a low of $368 · ha−1 from 0 to 5% of old growth to a high of $1145 · ha−1 for 21% to 22% of old growth.

Figure 4
figure 4

Maximum net present value (NPV) of harvested wood at increasing levels of the fraction of the forest landscape kept in an old-growth state. The old growth fraction is the constant annual equivalent of the present value of the old growth fraction over an infinite horizon.

As an example of application with multiple ecological objectives, Figure 5 shows the results of maximizing annual wood production (sawtimber and pulpwood), subject to constraints on the amount of CO2e sequestered in the growing stock. Both objectives were undiscounted and the constraint was as in equation (19). As indicated by Figure 5, there was no conflict between wood production and carbon sequestration up to approximately 100 t · ha−1 of CO2e. Beyond that point, wood harvest decreased rapidly at an increasing rate. Increasing carbon storage from 100 to 150 t · ha−1 decreased wood production by 1.2 m3 · ha−1 · y−1, and increasing it from 200 to 250 t · ha−1 decreased production by 3.3 m3 · ha−1 · y−1. As shown above, the maximum carbon storage achievable was about 245 t · ha−1, and at that point there was hardly any wood harvest.

Figure 5
figure 5

Maximum expected annual wood harvest at increasing levels of CO 2 e stored in growing stock. Both the annual harvest and the stored CO2e are undiscounted.

Discussion and conclusions

In a 1999 interview, G.B. Dantzig, the inventor of linear programming, remarked that “all problems that are solved under deterministic means have that fundamental weakness—they don’t properly take uncertainty into account” (Dantzig 1999). Forestry planning is no different. Risk, uncertainty, and stochastic behavior are critical parts of how forest ecosystems work. As illustrated above, simulations of mixed loblolly pine-hardwood stands reveal how ignoring random disturbances due to biological or catastrophic shocks in forest growth models can lead to wrong ecological predictions, such as the disappearance of pines in this context. To these natural sources of risk must be added the high risk that stems from price fluctuations which complicates decision making with financial objectives.

The methods suggested in this paper to handle risk in forest decision making follow the modeling approach outlined by Holling et al. (1986). The first step is to “bring the world to the laboratory” with possibly complex and non-linear stochastic models, and then simplify the models to allow for efficient optimization with Markov decision process models. MDPs are attractive for their simplicity while keeping the essence of planning problems under riskl. They recognize that the future is unknown and that it may be described only in terms of probabilities. And, they lead to best adaptive management policies whereby decisions depend entirely on the systems state, i.e. the state of knowledge of the decision maker at decision time.

The fact that with MDP models the future state of a forest ecosystem depends only on its current state should not be viewed as a shortcoming, but recognition of fact. Predictions can rely only on current knowledge. This current knowledge may include past behavior of the system (such as past tree growth) if necessary. Thus, Markov models are not necessarily “memory less” and apparent failure of Markov models to correctly represent forest growth (Binkley 1980; Roberts and Hruska 1986; Johnson et al. 1991) may be due to incomplete description of the current ecosystem state rather than to a shortcoming of the Markov model. Similarly, Markov models of price changes are very general, embracing random walk, rational expectations, autoregressive, and “any stochastic model in which the price is conditional on previous prices” (Taylor 1984). In sum, Markov chains allow the simplification of complex, multi-dimensional stochastic processes and make their optimization easier or at all possible (Holling et al. 1986; Insley and Rollins 2005).

There exists considerable potential and flexibility to enrich the models presented in this paper as they are relatively small compared with the size of the problems that can be solved with current linear programming software. Still, to be successful in forest ecosystem management the MDPs should adhere to the general principles of parsimony and simplicity. To this end, the number of system states must be kept as small as possible, with a few state variables (tree species and size categories, prices) and a few levels of these variables.

In the same spirit of simplicity, multiple objectives such as optimizing financial returns subject to ecological constraints, of vice-versa, have been treated here with linear models. This requires expressing both the objective function and the constraints in either discounted or undiscounted criteria. While discounting financial returns from forestry has long been a standard procedure (Faustmann 1849), it is less so for ecological criteria. Nevertheless, it is plausible to give more weight to the present ecological characteristics of a forest ecosystem than to their future values. Indeed, it seems preferable to maintain desirable current states (such as old-growth stands) rather than lose them, even if they could be restored later (albeit after a long delay). It has also been argued on theoretical grounds that “the future benefits provided by a generic public good--environmental quality--should be discounted at a rate that is close to the market rate of return for risk-free financial assets” (Howarth 2009), a principle that has been applied for example in discounting future carbon sequestered in forests (Boscolo et al. 1997)m.

The methods presented here can be expanded in many directions, in particular to deal with the new carbon markets and the highly stochastic prices for CO2e sequestered in forests. The treatment is parallel to the stochastic prices of timber. This would increase substantially the size of the problem, but still keep it well within the capabilities of current linear programming software.

Other issues may require deeper modifications, to deal with non-stationary processes. For example, “climate change can affect forests by altering the frequency, intensity, duration, and timing of fire, drought, introduced species, insect and pathogens outbreaks, hurricanes, windstorms, ice storms, or landslides” (Dale et al. 2001). This implies a change in the transition probabilities over time. Although non-stationary problems can in principle be converted to stationary ones by a reformulation (Bertsekas 1995, p. 167), efficient numerical methods are still elusive (Ghate and Smith 2013). Nevertheless, the consequences of changes in transition probabilities on the steady-state criteria can be readily explored with the methods described in this paper. In sum, given the conceptual generality and the well-developed theory of MDPs, coupled with the powerful solution techniques available, the MDP approach is well suited to deal with risk in forest ecosystem management, and to develop practical adaptive management policies with both economic and ecological objectives.

Endnotes

aIn this paper the terms “risk” or “uncertainty” are synonymous. Although Knight (1921) defines uncertainty as the absence of probabilities, it may be argued that probabilities, objective (data based) or subjective (opinion based) always exist to some degree.

bWhile modern linear programming software places almost no limit (computationally) on the number of states, it is best for all practical purposes (identifying states of nature, doing the computation, and laying out the recommendations) to use the minimum number of states needed for sufficient accuracy and realism of applications.

cFor other applications, any stochastic model of forest growth can be used to compute the transition probabilities used in the MDP approach.

dIf the price series is not stationary, the price data are first de-trended, and the trend is added to the discount rate in present value calculations.

eIn the particular application described here, prices vary by tree size and species, and all the prices change in parallel with the probabilistic changes of the price index described in Table 3.

fThat is: y ik  = z 0ik  + dz 1ik  + d 2 z 2ik  + d 3 z 3ik  + … where z tik is the probability of state i and decision k at time t and d = 1/(1 + ρ) is the discount factor (Hillier and Lieberman 2005, p. 921).

gAs applied to forestry, the MDP described by equations (10) and (11) is a generalization of (Faustmann 1849) formula recognizing that future stand states and prices are known only as probability distributions. The classical, deterministic Faustmann formula is a special case in which the transition probabilities are 0 or 1 (Buongiorno 2001).

hAccording to Timber Mart South (2014), the average price of softwood sawtimber stumpage in the South of the United States in the second quarter of 2014 was $25 per short ton or approximately $19 · m−3, placing it in the low range according to the definitions in Table 3.

iPrevious results also show that the optimum adaptive policy derived here is superior in terms of NPV to an optimum fixed policy that converts stands to a chosen state at fixed intervals (Zhou et al. 2008b).

jThe reduction in old-growth fraction from 29% to 15% was due to the best decision calling for a harvest when the stand was in the old-growth state #32 (011,111). Although the harvest occurred only every 86 years on average, this was equivalent to a low frequency natural catastrophe reducing substantially the fraction of the landscape in old-growth state.

kIn Figure 4, the old growth fraction is the right-hand side of equation (18) with O* replaced by its annual constant perpetual equivalent, ρO */(1 + ρ).

lSee Buongiorno and Gilless (2003, p. 337–371) for an introduction to Markov and MDP models in forestry.

mAlternatively, models (10)–(11) and (12–13) can be extended to deal with discounted objective functions and undiscounted constraints, or the reverse, by introducing non-linear constraints (e.g. Rollin et al. 2005; Zhou 2005), but at the cost of the attendant numerical difficulties.

References

  • AFC (Alabama Forestry Commission) (2014) How much carbon have your trees stored? Alabama Forestry Commission. http://www.forestry.state.al.us/HowMuchCarbonHaveYourTreesStored.aspx?bv=5&s=0. Accessed 26 Feb 2015

  • Bertsekas DP (1995) Dynamic programming and optimal control, vol 2. Athena Scientific, Belmont, p 292

    Google Scholar 

  • Binkley CS (1980) Is succession in hardwood forests a stationary Markov process? For Sci 26:566–570

    Google Scholar 

  • Boscolo M, Buongiorno J, Panayotou T (1997) Simulating options for carbon sequestration through improved management of a lowland tropical rainforest. Environ Dev Econ 2:241–263

    Article  Google Scholar 

  • Buongiorno J (2001) Generalization of Faustman’s formula for stochastic forest growth and prices with Markov decision process models. For Sci 47(4):466–474

    Google Scholar 

  • Buongiorno J, Gilless JK (2003) Decision methods for forest resource management. Academic Press/Elsevier, Amsterdam, p 439

    Google Scholar 

  • d’Epenoux F (1963) A probalistic production and inventory problem. Manage Sci 30(1):98–108

    Article  Google Scholar 

  • Dale VH, Joyce LA, McNulty S, Neilson RP, Ayres MP, Flannigan MD, Hanson PJ, Irland LC, Lugo A, Peterson CJ, Simberloff D, Swanson FJ, Stocks BJ, Wotton BM (2001) Climate change and forest disturbances. BioScience 51(9):723–734

    Article  Google Scholar 

  • Dantzig GB (1999) Interview of Georges B. Dantzig by P. Horner. First published in the October 1999 issue of OR/MS Today. http://www.orms-today.org/orms-6-05/dantzig.html. Accessed 26 Feb 2015

  • de Finetti B (1937) La prévision, ses lois logiques, ses sources subjectives. Annales de l’Institut Henri Poincaré 7(1):1–68

    Google Scholar 

  • Faustmann M (1849) Calculation of the value which forest land and immature stands possess for forestry. Allgemaine forst-und jagdzeitung 15:441–455

    Google Scholar 

  • Ghate A, Smith RL (2013) A linear programming approach to nonstationary infinite horizon Markov decision processes. Oper Res 61(2):413–425

    Article  Google Scholar 

  • U.S. Government (2014) Economic report of the President. U.S. Government Printing Office. http://www.gpo.gov/fdsys/browse/collection.action?collectionCode=ERP. Accessed 26 Feb 2015

  • Hillier FS, Lieberman GJ (2005) Introduction to operations research, 8th edn. Mc Graw Hill, Boston, p 1061

    Google Scholar 

  • Holling CS, Dantzig GB, Winkler C (1986) Determining optimal policies for ecosystems. In: Kallio M et al (eds) Systems analysis in forestry and forest industries. TIMS studies in management science 21, North Holland, Amsterdam, pp 463–473

    Google Scholar 

  • Howarth RB (2009) Discounting, uncertainty, and revealed time preferences. Land Econ 85(1):24–40

    Google Scholar 

  • Hummel S, Calkin DE (2005) Costs of landscape silviculture for fire and habitat management. For Ecol Manage 207(3):385–404

    Article  Google Scholar 

  • Insley M, Rollins K (2005) On solving the multirotational timber harvesting problem with stochastic prices: a linear complementary formulation. Amer J Agric Econ 87:735–755

    Article  Google Scholar 

  • Johnson SE, Ferguson IS, Li RW (1991) Evaluation of a stochastic diameter growth model for mountain ash. For Sci 37:1671–1681

    Google Scholar 

  • Kant S, Alavalapati RR (2014) Handbook of forest resource economics. Part 4: Economics of risk, uncertainty and natural disturbances. Routledge, New York, p 559

    Google Scholar 

  • Kaya I, Buongiorno J (1987) Economic harvesting of uneven-aged northern hardwood stands under risk: A markovian decision model. For Sci 33(4):889–907

    Google Scholar 

  • Knight FH (1921) Risk, uncertainty, and profit. Houghton Mifflin, Boston

    Google Scholar 

  • Lembersky MR (1976) Maximum average annual volume for managed stands. For Sci 22:69–81

    Google Scholar 

  • Lembersky MR, Johnson KN (1975) Optimal policies for managed stands: An infinite Markov decision process approach. For Sci 21:109–122

    Google Scholar 

  • Lin CR, Buongiorno J (1998) Tree diversity, landscape diversity, and economics of maple-birch forests: Implications of markovian models. Manage Sci 44(10):1351–1366

    Article  Google Scholar 

  • Magurran AE (1988) Ecological diversity and its measurement. Princeton University Press, Princeton, NJ, p 179

    Book  Google Scholar 

  • Manne AS (1960) Linear programming and sequential decisions. Manage Sci 6:259–267

    Article  Google Scholar 

  • Oliver CD, Larson BC (1990) Forest stand dynamics. McGraw-Hill, Chicago, p 467

    Google Scholar 

  • Roberts MR, Hruska AJ (1986) Predicting diameter distributions: a test of the stationary Markov model. Can J For Res 16:130–135

    Article  Google Scholar 

  • Rollin F, Buongiorno J, Zhou M, Peyron JL (2005) Management of mixed-species, uneven-aged forests in the French Jura: From stochastic growth and price models to decision tables. For Sci 51(1):64–75

    Google Scholar 

  • Schulte B, Buongiorno J, Lin CR, Skog K (1998) SouthPro: A computer program for uneven-aged management of loblolly pine stands. USDA Forest Service, Forest Products Laboratory, General technical report FPL-GRT-112, Madison, p 47

    Google Scholar 

  • Shannon CE (1948) A mathematical theory of communication. Bell Syst Tech J 27(379–423):623–656

    Article  Google Scholar 

  • Taylor CR (1984) Stochastic dynamic duality: Theory and empirical applicability. Amer J Agric Econ 66:351–357

    Article  Google Scholar 

  • Timber Mart-South (2014) South wide pine stumpage prices. http://www.timbermart-south.com/prices.html. Accessed 26 Feb 2015

  • White DJ (1993) A survey of applications of Markov decision processes. J Oper Res Soc 44(11):1073–1096

    Article  Google Scholar 

  • Zhou M (2005) Managing resources for multiple purposes under risk: Markov models of southern mixed loblolly pine-hardwood forests. University of Wisconsin-Madison, Dissertation, p 135

    Google Scholar 

  • Zhou M, Buongiorno J (2004) Nonlinearity and noise interaction in a model of forest growth. Ecol Model 180:291–304

    Article  Google Scholar 

  • Zhou M, Buongiorno J (2006) Forest landscape management in a stochastic environment, with an application to mixed loblolly pine-hardwood forests. For Ecol Manage 223:170–182

    Article  Google Scholar 

  • Zhou M, Buongiorno J, Liang J (2008a) Economic and ecological effects of diameter caps: A Markov decision model for Douglas-fir/Western hemlock forests. For Sci 54(4):397–407

    Google Scholar 

  • Zhou M, Liang J, Buongiorno J (2008b) Adaptive versus fixed policies for economic and ecological objectives in forest management. For Ecol Manage 254:178–187

    Article  Google Scholar 

Download references

Acknowledgements

The preparation of this paper was supported in part by the USDA Forest Service, Southern Research Station, through a cooperative research agreement with Joseph Buongiorno, directed by Jeff Prestemon. It was also supported in part by the USDA McIntire-Stennis fund WVA00105. We are also grateful to the organizers of the workshop on Risk and Uncertainty in Ecosystems Dynamics, Beijing Forestry University, October 13–17, 2014, for motivating the study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Joseph Buongiorno.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

This paper was based on research carried jointly by JB and MZ. JB wrote the first draft of the paper which was revised and approved by MZ. Both authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, 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 license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Buongiorno, J., Zhou, M. Adaptive economic and ecological forest management under risk. For. Ecosyst. 2, 4 (2015). https://doi.org/10.1186/s40663-015-0030-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40663-015-0030-y

Keywords