Adaptive Procurement Guidelines for Automatic Selection of Renewable Forest Energy Sources within a Sustainable Energy Production System

An automatic forest-stand selection method was developed that integrates the procurement of profitable energy sources within a sustainable energy production system. We tested the method using a forest harvester simulator. We found that site-specific estimates of forest characteristics are important when predicting the potential of an energy-wood stand as a renewable energy source. Further, tree parameters can be combined with automatic stem measurements from a multiple-tree harvester to predict the energy-wood biomass in the stands. The selection process uses data from profitability studies. The selection process also uses environmental criteria to ensure that sufficient soil organic matter is left behind and to protect the soil against erosion. The integrated system of the harvester automatically adapts the system’s models and stand-selection rules to account for various site-specific stand parameters. Predicting the profitable and environmentally acceptable yield of stand biomass has great potential in sustainable forest resource management, but managers must decide whether the operational procurement guidelines provided by the stand-selection method is acceptable under their local real-world wood procurement conditions.


Introduction
In most countries in Europe, biomass resources are small compared to the European energy balance [1]. However, in Finland, forest biomass provides the potential of a 50% increase of renewable fuels for distributed bio-based economy. Along with economic sustainability, based on the idea of continuous forest yield, environmental sustainability has become a serious issue in Finland's forest management debate [2]. Since the mid-1990s, commitments to sustainable energy-wood cutting have strengthened in Finnish forestry. Therefore, the development of future district heating systems and technologies involves forest fuel procurement measures as an important part of such technology [3].
The design of sustainable energy systems, including regional resource management systems, are described in a number of recent studies [4][5][6][7][8][9]. Typically, Geographic Information Systems (GISs) are used in energy planning. The concepts of sustainable regional source and land use management have been presented by Yamamoto et al. [4]. Their system analyzes the overall biomass flows and geographic region resource management is strongly linked to the land management strategy for sustainable development, in which the biomass flow has the characteristics of spatial planning [7]. The strategic energy planning is usually focused on social and economic needs with integral considerations of the environment to be sustainable in the long-term environment.
In the short-term planning, the local and distributed nature of the forest biomass resource and it's usually low energy density require energy supply chains with large harvesting (cutting and ‚ from which wood harvested is only energy wood and/or pulpwood, and ‚ where the mean diameter at breast height is less than 16 cm, but more than 3 cm, and ‚ where the height of dominant trees is more than 7 m in Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) dominated stands, and more than 9 m in broadleaf stands.
In practice, young stands also contains larger trees, which must be cut in mechanized first thinning. Around two thirds of the wood harvested from young stands come from Scots pine-dominated stands in Finland. During the 2000 s, the average stem size of young stands has been 82 dm 3 , the average industrial pulpwood removal 47 m 3¨h a´1, and the forest haulage distance 295 m. In 2014, the average harvesting (cutting and forest haulage) costs of pulpwood in first thinnings carried out by the Finnish forest industries was 17.7 €¨m´3. When producing whole-tree chips from young stands for energy generation, the total supply chain costs are approximately 18-24 €¨MWh´1 of which the harvesting costs were around 40%-50% [13,19].
In Finnish forests, energy-wood cutting is carried out only on suitable sites where stands meet the criteria for short-term environmental sustainability [20]. For example, stands with a diversity of less than required are not harvested using energy wood assortments. Suitable stands must also meet economic criteria set by the harvesting entrepreneur and the wood procurement company for profitable cutting, such as having acceptable logistics costs for transport of the energy wood to energy plants. Typically, the harvested areas must be located within 100 km from the energy plants [5]. In this study, profitability of energy-wood cutting will be integrated with selection of trees in stand harvesting planning. To predict reliable yield of stand for an energy-wood entrepreneur, selection of trees will also be integrated with biomass calculations by model of Repola et al. [21]. Obtaining reliable predictions of stand biomass requires local site-specific guidelines that can be applied to guide the management of stand harvesting [22]. Nowadays, energy-wood cutting is a common practice in Finland, and efficient methods for forest silviculture have been developed to make this form of cutting sustainable [20]. Energy-wood cutting can also be combined with pulpwood cutting [23], although this increases the cost of conventional cutting. After energy-wood cutting, there are short-term effects on sustainability, and these effects are soil-, site-, and practice-specific [16][17][18]. The effects can be managed and controlled by a harvester operator. In practice, harvester operators are charged with deciding how many trees to leave per hectare, and they select trees conventionally without assistance from decision-support and control software. Therefore, sustainability problems arise as a result of the selection of trees by harvester operators.
Recently, researchers have developed operational stand-selection systems and suggested that the tree-selection rules that are in current use need improvement. For instance, control of cutting could be based on site-specific measurements of each tree in a stand [23,24]. In addition to the energy wood volume, the exact stem diameter and location (using an onboard GPS receiver and machine vision) could be determined based on these measurements. These measurements can be obtained and analyzed automatically using work study software, such as John Deere's TimberLink [25,26]. In the present study, we developed software to identify data from each stem that could be used to develop potential operating rules to increase the efficiency and environmental sustainability of energy-wood cutting. We then tested the potential of this software (system) in simulation of stand selection using real production data collected by a multiple-tree harvester's computer.
In Finland, wood procurement organizations use automatic logistics systems to achieve their procurement targets while minimizing operating costs. Therefore, managers need accurate estimates of stand characteristics for their management information systems so they can accurately quantify the amount of energy wood that is available as a profitable energy source. To meet their needs, the present study was designed to construct an automatic adaptive system to select profitable stands from suitable stands that would permit environmentally sustainable energy-wood cutting within an operational stand-selection system ( Figure 1). We solved the research problem by integrating site-specific stem-related information, recorded by a John Deere multiple-tree harvester with profitability information obtained for the same harvester. The solution also provides practical tree-selection rules used to cut the young stands by the harvester operator.
Energies 2016, 9,155 3 these effects are soil-, site-, and practice-specific [16][17][18]. The effects can be managed and controlled by a harvester operator. In practice, harvester operators are charged with deciding how many trees to leave per hectare, and they select trees conventionally without assistance from decision-support and control software. Therefore, sustainability problems arise as a result of the selection of trees by harvester operators.
Recently, researchers have developed operational stand-selection systems and suggested that the tree-selection rules that are in current use need improvement. For instance, control of cutting could be based on site-specific measurements of each tree in a stand [23,24]. In addition to the energy wood volume, the exact stem diameter and location (using an onboard GPS receiver and machine vision) could be determined based on these measurements. These measurements can be obtained and analyzed automatically using work study software, such as John Deere's TimberLink [25,26]. In the present study, we developed software to identify data from each stem that could be used to develop potential operating rules to increase the efficiency and environmental sustainability of energy-wood cutting. We then tested the potential of this software (system) in simulation of stand selection using real production data collected by a multiple-tree harvester's computer.
In Finland, wood procurement organizations use automatic logistics systems to achieve their procurement targets while minimizing operating costs. Therefore, managers need accurate estimates of stand characteristics for their management information systems so they can accurately quantify the amount of energy wood that is available as a profitable energy source. To meet their needs, the present study was designed to construct an automatic adaptive system to select profitable stands from suitable stands that would permit environmentally sustainable energy-wood cutting within an operational stand-selection system ( Figure 1). We solved the research problem by integrating site-specific stem-related information, recorded by a John Deere multiple-tree harvester with profitability information obtained for the same harvester. The solution also provides practical tree-selection rules used to cut the young stands by the harvester operator.

