Adaptive tree multigrids and simplified spherical harmonics approximation in deterministic neutral and charged particle transport

A new deterministic three-dimensional neutral and charged particle transport code, MultiTrans, has been developed. In the novel approach, the adaptive tree multigrid technique is used in conjunction with simplified spherical harmonics approximation of the Boltzmann transport equation. The development of the new radiation transport code started in the framework of the Finnish boron neutron capture therapy (BNCT) project. Since the application of the MultiTrans code to BNCT dose planning problems, the testing and development of the MultiTrans code has continued in conventional radiotherapy and reactor physics applications. In this thesis, an overview of different numerical radiation transport methods is first given. Special features of the simplified spherical harmonics method and the adaptive tree multigrid technique are then reviewed. The usefulness of the new MultiTrans code has been indicated by verifying and validating the code performance for different types of neutral and charged particle transport problems, reported in separate publications.

Juvenile diversity of tree species in Theobroma cacao L. agroforestry and primary forest – Central Sulawesi, Indonesia

This thesis studies the tree species’ juvenile diversity in cacao (Theobroma cacao L.) based agroforestry and in primary forest in a natural conservation forest environment of Lore Lindu National Park, Sulawesi, Indonesia. Species’ adult composition in Lore Lindu National Park is relatively well studied, less is known about tree species’ diversity in seedling communities particularly in frequently disturbed cacao agroforestry field environment. Cacao production forms a potentially serious thread for maintaining the conservation areas pristine and forested in Sulawesi. The impacts of cacao production on natural environment are directly linked to the diversity and abundance of shade tree usage. The study aims at comparing differences between cacao agroforestry and natural forest in the surrounding area in their species composition in seedling and sapling size categories. The study was carried out in two parts. Biodiversity inventory of seedlings and saplings was combined with social survey with farmer interviews. Aim of the survey was to gain knowledge of the cacao fields, and farmers’ observations and choices regarding tree species associated with cacao. Data was collected in summer 2008. The assessment of the impact of environmental factors of solar radiation, weeding frequency, cacao tree planting density, distance to forest and distance to main park road, and type of habitat on seedling and sapling compositions was done with Non-metric Multidimensional Scaling (NMS). Outlier analysis was used to assess distorting variables for NMS, and Multi-Response Permutation Procedures (MRPP) analysis to differentiate the impact of categorical variables. Sampling success was estimated with rarefaction curves and jackknife estimate of species richness. In the inventory 135 species of trees and shrubs were found. Only some agroforestry related species were dominating. The most species rich were sapling communities in forest habitat. NMS was showing generally low linear correlation between variation of species composition and environmental variables. Solar radiation was having most significance as explaining variable. The most clearly separated in ordination were cacao and forest habitats. The results of seedling and sapling inventory were only partly coinciding with farmers’ knowledge of the tree species occurring on their fields. More research with frequent assessment of seedling cohorts is needed due to natural variability of cohorts and high mortality rate of seedlings.

Environmental and climatic dependences of stable isotope ratios in tree rings on different temporal scales

This work examines stable isotope ratios of carbon, oxygen and hydrogen in annual growth rings of trees. Isotopic composition in wood cellulose is used as a tool to study past climate. The method benefits from the accurate and precise dating provided by dendrochronology. In this study the origin, nature and the strength of climatic correlations are studied on different temporal scales and at different sites in Finland. The origin of carbon isotopic signal is in photosynthetic fractionation. The basic physical and chemical fractionations involved are reasonably well understood. This was confirmed by measuring instantaneous photosynthetic discrimination on Scots pine (Pinus sylvestris L.). The internal conductance of CO2 was recognized to have a significant impact on the observed fractionation, and further investigations are suggested to quantify its role in controlling the isotopic signal of photosynthates. Isotopic composition of the produced biomass can potentially be affected by variety of external factors that induce physiological changes in trees. Response of carbon isotopic signal in tree ring cellulose to changes in resource availability was assessed in a manipulation experiment. It showed that the signal was relatively stable despite of changes in water and nitrogen availability to the tree. Palaeoclimatic reconstructions are typically based on functions describing empirical relationship between isotopic and climatic parameters. These empirical relationships may change depending on the site conditions, species and timeframe studied. Annual variation in Scots pine tree ring carbon and oxygen isotopic composition was studied in northern and in central eastern Finland and annual variation in tree ring latewood carbon, oxygen and hydrogen isotopic ratio in Oak (Quercus robur L.) was studied in southern Finland. In all of the studied sites at least one of the studied isotope ratios was shown to record climate strongly enough to be used in climatic reconstructions. Using the observed relationships, four-century-long climate reconstructions from living Scots pine were created for northern and central eastern Finland. Also temporal stability of the relationships between three proxy indicators, tree ring growth and carbon and oxygen isotopic composition was studied during the four-hundred-year period. Isotope ratios measured from tree rings in Finland were shown to be sensitive indicators of climate. Increasing understanding of environmental controls and physiological mechanisms affecting tree ring isotopic composition will make possible more accurate interpretation of isotope data. This study also demonstrated that by measuring multiple isotopes and physical proxies from the same tree rings, additional information on tree physiology can be obtained. Thus isotopic ratios measured from tree ring cellulose provide means to improve the reliability of climate reconstructions.

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.

