Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

Understanding solid-state transformations during dehydration: new insights using vibrational spectroscopy and multivariate modelling

Many active pharmaceutical ingredients (APIs) have both anhydrate and hydrate forms. Due to the different physicochemical properties of solid forms, the changes in solid-state may result in therapeutic, pharmaceutical, legal and commercial problems. In order to obtain good solid dosage form quality and performance, there is a constant need to understand and control these phase transitions during manufacturing and storage. Thus it is important to detect and also quantify the possible transitions between the different forms. In recent years, vibrational spectroscopy has become an increasingly popular tool to characterise the solid-state forms and their phase transitions. It offers several advantages over other characterisation techniques including an ability to obtain molecular level information, minimal sample preparation, and the possibility of monitoring changes non-destructively in-line. Dehydration is the phase transition of hydrates which is frequently encountered during the dosage form production and storage. The aim of the present thesis was to investigate the dehydration behaviour of diverse pharmaceutical hydrates by near infrared (NIR), Raman and terahertz pulsed spectroscopic (TPS) monitoring together with multivariate data analysis. The goal was to reveal new perspectives for investigation of the dehydration at the molecular level. Solid-state transformations were monitored during dehydration of diverse hydrates on hot-stage. The results obtained from qualitative experiments were used to develop a method and perform the quantification of the solid-state forms during process induced dehydration in a fluidised bed dryer. Both in situ and in-line process monitoring and quantification was performed. This thesis demonstrated the utility of vibrational spectroscopy techniques and multivariate modelling to monitor and investigate dehydration behaviour in situ and during fluidised bed drying. All three spectroscopic methods proved complementary in the study of dehydration. NIR spectroscopy models could quantify the solid-state forms in the binary system, but were unable to quantify all the forms in the quaternary system. Raman spectroscopy models on the other hand could quantify all four solid-state forms that appeared upon isothermal dehydration. The speed of spectroscopic methods makes them applicable for monitoring dehydration and the quantification of multiple forms was performed during phase transition. Thus the solid-state structure information at the molecular level was directly obtained. TPS detected the intermolecular phonon modes and Raman spectroscopy detected mostly the changes in intramolecular vibrations. Both techniques revealed information about the crystal structure changes. NIR spectroscopy, on the other hand was more sensitive to water content and hydrogen bonding environment of water molecules. This study provides a basis for real time process monitoring using vibrational spectroscopy during pharmaceutical manufacturing.

Landscape-scale soil erosion modelling and ecological restoration for a mountainous watershed in Sichuan, China

Sichuanissa Tiibetin ylängön metsäkato on pysähtynyt mutta eroosio-ongelmat jatkuvat Viikin tropiikki-instituutin tutkija Ping ZHOU kartoitti trooppisen metsänhoidon alaan kuuluvassa väitöskirjatyössään maaperän eroosioalttiutta ja sen riippuvuutta metsäkasvillisuudesta Jangtsen tärkeää sivuhaaraa Min-jokea ympäröivällä n. 7400 neliökilometrin suuruisella valuma-alueella Sichuanin Aba-piirikunnassa. Aineistonaan hän käytti muun muassa satelliittikartoitustietoja ja mittaustuloksia yli 600 maastokoealalta. Tutkimuksen nimi suomeksi on "Maaperän eroosion mallinnus ja vuoristoisen valuma-alueen ekologinen ennallistaminen Sichuanissa Kiinassa". Aikaisempien tutkimusten perusteella oli tiedossa että metsien häviäminen tällä alueella pysähtyi jo 1980-luvun alussa. Sen jälkeen on metsien pinta-ala hitaasti kasvanut etupäässä sen vuoksi, että teollinen puunhakkuu luonnonmetsissä kiellettiin kokonaan v. 1998 ja 25 astetta jyrkemmillä rinteillä myös maatalouden harjoittaminen on saatu lopetetuksi viljelijöille tarjottujen taloudellisten houkuttimien avulla. Täten myös pelto- ja laidunmaata on voitu ennallistaa metsäksi. Ping Zhou pystyi jakamaan 5700 metrin korkeuteen saakka kohoavan vuoristoalueen eroosioalttiudeltaan erilaisiin vyöhykkeisiin rinteen kaltevuuden, sademäärän, kasvipeitteen ja maalajin perusteella. Noin 15 prosentilla tutkitun valuma-alueen pinta-alasta, lähinnä Min-joen pääuomaa ympäröivillä jyrkillä rinteillä, eroosioriski oli suuri tai erittäin suuri. Eri tyyppisellä kasvillisuudella oli hyvin erilainen vaikutus eroosioalttiuteen, ja myös alueen sijainti vuoriston eri korkeuksilla vaikutti eroosioon. Säästyneet lähes luonnontilaiset havumetsät, joita on etupäässä vuoriston ylimmissä osissa 2600-4000 metrin korkeudella, edistävät tehokkaasti metsän luontaista uudistumista ja levittäytymistä vaurioituneille alueille. Säilyneiden metsien puulajikoostumus antoi tutkimuksessa mahdollisuuden ennustaa metsien tulevaa kehitystä koko tutkitulla valuma-alueella sen eri korkeusvyöhykkeissä ja eri maaperätyypeillä. Ennallistamisen kannalta ongelmallisimpia olivat alueet joilta metsäpeite oli lähinnä puiden teollisen hakkuun vuoksi kokonaan hävinnyt ja joilla maaperä yleisesti oli eroosion pahoin kuluttama. Näillä alueilla ei ole tehty juuri mitään uudistamis- tai ennallistamistoimenpiteitä. Niillä metsien ennallistaminen vaatii myös puiden tai pensaiden istuttamista. Tähän sopivia ovat erityisesti ilmakehän typpeä sitovat lajit, joista alueella kasvaa luontaisena mm. sama tyrnilaji joka esiintyy myös Suomessa. Työssä tutkittiin yli kahdeksankymmenen paikallisen luontaisen puulajin (joista peräti noin kolmannes on havupuulajeja) ekologisia ominaisuuksia ja soveltuvuutta metsien ennallistamiseen. Avainasemassa työn onnistumisen kannalta ovat nyt paikalliset asukkaat, joiden maankäytön muutokset ovat jo selvästi edistänet luonnonmetsän ennalleen palautumista. Suomen Akatemia rahoitti vuosina 2004-2006 VITRI:n tutkimushanketta, josta Ping Zhou'n väitöskirjatyö muodosti keskeisen osan. Kenttätyö Sichuanissa avasi mahdollisuuden hedelmälliseen monitieteiseen yhteistyöhön ja tutkijavaihtoon Kiinan tiedeakatemian alaisen Chengdun biologiainstituutin (CIB) kanssa; tämä tieteellinen kanssakäyminen jatkuu edelleen.

