Silvicultural decisions based on simulation-optimization systems

Forest management is facing new challenges under climate change. By adjusting thinning regimes, conventional forest management can be adapted to various objectives of utilization of forest resources, such as wood quality, forest bioenergy, and carbon sequestration. This thesis aims to develop and apply a simulation-optimization system as a tool for an interdisciplinary understanding of the interactions between wood science, forest ecology, and forest economics. In this thesis, the OptiFor software was developed for forest resources management. The OptiFor simulation-optimization system integrated the process-based growth model PipeQual, wood quality models, biomass production and carbon emission models, as well as energy wood and commercial logging models into a single optimization model. Osyczka s direct and random search algorithm was employed to identify optimal values for a set of decision variables. The numerical studies in this thesis broadened our current knowledge and understanding of the relationships between wood science, forest ecology, and forest economics. The results for timber production show that optimal thinning regimes depend on site quality and initial stand characteristics. Taking wood properties into account, our results show that increasing the intensity of thinning resulted in lower wood density and shorter fibers. The addition of nutrients accelerated volume growth, but lowered wood quality for Norway spruce. Integrating energy wood harvesting into conventional forest management showed that conventional forest management without energy wood harvesting was still superior in sparse stands of Scots pine. Energy wood from pre-commercial thinning turned out to be optimal for dense stands. When carbon balance is taken into account, our results show that changing carbon assessment methods leads to very different optimal thinning regimes and average carbon stocks. Raising the carbon price resulted in longer rotations and a higher mean annual increment, as well as a significantly higher average carbon stock over the rotation.

Energiapuuharvennus osana metsänkasvatusta ja ilmastonmuutoksen hillitsemistä

The aim of this study was to compare the differences between forest management incorporating energy wood thinning and forest management based on silvicultural recommendations (baseline). Energy wood thinning was substituted for young stand thinning and the first commercial thinning of industrial wood. The study was based on the forest stand data from Southern Finland, which were simulated by the MOTTI-simulator. The main interest was to find out the climatic benefits resulting from carbon sequestration and energy substitution. The value of energy wood was set to substitute it for coal as an alternative energy fuel (emission trade). Other political instruments (Kemera subsidies) were also analysed. The largest carbon dioxide emission reductions were achieved as a combination of carbon sequestration and energy substitution (on average, a 26-90 % increase in discounted present value in the beginning of rotation) compared to the baseline. Energy substitution increased emission reductions more effectively than carbon sequestration, when maintaining dense young stands. According to the study, energy wood thinning as a part of forest management was more profitable than the baseline when the value of carbon dioxide averaged more than 15 €/CO2 and other political subsidies were unchanged. Alternatively, the price of energy wood should on average exceed 21 €/m3 on the roadside in order to be profitable in the absence of political instruments. The most cost-efficient employment of energy wood thinning occured when the dominant height was 12 meters, when energy substitution was taken into account. According to alternative forest management, thinning of sapling stands could be done earlier or less intensely than thinning based on silvicultural recommendations and the present criteria of subsidies. Consequently, the first commercial thinning could be profitable to carry out either as harvesting of industrial wood or energy wood, or as integrated harvesting depending on the costs of the harvesting methods available and the price level of small-size industrial wood compared to energy wood.

Nitrogen dynamics of organic farming in a crop rotation based on red clover (Trifolium pratense) leys

In agricultural systems which rely on organic sources of nitrogen (N), of which the primary source is biological N fixation (BNF), it is extremely important to use N as efficiently as possible with minimal losses to the environment. The amount of N through BNF should be maximised and the availability of the residual N after legumes should be synchronised to the subsequent plant needs in the crop rotation. Six field experiments in three locations in Finland were conducted in 1994-2006 to determine the productivity and amount of BNF in red clover-grass leys of different ages. The residual effects of the leys for subsequent cereals as well as the N leaching risk were studied by field measurements and by simulation using the CoupModel. N use efficiency (NUE) and N balances were also calculated. The yields of red clover-grass leys were highest in the two-year-old leys (6 700 kg ha-1) under study, but the differences between 2- and 3-year old leys were not high in most cases. BNF (90 kg ha-1 in harvested biomass) correlated strongly with red clover dry matter yield, as the proportion of red clover N derived from the atmosphere (> 85%) was high in our conditions of organically farmed field with low soil mineral N. A red clover content of over 40% in dry matter is targeted to avoid negative N-balances and to gain N for the subsequent crop. Surprisingly, the leys had no significant effect on the yields and N uptake of the two subsequent cereals (winter rye or spring wheat, followed by spring oats). On the other hand, yield and C:N of leys, as well as BNF-N and total-N incorporated into the soil influenced on subsequent cereal yields. NUE of cereals from incorporated ley crop residues was rather high, varying from 30% to 80% (mean 48%). The mineral N content of soil in the profile of 0-90 cm was low, mainly 15-30 kg ha-1. Simulation of N dynamics by CoupModel functioned satisfactorily and is considered a useful tool to estimate N flows in cropping systems relying on organic N sources. Understanding the long-term influence of cultivation history and soil properties on N dynamics remains to be a challenge to further research.

