Short-term effects of variable retention on epigaeic spiders and carabid beetles in Finland

Forestry has influenced forest dwelling organisms for centuries in Fennoscandia. For example, in Finland ca. 30% of the threatened species are threatened because of forestry. Nowadays forest management recommendations include practices aimed at maintaining biodiversity in harvesting, such as green-tree retention. However, the effects of these practices have been little studied. In variable retention, different numbers of trees are retained, varying from green-tree retention (at least a few live standing trees in clear-cuts) to thinning (only individual trees removed). I examined the responses of ground-dwelling spiders and carabid beetles to green-tree retention (with small and large tree groups), gap felling and thinning aimed at an uneven age structure of trees. The impacts of these harvesting methods were compared to those of clear-cutting and uncut controls. I aimed to test the hypothesis that retaining more trees positively affects populations of those species of spiders and carabids that were present before harvesting. The data come from two studies. First, spiders were collected with pitfall traps in south-central Finland in 1995 (pre-treatment) and 1998 (after-treatment) in order to examine the effects of clear-cutting, green-tree retention (with 0.01-0.02-ha sized tree groups), gap felling (with three 0.16-ha sized openings in a 1-ha stand), thinning aiming at an uneven age structure of trees and uncut control. Second, spiders and carabids were caught with pitfall traps in eastern Finland in 1998-2001 (pre-treatment and three post-treatment years) in eleven 0.09-0.55-ha sized retention-tree groups and clear-cuts adjacent to them. Original spider and carabid assemblages were better maintained after harvests that retained more trees. Thinning maintained forest spiders well. However, gap felling and large retention-tree groups maintained some forest spider and carabid species in the short-term, but negatively affected some species over time. However, use of small retention-tree groups was associated with negative effects on forest spider populations. Studies are needed on the long-term effects of variable retention on terrestrial invertebrates; especially those directed at defining appropriate retention patch size and on the importance of structural diversity provided by variable retention for invertebrate populations. However, the aims of variable retention should be specified first. For example, are retention-tree groups planned to constitute life-boats , stepping-stones or to create structural diversity? Does it suffice that some species are maintained, or do we want to preserve the most sensitive ones, and how are these best defined? Moreover, the ecological benefits and economic costs of modified logging methods should be compared to other approaches aimed at maintaining biodiversity.

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.

Blue mussel beds as biodiversity hotspots on the rocky shores of the northern Baltic Sea

ABSTRACT The Baltic Sea is a vulnerable ecosystem currently undergoing a number of changes, both natural and human induced. The changes are likely to affect the species found on these shores, e.g. their distribution and interactions with other species. Blue mussels (Mytilus trossulus x Mytilus edulis) provide one of the main biogenic hard structures on the shallow shores of the Baltic Sea where they aggregate into dense beds and provide a number of resources for over 40 associated macrofaunal species, thus functioning as ecosystem engineers. The blue mussel, being a marine species, is highly likely to be affected by any changes in sea water salinity, circulation and/or water balance. These changes could trickle down also to affect the associated macrofaunal communities. The aims of this thesis were three-fold: first, I examined and described the macrofaunal communities found within blue mussel patches since the fauna associated with mussel patches had never been described in the study area prior to this thesis. Second, I explored how changes in mussel density, size as well as patch size and shape would affect the mussel communities. Finally, I tested how general landscape theories derived from terrestrial studies function in blue mussel systems. Theories included the structural heterogeneity hypothesis, species-area relationships, edge effects and patch isolation effects. The work shows that blue mussels in the northern Baltic Sea have an indisputable function as diversity hotspots and that the faunal assemblages found in mussel patches are extremely rich and unique. Further on, it shows that changes in mussel biomass, size, patch size and amount of edge have the potential to alter the faunal assemblages and diversity within patches. Finally, it shows that although some landscape theories, such as the structural heterogeneity hypothesis, seem to apply also in blue mussel communities, others cannot be directly applied due to the different prevailing conditions in the study system. This is a pioneering work looking at diversity shaping processes on the rocky shores of the Gulf of Finland, making up over 40% of the total water basin. A focus on niche construction, positive facilitation effects and ecosystem engineering could provide new insights and methods for conservation biology, but before this can be done, we need to fully understand the circumstances under which a species becomes an ecosystem engineer and recognize the systems in which it functions.

