Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Experimental studies on new particle formation and ions

Atmospheric aerosol particles have significant climatic effects. Secondary new particle formation is a globally important source of these particles. Currently, the mechanisms of particle formation and the vapours participating in this process are, however, not truly understood. The recently developed Neutral cluster and Air Ion Spectrometer (NAIS) was widely used in field studies of atmospheric particle formation. The NAIS was calibrated and found to be in adequate agreement with the reference instruments. It was concluded that NAIS can be reliably used to measure ions and particles near the sizes where the atmospheric particle formation begins. The main focus of this thesis was to study new particle formation and participation of ions in this process. To attain this objective, particle and ion formation and growth rates were studied in various environments - at several field sites in Europe, in previously rarely studied sites in Antarctica and Siberia and also in an indoor environment. New particle formation was observed at all sites were studied and the observations were used as indicatives of the particle formation mechanisms. Particle size-dependent growth rates and nucleation mode hygroscopic growth factors were examined to obtain information on the particle growth. It was found that the atmospheric ions participate in the initial steps of new particle formation, although their contribution was minor in the boundary layer. The highest atmospheric particle formation rates were observed at the most polluted sites where the role of ions was the least pronounced. Furthermore, the increase of particle growth rate with size suggested that enhancement of the growth by ions was negligible. Participation of organic vapours in the particle growth was supported by laboratory and field observations. It was addressed that secondary new particle formation can also be a significant source of indoor air particles. These results, extending over a wide variety of environments, give support to previous observations and increase understanding on new particle formation on a global scale.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Kymmenen koirarodun perinnöllinen monimuotoisuus

Tämän tutkimuksen tavoitteena oli selvittää koirarotujen sisäistä ja välistä perinnöllistä muuntelua ja populaatioiden rakenteita. Tutkimuksessa käytettiin Finnzymes Oy:ltä saatua koirien mikrosatelliittimerkkeihin perustuvaa genotyypitysaineistoa. Lopullisessa aineistossa oli 395 koiraa kymmenestä keskenään varsin erilaisesta rodusta. Koirien määrä rotua kohti vaihteli 31:stä 53:een. Tutkimuksessa käytettiin 18 mikrosatelliittilokusta. Alleelirikkaus vaihteli mikrosatelliittilokuksissa välillä 2,0 – 9,9. Kaikkein muuntelevin lokus oli AHT137 ja vähiten muunteleva AHTk211. Jackrussellinterrierin alleelirikkaus oli yli kaikkien lokusten tarkasteltuna suurinta ja cavalier kingcharlesinspanielin pienintä. Eniten Hardy-Weinbergin tasapainosta poikkeavia mikrosatelliittilokuksia oli schipperke- rodulla. Coton de tulearin, saksanpaimenkoiran ja suomenlapinkoiran kaikki mikrosatelliittilokukset olivat Hardy-Weinbergin tasapainossa. Cavalier kingcharlesinspanielin havaittu heterotsygotia-aste oli matalin kaikkien lokusten yli tarkasteltuna (0,50) ja suomenlapinkoiran korkein (0,73). Ainoat tilastollisesti merkitsevät FIS-arvot olivat schipperken lokuksessa INU030 (0,39) ja kaikkien lokusten yli tarkasteltuna (0,11). Eniten populaatioiden välisiin eroihin perustuvaa muuntelua oli cavalier kingcharlesinspanielin ja pitkäkarvaisen collien välillä (FST = 0,34) ja vähiten chihuahuan ja coton de tulearin välillä (FST = 0,07). Koko aineistossa noin 17,7 % populaatioiden välisestä geneettisestä muuntelusta johtui populaatioiden välisistä eroista. Rodut ovat tulosten perusteella selvästi erillisiä populaatioita. Coton de tulearin alleeliparit olivat selvästi eniten kytkentäepätasapainossa keskenään (94) ja tiibetinspanielin vähiten (15). Pitkäkarvaisen collien tehollinen populaatiokoko oli pienin (35) ja chihuahuan suurin (86). U seiden populaatiogeneettisten tunnuslukujen perusteella nousivat esiin cavalier kingcharlesinspanieli, pitkäkarvainen collie ja schipperke perinnöllisen muuntelun vähäisyyden perusteella ja chihuahua, jackrussellinterrieri ja suomenlapinkoira keskimääräistä suuremman perinnöllisen muuntelun perusteella. Selityksiä geneettisen monimuotoisuuden vaihteluun näillä roduilla löytyy rotujen historiasta.

