TIMBER SUPPLY UNDER DEMAND UNCERTAINTY: WELFARE GAINS FROM PERFECT COMPETITION WITH RATIONAL EXPECTATIONS

ABSTRACT. This paper presents a noneconometric approach to estimating the short‐run timber supply function based on optimal harvest decisions. Determination of optimal harvest levels and estimation of supply function coefficients are integrated into one step by incorporating a parametric short‐run timber supply function into the harvest decision model. In this manner we convert the original harvest decision model into a new optimization problem with the supply function coefficients functioning as “decision variables.” Optimal solution to the new decision model gives the coefficients of the short‐run supply function and, indirectly, the optimal harvest levels. This approach enables us to develop stochastic models of the timber market that are particularly useful for forest sector analysis involving comparison of alternative institutional regimes or policy proposals and when the timber market is affected by stochastic variables. For demonstration purposes, we apply this method to compare the performances of two timber market regimes (perfect competition and monopoly) under demand uncertainty, using the Swedish data. The results show that the expected timber price is 22 percent lower and the expected annual timber supply is 43 percent higher in the competitive market than in the monopoly market. This confirms the theoretical result that monopoly reduces supply and increases price. The expected social welfare gain from perfect competition over monopoly is about 24 percent.     

AN OPTIMIZATION APPROACH TO IDENTIFYING TIMBER SUPPLY FUNCTION COEFFICIENTS

This paper examines an optimization approach to identifying short-run timber supply function coefficients when the form of the supply function is known. By definition, a short-run timber supply function is a functional relationship between the optimal harvest level in each period (e.g., each year) and the actual forest-market state in the same period. The short-run timber supply function represents the optimal harvest decision policy, and therefore, the problem of optimal harvesting can be formulated as a problem of determining this function. When the form of the supply function is known, the problem becomes one of identifying the coefficients of the supply function. If the management objective is to maximize the expected present value of net revenues from timber harvesting over an infinite time horizon, and the timber price process is, in a sense, stationary, the supply function coefficients correspond to the optimal solution to an anticipative optimization problem. In this case, the supply function coefficients can be determined by maximizing the expected present value of the net revenues from timber harvesting, where periodic harvest levels are determined using the supply function. Numerical results show that the short-run supply functions determined using this approach gives good approximations of the true supply function.     

TEMPORAL AGGREGATION AND TESTING FOR TIMBER PRICE BEHAVIOR

ABSTRACT. Different harvest timing models make different assumptions about timber price behavior. Those seeking to optimize harvest timing are thus first faced with a decision regarding which assumption of price behavior is appropriate for their market, particularly regarding the presence of a unit root in the timber price time series. Unfortunately for landowners and investors, the literature provides conflicting guidance on this subject. One source for the ambiguous results of unit root tests of timber prices may involve data problems. We used Monte Carlo simulations to show that aggregating observations below their observed rate resulted in similar power reductions and empirical size distortions across three classes of unit root tests. Moving‐average error structures can also affect power and sizes of tests on period‐averaged data. Such error structures can also be created by the kind of temporal averaging common in reported timber prices. If we take timber prices at their face value and therefore ignore these sampling error and temporal aggregation complications, we find that unit root tests on southern timber prices support a unit root in 158 out of 208 product‐deflation combinations tested, random walks in 38 of the series found to be nonsta‐tionary, and stationarity in none. However, if we recognize temporal aggregation errors, unit root tests more commonly favor stationarity, especially for pulpwood stumpage. Because price trends for sawtimber and pulpwood products may behave differently even in the same region, stochastic harvest timing models must be developed that allow their multiple products to follow different price paths.     

The optimal harvesting problem under price uncertainty

In this paper we study the exploitation of a one species forest plantation when timber price is governed by a stochastic process. The work focuses on providing closed expressions for the optimal harvesting policy in terms of the parameters of the price process and the discount factor, with finite and infinite time horizon. We assume that harvest is restricted to mature trees older than a certain age and that growth and natural mortality after maturity are neglected. We use stochastic dynamic programming techniques to characterize the optimal policy and we model price using a geometric Brownian motion and an Ornstein–Uhlenbeck process. In the first case we completely characterize the optimal policy for all possible choices of the parameters. In the second case we provide sufficient conditions, based on explicit expressions for reservation prices, assuring that harvesting everything available is optimal. In addition, for the Ornstein–Uhlenbeck case we propose a policy based on a reservation price that performs well in numerical simulations. In both cases we solve the problem for every initial condition and the best policy is obtained endogenously, that is, without imposing any ad hoc restrictions such as maximum sustained yield or convergence to a predefined final state.     

