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.

Temporomandibular disorders and related psychosocial factors in non-patients : A survey and a clinical follow-up study based on the RDC/TMD

Temporomandibular disorders (TMD) and psychosocial factors reportedly associate. The underlying factors remain partially obscure, however, and further studies are required to clarify the relationships. The aims of this study were thus to assess in a non-patient working population the prevalence of TMD and related symptoms, and to clinically diagnose and follow the natural courses of TMD over a one-year period. In addition, possible comorbidity of temporomandibular and/or neck muscle pain and perceived stress and their impact on work performance were investigated, as well as how various psychosocial aspects relate to TMD. A questionnaire was mailed to all 30- to 55-year-old employees of the Finnish Broadcasting Company Ltd. whose employment in the Helsinki area had lasted at least five years (n = 1784). Of the 1339 subjects, who returned the questionnaire, 241 were examined according to the RDC/TMD and standard neck muscle palpation methods. Clinical signs of temporomandibular and/or neck muscle pain were found in 118 subjects. One-year follow-up TMD examinations were conducted on 211 subjects. The prevalence of frequent painless TMJ-related symptoms was 10%, orofacial pain 7%, neck pain 38%, and headache 15%. TMD diagnoses were: myofascial pain (13%), disc displacements (16%), and arthralgia, osteoarthritis, osteoarthrosis (4%). Chronic myofascial pain was present in 7% and chronic disc displacement with reduction in 11% of the subjects. Symptoms were significantly associated with almost all the studied psychosocial symptoms. Reduced work performance was significantly positively associated with continuous pain, severity of pain, and health stress perception, and according to logistic regression, somatization with the probability of having chronic myofascial pain. It could be concluded based on the results of this study among a non-patient working population that TMD and related symptoms are common and associated with psychosocial factors. Moreover, myofascial pain and disc displacement with reduction are the most common diagnoses of TMD. In addition, self-reported health related stress, and continuous pain in temporomandibular and/or neck muscles are associated with reduced work performance, and somatization is significantly associated with chronic myofascial pain.

Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Tropical dryland agroforestry on clay soils: : Analysis of systems based on Acacia senegal in the Blue Nile region, Sudan

Acacia senegal, the gum arabic producing tree, is the most important component in traditional dryland agroforestry systems in the Blue Nile region, Sudan. The aim of the present study was to provide new knowledge on the potential use of A. senegal in dryland agroforestry systems on clay soils, as well as information on tree/crop interaction, and on silvicultural and management tools, with consideration on system productivity, nutrient cycling and sustainability. Moreover, the aim was also to clarify the intra-specific variation in the performance of A. senegal and, specifically, the adaptation of trees of different origin to the clay soils of the Blue Nile region. In agroforestry systems established at the beginning of the study, tree and crop growth, water use, gum and crop yields, nutrient cycling and system performance were investigated for a period of four years (1999 to 2002). Trees were grown at 5 x 5 m and 10 x 10 m spacing alone or in mixture with sorghum or sesame; crops were also grown in sole culture. The symbiotic biological N2 fixation by A. senegal was estimated using the 15N natural abundance (δ15N) procedure in eight provenances collected from different environments and soil types of the gum arabic belt and grown in clay soil in the Blue Nile region. Balanites aegyptiaca (a non-legume) was used as a non-N-fixing reference tree species, so as to allow 15N-based estimates of the proportion of the nitrogen in trees derived from the atmosphere. In the planted acacia trees, measurements were made on shoot growth, water-use efficiency (as assessed by the δ13C method) and (starting from the third year) gum production. Carbon isotope ratios were obtained from the leaves and branch wood samples. The agroforestry system design caused no statistically significant variation in water use, but the variation was highly significant between years, and the highest water use occurred in the years with high rainfall. No statistically significant differences were found in sorghum or sesame yields when intercropping and sole crop systems were compared (yield averages were 1.54 and 1.54 ha-1 for sorghum and 0.36 and 0.42 t ha-1 for sesame in the intercropped and mono-crop plots, respectively). Thus, at an early stage of agroforestry system management, A. senegal had no detrimental effect on crop yield, but the pattern of resource capture by trees and crops may change as the system matures. Intercropping resulted in taller trees and larger basal and crown diameters as compared to the development of sole trees. It also resulted in a higher land equivalent ratio. When gum yields were analysed it was found that a significant positive relationship existed between the second gum picking and the total gum yield. The second gum picking seems to be a decisive factor in gum production and could be used as an indicator for the total gum yield in a particular year. In trees, the concentrations of N and P were higher in leaves and roots, whereas the levels of K were higher in stems, branches and roots. Soil organic matter, N, P and K contents were highest in the upper soil stratum. There was some indication that the P content slightly increased in the topsoil as the agroforestry plantations aged. At a stocking of 400 trees ha-1 (5 x 5 m spacing), A. senegal accumulated in the biomass a total of 18, 1.21, 7.8 and 972 kg ha-1of N, P, K and OC, respectively. Trees contributed ca. 217 and 1500 kg ha-1 of K and OC, respectively, to the top 25-cm of soil over the first four years of intercropping. Acacia provenances of clay plain origin showed considerable variation in seed weight. They also had the lowest average seed weight as compared to the sandy soil (western) provenances. At the experimental site in the clay soil region, the clay provenances were distinctly superior to the sand provenances in all traits studied but especially in basal diameter and crown width, thus reflecting their adaptation to the environment. Values of δ13C, indicating water use efficiency, were higher in the sand soil group as compared to the clay one, both in leaves and in branch wood. This suggests that the sand provenances (with an average value of -28.07 ) displayed conservative water use and high drought tolerance. Of the clay provenances, the local one (Bout) displayed a highly negative (-29.31 ) value, which indicates less conservative water use that resulted in high productivity at this particular clay-soil site. Water use thus appeared to correspond to the environmental conditions prevailing at the original locations for these provenances. Results suggest that A. senegal provenances from the clay part of the gum belt are adapted for a faster growth rate and higher biomass and gum productivity as compared to provenances from sand regions. A strong negative relationship was found between the per-tree gum yield and water use efficiency, as indicated by δ13C. The differences in water use and gum production were greater among provenance groups than within them, suggesting that selection among rather than within provenances would result in distinct genetic gain in gum yield. The relative δ15N values ( ) were higher in B. aegyptiaca than in the N2-fixing acacia provenances. The amount of Ndfa increased significantly with age in all provenances, indicating that A. senegal is a potentially efficient nitrogen fixer and has an important role in t agroforestry development. The total above-ground contribution of fixed N to foliage growth in 4-year-old A. senegal trees was highest in the Rahad sand-soil provenance (46.7 kg N ha-1) and lowest in the Mazmoom clay-soil provenance (28.7 kg N ha-1). This study represents the first use of the δ15N method for estimating the N input by A. senegal in the gum belt of Sudan. Key words: Acacia senegal, agroforestry, clay plain, δ13C, δ15N, gum arabic, nutrient cycling, Ndfa, Sorghum bicolor, Sesamum indicum