Language universals and linguistic complexity : Three case studies in core argument marking

In this dissertation I study language complexity from a typological perspective. Since the structuralist era, it has been assumed that local complexity differences in languages are balanced out in cross-linguistic comparisons and that complexity is not affected by the geopolitical or sociocultural aspects of the speech community. However, these assumptions have seldom been studied systematically from a typological point of view. My objective is to define complexity so that it is possible to compare it across languages and to approach its variation with the methods of quantitative typology. My main empirical research questions are: i) does language complexity vary in any systematic way in local domains, and ii) can language complexity be affected by the geographical or social environment? These questions are studied in three articles, whose findings are summarized in the introduction to the dissertation. In order to enable cross-language comparison, I measure complexity as the description length of the regularities in an entity; I separate it from difficulty, focus on local instead of global complexity, and break it up into different types. This approach helps avoid the problems that plagued earlier metrics of language complexity. My approach to grammar is functional-typological in nature, and the theoretical framework is basic linguistic theory. I delimit the empirical research functionally to the marking of core arguments (the basic participants in the sentence). I assess the distributions of complexity in this domain with multifactorial statistical methods and use different sampling strategies, implementing, for instance, the Greenbergian view of universals as diachronic laws of type preference. My data come from large and balanced samples (up to approximately 850 languages), drawn mainly from reference grammars. The results suggest that various significant trends occur in the marking of core arguments in regard to complexity and that complexity in this domain correlates with population size. These results provide evidence that linguistic patterns interact among themselves in terms of complexity, that language structure adapts to the social environment, and that there may be cognitive mechanisms that limit complexity locally. My approach to complexity and language universals can therefore be successfully applied to empirical data and may serve as a model for further research in these areas.

Tuning the interfacial properties of supported metal nanoclusters

Nanoclusters are objects made up of several to thousands of atoms and form a transitional state of matter between single atoms and bulk materials. Due to their large surface-to-volume ratio, nanoclusters exhibit exciting and yet poorly studied size dependent properties. When deposited directly on bare metal surfaces, the interaction of the cluster with the substrate leads to alteration of the cluster properties, making it less or even non-functional. Surfaces modified with self-assembled monolayers (SAMs) were shown to form an interesting alternative platform, because of the possibility to control wettability by decreasing the surface reactivity and to add functionalities to pre-formed nanoclusters. In this thesis, the underlying size effects and the influence of the nanocluster environment are investigated. The emphasis is on the structural and magnetic properties of nanoclusters and their interaction with thiol SAMs. We report, for the first time, a ferromagnetic-like spin-glass behaviour of uncapped nanosized Au islands tens of nanometres in size. The flattening kinetics of the nanocluster deposition on thiol SAMs are shown to be mediated mainly by the thiol terminal group, as well as the deposition energy and the particle size distribution. On the other hand, a new mechanism for the penetration of the deposited nanoclusters through the monolayers is presented, which is fundamentally different from those reported for atom deposition on alkanethiols. The impinging cluster is shown to compress the thiol layer against the Au surface and subsequently intercalate at the thiol-Au interface. The compressed thiols try then to straighten and push the cluster away from the surface. Depending on the cluster size, this restoring force may or may not enable a covalent cluster-surface bond formation, giving rise to various cluster-surface binding patterns. Compression and straightening of the thiol molecules pinpoint the elastic nature of the SAMs, which has been investigated in this thesis using nanoindentation. The nanoindenation method has been applied to SAMs of varied tail groups, giving insight into the mechanical properties of thiol modified metal surfaces.

Endemic malaria: an 'indoor' disease in northern Europe

Background: Endemic northern malaria reached 68°N latitude in Europe during the 19th century, where the summer mean temperature only irregularly exceeded 16°C, the lower limit needed for sporogony of Plasmodium vivax. Because of the available historical material and little use of quinine, Finland was suitable for an analysis of endemic malaria and temperature. Methods: Annual malaria death frequencies during 1800–1870 extracted from parish records were analysed against long-term temperature records in Finland, Russia and Sweden. Supporting data from 1750–1799 were used in the interpretation of the results. The life cycle and behaviour of the anopheline mosquitoes were interpreted according to the literature. Results: Malaria frequencies correlated strongly with the mean temperature of June and July of the preceding summer, corresponding to larval development of the vector. Hatching of imagoes peaks in the middle of August, when the temperature most years is too low for the sporogony of Plasmodium. After mating some of the females hibernate in human dwellings. If the female gets gametocytes from infective humans, the development of Plasmodium can only continue indoors, in heated buildings. Conclusion: Northern malaria existed in a cold climate by means of summer dormancy of hypnozoites in humans and indoor transmission of sporozoites throughout the winter by semiactive hibernating mosquitoes. Variable climatic conditions did not affect this relationship. The epidemics, however, were regulated by the population size of the mosquitoes which, in turn, ultimately was controlled by the temperatures of the preceding summer.