Harvester Stem Data
The automatic stem data was collected from one stand (i.e., logging site) in Orimattila in Finland (60°48'P, 25°44'E). Thinning of the stand was conducted using a John-Deere multiple-tree harvester, and the area of the thinning site was 8.8 ha. Cutting instructions and intensities were equal on the site, but site volume (m 3 ·ha −1 ) varied inside the site before cutting. Therefore, the total removal of Figure 1. The stand-selection system. u, stem diameter; w, environmental parameter; BT, volume from the biomass equation; CBT, volume from the calibrated biomass equation; B, the stand biomass (m 3 per stand); y, cutting times; c, costs; r, revenues; β 0 and β 1 , regression parameters.

Harvester Stem Data
The automatic stem data was collected from one stand (i.e., logging site) in Orimattila in Finland (60˝48'P, 25˝44'E). Thinning of the stand was conducted using a John-Deere multiple-tree harvester, and the area of the thinning site was 8.8 ha. Cutting instructions and intensities were equal on the site, but site volume (m 3¨h a´1) varied inside the site before cutting. Therefore, the total removal of energy and pulp wood from the stand ranged from 35 m 3¨h a´1 to 74 m 3¨h a´1 and averaged 59.9 m 3¨h a´1. The harvester data consisted of 5125 stems: Norway spruce (9%), Scots pine (73%), and Silver birch (Betula pendula Roth., 18%). The measurements for each stem were stored in files that were transmitted daily to the entrepreneur's enterprise-resource-planning system via a mobile phone network connection. The harvester's sensor recorded stem diameters along the pulp wood stem at 10-cm intervals, starting from the lowest cutting point and ending at the last cutting height. The harvester recorded one diameter along the energy-wood stem at the cutting point. In the calculation system used in this study [24], the stem data were combined into a single file (the stem database), which contained information on the tree species, stem diameters and stem height.
For pulp wood stems we calculated the diameter at breast height (dbh, 130 cm) by linear interpolation between the measurement diameters above and below this height, with the distance between the measurements on the stem used to estimate the position of dbh between these two diameters. For energy-wood stems we also used the time and motion study data in the calculation of these stem characteristics. Tree height included parts above the highest cut that were not harvested. To derive the relative position of each diameter measurement along the stem, the total tree height must be estimated. To do so, we estimated the length of the top part of the stem (i.e., the part above the highest cut position) for each tree using the model published by Varjo [27]. The total tree height was computed as the sum of the height at the top cut and the estimated length of the top part of the stem above this cut based on the harvester measurements and Varjo's model.

