Effects of environmental variation on ecological and evolutionary dynamics

Environmental variation is a fact of life for all the species on earth: for any population of any particular species, the local environmental conditions are liable to vary in both time and space. In today's world, anthropogenic activity is causing habitat loss and fragmentation for many species, which may profoundly alter the characteristics of environmental variation in remaining habitat. Previous research indicates that, as habitat is lost, the spatial configuration of remaining habitat will increasingly affect the dynamics by which populations are governed. Through the use of mathematical models, this thesis asks how environmental variation interacts with species properties to influence population dynamics, local adaptation, and dispersal evolution. More specifically, we couple continuous-time continuous-space stochastic population dynamic models to landscape models. We manipulate environmental variation via parameters such as mean patch size, patch density, and patch longevity. Among other findings, we show that a mixture of high and low quality habitat is commonly better for a population than uniformly mediocre habitat. This conclusion is justified by purely ecological arguments, yet the positive effects of landscape heterogeneity may be enhanced further by local adaptation, and by the evolution of short-ranged dispersal. The predicted evolutionary responses to environmental variation are complex, however, since they involve numerous conflicting factors. We discuss why the species that have high levels of local adaptation within their ranges may not be the same species that benefit from local adaptation during range expansion. We show how habitat loss can lead to either increased or decreased selection for dispersal depending on the type of habitat and the manner in which it is lost. To study the models, we develop a recent analytical method, Perturbation expansion, to enable the incorporation of environmental variation. Within this context, we use two methods to address evolutionary dynamics: Adaptive dynamics, which assumes mutations occur infrequently so that the ecological and evolutionary timescales can be separated, and via Genotype distributions, which assume mutations are more frequent. The two approaches generally lead to similar predictions yet, exceptionally, we show how the evolutionary response of dispersal behaviour to habitat turnover may qualitatively depend on the mutation rate.

Ecological communities in variable environments : dynamics and diversity under coloured environmental stochasticity

While environmental variation is an ubiquitous phenomenon in the natural world which has for long been appreciated by the scientific community recent changes in global climatic conditions have begun to raise consciousness about the economical, political and sociological ramifications of global climate change. Climate warming has already resulted in documented changes in ecosystem functioning, with direct repercussions on ecosystem services. While predicting the influence of ecosystem changes on vital ecosystem services can be extremely difficult, knowledge of the organisation of ecological interactions within natural communities can help us better understand climate driven changes in ecosystems. The role of environmental variation as an agent mediating population extinctions is likely to become increasingly important in the future. In previous studies population extinction risk in stochastic environmental conditions has been tied to an interaction between population density dependence and the temporal autocorrelation of environmental fluctuations. When populations interact with each other, forming ecological communities, the response of such species assemblages to environmental stochasticity can depend, e.g., on trophic structure in the food web and the similarity in species-specific responses to environmental conditions. The results presented in this thesis indicate that variation in the correlation structure between species-specific environmental responses (environmental correlation) can have important qualitative and quantitative effects on community persistence and biomass stability in autocorrelated (coloured) environments. In addition, reddened environmental stochasticity and ecological drift processes (such as demographic stochasticity and dispersal limitation) have important implications for patterns in species relative abundances and community dynamics over time and space. Our understanding of patterns in biodiversity at local and global scale can be enhanced by considering the relevance of different drift processes for community organisation and dynamics. Although the results laid out in this thesis are based on mathematical simulation models, they can be valuable in planning effective empirical studies as well as in interpreting existing empirical results. Most of the metrics considered here are directly applicable to empirical data.

Dynamics of contour, object and face processing in the human visual cortex

The neural basis of visual perception can be understood only when the sequence of cortical activity underlying successful recognition is known. The early steps in this processing chain, from retina to the primary visual cortex, are highly local, and the perception of more complex shapes requires integration of the local information. In Study I of this thesis, the progression from local to global visual analysis was assessed by recording cortical magnetoencephalographic (MEG) responses to arrays of elements that either did or did not form global contours. The results demonstrated two spatially and temporally distinct stages of processing: The first, emerging 70 ms after stimulus onset around the calcarine sulcus, was sensitive to local features only, whereas the second, starting at 130 ms across the occipital and posterior parietal cortices, reflected the global configuration. To explore the links between cortical activity and visual recognition, Studies II III presented subjects with recognition tasks of varying levels of difficulty. The occipito-temporal responses from 150 ms onwards were closely linked to recognition performance, in contrast to the 100-ms mid-occipital responses. The averaged responses increased gradually as a function of recognition performance, and further analysis (Study III) showed the single response strengths to be graded as well. Study IV addressed the attention dependence of the different processing stages: Occipito-temporal responses peaking around 150 ms depended on the content of the visual field (faces vs. houses), whereas the later and more sustained activity was strongly modulated by the observers attention. Hemodynamic responses paralleled the pattern of the more sustained electrophysiological responses. Study V assessed the temporal processing capacity of the human object recognition system. Above sufficient luminance, contrast and size of the object, the processing speed was not limited by such low-level factors. Taken together, these studies demonstrate several distinct stages in the cortical activation sequence underlying the object recognition chain, reflecting the level of feature integration, difficulty of recognition, and direction of attention.