Population dynamics of blue mussels in a variable environment at the edge of their range

Spatial and temporal variation in the abundance of species can often be ascribed to spatial and temporal variation in the surrounding environment. Knowledge of how biotic and abiotic factors operate over different spatial and temporal scales in determining distribution, abundance, and structure of populations lies at the heart of ecology. The major part of the current ecological theory stems from studies carried out in central parts of the distributional range of species, whereas knowledge of how marginal populations function is inadequate. Understanding how marginal populations, living at the edge of their range, function is however in a key position to advance ecology and evolutionary biology as scientific disciplines. My thesis focuses on the factors affecting dynamics of marginal populations of blue mussels (Mytilus edulis) living close to their tolerance limits with regard to salinity. The thesis aims to highlight the dynamics at the edge of the range and contrast these with dynamics in more central parts of the range in order to understand the potential interplay between the central and the marginal part in the focal system. The objectives of the thesis are approached by studies on: (1) factors affecting regional patterns of the species, (2) long-term temporal dynamics of the focal species spaced along a regional salinity gradient, (3) selective predation by increasing populations of roach (Rutilus rutilus) when feeding on their main food item, the blue mussel, (4) the primary and secondary effects of local wave exposure gradients and (5) the role of small-scale habitat heterogeneity as determinants of large-scale pattern. The thesis shows that populations of blue mussels are largely determined by large scale changes in sea water salinity, affecting mainly recruitment success and longevity of local populations. In opposite to the traditional view, the thesis strongly indicate that vertebrate predators strongly affect abundance and size structure of blue mussel populations, and that the role of these predators increases towards the margin where populations are increasingly top-down controlled. The thesis also indicates that the positive role of biogenic habitat modifiers increases towards the marginal areas, where populations of blue mussels are largely recruitment limited. Finally, the thesis shows that local blue mussel populations are strongly dependent on high water turbulence, and therefore, dense populations are constrained to offshore habitats. Finally, the thesis suggests that ongoing sedimentation of rocky shores is detrimental for the species, affecting recruitment success and post-recruit survival, pushing stable mussel beds towards offshore areas. Ongoing large scale changes in the Baltic Sea, especially dilution processes with attendant effects, are predicted to substantially contract the distributional range of the mussel, but also affect more central populations. The thesis shows that in order to understand the functioning of marginal populations, research should (1) strive for multi-scale approaches in order to link ecosystem patterns with ecosystem processes, and (2) challenge the prevailing tenets that origin from research carried out in central areas that may not be valid at the edge.

Mutations as molecular tools : The metabolic-rate dependent molecular clock and DNA barcoding of allied species