A cutting plane method for solving harvest scheduling models with area restrictions

We describe a cutting plane algorithm for an integer programming problem that arises in forest harvest scheduling. Spatial harvest scheduling models optimize the binary decisions of cutting or not cutting forest management units in different time period subject to logistical, economic and environmental restrictions. One of the most common constraints requires that the contiguous size of harvest openings (i.e., clear-cuts) cannot exceed an area threshold in any given time period or over a set of periods called green-up. These so-called adjacency or green-up constraints make the harvest scheduling problem combinatorial in nature and very hard to solve. Our proposed cutting plane algorithm starts with a model without area restrictions and adds constraints only if a violation occurs during optimization. Since violations are less likely if the threshold area is large, the number of constraints is kept to a minimum. The utility of the approach is illustrated by an application, where the landowner needs to assess the cost of forest certification that involves clear-cut size restrictions stricter than what is required by law. We run empirical tests and find that the new method performs best when existing models fail: when the number of units is high or the allowable clear-cut size is large relative to average unit size. Since this scenario is the norm rather than the exception in forestry, we suggest that timber industries would greatly benefit from the method. In conclusion, we describe a series of potential applications beyond forestry.     

WHEN TO CUT A STAND OF TREES?

This paper revisits the debate over the economic optimality of different timber harvest rules. The traditional Faustmann and Maximum Sustained Yield rotation determinations are confined to a deterministic world. Once stochasticity is introduced into the model formulation and the additional rental and management costs due to postponement of harvest are taken into account, we find that the optimal stopping time becomes random and varies in response to changes in the underlying price and growth processes. As a result, this stochastic optimal stopping time is bounded by the Faustmann cutting age from below if variabilities of the stochastic processes diminish to zero, but not necessarily by the MSY rotation from above.     

A transient stochastic simulation–optimization model for operational fuel planning in-theater

Army fuel planners are responsible for developing daily loading plans that specify which tankers to load, with what fuel, and where to send the loaded tankers. The tools used to accomplish this task are custom built spreadsheets which require large amounts of time and effort to use, update, and keep free of errors. This research presents a transient stochastic simulation–optimization model of the in-theater bulk fuel supply chain, where the simulation model is used to simulate the performance of the fuel supply chain under a particular fuel distribution policy and the optimization portion is used to update the policy so that it results in the performance desired by the Army fuel planner. The fuel distribution policy can then be used to derive the daily loading plan. Due to the multi-objective nature of the problem, the set of policies that form the efficient frontier are all candidate policies for the Army fuel planner to select from. Results of experimentation with a wide variety of supply chain scenarios indicate that, for a given supply chain scenario, the optimization portion of the model identifies a set of fuel distribution policies that address the objectives of the Army fuel planner. In addition, the simulation–optimization model comfortably solves the largest supply chain scenarios the Army fuel planner would reasonably be expected to encounter.     

Optimal harvest strategy for even‐aged stands with price uncertainty and risk of natural disturbances

We present a reservation price model to examine the joint impacts of natural disturbances and stumpage price uncertainty on the optimal harvesting decision for even‐aged forest stands. We consider a landowner who manages a loblolly pine stand to produce timber and amenities, under age‐dependent risk of wildfires and uncertainty in future timber prices. We show that the incorporation of risk of wildfires decreases the optimal reservation prices. The inclusion of risk of wildfires leads to lower land values and reduces the mean harvest age compared with the case of no‐risk of wildfires. Higher economic gains are obtained with the reservation price strategy compared with the deterministic rotation age model—a difference in the land value of $2,326 ha^?1 (21%) between the two approaches. Recommendations for Resource Managers

• Our adaptive harvest strategy shows that the incorporation of risk of wildfires decreases the optimal reservation prices compared with the case of no‐risk of wildfires.
• Low reservation prices—a price that makes the landowner indifferent between harvesting or waiting longer—result in lower economic benefits for landowners and potential conversions of lands to nonforest use.
• Forest management practices oriented to reduce the effects of catastrophic disturbances, for example, creating a more complex forest structure with different stand densities, become imperative to ensure the sustainability of forestlands in the US South.
• Our analysis also suggests that the valuation of forestry investments should consider not only the risk of catastrophic events but also uncertainty in future timber prices. Higher appraisals of land value are obtained when timber price uncertainty is explicitly recognized, providing financial incentives for landowners to invest in forestry.

     

