**Open Journal of Forestry**

Vol.05 No.06(2015), Article ID:58636,13
pages

10.4236/ojf.2015.56053

Optimal Stochastic Pine Stands Harvest Rotation Policies

Eduardo Navarrete

Department of System Engineering, Universidad de la Frontera, Temuco, Chile

Email: eduardo.navarrete@ufrontera.cl

Copyright © 2015 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 15 May 2015; accepted 1 August 2015; published 6 August 2015

ABSTRACT

A new Faustmann optimal rotation harvesting stands’ problem under Brown geometric price and Logistic and Gompertz wood stock, diffusions is presented. The optimal cut policies for the stochastic Faustmann model and the single harvest rotation or Vicksell model are evaluated in the case of a Chilean Radiata pine forest company. The company cut policy validates the Vicksell model, its optimal cut policies overestimate the company policy cut in 1.2%, in the Gompertz case, and underestimate it in 2.3%, in the Logistic case. The Faustmann optimal cut policies present a larger underestimation of the company cut policy in 10.1%, in the Gompertz case, and in 21.5%, in the Logistic case. The preference for shorter evaluation period that the company shows is due to the organizational risk that the forest economic sectors has in Chile.

**Keywords:**

Optimal Tree Cutting, Faustmann Stochastic Formula Component, Optimal Stopping Problem

1. Introduction

The optimal forest harvest models development is marked by the emergence of two controversial issues. The first issue was formulated in its early beginning by Faustmann in 1849 . He formulated a multiple rotation version considering the dynamic effect that the future renovations had in the rotation period of the cut and planting of the trees. The issue was not considered relevant by researchers and forester management, who preferred a simple rotation or Vicksell model. The Faustmann model was finally rediscovered, first by Gaffney in 1957 , and finally by Samuelson in 1976 , who validated Faustmann’s deterministic formula as the correct one, since it was the only one to consider the cost of land rent. An increasing number of researchers continued extending this model (Brazee, 2001; Chang, 2001; Amacher et al., 2011) . The second issue was strongly formulated by Samuelson (1976) , calling to replace the simple notion of stationary equilibrium by the notion of perpetual Brownian motion of wood price. Many researchers followed his recommendation and considered the single rotation problem as an American call option, taking into account wood price as a geometric Brown diffusion (see, Newman, 2002; Clarke & Reed, 1989 ; Thomson, 1992; Platinga et al., 1998). Other authors also considered the wood stock as a second geometric Brown diffusion (Morck & Schwartz, 1989; Insley, 2002) , or as a Logistic diffusion (Alvarez & Koskela, 2007; Navarrete, 2011) . Willassen (1998) introduced uncertainty to the Faustmann model, incorporating the forest growth as a stochastic Markov diffusion process and characterizing the properties of the optimal solutions. Sodal (2002) simplified Willassen’s approach in a closed-form rotation formula for the same state stochastic variable and Insley & Rollins (2005) proposed a numeric algorithm to solve its Hamilton-Jacobi- Bellman solution. Navarrete & Bustos (2013) extended the Faustmann model considering, price as a geometric Brown and wood stock as a Logistical or Gompertz, diffusions processes, transforming the model into an equivalent optimal stopping problem, in spite of its contribution. This paper requires some improvement which will be done in the present paper. Therefore this paper has two basic objectives. In section 2.3, the proof of the equivalence of the Faustmann rotation model under Brown Price and ITO wood stock independent diffusion processes and a one dimensional modified wood stock Optimal Stopping problem is formalized as a reformulation lemma. In section 4, a new algorithm is developed to obtain the Faustmann equivalent one dimensional optimal stopping solution, intersecting the optimal condition curve with the expected wood stock growth feasibility restriction. Finally new stands growth data for the middle age 11 and 12 years is incorporated in Section 3, in order to validate Vicksell and Faustmann stands cut models, under price and wood stock uncertainty.

2. Stochastic Harvest Models

2.1. Stochastic Harvest Stands Rotation Models

The basic model considers two independent stochastic processes: an ITO diffusion for the wood stock, Equation (1) and a geometric Brownian diffusion, Equation (2). The functional objective for the single rotation is given by Equation (3) and for the multiple rotations by Equation (4) see Johnson (2006) :

(1)

(2)

(3)

(4)

Model notation:

Deterministic state variable Economic Parameter

t = Wood stock age C = Stand regeneration cost

Stochastic sate variables P_{0} = Price at time 0

V_{t} = Wood stock c = C/P_{0}