Harvester's Profitability Data in Operational Planning
Profitability data were collected in 12 stands in 2009, when an experienced harvester operator carefully selected trees either for multiple-tree processing or single-tree processing. We used a general work cycle model of cutting that also facilitated the use of automatically recorded data [28]. The data were used in the cost calculation of the stand-selection system. The energy-wood cutting costs amounted to 65 €/E15-hour. The price of the harvester was 230,000 € (at 0% VAT). The total annual work cycle was 2240 operating hours (E15-hours, which include delay times shorter than 15 min). The depreciation period for the base machine and the harvesting head were both set to 5.0 years. The operator's salary was set to 14.7 €/E15-hour, including indirect salary costs (59%) and including an amount for compensated travel time. The fuel consumption for the machine was 12 L/E15-hour at 1.07 €/L [29]. Repair and service costs were estimated to be 3.0 €/E15-hour, administration, maintenance, and insurance costs were 4400 € per year.
The energy-wood cutting cost decreases significantly as stem diameter increases [30]. For example, the cost decreases by 30% as the stem diameter increases from 5 cm to 20 cm, representing a change from 35 €¨m´3 Sob to 25 €¨m´3 Sob (solid volume, scaled outside the bark), respectively. When the stem diameter is 30 cm, the energy-wood cutting costs are 22.5 €¨m´3 Sob. When the stand diameter is 40 cm, the energy-wood cutting costs are 20 €¨m´3 Sob. When energy-wood cutting is integrated with pulp wood cutting, this operation decreases the costs of cutting by 5% to 10%, compared to cutting without the integrated technique [23,31]. In practice, wood procurement companies pay entrepreneurs an average cutting revenue that was 8.9 €¨m´3 sob in the present study.

Stand-Selection System
In our simulation study, one young stand was cut stem by stem using the stand-selection system. The profitability of stem cutting (by the entrepreneur) and of stand procurement (by the energy producer) must both be considered in the stand-selection system, because the productivity and cost-efficiency of stem cutting increases with increasing stem diameter [13,24]. In addition, energy producers pay the entrepreneur. To test the system and calculate the entrepreneur's profitability, we developed a method from a recent automatic time and motion study [24]. The stem-cutting time (s) was calculated for the diameter at the dbh, and the resulting operational profitability data were analyzed. Yield of stand was calculated by model of Repola et al. [21]. In this study, to maintain entrepreneur's profitability, the average hourly operating cost must be less than 65 €/E15-hour for the harvester.
The stem-selection process mined stem data from the stand into energy wood and pulp wood files. The selection process identified energy-wood stem, following the company's selection criteria for sustainable energy production. In practice, due to environmental sustainability criteria and watershed protection guidelines, tops and branches from the largest 75% of the trees were left in the forest. We also left small trees (Norway spruce) to increase biodiversity and to avoid a decrease in soil C. The heuristic energy-wood rules used to select stems for multiple-tree processing were the following:

