Markov random fields and spatial smoothing in improving of forest inventory estimates

Markov random fields (MRF) are popular in image processing applications to describe spatial dependencies between image units. Here, we take a look at the theory and the models of MRFs with an application to improve forest inventory estimates. Typically, autocorrelation between study units is a nuisance in statistical inference, but we take an advantage of the dependencies to smooth noisy measurements by borrowing information from the neighbouring units. We build a stochastic spatial model, which we estimate with a Markov chain Monte Carlo simulation method. The smooth values are validated against another data set increasing our confidence that the estimates are more accurate than the originals.

Fast- and drift-ice communities in the Bothnian Bay and the impact of UVA radiation on the Baltic Sea ice ecology

The aim of this thesis was to study ecology of Baltic Sea ice from two perspectives. In the first two studies, sea-ice ecology from riverine-influenced fast ice to drift ice in the Bothnian Bay was investigated, whereas the last two studies focus on the sensitivity of sea-ice bacteria and algae to UVA examined in situ. The seasonal sea ice cover is one of the main characteristics of the Baltic Sea, and despite the brackish parental water, the ice structure is similar to polar ice with saline brine inclusions, the sea ice habitat. The decreasing seawater salinity from the northern Baltic Sea to the Bothnian Bay translates to decreasing brine volumes along the gradient, governing the size and community structure of the food webs in ice. However, the drift and fast ice in the Bothnian Bay may differ greatly in this sense, as drift ice may have been formed at more southern locations. Rafting and the formation of snow ice are common processes in the ice field of the Bothnian Bay. As evidenced in this thesis, rafting altered the vertical distribution of organisms and snow-ice formation provided habitable space in the better-illuminated, nitrogen-rich surface layer. The divergence between fast and drift ice became apparent at the more advanced stages, and chlorophyte biomass decreased from fast to drift ice, while the opposite held true for protozoan and metazoan biomass. The brine volumes affected the communities somewhat, and a higher percentage of flagellate species was generally linked to lower brine volumes, whereas chain-forming diatoms were mostly concentrated in layers with larger brine volumes. These results add to knowledge of the ecological significance of the ice cover lasting up to 7 months per year in this area. Sea-ice food webs are generally light-limited, but while increasing light irradiances typically enhance the primary production and further, the secondary production in sea ice, any increase in solar radiation also includes an increase in harmful UVA radiation. The Baltic Sea ice microbial communities were clearly sensitive to UVA and the responses were strongly linked to the earlier light history, as well as to the solar irradiances they were exposed to. The increased biomass of chlorophytes and pennate diatoms, when UVA was excluded, indicates that their normally minor contribution to the biomass in the upper layers of sea ice might be partly dictated by UVA. The effects of UVA on bacterial production in Baltic Sea ice mostly followed the responses in algal growth, but occasionally the exposure to UVA even enhanced the bacterial production. The dominant bacterial class, Flavobacteria, seemed to be UVA-tolerant, whereas all the Alpha-, Beta- and Gammaproteobacteria present in the surface layer showed UVA sensitivity. These results indicate that changes in the light field of ice may alter the community structure and affect the functioning of ice food webs, and are of importance when the effects of thinning of the ice cover are assessed.

CSR in forest industry – Case study of reporting and implementation of social responsibility in three international companies