P_{t} = Wood stumpage price at time t R, Q = Probabilistic metrics

Diffusion Parameters F, Z = Functional objective

W_{v} = Wiener wood stock i = Risky rate of return

μ(V) = Wood stock diffusion drift rate parameter r_{T} =, Capitalized rate of return

σ(V) = Wood stock volatility parameter Optimal Parameters

W_{p} = Wiener price T = Optimal cut time

α = Wood price diffusion drift rate = Expected optimal cut volume

β = Wood price volatility

2.2. Reformulation of the Simple Harvest Rotation Problem

The stochastic model (1) (2) and (3) is difficult to solve due to both diffusions. But it is equivalent to the following one dimensional diffusion Optimal Stopping problem, (5) and (6) see Navarrete (2011) .

(5)

(6)

The solution of this problem is given by the (HJB) Equation (7), see Navarrete (2011) .

(7)

Assuming the existence of a frontier V^{*} that divides the zone into a
continuation zone (no-cutting) and a stopping zone (immediate-cutting), the solution
to the equation HJB is given by (8) for the continuation zone (8) and stopping zone
(9). A solution of (8) is given by (10) (see Johnson,
2006 ).

(8)

(9)

(10)

where
(resp.,) is
strictly increasing (resp., decreasing), function since the payoff function are
bounded and small, V is positive and should remain bounded and positive a V→0, necessarily
then B→0. The solution also fulfils the smooth-pasting condition at the free boundary
point V^{*}, equations (11) and eliminating the constant A gives Equation
(12).

(11)

(12).

2.3. Reformulation of the Multiple Harvest Rotation Problems

The Faustmann optimal rotation model (1) (2) and (4) its equivalent to the following Optimal Stopping one dimensional diffusion problem (13) and (14), see Lemma 1.

Lemma 1. Reformulation Lemma:

The Faustmann optimal rotation problem is equivalent to the following Optimal Stopping one dimensional diffusion problem (13) and (14).

(13)

(14)

Proof:

Given the independence of both variables the expectation can be calculated as the product of two independent expectation “” that are independent of the deterministic parameters and and the constant C, so they can be taken outside its domain. Introducing the integration of the price Brown diffusion

with
and C = P_{0}c in the functional objective (4) results in the objective
(15) and reduce the Faustmann problem to the following one dimensional diffusion
problem (15) and (16).

(15)

(16)

Dividing the objective functional (15) by the constant P_{0} and applying
the Thijssen version of the Girsanov theorem to Equations (15) and (16), for the
Martingale M_{t} and the Radom-Nykodym derivative
(see Thijssen, 2010, Appendix A) permits to reformulate
the problem (13) (14) by its equivalent Q metric optimal stopping problem (16) and
(17).

The formulation of the HJB equation for this problem is given by Equation (17) for the capitalized risky rate of return see Navarrete & Bustos (2013) .

(17)

In this case the differential equation for the continuation region is given by the non homogenous differential Equation (18), and by Equation (19) for the stopping zone.

(18)

(19)

The solution to the ordinary differential equation, (18) under the initial condition
for a given risky rate of return r_{T} is given in Equation (20), with,
the positive increasing solution of the homogenous part and
the particular solution of Equation (18).

(20)

In this case the smooth pasting condition for each parameter r_{T} are given
by Equation (21):

(21)

So
must fulfill smooth-pasting condition (22) for each parameter r_{T}.

(22)

2.4. Wood Stock Sigmoid Diffusion Equations

The basic requirement of a pine stand growing diffusion is its sigmoid pattern (Garcia, 2005) . The logistic diffusion, Equation (23) is a special case of the sigmoid model given by and, where µ and γ are the drift and saturation parameters and σ is the volatility parameter.

(23)

The integration of the value of V is given by Equation (24) (Kloeden & Platen, 1992: p. 125) and its expected value is given by Equation (25)

(24)

(25)

Other important sigmoid diffusion is the Gompertz geometrical diffusion, which is given by Equation (26), which is integrated to the expression (27), and its expected values given by expression (28) (see Gutierrez, 2009).

(26)

(27)

(28)

The Geometric Brown Price diffusion is given by Equation (29), and integrates in Equation (30).

(29)

(30)

3. Experimental Data and Parameter Fitting

3.1. Logistic Diffusion Fitting