Tree mortality and deadwood dynamics in late-successional boreal forests

Here I aimed at quantifying the main components of deadwood dynamics, i.e. tree mortality, deadwood pools, and their decomposition, in late-successional boreal forests. I focused on standing dead trees in three stand types dominated by Picea mariana and Abies balsamea in eastern Canada, and on standing and down dead trees in Picea abies-dominated stands in three areas in Northern Europe. Dead and living trees were measured on five sample plots of 1.6-ha size in each study area and stand type. Stem disks from dead trees were sampled to determine wood density and year of death, using dendrochronological methods. The results were applied to reconstruct past tree mortality and to model deadwood decay class dynamics. Site productivity, stand developmental stage, and the occurrence of episodic tree mortality influenced deadwood volume and quality. In all study areas tree mortality was continuous, leading to continuity in deadwood decay stage distribution. Episodic tree mortality due to either autogenic or allogenic causes influenced deadwood volume and quality in all but one study area. However, regardless of productivity and disturbance history deadwood was abundant, accounting for 20 53% of total wood volume in European study areas, and 15 27% of total standing volume in eastern Canada. Deadwood was a persistent structural component, since its expected residence time in early- and midstages of decay was 18 yr even in the area with the most rapid decomposition. The results indicated that in the absence of episodic tree mortality, stands may eventually develop to a steady state, in which deadwood volume fluctuates around an equilibrium state. However, in many forests deadwood is naturally variable, due to recurrent moderate-severity disturbances. This variability, the continuous tree mortality, and variation in rates of wood decomposition determine the dynamics and availability of deadwood as a habitat and carbon storage medium in boreal coniferous forest ecosystems.

Fen ecosystem carbon gas dynamics in changing hydrological conditions

Northern peatlands are thought to store one third of all soil carbon (C). Besides the C sink function, peatlands are one of the largest natural sources of methane (CH4) to the atmosphere. Climate change may affect the C gas dynamics as well as the labile C pool. Because the peatland C sequestration and CH4 emissions are governed by high water levels, changes in hydrology are seen as the driving factor in peatland ecosystem change. This study aimed to quantify the carbon dioxide (CO2) and CH4 dynamics of a fen ecosystem at different spatial scales: plant community components scale, plant community scale and ecosystem scale, under hydrologically normal and water level drawdown conditions. C gas exchange was measured in two fens in southern Finland applying static chamber and eddy covariance techniques. During hydrologically normal conditions, the ecosystem was a CO2 sink and CH4 source to the atmosphere. Sphagnum mosses and sedges were the most important contributors to the community photosynthesis. The presence of sedges had a major positive impact on CH4 emissions while dwarf shrubs had a slightly attenuating impact. C fluxes varied considerably between the plant communities. Therefore, their proportions determined the ecosystem scale fluxes. An experimental water level drawdown markedly reduced the photosynthesis and respiration of sedges and Sphagnum mosses and benefited shrubs. Consequently, changes were smaller at the ecosystem scale than at the plant group scale. The decrease in photosynthesis and the increase in respiration, mostly peat respiration, made the fen a smaller CO2 sink. CH4 fluxes were significantly lowered, close to zero. The impact of natural droughts was similar to, although more modest than, the impact of the experimental water level drawdown. The results are applicable to the short term impacts of the water level drawdown and to climatic conditions in which droughts become more frequent.

Stand structural dynamics on pristine and managed boreal peatlands