面对成比例交易成本的情况下外汇期权的定价

The purpose of present work is to examine the financial problem of finding the universal reservation prices of a European call option written on exchange rate when there is proportional transaction costs of trading foreign currency in the market. An approach is suggested to compute the reservation bid-ask price of foreign currency call option based on maximizing the investor's expected utility. Option prices are determined from the investor's basic portfolio selection problem, without the need to solve a more complex optimization problem involving the insertion of the option payoffs into the terminal value function. Option prices are computed numerically in a Markov chain approximation for the case of exponential utility.Numerical results show that the option price bounds are almost independent of the alternative risk aversion parameter, but the bounds of NT region becomes narrower and the range of values of the initial holding for which the fair price lies within the bid-ask spread is shifted to a lower value when the risk aversion parameter increases.     

FOREIGN CURRENCY OPTION PRICING WITH PROPORTIONAL TRANSACTION COSTS

The purpose of present work is to examine the financial problem of finding the universal reservation prices of a European call option written on exchange rate when there is proportional transaction costs of trading foreign currency in the market. An approach is suggested to compute the reservation bid-ask price of foreign currency call option based on maximizing the investor＇s expected utility. Option prices are determined from the investor＇s basic portfolio selection problem, without the need to solve a more complex optimization problem involving the insertion of the option payoffs into the terminal value function. Option prices are computed numerically in a Markov chain approximation for the case of exponential utility. Numerical results show that the option price bounds are almost independent of the alternative risk aversion parameter, but the bounds of NT region becomes narrower and the range of values of the initial holding for which the fair price lies within the bid-ask spread is shifted to a lower value when the risk aversion parameter increases.     

Medium range optimization of copper extraction planning under uncertainty in future copper prices

Deterministic mine planning models along a time horizon have proved to be very effective in supporting decisions on sequencing the extraction of material in copper mines. Some of these models have been developed for, and used successfully by CODELCO, the Chilean state copper company. In this paper, we wish to consider the uncertainty in a very volatile parameter of the problem, namely, the copper price along a given time horizon. We represent the uncertainty by a multistage scenario tree. The resulting stochastic model is then converted into a mixed 0–1 Deterministic Equivalent Model using a compact representation. We first introduce the stochastic model that maximizes the expected profit along the time horizon over all scenarios (i.e., as in a risk neutral environment). We then present several approaches for risk management in a risk averse environment. Specifically, we consider the maximization of the Value-at-Risk and several variants of the Conditional Value-at-Risk (one of them is new), the maximization of the expected profit minus the weighted probability of having an undesirable scenario in the solution provided by the model, and the maximization of the expected profit subject to stochastic dominance constraints recourse-integer for a set of profiles given by the pairs of target profits and bounds on either the probability of failure or the expected profit shortfall. We present an extensive computational experience on the actual problem, by comparing the risk neutral approach, the tested risk averse strategies and the performance of the traditional deterministic approach that uses the expected value of the uncertain parameters. The results clearly show the advantage of using the risk neutral strategy over the traditional deterministic approach, as well as the advantage of using any risk averse strategy over the risk neutral one.     

ECONOMIC EVALUATION OF BIOTECHNOLOGICAL PROGRESS: THE EFFECT OF CHANGING MANAGEMENT BEHAVIOR

Abstract The paper assesses the welfare effects of biotechnological progress, as exemplified by tree improvements, using a partial equilibrium model. Timber demand is assumed to be stochastic and the distributions of its coefficients known. The coefficients of a log‐linear supply function are determined by maximizing the expected present value of the total surplus of timber production, both in the presence and in the absence of genetically improved regeneration materials. The supply functions are then used to estimate the expected present values of the total surplus in different cases through simulation. These estimates enable us to assess the direct effect and the effect of changing harvest behavior on the expected present value of the total surplus. The main results of the study are (i) the presence of genetically improved regeneration materials has significant impacts on the aggregate timber supply function; (ii) the application of genetically improved regeneration materials leads to a significant increase in the expected present value of the total surplus; and (iii) a considerable proportion of the welfare gain results from the change in harvest behavior. A conclusion we draw from this study is that ignoring the influences of technological and policy changes on behavior can lead to significantly biased welfare estimates. We view the model as a potential approach to conducting counterfactual policy comparisons in economics without forward‐looking data.     