Soil carbon modelling as a tool for carbon balance studies in forestry

Soils represent a remarkable stock of carbon, and forest soils are estimated to hold half of the global stock of soil carbon. Topical concern about the effects of climate change and forest management on soil carbon as well as practical reporting requirements set by climate conventions have created a need to assess soil carbon stock changes reliably and transparently. The large spatial variability of soil carbon commensurate with relatively slow changes in stocks hinders the assessment of soil carbon stocks and their changes by direct measurements. Models therefore widely serve to estimate carbon stocks and stock changes in soils. This dissertation aimed to develop the soil carbon model YASSO for upland forest soils. The model was aimed to take into account the most important processes controlling the decomposition in soils, yet remain simple enough to ensure its practical applicability in different applications. The model structure and assumptions were presented and the model parameters were defined with empirical measurements. The model was evaluated by studying the sensitivities of the model results to parameter values, by estimating the precision of the results with an uncertainty analysis, and by assessing the accuracy of the model by comparing the predictions against measured data and to the results of an alternative model. The model was applied to study the effects of intensified biomass extraction on the forest carbon balance and to estimate the effects of soil carbon deficit on net greenhouse gas emissions of energy use of forest residues. The model was also applied in an inventory based method to assess the national scale forest carbon balance for Finland’s forests from 1922 to 2004. YASSO managed to describe sufficiently the effects of both the variable litter and climatic conditions on decomposition. When combined with the stand models or other systems providing litter information, the dynamic approach of the model proved to be powerful for estimating changes in soil carbon stocks on different scales. The climate dependency of the model, the effects of nitrogen on decomposition and forest growth as well as the effects of soil texture on soil carbon stock dynamics are areas for development when considering the applicability of the model to different research questions, different land use types and wider geographic regions. Intensified biomass extraction affects soil carbon stocks, and these changes in stocks should be taken into account when considering the net effects of forest residue utilisation as energy. On a national scale, soil carbon stocks play an important role in forest carbon balances.

Multiplication of hybrid aspen (Populus tremula L. x P. tremuloides Michx.) from cuttings

The primary aim of the present study was to find an efficient and simple method of vegetative propagation for producing large numbers of hybrid aspen (Populus tremuloides L. x P. tremula Michx.) plants for forest plantations. The key objectives were to investigate the main physiological factors that affect the ability of cuttings to regenerate and to determine whether these factors could be manipulated by different growth conditions. In addition, clonal variation in traits related to propagation success was examined. According to our results, with the stem cutting method, depending on the clone, it is possible to obtain only 1−8 plants from one stock plant per year. With the root cutting method the corresponding values for two-year-old stock plants are 81−207 plants. The difference in number of cuttings between one- and two-year-old stock plants is so pronounced that it is economically feasible to grow stock plants for two years. There is no reason to use much older stock plants as a source of cuttings, as it has been observed that rooting ability diminishes as root diameter increases. Clonal variation is the most important individual factor in propagation of hybrid aspen. The fact that the efficiently sprouted clones also rooted best facilitates the selection of clones for large-scale propagation. In practice, root cuttings taken from all parts of the root system of hybrid aspen were capable of producing new shoots and roots. However, for efficient rooting it is important to use roots smaller than one centimeter in diameter. Both rooting and sprouting, as well as sprouting rate, were increased by high soil temperature; in our studies the highest temperature tested (30ºC) was the best. Light accelerated the sprouting of root cuttings, but they rooted best in dark conditions. Rooting is essential because without roots the sprouted cutting cannot survive long. For aspen the criteria for clone selection are primarily fiber qualities and growth rate, but ability to regenerate efficiently is also essential. For large-scale propagation it is very important to find clones from which many cuttings per stock plant can be obtained. In light of production costs, however, it is even more important that the regeneration ability of the produced cuttings be high.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.