Plant species diversity of buffer zones in agricultural landscapes: in search of determinants from the local to regional scale

Buffer zones are vegetated strip-edges of agricultural fields along watercourses. As linear habitats in agricultural ecosystems, buffer strips dominate and play a leading ecological role in many areas. This thesis focuses on the plant species diversity of the buffer zones in a Finnish agricultural landscape. The main objective of the present study is to identify the determinants of floral species diversity in arable buffer zones from local to regional levels. This study was conducted in a watershed area of a farmland landscape of southern Finland. The study area, Lepsämänjoki, is situated in the Nurmijärvi commune 30 km to the north of Helsinki, Finland. The biotope mosaics were mapped in GIS. A total of 59 buffer zones were surveyed, of which 29 buffer strips surveyed were also sampled by plot. Firstly, two diversity components (species richness and evenness) were investigated to determine whether the relationship between the two is equal and predictable. I found no correlation between species richness and evenness. The relationship between richness and evenness is unpredictable in a small-scale human-shaped ecosystem. Ordination and correlation analyses show that richness and evenness may result from different ecological processes, and thus should be considered separately. Species richness correlated negatively with phosphorus content, and species evenness correlated negatively with the ratio of organic carbon to total nitrogen in soil. The lack of a consistent pattern in the relationship between these two components may be due to site-specific variation in resource utilization by plant species. Within-habitat configuration (width, length, and area) were investigated to determine which is more effective for predicting species richness. More species per unit area increment could be obtained from widening the buffer strip than from lengthening it. The width of the strips is an effective determinant of plant species richness. The increase in species diversity with an increase in the width of buffer strips may be due to cross-sectional habitat gradients within the linear patches. This result can serve as a reference for policy makers, and has application value in agricultural management. In the framework of metacommunity theory, I found that both mass effect(connectivity) and species sorting (resource heterogeneity) were likely to explain species composition and diversity on a local and regional scale. The local and regional processes were interactively dominated by the degree to which dispersal perturbs local communities. In the lowly and intermediately connected regions, species sorting was of primary importance to explain species diversity, while the mass effect surpassed species sorting in the highly connected region. Increasing connectivity in communities containing high habitat heterogeneity can lead to the homogenization of local communities, and consequently, to lower regional diversity, while local species richness was unrelated to the habitat connectivity. Of all species found, Anthriscus sylvestris, Phalaris arundinacea, and Phleum pretense significantly responded to connectivity, and showed high abundance in the highly connected region. We suggest that these species may play a role in switching the force from local resources to regional connectivity shaping the community structure. On the landscape context level, the different responses of local species richness and evenness to landscape context were investigated. Seven landscape structural parameters served to indicate landscape context on five scales. On all scales but the smallest scales, the Shannon-Wiener diversity of land covers (H') correlated positively with the local richness. The factor (H') showed the highest correlation coefficients in species richness on the second largest scale. The edge density of arable field was the only predictor that correlated with species evenness on all scales, which showed the highest predictive power on the second smallest scale. The different predictive power of the factors on different scales showed a scaledependent relationship between the landscape context and local plant species diversity, and indicated that different ecological processes determine species richness and evenness. The local richness of species depends on a regional process on large scales, which may relate to the regional species pool, while species evenness depends on a fine- or coarse-grained farming system, which may relate to the patch quality of the habitats of field edges near the buffer strips. My results suggested some guidelines of species diversity conservation in the agricultural ecosystem. To maintain a high level of species diversity in the strips, a high level of phosphorus in strip soil should be avoided. Widening the strips is the most effective mean to improve species richness. Habitat connectivity is not always favorable to species diversity because increasing connectivity in communities containing high habitat heterogeneity can lead to the homogenization of local communities (beta diversity) and, consequently, to lower regional diversity. Overall, a synthesis of local and regional factors emerged as the model that best explain variations in plant species diversity. The studies also suggest that the effects of determinants on species diversity have a complex relationship with scale.