Management of bushland vegetation using rainwater harvesting in eastern Kenya.

Microcatchment water harvesting (MCWH) improved the survival and growth of planted trees on heavy soils in eastern Kenya five to six years after planting. In the best method, the cross-tied furrow microcatchments, the mean annual increments (MAI; based on the average biomass of living trees multiplied by tree density and survival) of the total and usable biomass in Prosopis juliflora were 2787 and 1610 kg ha-1 a-1 respectively, when the initial tree density was 500 to 1667 trees per hectare. Based on survival, the indigenous Acacia horrida, A. mellifera and A. zanzibarica were the most suitable species for planting using MCWH. When both survival and yield were considered, a local seed source of the introduced P. juliflora was superior to all other species. The MAI in MCWH was at best distinctly higher than that in the natural vegetation (163­307 and 66­111 kg ha-1 a-1 for total and usable biomass respectively); this cannot satisfy the fuelwood demand of concentrated populations, such as towns or irrigation schemes. The density of seeds of woody species in the topsoil was 40.1 seeds m-2 in the Acacia-Commiphora bushland and 12.6 seeds m-2 in the zone between the bushland and the Tana riverine forest. Rehabilitation of woody vegetation using the soil seed bank alone proved difficult due to the lack of seeds of desirable species. The regeneration and dynamics of woody vegetation were also studied both in cleared and undisturbed bushland. A sub-type of Acacia-Commiphora bushland was identified as Acacia reficiens bushland, in which the dominant Commiphora species is C. campestris. Most of the woody species did not have even-aged populations but cohort structures that were skewed towards young individuals. The woody vegetation and the status of soil nutrients were estimated to recover in 15­20 years on Vertic Natrargid soils after total removal of above-ground vegetation.

Alikasvoksen mittaus ja kartoitus laserkeilauksella