Mutation and recombination are the fundamental processes leading to genetic variation in natural populations. This variation forms the raw material for evolution through natural selection and drift. Therefore, studying mutation rates may reveal information about evolutionary histories as well as phylogenetic interrelationships of organisms. In this thesis two molecular tools, DNA barcoding and the molecular clock were examined. In the first part, the efficiency of mutations to delineate closely related species was tested and the implications for conservation practices were assessed. The second part investigated the proposition that a constant mutation rate exists within invertebrates, in form of a metabolic-rate dependent molecular clock, which can be applied to accurately date speciation events. DNA barcoding aspires to be an efficient technique to not only distinguish between species but also reveal population-level variation solely relying on mutations found on a short stretch of a single gene. In this thesis barcoding was applied to discriminate between Hylochares populations from Russian Karelia and new Hylochares findings from the greater Helsinki region in Finland. Although barcoding failed to delineate the two reproductively isolated groups, their distinct morphological features and differing life-history traits led to their classification as two closely related, although separate species. The lack of genetic differentiation appears to be due to a recent divergence event not yet reflected in the beetles molecular make-up. Thus, the Russian Hylochares was described as a new species. The Finnish species, previously considered as locally extinct, was recognized as endangered. Even if, due to their identical genetic make-up, the populations had been regarded as conspecific, conservation strategies based on prior knowledge from Russia would not have guaranteed the survival of the Finnish beetle. Therefore, new conservation actions based on detailed studies of the biology and life-history of the Finnish Hylochares were conducted to protect this endemic rarity in Finland. The idea behind the strict molecular clock is that mutation rates are constant over evolutionary time and may thus be used to infer species divergence dates. However, one of the most recent theories argues that a strict clock does not tick per unit of time but that it has a constant substitution rate per unit of mass-specific metabolic energy. Therefore, according to this hypothesis, molecular clocks have to be recalibrated taking body size and temperature into account. This thesis tested the temperature effect on mutation rates in equally sized invertebrates. For the first dataset (family Eucnemidae, Coleoptera) the phylogenetic interrelationships and evolutionary history of the genus Arrhipis had to be inferred before the influence of temperature on substitution rates could be studied. Further, a second, larger invertebrate dataset (family Syrphidae, Diptera) was employed. Several methodological approaches, a number of genes and multiple molecular clock models revealed that there was no consistent relationship between temperature and mutation rate for the taxa under study. Thus, the body size effect, observed in vertebrates but controversial for invertebrates, rather than temperature may be the underlying driving force behind the metabolic-rate dependent molecular clock. Therefore, the metabolic-rate dependent molecular clock does not hold for the here studied invertebrate groups. This thesis emphasizes that molecular techniques relying on mutation rates have to be applied with caution. Whereas they may work satisfactorily under certain conditions for specific taxa, they may fail for others. The molecular clock as well as DNA barcoding should incorporate all the information and data available to obtain comprehensive estimations of the existing biodiversity and its evolutionary history.

Intraspecific variation in brain size and architecture : population divergence and phenotypic plasticity

Brain size and architecture exhibit great evolutionary and ontogenetic variation. Yet, studies on population variation (within a single species) in brain size and architecture, or in brain plasticity induced by ecologically relevant biotic factors have been largely overlooked. Here, I address the following questions: (i) do locally adapted populations differ in brain size and architecture, (ii) can the biotic environment induce brain plasticity, and (iii) do locally adapted populations differ in levels of brain plasticity? In the first two chapters I report large variation in both absolute and relative brain size, as well as in the relative sizes of brain parts, among divergent nine-spined stickleback (Pungitius pungitius) populations. Some traits show habitat-dependent divergence, implying natural selection being responsible for the observed patterns. Namely, marine sticklebacks have relatively larger bulbi olfactorii (chemosensory centre) and telencephala (involved in learning) than pond sticklebacks. Further, I demonstrate the importance of common garden studies in drawing firm evolutionary conclusions. In the following three chapters I show how the social environment and perceived predation risk shapes brain development. In common frog (Rana temporaria) tadpoles, I demonstrate that under the highest per capita predation risk, tadpoles develop smaller brains than in less risky situations, while high tadpole density results in enlarged tectum opticum (visual brain centre). Visual contact with conspecifics induces enlarged tecta optica in nine-spined sticklebacks, whereas when only olfactory cues from conspecifics are available, bulbus olfactorius become enlarged.Perceived predation risk results in smaller hypothalami (complex function) in sticklebacks. Further, group-living has a negative effect on relative brain size in the competition-adapted pond sticklebacks, but not in the predation-adapted marine sticklebacks. Perceived predation risk induces enlargement of bulbus olfactorius in pond sticklebacks, but not in marine sticklebacks who have larger bulbi olfactorii than pond fish regardless of predation. In sum, my studies demonstrate how applying a microevolutionary approach can help us to understand the enormous variation observed in the brains of wild animals a point-of-view which I high-light in the closing review chapter of my thesis.

Living on the edge : Population genetics of Finno-Ugric-speaking humans in North Eurasia