Quality Assessment of a Forest Inventory Based on Airborne Laser Scanning

In this study, a quality assessment method based on sampling of primary laser inventory units (microsegments) was analysed. The accuracy of a laser inventory carried out in Kuhmo was analysed as a case study. Field sample plots were measured on the sampled microsegments in the Kuhmo inventory area. Two main questions were considered. Did the ALS based inventory meet the accuracy requirements set for the provider and how should a reliable, cost-efficient and independent quality assessment be undertaken. The agreement between control measurement and ALS based inventory was analysed in four ways: 1) The root mean squared errors (RMSEs) and bias were calculated. 2) Scatter plots with 95% confidence intervals were plotted and the placing of identity lines was checked. 3) Bland-Altman plots were drawn so that the mean difference of attributes between the control method and ALS-method was calculated and plotted against average value of attributes. 4) The tolerance limits were defined and combined with Bland-Altman plots. The RMSE values were compared to a reference study from which the accuracy requirements had been set to the service provider. The accuracy requirements in Kuhmo were achieved, however comparison of RMSE values proved to be difficult. Field control measurements are costly and time-consuming, but they are considered to be robust. However, control measurements might include errors, which are difficult to take into account. Using the Bland-Altman plots none of the compared methods are considered to be completely exact, so this offers a fair way to interpret results of assessment. The tolerance limits to be set on order combined with Bland-Altman plots were suggested to be taken in practise. In addition, bias should be calculated for total area. Some other approaches for quality control were briefly examined. No method was found to fulfil all the required demands of statistical reliability, cost-efficiency, time efficiency, simplicity and speed of implementation. Some benefits and shortcomings of the studied methods were discussed.

Influence of harvesting strategy on nutrient supply and production of dairy cows consuming diets based on grass and red clover silage