Tasaikäisen metsän alle muodostuvilla alikasvoksilla on merkitystä puunkorjuun, metsänuudistamisen, näkemä-ja maisema-analyysien sekä biodiversiteetin ja hiilitaseen arvioinnin kannalta. Ilma-aluksista tehtävä laserkeilaus on osoittautunut tehokkaaksi kaukokartoitusmenetelmäksi varttuneiden puustojen mittauksessa. Laserkeilauksen käyttöönotto operatiivisessa metsäsuunnittelussa mahdollistaa aiempaa tarkemman tiedon tuottamisen alikasvoksista, mikäli alikasvoksen ominaisuuksia voidaan tulkita laseraineistoista. Tässä työssä käytettiin tarkasti mitattuja maastokoealoja ja kaikulaserkeilausaineistoja (discrete return LiDAR) usealta vuodelta (1–2 km lentokorkeus, 0,9–9,7 pulssia m-2). Laserkeilausaineistot oli hankittu Optech ALTM3100 ja Leica ALS50-II sensoreilla. Koealat edustavat suomalaisia tasaikäisiä männiköitä eri kehitysvaiheissa. Tutkimuskysymykset olivat: 1) Minkälainen on alikasvoksesta saatu lasersignaali yksittäisen pulssin tasolla ja mitkä tekijät signaaliin vaikuttavat? 2) Mikä on käytännön sovelluksissa hyödynnettävien aluepohjaisten laserpiirteiden selitysvoima alikasvospuuston ominaisuuksien ennustamisessa? Erityisesti haluttiin selvittää, miten laserpulssin energiahäviöt ylempiin latvuskerroksiin vaikuttavat saatuun signaaliin, ja voidaanko laserkaikujen intensiteetille tehdä energiahäviöiden korjaus. Puulajien väliset erot laserkaiun intensiteetissä olivat pieniä ja vaihtelivat keilauksesta toiseen. Intensiteetin käyttömahdollisuudet alikasvoksen puulajin tulkinnassa ovat siten hyvin rajoittuneet. Energiahäviöt ylempiin latvuskerroksiin aiheuttivat alikasvoksesta saatuun lasersignaaliin kohinaa. Energiahäviöiden korjaus tehtiin alikasvoksesta saaduille laserpulssin 2. ja 3. kaiuille. Korjauksen avulla pystyttiin pienentämään kohteen sisäistä intensiteetin hajontaa ja parantamaan kohteiden luokittelutarkkuutta alikasvoskerroksessa. Käytettäessä 2. kaikuja oikeinluokitusprosentti luokituksessa maan ja yleisimmän puulajin välillä oli ennen korjausta 49,2–54,9 % ja korjauksen jälkeen 57,3–62,0 %. Vastaavat kappa-arvot olivat 0,03–0,13 ja 0,10–0,22. Tärkein energiahäviöitä selittävä tekijä oli pulssista saatujen aikaisempien kaikujen intensiteetti, mutta hieman merkitystä oli myös pulssin leikkausgeometrialla ylemmän latvuskerroksen puiden kanssa. Myös 3. kaiuilla luokitustarkkuus parani. Puulajien välillä havaittiin eroja siinä, kuinka herkästi ne tuottavat kaiun laserpulssin osuessa puuhun. Kuusi tuotti kaiun suuremmalla todennäköisyydellä kuin lehtipuut. Erityisen selvä tämä ero oli pulsseilla, joissa oli energiahäviöitä. Laserkaikujen korkeusjakaumapiirteet voivat siten olla riippuvaisia puulajista. Sensorien välillä havaittiin selviä eroja intensiteettijakaumissa, mikä vaikeuttaa eri sensoreilla hankittujen aineistojen yhdistämistä. Myös kaiun todennäköisyydet erosivat jonkin verran sensorien välillä, mikä aiheutti pieniä eroavaisuuksia kaikujen korkeusjakaumiin. Aluepohjaisista laserpiirteistä löydettiin alikasvoksen runkolukua ja keskipituutta hyvin selittäviä piirteitä, kun rajoitettiin tarkastelu yli 1 m pituisiin puihin. Piirteiden selitysvoima oli parempi runkoluvulle kuin keskipituudelle. Selitysvoima ei merkittävästi alentunut pulssitiheyden pienentyessä, mikä on hyvä asia käytännön sovelluksia ajatellen. Lehtipuun osuutta ei pystytty selittämään. Tulosten perusteella kaikulaserkeilausta voi olla mahdollista hyödyntää esimerkiksi ennakkoraivaustarpeen arvioinnissa. Sen sijaan alikasvoksen tarkempi luokittelu (esim. puulajitulkinta) voi olla vaikeaa. Kaikkein pienimpiä alikasvospuita ei pystytä havaitsemaan. Lisää tutkimuksia tarvitaan tulosten yleistämiseksi erilaisiin metsiköihin.

Metsän erirakenteisuutta edistävien hakkuiden kannattavuustarkastelu metsäaluetasolla