The experimental data belonged to 122 harvest pine stands of a Chilean Forest Company
in the Araucania Chile region, between years 1999 and 2005 (see Appendix A), and
came from different sample plots, with site indexes between 30 and 35 meters, representing
sites with high forest aptitude and a tree average initial volume of 32 m^{3}/ha
at the first 4 years after initial seed cultivation period. The additional 20 point
for years 11 and 12 were taken from Alvarez et al. (2012). The business harvest
cut data for a Logistic diffusion for the 95% confidence range is given by data
point plotted in Figure 1.

The logistic diffusion parameter cannot be adjusted by maximum verisimilitude, see
(Beskos et al., 2006) , so it was fitted using
a logistical nonlinear regression and a Monte Carlo/Bootstrap simulation sampling
method, implemented by Meyer et al. (Loglet Lab.1 software,
1999) . Choosing V_{0} = 1/2γ = half saturation volume, and T_{0}
= T_{m}, time to achieve that volume, Equation (25) is transformed into
the more conventional expression (31).

Figure 1. Wood volume per
hectare, (m^{3}/ha) versus years.

(31)

With; 1/g = saturation volume, m = growth rate parameter, V_{s} = Saturation
volume and t_{m} = time to achieve the midpoint of the saturation volume.
The standard deviation Sd(∞), at the saturation zone is constant, and given by Equation
(32).

(32)

Since the saturation volume V_{s} is also constant, sigma can easily be
estimated by Equation (33) and the summary of the parameter fitting is shown in
Table 1.

(33)

3.2. Gompertz Diffusion Fitting

The Gompertz model can be fitted by common statistic features, such as maximum verisimilitude (see Gutierrez et al., 2008 ), but the lack of equal time distribution of data make it difficult. So a quadratic fitting method was development using the SSPS software. By taking natural logarithm and arranging it we get Equation (34).

(34)

with (35)

Given a value for k, a quadratic fitting for
and
was done estimating the value of A, B and C until a common value for θ was obtained
from A and B, determining the estimation for θ, k and σ. The deterministic parameter
only requires a linear fitting with.
Both fittings were done for the initial value V_{0} = 32 (m^{3}/ha)
and the results are summarized in Table 2.

3.3. Wood Price Diffusion Fitting Equations and Regeneration Costs

The stumpage stand price Brownian diffusion parameters are estimated from saw logs and pulp logs exportation prices (see Appendix B). The summary of Brown diffusion parameters for the stumpage price and the regeneration cost of radiata Pine Stands are given in Table 3, see Navarrete & Bustos (2013) .

4. Stochastic Optimal Pine Harvestings Results

4.1. Logistic Wood Stock and Brown Stumpage Price models

The deterministic optimal solution given by, for the simple rotation model and by equation (36) for the multiple rotation model, see Navarrete (2013).

Table 1. Logistic fitting parameters.

Table 2. Gompertz diffusions parameters estimations.

Table 3. Radiata pine price diffusion parameters and capital and regeneration stands costs.

Source: Appendix B.

(36)

The solution of the Vicksell and Faustmann model under Brown price and Logistic
wood stock diffusion requires the solution of the differential Equations (8) and
(19) for
and.
The positive increasing function ψ(V) for the Vicksell model, is the solution of
the homogenous component (37) of the differential Equation (8) with r_{T}
= i in the Vicksell case.

The solution of Equation (37) is given by the Kummer’s confluent hyper geometric function, expression (38) with the positive root θ by Equation (39), ψ(V) also must fulfill the smooth pasting condition (12), which was programmed in Mapple 15, see Navarrete (2013).

(37)

(38)

(39)

The Faustmann stochastic optimum requires the solution of the differential Equation
(16), so the positive function ψ(V) is the solution of its homogenous Equation (37)
for the parameter r_{T} and also fulfill the smooth pasting condition, (22).
For a given value of the capitalized interest
Equation (37) can be solved using the same methodology of the Vicksell solution,
generating a family of optimal cuts for the different parameters r_{T}.
The optimum solution is obtained intersecting the optimal parametric solution for
different T with the corresponding expected value logistical modify diffusion, (14)
for,
see Figure 2. The summary of all optimal cuts results
for the aggregate 30/35 site index series of the multiple and simple rotation harvest
model under Brown price and Logistic wood stock diffusions is given in Table 4.

The simple optimal cut is a better explanation of the company cut policy, since it only underestimates it in a 2.5%. As expected the Faustmann optimal cut is lower and underestimate the company policy in 11.5%. The deterministic optimal cuts are even worse and underestimate the company policy in 23.4% in the simple model and in 25.6% in the multiple cases. Finally the Logistic wood stock underestimates the saturation Volume in 4.2%.