The main objective of this thesis was to elucidate the effects of regrowth grass silage and red clover silage on nutrient supply and milk production of dairy cows as compared with primary growth grass silages. In the first experiment (publication I), two primary growth and four regrowth grass silages were harvested at two stages of growth. These six silages were fed to 24 lactating dairy cows with two levels of concentrate allowance. Silage intake and energy corrected milk yield (ECM) responses, and the range in these response variables between the diets, were smaller when regrowth silages rather than primary growth silages were fed. Milk production of dairy cows reflected the intake of metabolizable energy (ME), and no differences in the ME utilization were found between the diets based on silages harvested from primary growth and regrowth. The ECM response to increased concentrate allowance was, on average, greater when regrowth rather than primary growth silages were fed. In the second experiment (publication II), two silages from primary growth and two from regrowth used in I were fed to rumen cannulated lactating dairy cows. Cows consumed less feed dry matter (DM), energy and protein, and produced less milk, when fed diets based on regrowth silages rather than primary growth silages. Lower milk production responses of regrowth grass silage diets were mainly due to the lower silage DM intake, and could not be accounted for by differences in energy or protein utilization. Regrowth grass silage intake was not limited due to neutral detergent fibre (NDF) digestion or rumen fill or passage kinetics. However, lower intake may be at least partly attributable to plant diseases such as leaf spot infections, dead deteriorating material or abundance of weeds, which are all higher in regrowth compared with primary growth, and increase with advancing regrowth. In the third experiment (publications III and IV), red clover silages and grass silages harvested at two stages of growth, and a mixed diet of red clover and grass silages, were fed to five rumen cannulated lactating dairy cows. In spite of the lower average ME intake for red clover diets, the ECM production remained unchanged suggesting more efficient utilisation of ME for red clover diets compared with grass diets. Intake of N, and omasal canal flows of total non-ammonia N (NAN), microbial and non-microbial NAN were higher for red clover than for grass silage diets, but were not affected by forage maturity. Delaying the harvest tended to decrease DM intake of grass silage and increase that of red clover silage. The digestion rate of potentially digestible NDF was faster for red clover diets than for grass silage diets. Delaying the harvest decreased the digestion rate for grass but increased it for red clover silage diets. The low intake of early-cut red clover silage could not be explained by silage digestibility, fermentation quality, or rumen fill but was most likely related to the nutritionally suboptimal diet composition because inclusion of moderate quality grass silage in mixed diet increased silage DM intake. Despite the higher total amino acid supply of cows fed red clover versus grass silage diets, further milk production responses on red clover diets were possibly compromised by an inadequate supply of methionine as evidenced by lower methionine concentration in the amino acid profile of omasal digesta and plasma. Increasing the maturity of ensiled red clover does not seem to affect silage DM intake as consistently as that of grasses. The efficiency of N utilization for milk protein synthesis was lower for red clover diets than for grass diets. It was negatively related to diet crude protein concentration similarly to grass silage diets.

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.

Topological Quantum Computation - An Analysis of an Anyon Model Based on Quantum Double Symmetries

There exists various suggestions for building a functional and a fault-tolerant large-scale quantum computer. Topological quantum computation is a more exotic suggestion, which makes use of the properties of quasiparticles manifest only in certain two-dimensional systems. These so called anyons exhibit topological degrees of freedom, which, in principle, can be used to execute quantum computation with intrinsic fault-tolerance. This feature is the main incentive to study topological quantum computation. The objective of this thesis is to provide an accessible introduction to the theory. In this thesis one has considered the theory of anyons arising in two-dimensional quantum mechanical systems, which are described by gauge theories based on so called quantum double symmetries. The quasiparticles are shown to exhibit interactions and carry quantum numbers, which are both of topological nature. Particularly, it is found that the addition of the quantum numbers is not unique, but that the fusion of the quasiparticles is described by a non-trivial fusion algebra. It is discussed how this property can be used to encode quantum information in a manner which is intrinsically protected from decoherence and how one could, in principle, perform quantum computation by braiding the quasiparticles. As an example of the presented general discussion, the particle spectrum and the fusion algebra of an anyon model based on the gauge group S_3 are explicitly derived. The fusion algebra is found to branch into multiple proper subalgebras and the simplest one of them is chosen as a model for an illustrative demonstration. The different steps of a topological quantum computation are outlined and the computational power of the model is assessed. It turns out that the chosen model is not universal for quantum computation. However, because the objective was a demonstration of the theory with explicit calculations, none of the other more complicated fusion subalgebras were considered. Studying their applicability for quantum computation could be a topic of further research.

Synthesis, Reactivity and Biological Activity of 17β-HSD1 Inhibitors Based on a Thieno[2,3-d]pryrimidin-4(3H)-one Core

