Abstract. In this study, wood production and hydrologic functions of forests were accommodated within a planning procedure for separate working circles (areas dedicated to certain forest functions) that were delineated according to an Ecosystem-Based Functional Planning approach. Mixed integer goal programming was used as the optimization technique. The timing and scheduling of a maintenance cutting (partial harvest) was the decision variable in the modeling effort, and an original formulation was developed as a multiobjective planning procedure. Four sample planning strategies were developed and model outputs were evaluated according to these strategies. Spatial characteristics of stands were considered, and used to prohibit the regeneration of adjacent stands during the same time period. Because of the positive relationship between qualified water production and standing timber volume in the forest, the model attempts to maximize qualified water * Corresponding author. Hayati Zengin,

© iForest -Biogeosciences and Forestry IntroductionInformation about extant forest communities is of considerable importance in forest management planning. The stand is the smallest production and silvicultural unit, and the most crucial part of this information. The location of planning units, their size and silvicultural features are used to select management objectives and determine the working cycle, management regulations and silvicultural goals. Thus, the precision of decisions to be implemented and the planning to be realized are both closely related to the accuracy of the stand maps. Moreover, map accuracy depends on the type of forest inventory, the sample size, and the techniques used to evaluate the inventory results.Stand maps have previously been produced subjectively (Feng et al. 2006) in conjunction with current forest management practices. In conventional management planning, stands are defined through the planning process. The type of of the management objectives is a component of this process (Gunnarsson et al. 1998). However, the use of computers has changed this conventional approach, and spatially and temporally dynamic description and treatment units can readily be produced (Holmgren & Thuresson 1997). Thus, dynamic forest planning can be achieved with the use of geopositioned field plot data and computers (Wallerman et al. 2002).The sampling methods used in forest inventories vary according to different countries' forestry objectives and forest structures. Systematic sampling methods are recommended for large and homogeneous forest areas and have been implemented in all production and conservation areas of Turkey since 1964. However, if this sampling method is implemented with an insufficient number of sample plots and/or in heterogeneous forests, the results will be highly questionable because they may fail to reflect the true forest composition. Although Sherman (1996), Aurbi & Debouzie (2000), Flores et al. (2003) andD'Orazio (2003) used remote sensing data to improve the results of such inventories, the desired outcome was never achieved in practice.Spatial interpolation methods, which have been classified by Li & Heap (2008) as nongeostatistical, geostatistical and combined, came into use at the end of the 1960s and have been investigated for use in forest management. As described in Akhavan et al. (2010), Guibal (1973 was the first study to use kriging in the forest inventory process. Jost (1993) also used geostatistical methods to compare conventional inventory results based on systematic sampling with the results of the kriging method. Geostatistically based methods that utilize textural information are frequently used to analyze remote-sensing (RS) images (Zawadzki et al. 2005). The quality and quantity of these methods have increased through advances in the computer sciences and in the science of geographical information systems. These methods allow the use of a variety of ecological and technical parameters in the assessment of inventory results. Spatial interpolation and multi-...