4.2. Gompertz Wood Stock and Brown Stumpage Price Case

The deterministic optimal solution, is given by for the simple rotation model and for Equation (41) for the multiple rotation case the optimal is obtained intercepting condition (40) with Equation (27). See Navarrete (2013).

(40)

Figure 2. Faustmann logistic optimal cut.

Table 4. Optimal harvest rotation optimal results, Brown price logistic wood stock diffusion.

The solution of the Vicksell and Faustmann model under Brown price and Gompertz
wood stock diffusion requires the solution of the differential Equations, (8) and
(16) for
and.
The stochastic increasing function ψ(V), for these cases is given by the solution
of the homogenous Equation (41) with r_{T} = i for the simple rotation case
and the capitalized rate of interest r_{T} for the Faustmann case.

(41)

Replacing the parameters given in Equations (42) in (41), gives the differential Equation (43) of the Exponential, Ornstein-Uhlembeck diffusion whose positive increasing solution ψ(V) is given by Equation (44) (see Johnson, 2005), with:, b = 0.5 and

(42)

(43)

(44)

The optimal solutions also requires that ψ(V) fulfills the smooth pasting condition (38), which is programmed in Mapple 15 generating one optimal cut solution for the Vicksell model and a family of optimal cut condition parametrized by T for the Faustmann case. In this late case the optimal solution is obtain intercepting the family with the expected volume of the Gompertz modify diffusion (20), see Figure 3.

The deterministic and stochastic optimum were programed in Maple 15, using in this case the KummerU function of the program, the results are summarized in Table 5. These results show that the simple stochastic rotation model under Brown Price and Gompertz wood stock diffusion is a better explanation of the company cut policy, since it only overestimates it in 1.2%. As expected, the Faustmann model has lower optimal cut and underestimate it in 9.2% and the deterministic optimal cut is even lower, and underestimate the company policy in 27.9% in the simple case and in 33.1% in the multiple cases. Finally, the Gompertz wood stock overestimated the Saturation volume in 8.6%.

5. Conclusions

1) The optimal cut company policy validates the use of the simple stochastic rotations model under Bown price and Logistic or Gompertz wood stock diffusion.

2) The discrepancy in the theoretical and practical cut policy can be explained by the preference that the business policy gives to shorter rotations periods 25 or less years due to the high organizational risk of the industrial sector in Chile.

3) The Gompertz and Logistic diffusion models present small estimation differences in the growing phase of wood stock, but significant differences in the saturation volume of the wood stock, which should be crucial in the model diffusion selection.

Acknowledgements

The author acknowledge the collaboration of P. Santibañez of MININCO (2006) forest company for providing with the Radiata pine stands harvest near Temuco and to “Modelo Nacional de Simulación de Pino Radiata” for the complementary stands data of the region.

Figure 3. Faustmann gompertz optimal cut.

Table 5. Optimal harvest rotation optimal results, Brown price Gompertz wood stock diffusion.

Cite this paper

EduardoNavarrete, (2015) Optimal Stochastic Pine Stands Harvest Rotation Policies.
*Open Journal of Forestry*,**05**,593-606. doi: 10.4236/ojf.2015.56053

References

- 1. Alvarez, L. R., & Koskela, E. (2007). Optimal Harvesting under Resource Stock and Price Uncertainty. Journal of Economic Dynamics and Control, 31, 2461-2485.
- 2. Amacher, G. S., Brazee, R. J., & Deegan, P. (2011). Faustmann Continues to Yield. Journal of Forest Economics, 17, 231-234. http://dx.doi.org/10.1016/j.jfe.2011.06.001
- 3. Brazee, R.J. (2001). The Faustmann Formula. Forest Science, 47, 44-49.
- 4. Beskos, A., Papaspliopoulos, O., Roberts, G., & Fearnhead, P. (2006). Exact and Computationally Efficient Likelihood-Based Estimation for Discretely Observed Diffusion Processes (with Discussion). Journal of the Royal Statistical Society: Series B (Statistical Methodology), 68, 333-382.
- 5. Chang, S. J. (2001). One Formula Myriad Conclusions, 150
Years of Practicing the Faustmann Formula in Central Europe and the USA. Forest
Policy and Economics, 2, 97-99.