Breast cancer is the most common cancer in women in Western countries. In the early stages of development most breast cancers are hormone-dependent, and estrogens, especially estradiol, have a pivotal role in their development and progression. One approach to the treatment of hormone-dependent breast cancers is to block the formation of the active estrogens by inhibiting the action of the steroid metabolising enzymes. 17beta-Hydroxysteroid dehydrogenase type 1 (17beta-HSD1) is a key enzyme in the biosynthesis of estradiol, the most potent female sex hormone. The 17beta-HSD1 enzyme catalyses the final step and converts estrone into the biologically active estradiol. Blocking 17beta-HSD1 activity with a specific enzyme inhibitor could provide a means to reduce circulating and tumour estradiol levels and thus promote tumour regression. In recent years 17beta-HSD1 has been recognised as an important drug target. Some inhibitors of 17beta-HSD1 have been reported, however, there are no inhibitors on the market nor have clinical trials been announced. The majority of known 17beta-HSD1 inhibitors are based on steroidal structures, while relatively little has been reported on non-steroidal inhibitors. As compared with 17beta-HSD1 inhibitors based on steroidal structures, non-steroidal compounds could have advantages of synthetic accessibility, drug-likeness, selectivity and non-estrogenicity. This study describes the synthesis of large group of novel 17beta-HSD1 inhibitors based on a non-steroidal thieno[2,3-d]pyrimidin-4(3H)-one core. An efficient synthesis route was developed for the lead compound and subsequently employed in the synthesis of thieno[2,3-d]pyrimidin-4(3H)-one based molecule library. The biological activities and binding of these inhibitors to 17beta-HSD1 and, finally, the quantitative structure activity relationship (QSAR) model are also reported. In this study, several potent and selective 17beta-HSD1 inhibitors without estrogenic activity were identified. This establishment of a novel class of inhibitors is a progressive achievement in 17beta-HSD1 inhibitor development. Furthermore, the 3D-QSAR model, constructed on the basis of this study, offers a powerful tool for future 17beta-HSD1 inhibitor development. As part of the fundamental science underpinning this research, the chemical reactivity of fused (di)cycloalkeno thieno[2,3-d]pyrimidin-4(3H)-ones with electrophilic reagents, i.e. Vilsmeier reagent and dimethylformamide dimethylacetal, was investigated. These findings resulted in a revision of the reaction mechanism of Vilsmeier haloformylation and further contributed to understanding the chemical reactivity of this compound class. This study revealed that the reactivity is dependent upon a stereoelectronic effect arising from different ring conformations.

Detection and progression of chronic allograft nephropathy : A study based on protocol biopsies

Background. Kidney transplantation (KTX) is considered to be the best treatment of terminal uremia. Despite improvements in short-term graft survival, a considerable number of kidney allografts are lost due to the premature death of patients with a functional kidney and to chronic allograft nephropathy (CAN). Aim. To investigate the risk factors involved in the progression of CAN and to analyze diagnostic methods for this entity. Materials and methods. Altogether, 153 implant and 364 protocol biopsies obtained between June 1996 and April 2008 were analyzed. The biopsies were classified according to Banff ’97 and chronic allograft damage index (CADI). Immunohistochemistry for TGF-β1 was performed in 49 biopsies. Kidney function was evaluated by creatinine and/or cystatin C measurement and by various estimates of glomerular filtration rate (GFR). Demographic data of the donors and recipients were recorded after 2 years’ follow-up. Results. Most of the 3-month biopsies (73%) were nearly normal. The mean CADI score in the 6-month biopsies decreased significantly after 2001. Diastolic hypertension correlated with ΔCADI. Serum creatinine concentration at hospital discharge and glomerulosclerosis were risk factors for ΔCADI. High total and LDL cholesterol, low HDL and hypertension correlated with chronic histological changes. The mean age of the donors increased from 41 -52 years. Older donors were more often women who had died from an underlying disease. The prevalence of delayed graft function increased over the years, while acute rejections (AR) decreased significantly over the years. Sub-clinical AR was observed in 4% and it did not affect long-term allograft function or CADI. Recipients´ drug treatment was modified along the Studies, being mycophenolate mophetil, tacrolimus, statins and blockers of the renine-angiotensin-system more frequently prescribed after 2001. Patients with a higher ΔCADI had lower GFR during follow-up. CADI over 2 was best predicted by creatinine, although with modest sensitivity and specificity. Neither cystatin C nor other estimates of GFR were superior to creatinine for CADI prediction. Cyclosporine A toxicity was seldom seen. Low cyclosporin A concentration after 2 h correlated with TGF- β1 expression in interstitial inflammatory cells, and this predicted worse graft function. Conclusions. The progression of CAN has been affected by two major factors: the donors’ characteristics and the recipients’ hypertension. The increased prevalence of DGF might be a consequence of the acceptance of older donors who had died from an underlying disease. Implant biopsies proved to be of prognostic value, and they are essential for comparison with subsequent biopsies. The progression of histological damage was associated with hypertension and dyslipidemia. The augmented expression of TGF-β1 in inflammatory cells is unclear, but it may be related to low immunosuppression. Serum creatinine is the most suitable tool for monitoring kidney allograft function on every-day basis. However, protocol biopsies at 6 and 12 months predicted late kidney allograft dysfunction and affected the clinical management of the patients. Protocol biopsies are thus a suitable surrogate to be used in clinical trials and for monitoring kidney allografts.