The major aim of this thesis was to examine the origins and distribution of uniparental and autosomal genetic variation among the Finno-Ugric-speaking human populations living in Boreal and Arctic regions of North Eurasia. In more detail, I aimed to disentangle the underlying molecular and population genetic factors which have produced the patterns of uniparental and autosomal genetic diversity in these populations. Among Finno-Ugrics the genetic amalgamation and clinal distribution of West and East Eurasian gene pools were observed within uniparental markers. This admixture indicates that North Eurasia was colonized through Central Asia/ South Siberia by human groups already carrying both West and East Eurasian lineages. The complex combination of founder effects, gene flow and genetic drift underlying the genetic diversity of the Finno-Ugric- speaking populations were emphasized by low haplotype diversity within and among uniparental and biparental markers. A high prevalence of lactase persistence allele among the North Eurasian Finno- Ugric agriculturalist populations was also shown indicating a local adaptation to subsistence change with lactose rich diet. Moreover, the haplotype background of lactase persistence allele among the Finno- Ugric-speakers strongly suggested that the lactase persistence T-13910 mutation was introduced independently more than once to the North Eurasian gene pool. A significant difference in genetic diversity, haplotype structure and LD distribution within the cytochrome P450 CYP2C and CYP2D regions revealed the unique gene pool of the Finno-Ugric Saami created mainly by population genetic processes compared to other Europeans and sub-Saharan Mandenka population. From all studied populations the Saami showed also significantly the highest allele frequency of a CYP2C19 gene mutation causing variable drug reactions. The diversity patterns observed within CYP2C and CYP2D regions emphasize the strong effect of demographic history shaping genetic diversity and LD especially among such small and constant size populations as the Finno-Ugric-speaking Saami. Moreover, the increased LD in Saami due to genetic drift and/or admixture was shown to offer an advantage for further attempts to identify alleles associated to common complex pharmacogenetic traits.

Plant secondary compounds and soil microbial processes in carbon and nitrogen cycling in relation to tree species

The aim of this study was to explore soil microbial activities related to C and N cycling and the occurrence and concentrations of two important groups of plant secondary compounds, terpenes and phenolic compounds, under silver birch (Betula pendula Roth), Norway spruce (Picea abies (L.) Karst) and Scots pine (Pinus sylvestris L.) as well as to study the effects of volatile monoterpenes and tannins on soil microbial activities. The study site, located in Kivalo, northern Finland, included ca. 70-year-old adjacent stands dominated by silver birch, Norway spruce and Scots pine. Originally the soil was very probably similar in all three stands. All forest floor layers (litter (L), fermentation layer (F) and humified layer (H)) under birch and spruce showed higher rates of CO2 production, greater net mineralisation of nitrogen and higher amounts of carbon and nitrogen in microbial biomass than did the forest floor layers under pine. Concentrations of mono-, sesqui-, di- and triterpenes were higher under both conifers than under birch, while the concentration of total water-soluble phenolic compounds as well as the concentration of condensed tannins tended to be higher or at least as high under spruce as under birch or pine. In general, differences between tree species in soil microbial activities and in concentrations of secondary compounds were smaller in the H layer than in the upper layers. The rate of CO2 production and the amount of carbon in the microbial biomass correlated highly positively with the concentration of total water-soluble phenolic compounds and positively with the concentration of condensed tannins. Exposure of soil to volatile monoterpenes and tannins extracted and fractionated from spruce and pine needles affected carbon and nitrogen transformations in soil, but the effects were dependent on the compound and its molecular structure. Monoterpenes decreased net mineralisation of nitrogen and probably had a toxic effect on part of the microbial population in soil, while another part of the microbes seemed to be able to use monoterpenes as a carbon source. With tannins, low-molecular-weight compounds (also compounds other than tannins) increased soil CO2 production and nitrogen immobilisation by soil microbes while the higher-molecular-weight condensed tannins had inhibitory effects. In conclusion, plant secondary compounds may have a great potential in regulation of C and N transformations in forest soils, but the real magnitude of their significance in soil processes is impossible to estimate.

Context-dependent mate choice in the sand goby, Pomatoschistus minutus

The results presented in this thesis show that all females of a given population do not necessarily choose similar mating partners. Specifically, partner preferences of a fish, the sand goby (Pomatoschistus minutus), varied among individual females and depended on the social context at the time of choice. I also show that females assess multiple mate choice cues simultaneously; partner preferences were based more strongly on an interaction effect between different choice cues than on any individual cue. Furthermore, I found that preferred matings involved fitness benefits in the form of increased offspring success, but these benefits were not significantly affected by mate compatibility. Hence, mate choice for partner compatibility does not appear to be an important determinant of the observed variation in female mate preferences in this species. The context-dependency of female mating preferences revealed is relevant to how genetic variation in sexually selected traits might be maintained: as the mating success of a certain male type varies according to the choice context, directional sexual selection on male traits is shown to be less intense than generally thought making for a slower loss of genetic variation in these traits. Mating preferences of sand gobies were assessed by giving females a binary choice between males that differed in body size and/or other focus traits. These association preferences were found to be sexually motivated, repeatable and to correspond to actual mating decisions.