Corporate Social Responsibility (CSR) has become increasingly important topic in forest industries, and other global companies, in recent years. Globalisation, faster information delivery and demand for sustainable development have set new challenges for global companies in their business operations. Also the importance of stakeholder relations, and pressure to become more transparent has increased in the forest industries. Three dimensions of corporate responsibility economic, environmental and social, are often included in the concept of CSR. Global companies mostly claim that these dimensions are equally important. This study analyses CSR in forest industry and has focus on reporting and implementation of social responsibility in three international companies. These case-companies are Stora Enso, SCA and Sappi, and they have different geographical base, product portfolios and therefore present interesting differences about forest industry strategy and CSR. Global Reporting Initiative (GRI) has created the most known and used reporting framework in CSR reporting. GRI Guidelines have made CSR reporting a uniform function, which can also be measured between companies and different sectors. GRI Guidelines have also made it possible to record and control CSR data in the companies. In recent years the use of GRI Guidelines has increased substantially. Typically CSR reporting on economic and environmental responsibility have been systematic in the global companies and often driven by legistlation and other regulations. However the social responsibility has been less regulated and more difficult to compare. Therefore it has previously been often less focused in the CSR reporting of the global companies. The implementation and use of GRI Guidelines have also increased dialogue on social responsibility issues and stakeholder management in global companies. This study analyses the use of GRI´s framework in the forest industry companies´ CSR reporting. This is a qualitative study and the disclosure of data is empricially analysed using content analysis. Content analysis has been selected as a method for this study because it makes it possible to use different sources of information. The data of this study consists of existing academic literature of CSR, sustainability reports of thecase-companies during 2005-2009, and the semi-structured interviews with company representatives. Different sources provide the possibility to look at specific subject from more than one viewpoint. The results of the study show that all case-companies have relatively common themes in their CSR disclosure, and the differences rise mainly from their product-portfolios, and geographic base. Social impacts to local communities, in the CSR of the companies, were mainly dominated by issues concerning creating wealth to the society and impacting communities through creation of work. The comparability of the CSR reporting, and especially social indicators increased significally from 2007 onwards in all case-companies. Even though the companies claim that three dimensions of CSR economic, environmental and social are equally important economic issues and profit improvement still seem to drive most of the operations in the global companies. Many issues that are covered by laws and regulations are still essentially presented as social responsibility in CSR. However often the unwelcome issues in companies like closing operations are covered just briefly, and without adequate explanation. To make social responsibility equally important in the CSR it would demand more emphasis from all the case-companies. A lot of emphasis should be put especially on the detail and extensiveness of the social reponsibility content in the CSR.

Effect of soil preparation of boreal spruce forest on air and soil temperature conditions in forest regeneration areas.

The effect of scarification, ploughing and cross-directional plouhing on temperature conditions in the soil and adjacent air layer have been studied during 11 consecutive growth periods by using an unprepared clear-cut area as a control site. The maximum and minimum temperatures were measured daily in the summer months, and other temperature observations were made at four-hour intervals by means of a Grant measuring instrument. The development of the seedling stand was also followed in order to determine its shading effect on the soil surface. Soil preparation decreased the daily temperature amplitude of the air at the height of 10 cm. The maximum temperatures on sunny days were lower in the tilts of the ploughed and in the humps of the cross-directional ploughed sites compared with the unprepared area. Correspondingly, the night temperatures were higher and so the soil preparation considerably reduced the risk of night frost. In the soil at the depth of 5 cm, soil preparation increased daytime temperatures and reduced night temperatures compared with unprepared area. The maximum increase in monthly mean temperatures was almost 5 °C, and the daily variation in the surface parts of the tilts and humps increased so that excessively high temperatures for the optimal growth of the root system were measured from time to time. The temperature also rose at the depths of 50 and 100 cm. Soil preparation also increased the cumulative temperature sum. The highest sums accumulated during the summer months were recorded at the depth of 5 cm in the humps of cross-directional ploughed area (1127 dd.) and in the tilts of the ploughed area (1106 dd.), while the corresponding figure in the unprepared soil was 718 dd. At the height of 10 cm the highest temperature sum was 1020 dd. in the hump, the corresponding figure in the unprepared area being 925 dd. The incidence of high temperature amplitudes and percentage of high temperatures at the depth of 5 cm decreased most rapidly in the humps of cross-directional ploughed area and in the ploughing tilts towards the end of the measurement period. The decrease was attributed principally to the compressing of tilts, the ground vegetation succession and the growth of seedlings. The mean summer temperature in the unprepared area was lower than in the prepared area and the difference did not diminish during the period studied. The increase in temperature brought about by soil preparation thus lasts at least more than 10 years.

Understorey vegetation as an indicator of forest site potential in Southern Finland.

The relationship between site characteristics and understorey vegetation composition was analysed with quantitative methods, especially from the viewpoint of site quality estimation. Theoretical models were applied to an empirical data set collected from the upland forests of southern Finland comprising 104 sites dominated by Scots pine (Pinus sylvestris L.), and 165 sites dominated by Norway spruce (Picea abies (L.) Karsten). Site index H100 was used as an independent measure of site quality. A new model for the estimation of site quality at sites with a known understorey vegetation composition was introduced. It is based on the application of Bayes' theorem to the density function of site quality within the study area combined with the species-specific presence-absence response curves. The resulting posterior probability density function may be used for calculating an estimate for the site variable. Using this method, a jackknife estimate of site index H100 was calculated separately for pine- and spruce-dominated sites. The results indicated that the cross-validation root mean squared error (RMSEcv) of the estimates improved from 2.98 m down to 2.34 m relative to the "null" model (standard deviation of the sample distribution) in pine-dominated forests. In spruce-dominated forests RMSEcv decreased from 3.94 m down to 3.16 m. In order to assess these results, four other estimation methods based on understorey vegetation composition were applied to the same data set. The results showed that none of the methods was clearly superior to the others. In pine-dominated forests, RMSEcv varied between 2.34 and 2.47 m, and the corresponding range for spruce-dominated forests was from 3.13 to 3.57 m.