Tutkimuksessa vertailtiin metsän erirakenteisuutta edistävien poimintahakkuiden ja pienaukkohakkuiden kannattavuutta metsänhoitosuositusten mukaiseen metsänkasvatukseen Keski-Suomessa. Poimintahakkuut ja pienaukkohakkuut ovat menetelmiä, joilla voidaan lisätä luonnonmetsän häiriödynamiikan mukaista pienipiirteistä elinympäristöjen vaihtelua ja siksi ne sopivat etenkin erityiskohteisiin monimuotoisuuden, maiseman tai metsien monikäytön vuoksi. Ne johtavat yleensä vähitellen eri-ikäisrakenteiseen metsään, jossa puuston läpimittaluokkajakauma muistuttaa käänteistä J-kirjainta. Eri-ikäisrakenteisen metsänkäsittelyn taloudellista kannattavuutta puoltavat uudistumiskustannusten poisjäänti ja tukkipuihin painottuvat säännöllisin väliajoin toteutuvat hakkuut. Menetelmän soveltumista Suomen olosuhteisiin pidetään kuitenkin epävarmana. Tässä tutkimuksessa tarkasteltiin tasaikäisrakenteisen metsän muuttamista eri-ikäisrakenteiseksi 40 vuoden siirtymäaikana Metsähallituksen hallinnoimassa Isojäven ympäristöarvometsässä Kuhmoisissa. Tutkimusaineisto koostui 405 kuusivaltaisesta tasaikäisestä kuviosta, joiden pinta-alasta metsämaata on 636 hehtaaria. Metsän kehitystä simuloitiin puutason kasvumalleja käyttäen ja käsittelytoimenpiteet simuloitiin viisivuotiskausittain SIMO-metsäsuunnitteluohjelmistolla. Simulointien avulla selvitettiin jokaisen käsittelyskenaarion hakkuumäärät puutavaralajeittain, diskontatut kassavirrat ja puustopääoman muutos tarkasteluajanjakson aikana. Puunkorjuun yksikkökustannusten laskennan apuna käytettiin automatisoitua seurantajärjestelmää, jossa metsäkoneisiin asennettuilla matkapuhelimilla kerättiin MobiDoc2-sovelluksella metsäkoneiden käytöstä kiihtyvyystiedot, GPS-paikkatiedot ja syötetiedot. Lopulta jokaiselle käsittelyskenaariolle laskettiin metsän puuntuotannollista arvoa kuvaavalla tuottoarvon yhtälöllä nettonykyarvot, josta vähennettiin diskontatut puunkorjuun kustannukset. Tutkimuksen tulosten mukaan poimintahakkuun NPV oli 3 prosentin korkokannalla noin 91 % (7420 €/ha) ja pienaukkohakkuiden noin 99 % (8076 €/ha) metsänhoitosuositusten mukaisesta käsittelystä (8176 €/ha). Komparatiivinen statiikka osoitti, että korkokannan kasvattaminen 5 prosenttiin ei olennaisesti lisännyt nettonykyarvojen eroja. Poimintahakkuiden puunkorjuun yksikkökustannukset olivat 0,8 €/m3 harvennushakkuita pienemmät ja 7,2 €/m3 uudistushakkuita suuremmat. Pienaukkohakkuiden yksikkökustannukset olivat 0,7 €/m3 uudistushakkuita suuremmat.Tulosten perusteella on väistämätöntä että siirtymävaihe tasaikäisrakenteisesta eri-ikäisrakenteiseksi metsäksi aiheuttaa taloudellisia tappioita siitäkin huolimatta, että hakkuut ovat voimakkaita ja tehdään varttuneeseen kasvatusmetsään. Tappion määrä on metsän peitteisyyden ylläpidosta aiheutuva vaihtoehtoiskustannus.

The Biocontrol Agent Phlebiopsis gigantea: Efficacy and Impacts on the Stump Bacterial Biota and Conifer tree Defences

Phlebiopsis gigantea has been for a long time known as a strong competitor against Heterobasidion annosum and intensively applied as a biological control agent on stump surfaces of Picea abies in Fennoscandia. However, the mechanism underlying its antagonistic activity is still unknown. A primary concern is the possible impact of P. gigantea treatment on resident non-target microbial biota of conifer stumps. Additional risk factor is the potential of P. gigantea to acquire a necrotrophic habit through adaptation to living wood tissues. This study focused on the differential screening of several P. gigantea isolates from diverse geographical sources as well as the use of breeding approach to enhance the biocontrol efficacy against H. annosum infection. The results showed a significant positive correlation between growth rate in wood and high biocontrol efficacy. Furthermore, with aid of breeding approach, several progeny strains were obtained that had better growth rate and control efficacy than parental isolates. To address the issue of the potential of P. gigantea to acquire necrotrophic capability, a combination of histochemical, molecular and transcript profiling (454 sequencing) were used to investigate the interactions between these two fungi and ten year old P. sylvestris seedlings. The results revealed that both P. gigantea and H. annosum provoked strong necrotic lesions, but after prolonged incubation, P. gigantea lesions shrank and ceased to expand further. Tree seedlings pre-treated with P. gigantea further restricted H. annosum-induced necrosis and had elevated transcript levels of genes important for lignification, cell death regulation and jasmonic acid signalling. These suggest that induced localized resistance is a contributory factor for the biocontrol efficacy of P.gigantea, and it has a comparatively limited necrotrophic capability than H. annosum. Finally, to investigate the potential impact of P. gigantea on the stump bacterial biota, 16S rDNA isolated from tissue samples from stumps of P. abies after 1-, 6- and 13-year post treatment was sequenced using bar-coded 454 Titanium pyrosequencing. Proteobacteria were found to be the most abundant at the initial stages of stump decay but were selectively replaced by Acidobacteria at advanced stages of the decay. Moreover, P. gigantea treatment significantly decreased the bacterial richness at initial decay stage in the stumps. Over time, the bacterial community in the stumps gradually recovered and the negative effects of P. gigantea was attenuated.