Chemical and source characterisation of size-segregated urban air particulate matter

Recent epidemiological studies have shown a consistent association of the mass concentration of urban air thoracic (PM10) and fine (PM2.5) particles with mortality and morbidity among cardiorespiratory patients. However, the chemical characteristics of different particulate size ranges and the biological mechanisms responsible for these adverse health effects are not well known. The principal aims of this thesis were to validate a high volume cascade impactor (HVCI) for the collection of particulate matter for physicochemical and toxicological studies, and to make an in-depth chemical and source characterisation of samples collected during different pollution situations. The particulate samples were collected with the HVCI, virtual impactors and a Berner low pressure impactor in six European cities: Helsinki, Duisburg, Prague, Amsterdam, Barcelona and Athens. The samples were analysed for particle mass, common ions, total and water-soluble elements as well as elemental and organic carbon. Laboratory calibration and field comparisons indicated that the HVCI can provide a unique large capacity, high efficiency sampling of size-segregated aerosol particles. The cutoff sizes of the recommended HVCI configuration were 2.4, 0.9 and 0.2 μm. The HVCI mass concentrations were in a good agreement with the reference methods, but the chemical composition of especially the fine particulate samples showed some differences. This implies that the chemical characterization of the exposure variable in toxicological studies needs to be done from the same HVCI samples as used in cell and animal studies. The data from parallel, low volume reference samplers provide valuable additional information for chemical mass closure and source assessment. The major components of PM2.5 in the virtual impactor samples were carbonaceous compounds, secondary inorganic ions and sea salt, whereas those of coarse particles (PM2.5-10) were soil-derived compounds, carbonaceous compounds, sea salt and nitrate. The major and minor components together accounted for 77-106% and 77-96% of the gravimetrically-measured masses of fine and coarse particles, respectively. Relatively large differences between sampling campaigns were observed in the organic carbon content of the PM2.5 samples as well as the mineral composition of the PM2.5-10 samples. A source assessment based on chemical tracers suggested clear differences in the dominant sources (e.g. traffic, residential heating with solid fuels, metal industry plants, regional or long-range transport) between the sampling campaigns. In summary, the field campaigns exhibited different profiles with regard to particulate sources, size distribution and chemical composition, thus, providing a highly useful setup for toxicological studies on the size-segregated HVCI samples.

Topographical, structural and geophysical characterization of fracture zones: implications for groundwater flow and vulnerability