Non-industrial private forest landowners’ choices of timber management strategies and potential allowable cut.

In the study, the potential allowable cut in the district of Pohjois-Savo - based on the non-industrial private forest landowners' (NIPF) choices of timber management strategies - was clarified. Alternative timber management strategies were generated, and the choices and factors affecting the choices of timber management strategies by NIPF landowners were studied. The choices of timber management strategies were solved by maximizing the utility functions of the NIPF landowners. The parameters of the utility functions were estimated using the Analytic Hierarchy Process (AHP). The level of the potential allowable cut was compared to the cutting budgets based on the 7th and 8th National Forest Inventories (NFI7 and NFI8), to the combining of private forestry plans, and to the realized drain from non-industrial private forests. The potential allowable cut was calculated using the same MELA system as has been used in the calculation of the national cutting budget. The data consisted of the NIPF holdings (from the TASO planning system) that had been inventoried compartmentwise and had forestry plans made during the years 1984-1992. The NIPF landowners' choices of timber management strategies were clarified by a two-phase mail inquiry. The most preferred strategy obtained was "sustainability" (chosen by 62 % of landowners). The second in order of preference was "finance" (17 %) and the third was "saving" (11 %). "No cuttings", and "maximum cuttings" were the least preferred (9 % and 1 %, resp.). The factors promoting the choices of strategies with intensive cuttings were a) "farmer as forest owner" and "owning fields", b) "increase in the size of the forest holding", c) agriculture and forestry orientation in production, d) "decreasing short term stumpage earning expectations", e) "increasing intensity of future cuttings", and f) "choice of forest taxation system based on site productivity". The potential allowable cut defined in the study was 20 % higher than the average of the realized drain during the years 1988-1993, which in turn, was at the same level as the cutting budget based on the combining of forestry plans in eastern Finland. Respectively, the potential allowable cut defined in the study was 12 % lower than the NFI8-based greatest sustained allowable cut for the 1990s. Using the method presented in this study, timber management strategies can be clarified for non-industrial private forest landowners in different parts of Finland. Based on the choices of timber managemet strategies, regular cutting budgets can be calculated more realistically than before.

Non-industrial private forest landowners’ choices of timber management strategies.

The factors affecting the non-industrial, private forest landowners' (hereafter referred to using the acronym NIPF) strategic decisions in management planning are studied. A genetic algorithm is used to induce a set of rules predicting potential cut of the landowners' choices of preferred timber management strategies. The rules are based on variables describing the characteristics of the landowners and their forest holdings. The predictive ability of a genetic algorithm is compared to linear regression analysis using identical data sets. The data are cross-validated seven times applying both genetic algorithm and regression analyses in order to examine the data-sensitivity and robustness of the generated models. The optimal rule set derived from genetic algorithm analyses included the following variables: mean initial volume, landowner's positive price expectations for the next eight years, landowner being classified as farmer, and preference for the recreational use of forest property. When tested with previously unseen test data, the optimal rule set resulted in a relative root mean square error of 0.40. In the regression analyses, the optimal regression equation consisted of the following variables: mean initial volume, proportion of forestry income, intention to cut extensively in future, and positive price expectations for the next two years. The R2 of the optimal regression equation was 0.34 and the relative root mean square error obtained from the test data was 0.38. In both models, mean initial volume and positive stumpage price expectations were entered as significant predictors of potential cut of preferred timber management strategy. When tested with the complete data set of 201 observations, both the optimal rule set and the optimal regression model achieved the same level of accuracy.

A participatory approach to tactical forest planning.