Dispersal and related life history traits in the Glanville fritillary butterfly

Most studies of life history evolution are based on the assumption that species exist at equilibrium and spatially distinct separated populations. In reality, this is rarely the case, as populations are often spatially structured with ephemeral local populations. Therefore, the characteristics of metapopulations should be considered while studying factors affecting life history evolution. Theoretical studies have examined spatial processes shaping the evolution of life history traits to some extent, but there is little empirical data and evidence to investigate model predictions. In my thesis I have tried to bridge the gap between theoretical and empirical studies by using the well-known Glanville fritillary (Melitaea cinxia) metapopulation as a model system. The long-term persistence of classic metapopulations requires sufficient dispersal to establish new local populations to compensate for local extinctions. Previous studies on the Glanville fritillary have shown that females establishing new populations are not a random sample from the metapopulation, but they are in fact more dispersive than females in old populations. Many other life-history traits, such as body size, fecundity and lifespan, may be related to dispersal rate. Therefore, I examined a range of correlated traits for their evolutionary and ecological consequences. I was particularly interested in how the traits vary under natural environmental conditions, hence all studies were conducted in a large (32 x 26 m) outdoor population cage built upon a natural habitat patch. Individuals for the experiments were sampled from newly-established and old populations within a large metapopulation. Results show that females originating from newly-established populations had higher within-habitat patch mobility than females from old populations. I showed that dispersal rate is heritable and that flight activity is related to variation in a gene encoding the glycolytic enzyme phosphoglucose isomerase. Both among-individual and among-population variation in dispersal are correlated with the reproductive performance of females, though I found no evidence for a trade-off between dispersal and fecundity in terms of lifetime egg production or clutch size. Instead, the results suggest that highly dispersive females from newly-established populations have a shorter lifespan than females from old populations, and that dispersive females may pay a cost in terms of reduced lifetime reproductive success due to increased time spent outside habitat patches. In summary, the results of this thesis show that genotype-dependent dispersal rate correlates with other life history traits in the Glanville fritillary, and that the rapid turnover of local populations (extinctions and re-colonisations) is likely to be the mechanism that maintains phenotypic variation in many life history traits at the metapopulation level.

Conceptual and statistical modelling of environmental effects in population dynamics

Population dynamics are generally viewed as the result of intrinsic (purely density dependent) and extrinsic (environmental) processes. Both components, and potential interactions between those two, have to be modelled in order to understand and predict dynamics of natural populations; a topic that is of great importance in population management and conservation. This thesis focuses on modelling environmental effects in population dynamics and how effects of potentially relevant environmental variables can be statistically identified and quantified from time series data. Chapter I presents some useful models of multiplicative environmental effects for unstructured density dependent populations. The presented models can be written as standard multiple regression models that are easy to fit to data. Chapters II IV constitute empirical studies that statistically model environmental effects on population dynamics of several migratory bird species with different life history characteristics and migration strategies. In Chapter II, spruce cone crops are found to have a strong positive effect on the population growth of the great spotted woodpecker (Dendrocopos major), while cone crops of pine another important food resource for the species do not effectively explain population growth. The study compares rate- and ratio-dependent effects of cone availability, using state-space models that distinguish between process and observation error in the time series data. Chapter III shows how drought, in combination with settling behaviour during migration, produces asymmetric spatially synchronous patterns of population dynamics in North American ducks (genus Anas). Chapter IV investigates the dynamics of a Finnish population of skylark (Alauda arvensis), and point out effects of rainfall and habitat quality on population growth. Because the skylark time series and some of the environmental variables included show strong positive autocorrelation, the statistical significances are calculated using a Monte Carlo method, where random autocorrelated time series are generated. Chapter V is a simulation-based study, showing that ignoring observation error in analyses of population time series data can bias the estimated effects and measures of uncertainty, if the environmental variables are autocorrelated. It is concluded that the use of state-space models is an effective way to reach more accurate results. In summary, there are several biological assumptions and methodological issues that can affect the inferential outcome when estimating environmental effects from time series data, and that therefore need special attention. The functional form of the environmental effects and potential interactions between environment and population density are important to deal with. Other issues that should be considered are assumptions about density dependent regulation, modelling potential observation error, and when needed, accounting for spatial and/or temporal autocorrelation.