‚
Step 1: The tree species is not spruce and stem belongs to the energy-wood diameter (<16 cm) frequency distribution.

‚
Step 2: The diameter at dbh height is less than or equal to 8 cm.

‚
Step 3: The diameter of next tree at dbh height is less than or equal to 8 cm.
In this process, the tree species and diameter can be considered as a kind of economic criterion for energy-wood procurement by means of stem cutting. To mine the stem data more efficiently for stand selection, we developed a process that integrated the profitability criterion with the criteria for environmentally sustainable energy-wood cutting (Figure 1). We computed profitability using a harvester simulator's onboard computer [32]. The onboard computer was similar to that used by machines in the forests, and the software was programmed accordingly.
We developed software that let us automatically obtain detailed data for each stem. The software identified the tree characteristics used for stand selection based on the protocols of the adaptive system controller (Figure 2). During the energy-wood cutting experiment, the software automatically monitored the diameter and compared it with the stand-selection rules for the stand. The resulting data were also automatically downloaded to a database after the experiment. To facilitate our analysis the software displayed work reports that let us compare stand-cutting expenses and revenues across a range of diameters.
Energies 2016, 9,155 The stem-selection process mined stem data from the stand into energy wood and pulp wood files. The selection process identified energy-wood stem, following the company's selection criteria for sustainable energy production. In practice, due to environmental sustainability criteria and watershed protection guidelines, tops and branches from the largest 75% of the trees were left in the forest. We also left small trees (Norway spruce) to increase biodiversity and to avoid a decrease in soil C. The heuristic energy-wood rules used to select stems for multiple-tree processing were the following:


Step 1: The tree species is not spruce and stem belongs to the energy-wood diameter (<16 cm) frequency distribution.  Step 2: The diameter at dbh height is less than or equal to 8 cm.  Step 3: The diameter of next tree at dbh height is less than or equal to 8 cm.
In this process, the tree species and diameter can be considered as a kind of economic criterion for energy-wood procurement by means of stem cutting. To mine the stem data more efficiently for stand selection, we developed a process that integrated the profitability criterion with the criteria for environmentally sustainable energy-wood cutting ( Figure 1). We computed profitability using a harvester simulator's onboard computer [32]. The onboard computer was similar to that used by machines in the forests, and the software was programmed accordingly.
We developed software that let us automatically obtain detailed data for each stem. The software identified the tree characteristics used for stand selection based on the protocols of the adaptive system controller (Figure 2). During the energy-wood cutting experiment, the software automatically monitored the diameter and compared it with the stand-selection rules for the stand. The resulting data were also automatically downloaded to a database after the experiment. To facilitate our analysis the software displayed work reports that let us compare stand-cutting expenses and revenues across a range of diameters.

Results
Several stand-selection processes were analyzed based on tree characteristics. The best process classified trees into two groups: trees of which stems, top parts and branches were probably cut for energy, and trees (pulpwood) of which top parts and branches were probably left in the forest. The