The paper examines the needs, premises and criteria for effective public participation in tactical forest planning. A method for participatory forest planning utilizing the techniques of preference analysis, professional expertise and heuristic optimization is introduced. The techniques do not cover the whole process of participatory planning, but are applied as a tool constituting the numerical core for decision support. The complexity of multi-resource management is addressed by hierarchical decision analysis which assesses the public values, preferences and decision criteria toward the planning situation. An optimal management plan is sought using heuristic optimization. The plan can further be improved through mutual negotiations, if necessary. The use of the approach is demonstrated with an illustrative example, it's merits and challenges for participatory forest planning and decision making are discussed and a model for applying it in general forest planning context is depicted. By using the approach, valuable information can be obtained about public preferences and the effects of taking them into consideration on the choice of the combination of standwise treatment proposals for a forest area. Participatory forest planning calculations, carried out by the approach presented in the paper, can be utilized in conflict management and in developing compromises between competing interests.

Factors affecting partitioning of private forest holdings in Finland.

Questions of the small size of non-industrial private forest (NIPF) holdings in Finland are considered and factors affecting their partitioning are analyzed. This work arises out of Finnish forest policy statements in which the small average size of holdings has been seen to have a negative influence on the economics of forestry. A survey of the literature indicates that the size of holdings is an important factor determining the costs of logging and silvicultural operations, while its influence on the timber supply is slight. The empirical data are based on a sample of 314 holdings collected by interviewing forest owners in the years 1980-86. In 1990-91 the same holdings were resurveyed by means of a postal inquiry and partly by interviewing forest owners. The principal objective in compiling the data is to assist in quantifying ownership factors that influence partitioning among different kinds of NIPF holdings. Thus the mechanism of partitioning were described and a maximum likelihood logistic regression model was constructed using seven independent holding and ownership variables. One out of four holdings had undergone partitioning in conjunction with a change in ownership, one fifth among family owned holdings and nearly a half among jointly owned holdings. The results of the logistic regression model indicate, for instance, that the odds on partitioning is about three times greater for jointly owned holdings than for family owned ones. Also, the probabilities of partitioning were estimated and the impact of independent dichotomous variables on the probability of partitioning ranged between 0.02 and 0.10. The low value of the Hosmer-Lemeshow test statistic indicates a good fit of the model and the rate of correct classification was estimated to be 88 per cent with a cutoff point of 0.5. The average size of holdings undergoing ownership changes decreased from 29.9 ha to 28.7 ha over the approximate interval 1983-90. In addition, the transition probability matrix showed that the trends towards smaller size categories mostly involved in the small size categories, less than 20 ha. The results of the study can be used in considering the effects of the small size of holdings for forestry and if the purpose is to influence partitioning through forest or rural policy.

Analysis of forest soil chemistry and hydrology with a dynamic model ACIDIC.

In this study we analyze how the ion concentrations in forest soil solution are determined by hydrological and biogeochemical processes. A dynamic model ACIDIC was developed, including processes common to dynamic soil acidification models. The model treats up to eight interacting layers and simulates soil hydrology, transpiration, root water and nutrient uptake, cation exchange, dissolution and reactions of Al hydroxides in solution, and the formation of carbonic acid and its dissociation products. It includes also a possibility to a simultaneous use of preferential and matrix flow paths, enabling the throughfall water to enter the deeper soil layers in macropores without first reacting with the upper layers. Three different combinations of routing the throughfall water via macro- and micropores through the soil profile is presented. The large vertical gradient in the observed total charge was simulated succesfully. According to the simulations, gradient is mostly caused by differences in the intensity of water uptake, sulfate adsorption and organic anion retention at the various depths. The temporal variations in Ca and Mg concentrations were simulated fairly well in all soil layers. For H+, Al and K there were much more variation in the observed than in the simulated concentrations. Flow in macropores is a possible explanation for the apparent disequilibrium of the cation exchange for H+ and K, as the solution H+ and K concentrations have great vertical gradients in soil. The amount of exchangeable H+ increased in the O and E horizons and decreased in the Bs1 and Bs2 horizons, the net change in whole soil profile being a decrease. A large part of the decrease of the exchangeable H+ in the illuvial B horizon was caused by sulfate adsorption. The model produces soil water amounts and solution ion concentrations which are comparable to the measured values, and it can be used in both hydrological and chemical studies of soils.

Differences in the climatic adaptation of silver birch (Betula pendula) and downy birch (B. pubescens) in Finland based on male flowering phenology.