VALUING A TIMBER HARVEST CONTRACT AS A HIGH‐DIMENSIONAL AMERICAN CALL OPTION VIA LEAST‐SQUARES MONTE CARLO SIMULATION

Abstract Industrial timberland ownership in the United States has shifted substantially in the last 20 years. Having sold their fee‐owned timberlands, forest products companies relied heavily on the open market for raw timber. To reduce their exposure to market risks, however, forest products companies have been using a number of supply chain instruments, such as timber harvest contracts. As these vehicles become increasingly important to the forest industry, it is necessary and important to determine their economic values. In this study, we treated a 3‐year timber harvest contract on a 30‐year‐old loblolly pine plantation as a high‐dimensional American call option and calculated its value by the least‐squares Monte Carlo simulation technique. The estimated values of such a contract ranged from $1,693/ac to $1,984/ac under two timber price assumptions. With reasonable starting timber prices and strike price in the simulation, random timber prices led to higher contract values. Results from this study can help private landowners, timber brokers, and forest products companies better manage their business risks.     

Timber harvest scheduling with price uncertainty using Markowitz portfolio optimization

Harvest scheduling models need to account for uncertain revenue predictions when minimizing risk of financial loss is an important management objective. In this paper, we present methods for estimating the means and covariances of stumpage prices and incorporating them in harvest scheduling models. We approached the estimation problem by fitting time-series models to loblolly pine sawtimber and pulpwood stumpage prices in Georgia, USA, and deriving formulas for means and covariances of price predictions. Statistical evidence supported integrated autoregressive models, which caused covariances of price predictions to increase with time. The means and covariances of price predictions were combined with timber yield and land value predictions to give exact formulas for the revenue means and covariances of timber management activities. Sawtimber regimes dominated pulpwood regimes by providing higher mean revenues across a wide range of revenue variances. Harvest scheduling results for a hypothetical forest of pine plantations showed that the forest plan that maximized mean income without concern for risk (expressed as the standard deviation of income) involved sawtimber production with a 35-year rotation age. Risk was reduced 30% with little effect on mean income by using shorter-rotation sawtimber regimes. Risk was reduced 80% by using a mix of short-rotation sawtimber and pulpwood regimes because pulpwood price was only weakly correlated with sawtimber price. The latter risk-reduction came at the expense of mean income, which was reduced by as much as 50%. The risks and compositions of optimal forest plans were extremely sensitive to assumptions about the range of future prices that were inherent in different prediction models. This sensitivity emphasizes the importance of carefully determining the decision makers beliefs about stumpage price behavior.     

The Sample Average Approximation Method Applied to Stochastic Routing Problems: A Computational Study

The sample average approximation (SAA) method is an approach for solving stochastic optimization problems by using Monte Carlo simulation. In this technique the expected objective function of the stochastic problem is approximated by a sample average estimate derived from a random sample. The resulting sample average approximating problem is then solved by deterministic optimization techniques. The process is repeated with different samples to obtain candidate solutions along with statistical estimates of their optimality gaps.We present a detailed computational study of the application of the SAA method to solve three classes of stochastic routing problems. These stochastic problems involve an extremely large number of scenarios and first-stage integer variables. For each of the three problem classes, we use decomposition and branch-and-cut to solve the approximating problem within the SAA scheme. Our computational results indicate that the proposed method is successful in solving problems with up to 2¹⁶⁹⁴ scenarios to within an estimated 1.0% of optimality. Furthermore, a surprising observation is that the number of optimality cuts required to solve the approximating problem to optimality does not significantly increase with the size of the sample. Therefore, the observed computation times needed to find optimal solutions to the approximating problems grow only linearly with the sample size. As a result, we are able to find provably near-optimal solutions to these difficult stochastic programs using only a moderate amount of computation time.     

CORO,a modeling and an algorithmic framework for oil supply,transformation and distribution optimization under uncertainty