Asteroid identification using statistical orbital inversion methods

An efficient and statistically robust solution for the identification of asteroids among numerous sets of astrometry is presented. In particular, numerical methods have been developed for the short-term identification of asteroids at discovery, and for the long-term identification of scarcely observed asteroids over apparitions, a task which has been lacking a robust method until now. The methods are based on the solid foundation of statistical orbital inversion properly taking into account the observational uncertainties, which allows for the detection of practically all correct identifications. Through the use of dimensionality-reduction techniques and efficient data structures, the exact methods have a loglinear, that is, O(nlog(n)), computational complexity, where n is the number of included observation sets. The methods developed are thus suitable for future large-scale surveys which anticipate a substantial increase in the astrometric data rate. Due to the discontinuous nature of asteroid astrometry, separate sets of astrometry must be linked to a common asteroid from the very first discovery detections onwards. The reason for the discontinuity in the observed positions is the rotation of the observer with the Earth as well as the motion of the asteroid and the observer about the Sun. Therefore, the aim of identification is to find a set of orbital elements that reproduce the observed positions with residuals similar to the inevitable observational uncertainty. Unless the astrometric observation sets are linked, the corresponding asteroid is eventually lost as the uncertainty of the predicted positions grows too large to allow successful follow-up. Whereas the presented identification theory and the numerical comparison algorithm are generally applicable, that is, also in fields other than astronomy (e.g., in the identification of space debris), the numerical methods developed for asteroid identification can immediately be applied to all objects on heliocentric orbits with negligible effects due to non-gravitational forces in the time frame of the analysis. The methods developed have been successfully applied to various identification problems. Simulations have shown that the methods developed are able to find virtually all correct linkages despite challenges such as numerous scarce observation sets, astrometric uncertainty, numerous objects confined to a limited region on the celestial sphere, long linking intervals, and substantial parallaxes. Tens of previously unknown main-belt asteroids have been identified with the short-term method in a preliminary study to locate asteroids among numerous unidentified sets of single-night astrometry of moving objects, and scarce astrometry obtained nearly simultaneously with Earth-based and space-based telescopes has been successfully linked despite a substantial parallax. Using the long-term method, thousands of realistic 3-linkages typically spanning several apparitions have so far been found among designated observation sets each spanning less than 48 hours.

Origins of the sex differences in handedness and mental rotation ability : genetic, environmental and hormonal effects

In humans, well-replicated and robust sex differences in cognitive functions exist for handedness and mental rotation ability. A common characteristic in human cognitive functions is the lateralization of language functions. Handedness is a common measure of laterality and is related to language lateralization. The prevalence of left-handedness is higher in males than in females, the male to female ratio being about 1.2. Among cognitive abilities, the largest sex difference is evident in the Vandenberg and Kuse Mental Rotation Test (MRT), which requires the ability to rotate objects in mental space. On average, males achieve scores one standard deviation higher than females in the MRT. The present thesis investigated the origins of the sex differences in laterality and spatial ability as represented by handedness and mental rotation ability, respectively. Two population-based Finnish twin cohorts were utilized in this study. Handedness was studied in 25 810 twins and 4068 singletons born before 1958 from the Older Finnish Twin Cohort, and in 4736 twins born in 1983-87 from the FinnTwin12. MRT was studied in a sub-sample of 804 young adult participants from the FinnTwin12 sample. The main findings of this study were: 1) the prevalence of left-handedness was higher among males than among females in both singletons and in twins; 2) males had significantly higher scores than females in MRT; 3) about one quarter of the variance in handedness and about half of the variance in MRT was explained by genetic effects, whereas the remainder of the variance in these traits was explained by environmental effects unique to each individual. The magnitude of the genetic effects was similar in both sexes; 4) left-handedness was significantly less common in female co-twins of a male than in female co-twins of a female, and female co-twins of a male scored significantly higher than did female co-twins of a female in the Mental Rotation Test. This dissertation discusses whether these differences between females from opposite- and same-sex twin pairs are due to the prenatal transfer of testosterone from the male fetus in females with male co-twins or whether they arise from postnatal socialization effects.

Fast distributed approximation algorithms for vertex cover and set cover in anonymous networks

We present a distributed algorithm that finds a maximal edge packing in O(Δ + log* W) synchronous communication rounds in a weighted graph, independent of the number of nodes in the network; here Δ is the maximum degree of the graph and W is the maximum weight. As a direct application, we have a distributed 2-approximation algorithm for minimum-weight vertex cover, with the same running time. We also show how to find an f-approximation of minimum-weight set cover in O(f2k2 + fk log* W) rounds; here k is the maximum size of a subset in the set cover instance, f is the maximum frequency of an element, and W is the maximum weight of a subset. The algorithms are deterministic, and they can be applied in anonymous networks.