The objectives of this study were to investigate the stand structure and succession dynamics in Scots pine (Pinus sylvestris L.) stands on pristine peatlands and in Scots pine and Norway spruce (Picea abies (L.) Karst.) dominated stands on drained peatlands. Furthermore, my focus was on characterising how the inherent and environmental factors and the intermediate thinnings modify the stand structure and succession. For pristine peatlands, the study was based on inventorial stand data, while for drained peatlands, longitudinal data from repeatedly measured stands were utilised. The studied sites covered the most common peatland site types in Finland. They were classified into two categories according to the ecohydrological properties related to microsite variation and nutrient levels within sites. Tree DBH and age distributions in relation to climate and site type were used to study the stand dynamics on pristine sites. On drained sites, the Weibull function was used to parameterise the DBH distributions and mixed linear models were constructed to characterise the impacts of different ecological factors on stand dynamics. On pristine peatlands, both climate and the ecohydrology of the site proved to be crucial factors determining the stand structure and its dynamics. Irrespective of the vegetation succession, enhanced site productivity and increased stand stocking they significantly affected the stand dynamics also on drained sites. On the most stocked sites on pristine peatlands the inter-tree competition seemed to also be a significant factor modifying stand dynamics. Tree age and size diversity increased with stand age, but levelled out in the long term. After drainage, the stand structural unevenness increased due to the regeneration and/or ingrowth of the trees. This increase was more pronounced on sparsely forested composite sites than on more fully stocked genuine forested sites in Scots pine stands, which further undergo the formation of birch and spruce undergrowth beneath the overstory as succession proceeds. At 20-30 years after drainage the structural heterogeneity started to decrease, indicating increased inter-tree competition, which increased the mortality of suppressed trees within stand. Peatland stands are more dynamic than anticipated and are generally not characterized by a balanced, self-perpetuating structure. On pristine sites, various successional pathways are possible, whereas on drained sites the succession has more uniform trend. Typically, stand succession proceeds without any distinct developmental stages on pristine peatlands, whereas on drained peatlands, at least three distinct stages could be identified. Thinnings had only little impact on the stand succession. The new information on stand dynamics may be utilised, e.g. in forest management planning to facilitate the allocation of the growth resources to the desired crop component by appropriate silvicultural treatments, as well as assist in assessing the effects of the climate change on the forested boreal peatlands.

The decline of northern malaria and population dynamics of Plamodium vivax

Europe was declared malaria free in 1975. The disappearance of malaria has traditionally been attributed to numerous deliberate actions like vector control, the screening of houses, more efficient medication etc. Malaria, however, disappeared from many countries like Finland before any counter measures had even started. The aim of this thesis is to study the population ecology of P. vivax and its interaction with the human host and the vector. By finding the factors that attributed to the extinction of vivax malaria it might be possible to improve the modern strategy against P. vivax. The parasite was studied with data from Finland, which provides the longest time series (1749-2008) of malaria statistics in the world. The malaria vectors, Anopheles messeae and A. beklemishevi are still common species in the country. The eradication of vivax malaria is difficult because the parasite has a dormant stage that can cause a relapse long after a primary infection. It was now shown that P. vivax is able to detect the presence of a potential vector. A dormant stage is triggered even from a bite of an uninfected Anopheles mosquito. This optimizes the chances for the Plasmodium to reach a mosquito vector for sexual reproduction. The longevity of the dormant stage could be shown to be at least nine years. The parasite spends several years in its human host and the behaviour of the human carrier had a profound impact on the decline of the disease in Finland. Malaria spring epidemics could be explained by a previous warm summer. Neither annual nor summer mean temperature had any impact on the long term malaria trend. Malaria disappeared slowly from Finland without mosquito control. The sociological change from extended families to nuclear families led to decreased household size. The decreased household size correlated strongly with the decline of malaria. That led to an increased isolation of the subpopulations of P. vivax. Their habitat consisted of the bedrooms in which human carriers slept together with the overwintering vectors. The isolation of the parasite ultimately led to the extinction of vivax malaria. Metapopulation models adapted to local conditions should therefore be implemented as a tool for settlement planning and socio-economic development and become an integrated part of the fight against malaria.

Nitrogen dynamics and nitrogen use efficiency of spring cereals under Finnish growing conditions