We present a modeling framework for the optimization of a multiperiod Supply, Transformation and Distribution (STD) scheduling problem under uncertainty on the product demand, spot supply cost and spot selling price. The Hydrocarbon and Chemical sector has been chosen as the pilot area, but the approach has a far more reaching application. A deterministic treatment of the problem provides unsatisfactory results. We use a 2-stage scenario analysis based on a partial recourse approach, where the STD policy can be implemented for a given set of initial time periods, such that the solution for the other periods does not need to be anticipated and, then, it depends on the scenario to occur. In any case, it takes into consideration all the given scenarios. Novel schemes are presented for modeling multiperiod linking constraints, such that they are satisfied through the scenario tree; they are modeled by using a splitting variable scheme, via a reduntant circular linking representation.     

Vehicle routing with stochastic time-dependent travel times

Assigning and scheduling vehicle routes in a stochastic time-dependent environment is a crucial management problem. The assumption that in a real-life environment everything goes according to an a priori determined static schedule is unrealistic. Our methodology builds on earlier work in which the traffic congestion is captured in an analytical way using queueing theory. The congestion is then applied to the VRP problem. In this paper, we introduce the variability in traffic flows into the model. This allows for an evaluation of the routes based on the uncertainty involved. Different experiments show that the risk taking behavior of the planner can be taken into account during optimization. As more weight is given to the variability component, the resulting optimal route will take a slightly longer travel time, but will be more reliable. We propose a powerful objective function that is easily implemented and that captures the trade-off between the average travel time and its variance. The evaluation of the solution is done in terms of the 95th-percentile of the travel time distribution (assumed to be lognormal), which reflects well the quality of the solution in this stochastic time-dependent environment.     

A system for the design of short term harvesting strategy

Short term harvesting requires decisions on which stands to harvest, what timber volume to cut, what bucking patterns (how to cut up the logs) to apply to logs in order to obtain products that satisfy demand and which harvesting machinery to use. This is an important problem in forest management and difficult to solve well in satisfying demand, while maximizing net profits. Traditionally, foresters have used manual approaches to find adequate solutions, which has shortcomings both in time spent by analysts and the quality of solutions. Since demand for timber products is defined by length, diameter and quality of each piece, this leads to a complex combinatorial problem in matching supply (standing trees) and demand. We developed one of the few reported approaches for solving the short term harvesting problem based on a computerized system, using a linear programming approach. Determining adequate bucking patterns is not trivial. We develop a column generation approach to generate such patterns. The subproblem is a specially designed branch and bound scheme. The generation of bucking patterns implemented within the LP formulation led to a significant improvement of solutions. We believe this is the first system implemented with this level of detail. This system has been advantageously implemented in several forest companies. The results obtained show improvements obtained by the firms of 5–8% in net revenues over traditional manual approaches.     

Design of insurance contracts using stochastic programming in forestry planning

This work addresses a tactical planning problem faced by a forestry firm, deciding which timber units to harvest and what roads to build to obtain the greatest possible benefits. We include uncertainty in prices by means of utility theory. This enables solutions to be found that the firm finds preferable to those obtained when risk aversion is ignored and makes it possible to design insurance contracts that benefit the firm while also being attractive to an insurer. Two types of contract are designed; one dependent on the firm’s operating result and the other independent of it. Metrics are then developed to quantify the benefits conferred by a contract, demonstrating that the latter contract type dominates the former. These results are then illustrated by applying them to a simplified planning problem of a forest owned by the Chilean forestry operator Millalemu.     

Modeling target volume flows in forest harvest scheduling subject to maximum area restrictions

In forest harvest scheduling problems, one must decide which stands to harvest in each period during a planning horizon. A typical requirement in these problems is a steady flow of harvested timber, mainly to ensure that the industry is able to continue operating with similar levels of machine and labor utilizations. The integer programming approaches described use the so-called volume constraints to impose such a steady yield. These constraints do not directly impose a limit on the global deviation of the volume harvested over the planning horizon or use pre-defined target harvest levels. Addressing volume constraints generally increases the difficulty of solving the integer programming formulations, in particular those proposed for the area restriction model approach. In this paper, we present a new type of volume constraint as well as a multi-objective programming approach to achieve an even flow of timber. We compare the main basic approaches from a computational perspective. The new volume constraints seem to more explicitly control the global deviation of the harvested volume, while the multi-objective approach tends to provide the best profits for a given dispersion of the timber flow. Neither approach substantially changed the computational times involved.