The main objective of this study is to evaluate selected geophysical, structural and topographic methods on regional, local, and tunnel and borehole scales, as indicators of the properties of fracture zones or fractures relevant to groundwater flow. Such information serves, for example, groundwater exploration and prediction of the risk of groundwater inflow in underground construction. This study aims to address how the features detected by these methods link to groundwater flow in qualitative and semi-quantitative terms and how well the methods reveal properties of fracturing affecting groundwater flow in the studied sites. The investigated areas are: (1) the Päijänne Tunnel for water-conveyance whose study serves as a verification of structures identified on regional and local scales; (2) the Oitti fuel spill site, to telescope across scales and compare geometries of structural assessment; and (3) Leppävirta, where fracturing and hydrogeological environment have been studied on the scale of a drilled well. The methods applied in this study include: the interpretation of lineaments from topographic data and their comparison with aeromagnetic data; the analysis of geological structures mapped in the Päijänne Tunnel; borehole video surveying; groundwater inflow measurements; groundwater level observations; and information on the tunnel s deterioration as demonstrated by block falls. The study combined geological and geotechnical information on relevant factors governing groundwater inflow into a tunnel and indicators of fracturing, as well as environmental datasets as overlays for spatial analysis using GIS. Geophysical borehole logging and fluid logging were used in Leppävirta to compare the responses of different methods to fracturing and other geological features on the scale of a drilled well. Results from some of the geophysical measurements of boreholes were affected by the large diameter (gamma radiation) or uneven surface (caliper) of these structures. However, different anomalies indicating more fractured upper part of the bedrock traversed by well HN4 in Leppävirta suggest that several methods can be used for detecting fracturing. Fracture trends appear to align similarly on different scales in the zone of the Päijänne Tunnel. For example, similarities of patterns were found between the regional magnetic trends, correlating with orientations of topographic lineaments interpreted as expressions of fracture zones. The same structural orientations as those of the larger structures on local or regional scales were observed in the tunnel, even though a match could not be made in every case. The size and orientation of the observation space (patch of terrain at the surface, tunnel section, or borehole), the characterization method, with its typical sensitivity, and the characteristics of the location, influence the identification of the fracture pattern. Through due consideration of the influence of the sampling geometry and by utilizing complementary fracture characterization methods in tandem, some of the complexities of the relationship between fracturing and groundwater flow can be addressed. The flow connections demonstrated by the response of the groundwater level in monitoring wells to pressure decrease in the tunnel and the transport of MTBE through fractures in bedrock in Oitti, highlight the importance of protecting the tunnel water from a risk of contamination. In general, the largest values of drawdown occurred in monitoring wells closest to the tunnel and/or close to the topographically interpreted fracture zones. It seems that, to some degree, the rate of inflow shows a positive correlation with the level of reinforcement, as both are connected with the fracturing in the bedrock. The following geological features increased the vulnerability of tunnel sections to pollution, especially when several factors affected the same locations: (1) fractured bedrock, particularly with associated groundwater inflow; (2) thin or permeable overburden above fractured rock; (3) a hydraulically conductive layer underneath the surface soil; and (4) a relatively thin bedrock roof above the tunnel. The observed anisotropy of the geological media should ideally be taken into account in the assessment of vulnerability of tunnel sections and eventually for directing protective measures.

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.

Disturbance in boreal spruce forest - immediate dynamics from stand to understorey level

The immediate effects of two human-related vegetation disturbances, (1) green tree retention (GTR) patch felling and scarification by harrowing and (2) experimental understorey vegetation layer removal, were examined in boreal forest stands in Finland. Effects of GTR patch felling and scarification on tree uprootings, on coarse woody debris (CWD) and on epixylic plant community were followed in upland and in paludified forest types. Uprootings increased considerably during 2-3 years after the fellings and were more frequent (47%) in the paludified than in the upland forest (13%). Scarification reduced 68% of the CWD in the felling area. Cover and especially species richness of epixylics declined in the both areas during 1-2 years after the felling. The increasing size of GTR patch correlated positively with the species richness. Regeneration of understorey vegetation community and Vaccinium myrtillus and Vaccinium vitis-idaea after different removals of vegetation layers in an old-growth forest took four years. The regeneration occurred mainly by vegetative means and it was faster in the terms of species richness than in the cover. In the most severe treatment, recovery occurred merely by sexual reproduction. V. myrtillus recovered mainly by producing new shoots. V. vitis-idaea recovered faster than V. myrtillus, mainly by increasing length growth. For ecological reasons, use of larger GTR patches on paludified biotope would be recommendable. In felling areas, scarification by harrowing could be replaced with some other spot-wise method. After moderate intensity level disturbance, recovery occurs rapidly by vegetative regrowth of the dominating species. High level of intensity may prevent the recovery of vegetation community for years, while enabling also the genetic regeneration of the initial species. Local anthropogenic-related disturbances are currently increasing and they can interact during temporally short times, which should be taken in to account in the future forest management plans.