Male flowering was studied at the canopy level in 10 silver birch (Betula pendula Roth) stands from 8 localities and in 14 downy birch (B. pubescens Ehrh.) stands from 10 localities in Finland from 1963 to 1973. Distributions of cumulative pollen catches were compared to the normal Gaussian distribution. The basis for the timing of flowering was the 50 per cent point of the anthesis-fitted normal distribution. To eliminate effects of background pollen, only the central, normally distributed part of the cumulative distribution was used. Development up to the median point of the distribution was measured and tested in calendar days, in degree days (> 5 °C) and in period units. The count of each parameter began on and included March 19. Male flowering in silver birch occurred from late April to late June depending on latitude, and flowering in downy birch took place from early May to early July. The heat sums needed for male flowering varied in downy birch stands latitudinally but there was practically no latitudinal variation in heat sums needed for silver birch flowering. The amount of male flowering in stands of both birch species were found to have a large annual variation but without any clear periodicity. The between years pollen catch variation in stands of either birch species did not show any significant latitudinal correlation in contrast to Norway spruce stands. The period unit heat sum gave the most accurate forecast of the timing of flowering for 60 per cent of the silver birch stands and for 78.6 per cent of the for downy birch stands. Calendar days, however, gave the best forecast for silver birch in 25 per cent of the cases, while degree days gave the best forecast for downy birch in 21.4 per cent of the cases. Silver birch seems to have a local inclination for a more fixed flowering date compared to downy birch, which could mean a considerable photoperiodic influence on flowering time of silver birch. Silver birch and downy birch had different geographical correlations. Frequent hybridization of birch species occurs more often in northern Finland in than in more southern latitudes. The different timing in flowering caused increasing scatter in flowering times in the north, especially in the case of downy birch. The chance of simultaneous flowering of silver birch and downy birch so increased northwards due to a more variable climate and also higher altitudinal variations. Compared with conifers, the reproduction cycles of both birch species were found to be well protected from damage by frost.

Bioremediation of diesel oil contaminated soil and water

Diesel spills contaminate aquatic and terrestrial environments. To prevent the environmental and health risks, the remediation needs to be advanced. Bioremediation, i.e., degradation by microbes, is one of the suitable methods for cleaning diesel contamination. In monitored natural attenuation technique are natural processes in situ combined, including bioremediation, volatilization, sorption, dilution and dispersion. Soil bacteria are capable of adapting to degrade environmental pollutants, but in addition, some soil types may have indigenous bacteria that are naturally suitable for degradation. The objectives for this work were (1) to find a feasible and economical technique to remediate oil spilled into Baltic Sea water and (2) to bioremediate soil contaminated by diesel oil. Moreover, the aim was (3) to study the potential for natural attenuation and the indigenous bacteria in soil, and possible adaptation to degrade diesel hydrocarbons. In the aquatic environment, the study concentrated on diesel oil sorption to cotton grass fiber, a natural by-product of peat harvesting. The impact of diesel pollution was followed in bacteria, phytoplankton and mussels. In a terrestrial environment, the focus was to compare the methods of enhanced biodegradation (biostimulation and bioaugmentation), and to study natural attenuation of oil hydrocarbons in different soil types and the effect that a history of previous contamination may have on the bioremediation potential. (1) In the aquatic environment, rapid removal of diesel oil was significant for survival of tested species and thereby diversity maintained. Cotton grass not only absorbed the diesel but also benefited the bacterial growth by providing a large colonizable surface area and hence oil-microbe contact area. Therefore use of this method would enhance bioremediation of diesel spills. (2) Biostimulation enhances bioremediation, and (3) indigenous diesel-degrading bacteria are present in boreal environments, so microbial inocula are not always needed. In the terrestrial environment experiments, the combination of aeration and addition of slowly released nitrogen advanced the oil hydrocarbon degradation. Previous contamination of soil gives the bacterial community the potential for rapid adaptation and efficient degradation of the same type of contaminant. When the freshly contaminated site needs addition of diesel degraders, previously contaminated and remediated soil could be used as a bacterial inoculum. Another choice of inoculum could be conifer forest soil, which provides a plentiful population of degraders, and based on the present results, could be considered as a safe non-polluted inoculum. According to the findings in this thesis, bioremediation (microbial degradation) and monitored natural attenuation (microbial, physical and chemical degradation) are both suitable techniques for remediation of diesel-contaminated sites in Finland.

Silva Fennica. A quarterly journal of forest science. Contents Vol. 28, 1994.

Julkaistu Silva Fennica Vol. 28(4) -numeron liitteenä.