Results
Several stand-selection processes were analyzed based on tree characteristics. The best process classified trees into two groups: trees of which stems, top parts and branches were probably cut for energy, and trees (pulpwood) of which top parts and branches were probably left in the forest. The process identified trees by integrating the economic selection criteria with the criteria for environmentally sustainable energy-wood cutting (Figure 1).
Initially, trees were defined as energy wood if they were not spruce, had a diameter at dbh height ď8 cm and belonged to small trees of the diameter frequency distribution. In practice, the manager can examine the frequency distribution for diameter at dbh height that is generated by the harvester computer and can determine whether a given sample or stand met the criteria using the information in Tables 1 and 2. Our control system computed the acceptable parameters for harvestable stems in stand following the protocol in Figure 2. To study the procurement guidelines, we analyzed site-specific data. The best-fit rules classified trees using stand diameter classes, costs, revenues and profits. The entrepreneur's cutting revenue must equal or exceed the cutting cost to earn a profit, and thus must be at least 65 €/E15-hour. As an example of the stand selection process, Table 2 summarizes the cutting profits for the observed stem-diameter classes based on the average diameters in each class for the stems. In this particular sample of stand, the average revenues from the two smallest diameter classes were under the minimum profitability criterion and stems with these diameters might not be harvested for energy. However, these stems should be harvested for energy, because the stand's total profit (sum of the stand profits of the diameter classes) was positive (25 €). Therefore, stems were defined as harvestable, when they had a diameter at dbh height ď8 cm. This diameter also became the new maximum energy-wood diameter in the stand-selection process, if stems must belong to the energy-wood diameter frequency distribution.