Nitrogen (N) is one of the main inputs in cereal cultivation and as more than half of the arable land in Finland is used for cereal production, N has contributed substantially to agricultural pollution through fertilizer leaching and runoff. Based on this global phenomenon, the European Community has launched several directives to reduce agricultural emissions to the environment. Trough such measures, and by using economic incentives, it is expected that northern European agricultural practices will, in the future, include reduced N fertilizer application rates. Reduced use of N fertilizer is likely to decrease both production costs and pollution, but could also result in reduced yields and quality if crops experience temporary N deficiency. Therefore, more efficient N use in cereal production, to minimize pollution risks and maximize farmer income, represents a current challenge for agronomic research in the northern growing areas. The main objective of this study was to determine the differences in nitrogen use efficiency (NUE) among spring cereals grown in Finland. Additional aims were to characterize the multiple roles of NUE by analysing the extent of variation in NUE and its component traits among different cultivars, and to understand how other physiological traits, especially radiation use efficiency (RUE) and light interception, affect and interact with the main components of NUE and contribute to differences among cultivars. This study included cultivars of barley (Hordeum vulgare L.), oat (Avena sativa L.) and wheat (Triticum aestivum L.). Field experiments were conducted between 2001 and 2004 at Jokioinen, in Finland. To determine differences in NUE among cultivars and gauge the achievements of plant breeding in NUE, 17-18 cultivars of each of the three cereal species released between 1909 and 2002 were studied. Responses to nitrogen of landraces, old cultivars and modern cultivars of each cereal species were evaluated under two N regimes (0 and 90 kg N ha-1). Results of the study revealed that modern wheat, oat and barley cultivars had similar NUE values under Finnish growing conditions and only results from a wider range of cultivars indicated that wheat cultivars could have lower NUE than the other species. There was a clear relationship between nitrogen uptake efficiency (UPE) and NUE in all species whereas nitrogen utilization efficiency (UTE) had a strong positive relationship with NUE only for oat. UTE was clearly lower in wheat than in other species. Other traits related to N translocation indicated that wheat also had a lower harvest index, nitrogen harvest index and nitrogen remobilisation efficiency and therefore its N translocation efficiency was confirmed to be very low. On the basis of these results there appears to be potential and also a need for improvement in NUE. These results may help understand the underlying physiological differences in NUE and could help to identify alternative production options, such as the different roles that species can play in crop rotations designed to meet the demands of modern agricultural practices.

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.

Water-ethanol mixtures by molecular dynamics and x-ray Compton scattering

Water-ethanol mixtures are commonly used in industry and house holds. However, quite surprisingly their molecular-level structure is still not completely understood. In particular, there is evidence that the local intermolecular geometries depend significantly on the concentration. The aim of this study was to gain information on the molecular-level structures of water-ethanol mixtures by two computational methods. The methods are classical molecular dynamics (MD), where the movement of molecules can be studied, and x-ray Compton scattering, in which the scattering cross section is sensitive to the electron momentum density. Firstly, the water-ethanol mixtures were studied with MD simulations, with the mixture concentration ranging from 0 to 100%. For the simulations well-established force fields were used for the water and ethanol molecules (TIP4P and OPLS-AA, respectively). Moreover, two models were used for ethanol, rigid and non-rigid. In the rigid model the intramolecular bond lengths are fixed, whereas in the non-rigid model the lengths are determined by harmonic potentials. Secondly, mixtures with three different concentrations employing both ethanol models were studied by calculating the experimentally observable x-ray quantity, the Compton profile. In the MD simulations a slight underestimation in the density was observed as compared to experiment. Furthermore, a positive excess of hydrogen bonding with water molecules and a negative one with ethanol was quantified. Also, the mixture was found more structured when the ethanol concentration was higher. Negligible differences in the results were found between the two ethanol models. In contrast, in the Compton scattering results a notable difference between the ethanol models was observed. For the rigid model the Compton profiles were similar for all the concentrations, but for the non-rigid model they were distinct. This leads to two possibilities of how the mixing occurs. Either the mixing is similar in all concentrations (as suggested by the rigid model) or the mixing changes for different concentrations (as suggested by the non-rigid model). Either way, this study shows that the choice of the force field is essential in the microscopic structure formation in the MD simulations. When the sources of uncertainty in the calculated Compton profiles were analyzed, it was found that more statistics needs to be collected to reduce the statistical uncertainty in the final results. The obtained Compton scattering results can be considered somewhat preliminary, but clearly indicative of the behaviour of the water-ethanol mixtures when the force field is modified. The next step is to collect more statistics and compare the results with experimental data to decide which ethanol model describes the mixture better. This way, valuable information on the microscopic structure of water-ethanol mixtures can be found. In addition, information on the force fields in the MD simulations and on the ability of the MD simulations to reproduce the microscopic structure of binary liquids is obtained.