http://dx.doi.org/10.1016/S1389-9341(01)00053-3 - 6. Clarke, H. R., & Reed, W. J. (1989). The Tree Cutting Problem in a Stochastic Environment. Journal of Economics Dynamic and Control, 13, 569-595. http://dx.doi.org/10.1016/0165-1889(89)90004-3
- 7. Faustmann, M. (1995). (Originally 1849). Calculation of the Value which Forest Land and Immature Stands Processess for Forestry. Journal of Forest Economics, 1, 7-44.
- 8. Garcia, O. (2005). Unifying Sigmoid Univariate Growth Equations. Forest Biometry, Modelling and Information Sciences (FBMIS), 1, 63-68.
- 9. Gaffney, M. M. (1957). Concepts of Financial Maturity of Timber and Other Assets. Agricultural Economics Information Series. Raleigth, NC: North Caroline State College.
- 10. Gutierrez, R., Gutierrez-Sanchez, R., & Nafidi, A. (2008). Modelling and Forecasting Vehicle Stocks Using Trends of Stochastic Gompertz Diffusion Models. Applied Stochastic Models in Business and Industry, 25, 385-405.
- 11. Insley, M. (2002). A Real Option Approach to the Valuation of a Forestry on Investment. Journal of Environmental Economics and Management, 44, 471-492.
- 12. Insley, M., & Rollins, K. (2005). On Solving the Multi-Rotational Timber Harvesting Problem with Stochastic Prices: A Linear Complimentarily Formulation. American Journal of Agriculture Economics, 87, 735-755. http://dx.doi.org/10.1111/j.1467-8276.2005.00759.x
- 13. Johnson, T. C. (2006). The Optimal Timing of Investment Decisions. PhD Thesis, London: University of London.
- 14. Kloeden, P., & Platen, E. (1992). Numerical Solution of Stochastic Differential Equation (p. 125). Berlin: Springer-Verlag. http://dx.doi.org/10.1007/978-3-662-12616-5
- 15. Meyer, P., Yung, J., & Ausubel, J. (1999). A Primer on Logistic Growth and Substitution: The Mathematics of the Logolet Lab Software. Technological Foresting and Social Change.
- 16. MININCO (2006). Santibáñez P. Los Angeles.
- 17. Morck, R., & Schwartz, E. (1989). The Valuation of Forestry
Resources under Stochastic Prices and Inventories. Journal of Financial and Quantitative
Analysis, 24, 473-487.

http://dx.doi.org/10.2307/2330980 - 18. Navarrete, E. (2011). Modelling Optimal Pine Stands Harvest under Stochastic Wood Stock and Price in Chile. Forest Policy and Economics, 15, 54-59.
- 19. Navarrete, E., & Bustos, J. (2013). Faustmann Optimal Pine Stands Stochastic Rotation Problem. Forest Policy and Economics, 30, 39-45. http://dx.doi.org/10.1016/j.forpol.2013.02.007
- 20. Samuelson, P. (1976). Economics of Forestry in an Evolving
Economy. Economic Inquiry, 14, 466-491.

http://dx.doi.org/10.1111/j.1465-7295.1976.tb00437.x - 21. Sodal, S. (2002). The Stochastic Rotation Problem: A Comment. Journal of Economics & Control, 26, 509-515. http://dx.doi.org/10.1016/S0165-1889(00)00076-2
- 22. Thijssen, J. J. J. (2010). Irreversible Investment and Discounting: An Arbitrage Pricing Approach. Annals of Finance, 6, 295-315. http://dx.doi.org/10.1007/s10436-008-0108-4
- 23. Willassen, Y. (1998). The Stochastic Rotation Problem: A
Generalization of Faustmann’s Formula to a Stochastic Forest Growth. Journal of
Economic Dynamics and Control, 22, 573-596.

http://dx.doi.org/10.1016/S0165-1889(97)00071-7

Appendix A: Mininco Radiata Pine Stands Harvest Data

Notes: TVOL = Total Volume (mts3/ha); TVOL = POD + DEBO + INDUST + COMERC + PULPA
POD + DEBO = Wood thinning cuts; Logs composition: Saw logs = INDUST + COMERC =
83.9%; Pulp logs = PULP + POD + DEBO = 16.1%; Saturation volume Vs = 601.6 m^{3}/ha,
(3% stands biggest TVOL average).

Appendix B: Radiata Pine Log Prices

Table B1. Nominal Radiata Pine log exportation prices.

Source CONAF-INFOR Chile.

The Price diffusion parameters can very easily be calculated by making the following logarithmic transformation. The parameters are given by Table B2.

Table B2. Log summary.