Discussion
In the present study, we analyzed the energy-wood data recorded by the multiple-tree harvester that operated in a young stand. To predict the sustainable yield of energy-wood biomass, we developed a stand-selection method that integrates profitability for energy-wood cutting with sustainability criteria. Real stand data collected by harvesters were used in the stand-selection method. Based on this representative stand, we found that both the tree characteristics and site-specific estimates of stands were important factors when predicting the potential of a stand as a renewable energy source (Figure 1).
The stand-selection system automatically adapted the models and the stand-selection rules to account for site-specific parameters. For harvesters, such advanced monitoring systems have been available for several years, and the software can collect and process CAN-bus (Controller Area Network) data about the operating conditions and performance, including stem dimensions [23,25,27]. As in previous studies, we applied this monitoring system to analyze the performance of the stand-selection rules. Currently, managers must decide whether their stand cutting conditions meet the sustainability criteria. However, this kind of decision-making is time-consuming and subject to human error. Using the integrated system developed in the present study would eliminate these problems by allowing the harvester software to automatically apply manager-defined criteria without requiring the use of additional management information systems or the manager's time. In addition to the profitability of the energy-wood cutting, the profitability of other procurement operations (e.g., forwarding, road transportation) can be integrated with the sustainability criteria. For instance, forwarding and road transportation costs can be subtracted from energy production revenues.
The results suggest that, despite the large number of potential error factors during the stand-selection process, the selection rules can be adjusted during the stand identification phase to improve the efficiency of collection of stands for energy production. By integrating the selection rules within the system, our method revealed operational procurement guidelines that could help managers to select forest stands in which energy wood can be harvested economically and in an environmentally sustainable way ( Figure 1, Table 1). Based on the profitability calculations without cutting subsidies, the economic sustainability of pulpwood cutting was acceptable for stems larger than 12 cm in diameter (Table 2), which was con1sistent with expectations, because an increase in an average stand diameter increases both cutting productivity and cost-efficiency [13]. It must be noted that pulpwood is more expensive than energy wood. Further, the economic and environmental sustainability criteria resulted in the selection of stands for energy-wood cutting, when a stem diameter at dbh was smaller or equal than 8 cm and the profitability calculations included cutting subsidies (+7 €¨m´3). These stems should be harvested for energy, because the stand's total profit (sum of the stand profits of the diameter classes) was positive (25 €) ( Table 2). On the other hand, all diameter classes of stand could be defined as harvestable for energy by procurement manager or harvester entrepreneur. Our method could also support harvester operators, because heuristic energy-wood rules used to select stems for multiple-tree processing are adjusted to managers' needs.
Palander and Perho [23] suggest that the harvester operator uses multiple-tree processing as the average diameter and tree height were smaller than 10 cm and 15 m, respectively. The operator also used multiple-tree processing for stems suitable for pulp wood cutting. These results support the present results. Significant findings of this study are that the efficiency of multiple-tree processing can be promoted in energy-wood cutting by careful selection of stands. Using the developed system, the efficiency of the multiple-tree processing can be compared to single-tree wood cutting technology. We did not study the effects of stand density as a potential stand character, but Palander and Perho [23] found that the harvester operator selected multiple trees into the harvester head as stand density was more than 2000 trees¨ha´1, and multiple-tree processing was more efficient than single-tree processing as stand density was more than 4000 trees¨ha´1.
The present study suggests that the minimum dbh diameter of 7 cm specified in previous studies for harvested young stands without cutting subsidies is consistent [13]. From an environmental sustainability perspective, the diameter selection will also depend on fertility differences among stands [16][17][18]. Thus, the most efficient operational stand-selection rules are site-specific, although the present results show how the rules can be identified systematically from the harvester's stand data. For example, the tree species couldn't be spruce in energy-wood cutting of this study [18]. This rule considered the cutting practice on this specific site, which had positive effects on forest diversity [16][17][18]. The stand-selection rules that we evaluated could be expanded to consider other important environmental criteria. For instance, a recent short-term study showed that 25 years after cutting of similar young stands, soil carbon (C) stocks after energy-wood cutting did not differ from those after conventional cutting without energy-wood cutting [13]. Soil C stocks after cutting could be considered as an additional environmental parameter in our stand-selection system (Figure 1). In addition, it would be useful to consider the magnitude of the variation of tree characteristics and site-specific stand characteristics, such as the fertility of the soil in long-term criteria to avoid potential risks (a decrease in soil C or regeneration) in some stands. These measures could assist the managers in identifying sites where young stands can be harvested with respecting long-term sustainability.
The design of sustainable energy system including local resource management system is described in a study of Vainio et al. [6]. Their system (GIS) analyzes the biomass flows from a forest level to plants. The system is used for strategic energy planning in the Finnish wood-procurement environment. The present study developed an operational energy planning system that can be used in the same wood procurement environment. The system automatically calculates and presents the operational procurement guidelines of sustainable cutting of young stands. If compared to regional resource management systems [4,5,7,9], the developed system can only be used for consideration of stem data of few stands. Therefore, large regional analyses are impossible to perform using current computers. However, the operational interactions between the parameters of energy supply chains, such as stem diameter, log length, energy content, wood moisture content, dry matter, solid and bulk density and infrastructure constraints could be evaluated more properly in order to deliver local and distributed forest biomass resources productively, cheaply and efficiently. Several recent studies have considered these parameters as important characteristics of management and control in solution of the complex energy planning problem. [3,7,8,10,11].
Although stands can be selected efficiently by the present computers and management information systems, the selection of trees by the harvester operators during the cutting phase is a potential source of error. In the modelling process of the present study, we assumed that the stand-selection rules could be followed more closely than they would be in real life. The tree selection problems can be overcome by developing better working rules for the harvester operators and by training them to better understand the sustainability criteria. Cutting could even be based on semi-automatic selection of each tree (e.g., using an onboard GPS receiver and machine vision) by the harvester [24,26,[30][31][32][33]. The harvester could suggest to the operator trees that meet both economic and environmental sustainability criteria. Furthermore, the computational method can be easily embedded in the software responsible for post-processing of harvester data, as well as in the stand-harvesting planning systems of the wood-procurement organization [22]. The results presented in this study will be applicable in Finland, but the models are likely to be generally applicable in harvesters and cutting simulators that have been developed for other countries, though possibly requiring some calibration.

Conclusions
The objective of the study was to present the operational procurement guidelines for sustainable cutting of young stands, utilizing the stem-selection rules in the harvester simulator's system.
The results were used to develop an automatic stand-selection system that would meet both economic and environmental sustainability criteria to guide energy-wood cutting work. Our results showed that there is a significant potential for meeting both criteria during energy-wood cutting, and that the system developed in this study would be more efficient than previous stand-selection systems. The new stand-selection system is superior to the manual system because it can automatically adapt the stem-selection rules to account for site-specific stem parameters. This would also be an important advantage for future mechanized energy-wood cutting operations, when tree selection is automated or semiautomated in harvesters. Predicting the sustainable yield of energy wood biomass has great potential, but the stand-selection variations remain large, even with the improved system, and managers must decide whether the degree of accuracy provided by the integrated system will be acceptable under their energy-wood cutting conditions.