Isotopic fingerprints in surficial waters : Stable isotope methods applied in hydrogeological studies

The driving force behind this study has been the need to develop and apply methods for investigating the hydrogeochemical processes of significance to water management and artificial groundwater recharge. Isotope partitioning of elements in the course of physicochemical processes produces isotopic variations to their natural reservoirs. Tracer property of the stable isotope abundances of oxygen, hydrogen and carbon has been applied to investigate hydrogeological processes in Finland. The work described here has initiated the use of stable isotope methods to achieve a better understanding of these processes in the shallow glacigenic formations of Finland. In addition, the regional precipitation and groundwater records will supplement the data of global precipitation, but as importantly, provide primary background data for hydrological studies. The isotopic composition of oxygen and hydrogen in Finnish groundwaters and atmospheric precipitation was determined in water samples collected during 1995 2005. Prior to this study, no detailed records existed on the spatial or annual variability of the isotopic composition of precipitation or groundwaters in Finland. Groundwaters and precipitation in Finland display a distinct spatial distribution of the isotopic ratios of oxygen and hydrogen. The depletion of the heavier isotopes as a function of increasing latitude is closely related to the local mean surface temperature. No significant differences were observed between the mean annual isotope ratios of oxygen and hydrogen in precipitation and those in local groundwaters. These results suggest that the link between the spatial variability in the isotopic composition of precipitation and local temperature is preserved in groundwaters. Artificial groundwater recharge to glaciogenic sedimentary formations offers many possibilities to apply the isotopic ratios of oxygen, hydrogen and carbon as natural isotopic tracers. In this study the systematics of dissolved carbon have been investigated in two geochemically different glacigenic groundwater formations: a typical esker aquifer at Tuusula, in southern Finland and a carbonate-bearing aquifer with a complex internal structure at Virttaankangas, in southwest Finland. Reducing the concentration of dissolved organic carbon (DOC) in water is a primary challenge in the process of artificial groundwater recharge. The carbon isotope method was used to as a tool to trace the role of redox processes in the decomposition of DOC. At the Tuusula site, artificial recharge leads to a significant decrease in the organic matter content of the infiltrated water. In total, 81% of the initial DOC present in the infiltrated water was removed in three successive stages of subsurface processes. Three distinct processes in the reduction of the DOC content were traced: The decomposition of dissolved organic carbon in the first stage of subsurface flow appeared to be the most significant part in DOC removal, whereas further decrease in DOC has been attributed to adsorption and finally to dilution with local groundwater. Here, isotope methods were used for the first time to quantify the processes of DOC removal in an artificial groundwater recharge. Groundwaters in the Virttaankangas aquifer are characterized by high pH values exceeding 9, which are exceptional for shallow aquifers on glaciated crystalline bedrock. The Virttaankangas sediments were discovered to contain trace amounts of fine grained, dispersed calcite, which has a high tendency to increase the pH of local groundwaters. Understanding the origin of the unusual geochemistry of the Virttaankangas groundwaters is an important issue for constraining the operation of the future artificial groundwater plant. The isotope ratios of oxygen and carbon in sedimentary carbonate minerals have been successfully applied to constrain the origin of the dispersed calcite in the Virttaankangas sediments. The isotopic and chemical characteristics of the groundwater in the distinct units of aquifer were observed to vary depending on the aquifer mineralogy, groundwater residence time and the openness of the system to soil CO2. The high pH values of > 9 have been related to dissolution of calcite into groundwater under closed or nearly closed system conditions relative to soil CO2, at a low partial pressure of CO2.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

Locus Bonus : The Relationship of the Roman Villa to its Environment in the Vicinity of Rome

The aim of this study is to examine the relationship of the Roman villa to its environment. The villa was an important feature of the countryside intended both for agricultural production and for leisure. Manuals of Roman agriculture give instructions on how to select a location for an estate. The ideal location was a moderate slope facing east or south in a healthy area and good neighborhood, near good water resources and fertile soils. A road or a navigable river or the sea was needed for transportation of produce. A market for selling the produce, a town or a village, should have been nearby. The research area is the surroundings of the city of Rome, a key area for the development of the villa. The materials used consist of archaeological settlement sites, literary and epigraphical evidence as well as environmental data. The sites include all settlement sites from the 7th century BC to 5th century AD to examine changes in the tradition of site selection. Geographical Information Systems were used to analyze the data. Six aspects of location were examined: geology, soils, water resources, terrain, visibility/viewability and relationship to roads and habitation centers. Geology was important for finding building materials and the large villas from the 2nd century BC onwards are close to sources of building stones. Fertile soils were sought even in the period of the densest settlement. The area is rich in water, both rainfall and groundwater, and finding a water supply was fairly easy. A certain kind of terrain was sought over very long periods: a small spur or ridge shoulder facing preferably south with an open area in front of the site. The most popular villa resorts are located on the slopes visible from almost the entire Roman region. A visible villa served the social and political aspirations of the owner, whereas being in the villa created a sense of privacy. The area has a very dense road network ensuring good connectivity from almost anywhere in the region. The best visibility/viewability, dense settlement and most burials by roads coincide, creating a good neighborhood. The locations featuring the most qualities cover nearly a quarter of the area and more than half of the settlement sites are located in them. The ideal location was based on centuries of practical experience and rationalized by the literary tradition.

Phosphorus retention in forest soils and the functioning of buffer zones used in forestry

Phosphorus (P) retention properties of soils typical for boreal forest, i.e. podzolic soil and peat soils, vary significantly, but the range of this variation has not been sufficiently documented. To assess the usefulness of buffer zones used in forestry in removing P from the discharge by chemical sorption in soil, and to estimate the risk of P leaching after forestry operations, more data is needed on soil P retention properties. P retention properties of soils were studied at clear-cut areas, unharvested buffer zones adjoining the clear-cut and at peatland buffer zone areas. Desorption-sorption isotherms were determined for the humus layer, the mineral soil horizons E, B and C of the Podzol profile and for the surface layer peat (0-15 cm) and the subsurface layer peat (15-30 cm). The efficiency of buffer zones in retaining P was studied at six peatland buffer zone areas by adding P-containing solute in the inflow. A tracer study was conducted at one of the buffer zone areas to determine the allocation of the added P in soil and vegetation. Measured sorption or desorption rather than parameter values of fitted sorption equations described P desorption and sorption behaviour in soil. The highest P retention efficiency was in the B horizon and consequently, if contact occurred or was established between the soluble P in the water and the soil B horizon, the risk of P leaching was low. Humus layer was completely incapable of retaining P after clear-cutting. In the buffer zones, the decrease in P retention properties in the humus layer and the low amount of P sorbed by it indicated that the importance of the layer in the functioning of buffer zones is low. The peatland buffer zone areas were efficient in retaining soluble P from inflow. P sorption properties of the peat soil at the buffer zone areas varied largely but the contribution of P sorption in the peat was particularly important during high flow in spring, when the vegetation was not fully developed. Factors contributing to efficient P retention were large buffer size and low hydrological load whereas high hydrological load combined with the formation of preferential flow paths, especially during early spring or late autumn was disadvantageous. However, small buffer zone areas, too, may be efficient in reducing P load.

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.

Phosphorus in agricultural soils of Finland - characterization of reserves and retention in mineral soil profiles

An overwhelming majority of all the research on soil phosphorus (P) has been carried out with soil samples taken from the surface soils only, and our understanding of the forms and the reactions of P at a soil profile scale is based on few observations. In Finland, the interest in studying the P in complete soil profiles has been particularly small because of the lack of tradition in studying soil genesis, morphology, or classification. In this thesis, the P reserves and the retention of orthophosphate phosphorus (PO4-P) were examined in four cultivated mineral soil profiles in Finland (three Inceptisols and one Spodosol). The soils were classified according to the U.S. Soil Taxonomy and soil samples were taken from the genetic horizons in the profiles. The samples were analyzed for total P concentration, Chang and Jackson P fractions, P sorption properties, concentrations of water-extractable P, and for concentrations of oxalate-extractable Al and Fe. Theoretical P sorption capacities and degrees of P saturation were calculated with the data from the oxalate-extractions and the P fractionations. The studied profiles can be divided into sections with clearly differing P characteristics by their master horizons Ap, B and C. The C (or transitional BC) horizons below an approximate depth of 70 cm were dominated by, assumingly apatitic, H2SO4-soluble P. The concentration of total P in the C horizons ranged from 729 to 810 mg kg-1. In the B horizons between the depths of 30 and 70 cm, a significant part of the primary acid-soluble P has been weathered and transformed to secondary P forms. A mean weathering rate of the primary P in the soils was estimated to vary between 230 and 290 g ha-1 year-1. The degrees of P saturation in the B and C horizons were smaller than 7%, and the solubility of PO4-P was negligible. The P conditions in the Ap horizons differed drastically from those in the subsurface horizons. The high concentrations of total P (689-1870 mg kg-1) in the Ap horizons are most likely attributable to long-term cultivation with positive P balances. A significant proportion of the P in the Ap horizons occurred in the NH4F- and NaOH-extractable forms and as organic P. These three P pools, together with the concentrations of oxalate-extractable Al and Fe, seem to control the dynamics of PO4-P in the soils. The degrees of P saturation in the Ap horizons were greater (8-36%) than in the subsurface horizons. This was also reflected in the sorption experiments: Only the Ap horizons were able to maintain elevated PO4-P concentrations in the solution phase − all the subsoil horizons acted as sinks for PO4-P. Most of the available sorption capacity in the soils is located in the B horizons. The results suggest that this capacity could be utilized in reducing the losses of soluble P from excessively fertilized soils by mixing highly sorptive material from the B horizons with the P-enriched surface soil. The drastic differences in the P characteristics observed between adjoining horizons have to be taken into consideration when conducting soil sampling. Sampling of subsoils has to be made according to the genetic horizons or at small depth increments. Otherwise, contrasting materials are likely to be mixed in the same sample; and the results of such samples are not representative of any material present in the studied profile. Air-drying of soil samples was found to alter the results of the sorption experiments and the water extractions. This indicates that the studies on the most labile P forms in soil should be carried out with moist samples.

Organic contaminants - occurrence and biological effects in the Baltic Sea

The Baltic Sea was studied with respect to selected organic contaminants and their ecotoxicology. The research consisted of analyses of total hydrocarbons, polycyclic aromatic hydrocarbons, bile metabolites, hepatic ethoxyresorufin-O-deethylase (EROD) activity, polychlorinated biphenyls (PCBs) and organochlorine pesticides (OCPs). The contaminants were measured from various matrices, such as seawater, sediment and biota. The methods of analysis were evaluated and refined to comparability of the results. Polyaromatic hydrocarbons, originating from petroleum, are known to be among the most harmful substances to the marine environment. In Baltic subsurface water, seasonal dependence of the total hydrocarbon concentrations (THCs) was seen. Although concentrations of parent polycyclic aromatic hydrocarbons (PAHs) in sediment surface varied between 64 and 5161 ug kg-1 (dw), concentrations above 860 ug kg-1 (dw) were found in all the studied sub-basins of the Baltic Sea. Concentrations commonly considered to substantially increase the risk of liver disease and reproductive impairment in fish, as well as potential effects on growth (above 1000 ug kg-1 dw), were found in all the studied sub-basins of the Baltic Sea except Kattegat. Thus, considerable pollution in sediments was indicated. In bivalves, the sums of 12 PAHs varied on a wet weight basis between 44 and 298 ug kg-1 (ww). The predominant PAHs were high molecular weight and the PAH profiles of M. balthica differed from those found in sediment from the same area. The PAHs were both pyrolytic and petrogenic in origin, and a contribution from diesel engines was found, which indicates pollution of the Baltic Sea, most likely caused by the steadily increasing shipping in the area. The HPLC methods developed for hepatic EROD activity and bile metabolite measurements proved to be fast and suitable for the study of biological effects. A mixed function oxygenase enzyme system in Baltic Sea perch collected from the Gulf of Finland was induced slightly: EROD activity in perch varied from 0.30 14 pmol min-1 mg-1 protein. This range can be considered to be comparable to background values. Recent PAH exposure was also indicated by enhanced levels (213 and 1149 ug kg-1) of the bile metabolite 1-hydroxypyrene. No correlation was indicated between hepatic EROD activity and concentration of 1-hydroxypyrene in bile. PCBs and OCPs were observed in Baltic Sea sediment, bivalves and herring. Sums of seven CBs in surface sediment (0 5 cm) ranged from 0.04 to 6.2 ug kg-1 (dw) and sums of three DDTs from 0.13 to 5.0 ug kg-1 (dw). The highest levels of contaminants were found in the most eastern area of the Gulf of Finland where the highest total carbon and nitrogen content was found and where the lowest percentage proportion of p,p -DDT was found. The highest concentrations of CBs and the lowest concentration of DDTs were found in M. balthica from the Gulf of Finland. The highest levels of DDTs were found in M. balthica from the Hanö Bight, which is the outer part of the Bornholm Basin close to the Swedish mainland. In bivalves, the sums of seven CBs were 72 108 ug kg-1 (lw) and the sums of three DDTs were 66 139 ug kg-1 (lw). Results from temporal trend monitoring showed, that during the period 1985 2002, the concentrations of seven CBs in two-year-old female Baltic herring were clearly decreased, from 9 16 to 2 6 ug kg-1 (ww) in the northern Baltic Sea. At the same time, concentrations of three DDTs declined from 8 15 to 1 5 ug kg-1 (ww). The total concentration of the fat-soluble CBs and DDTs in Baltic herring muscle was shown to be age-dependent; the average concentrations in ten-year-old Baltic herring were three to five-fold higher than in two-year-old herring. In Baltic herring and bivalves, as well as in surface sediments, CB 138 and CB153 were predominant among CBs, whereas among DDTs p,p'-DDD predominated in sediment and p,p'-DDE in bivalves and Baltic herring muscle. Baltic Sea sediments are potential sources of contaminants that may become available for bioaccumulation. Based on ecotoxicological assessment criteria, cause for concern regarding CBs in sediments was indicated for the Gulf of Finland and the northern Baltic Proper, and for the northern Baltic Sea regarding CBs in Baltic herring more than two years old. Statistical classification of selected organic contaminants indicated high-level contamination for p,p'-DDT, p,p'-DDD, p,p'-DDE, total DDTs, HCB, CB118 and CB153 in muscle of Baltic herring in age groups two to ten years; in contrast, concentrations of a-HCH and g-HCH were found to be moderate. The concentrations of DDTs and CBs in bivalves is sufficient to cause biological effects, and demonstrates that long-term biological effects are still possible in the case of DDTs in the Hanö Bight.

Spatial and temporal variability in midge (Nematocera) assemblages in shallow Finnish lakes (60−70 °N) : community-based modelling of past environmental change

Multi- and intralake datasets of fossil midge assemblages in surface sediments of small shallow lakes in Finland were studied to determine the most important environmental factors explaining trends in midge distribution and abundance. The aim was to develop palaeoenvironmental calibration models for the most important environmental variables for the purpose of reconstructing past environmental conditions. The developed models were applied to three high-resolution fossil midge stratigraphies from southern and eastern Finland to interpret environmental variability over the past 2000 years, with special focus on the Medieval Climate Anomaly (MCA), the Little Ice Age (LIA) and recent anthropogenic changes. The midge-based results were compared with physical properties of the sediment, historical evidence and environmental reconstructions based on diatoms (Bacillariophyta), cladocerans (Crustacea: Cladocera) and tree rings. The results showed that the most important environmental factor controlling midge distribution and abundance along a latitudinal gradient in Finland was the mean July air temperature (TJul). However, when the dataset was environmentally screened to include only pristine lakes, water depth at the sampling site became more important. Furthermore, when the dataset was geographically scaled to southern Finland, hypolimnetic oxygen conditions became the dominant environmental factor. The results from an intralake dataset from eastern Finland showed that the most important environmental factors controlling midge distribution within a lake basin were river contribution, water depth and submerged vegetation patterns. In addition, the results of the intralake dataset showed that the fossil midge assemblages represent fauna that lived in close proximity to the sampling sites, thus enabling the exploration of within-lake gradients in midge assemblages. Importantly, this within-lake heterogeneity in midge assemblages may have effects on midge-based temperature estimations, because samples taken from the deepest point of a lake basin may infer considerably colder temperatures than expected, as shown by the present test results. Therefore, it is suggested here that the samples in fossil midge studies involving shallow boreal lakes should be taken from the sublittoral, where the assemblages are most representative of the whole lake fauna. Transfer functions between midge assemblages and the environmental forcing factors that were significantly related with the assemblages, including mean air TJul, water depth, hypolimnetic oxygen, stream flow and distance to littoral vegetation, were developed using weighted averaging (WA) and weighted averaging-partial least squares (WA-PLS) techniques, which outperformed all the other tested numerical approaches. Application of the models in downcore studies showed mostly consistent trends. Based on the present results, which agreed with previous studies and historical evidence, the Medieval Climate Anomaly between ca. 800 and 1300 AD in eastern Finland was characterized by warm temperature conditions and dry summers, but probably humid winters. The Little Ice Age (LIA) prevailed in southern Finland from ca. 1550 to 1850 AD, with the coldest conditions occurring at ca. 1700 AD, whereas in eastern Finland the cold conditions prevailed over a longer time period, from ca. 1300 until 1900 AD. The recent climatic warming was clearly represented in all of the temperature reconstructions. In the terms of long-term climatology, the present results provide support for the concept that the North Atlantic Oscillation (NAO) index has a positive correlation with winter precipitation and annual temperature and a negative correlation with summer precipitation in eastern Finland. In general, the results indicate a relatively warm climate with dry summers but snowy winters during the MCA and a cool climate with rainy summers and dry winters during the LIA. The results of the present reconstructions and the forthcoming applications of the models can be used in assessments of long-term environmental dynamics to refine the understanding of past environmental reference conditions and natural variability required by environmental scientists, ecologists and policy makers to make decisions concerning the presently occurring global, regional and local changes. The developed midge-based models for temperature, hypolimnetic oxygen, water depth, littoral vegetation shift and stream flow, presented in this thesis, are open for scientific use on request.

Eurasian Neogene large herbivorous mammals and climate

The development and changes in the distribution of herbivorous mammal communities during the Neogene is complex. The Eurasian scale environmental patterns reflect the large scale geographical and climatic patterns. The reorganization of these affect the biome distribution throughout the continent. The distribution of mammal taxa was closely associated with the distribution of biomes. In Eurasia the Neogene development of environments was twofold. The early and middle Miocene that seemed to have been advantageous for mammals was followed by drying of environments during the late Neogene. The mid-latitude drying was the main trend, and it is the combined result of the retreat of Paratethys, the uplift of Tibetan Plateau and changes in the ocean currents and temperatures. The common mammals were "driving" the evolution of mammalian communities. During the late Miocene we see the drying affecting more and more regions, and we see changes in the composition of mammalian communities.

Sexual reproduction in chydorid cladocerans (Anomopoda, Chydoridae) in southern Finland - implications for paleolimnology

The relationship between sexual reproduction of littoral chydorid cladocerans (Anomopoda, Chydoridae) and environmental factors in aquatic ecosystems has been rarely studied, although the sexual behavior of some planktonic cladocerans is well documented. Ecological monitoring was used to study the relationship between climate-related and non-climatic environmental factors and chydorid sexual reproduction patterns in nine environmentally different lakes that were closely situated to each other in southern Finland. Furthermore, paleolimnological ephippium analysis was used to clarify how current sexual reproduction is reflected in surface sediments of the same nine lakes. Additionally, short sediment cores from two of the lakes were studied with ephippium analysis to examine how recent climate-related and non-climatic environmental changes were reflected in chydorid sexual reproduction. Ephippium analysis uses the subfossil shells of asexual individuals to represent asexual reproduction and the shells of sexual females, i.e. ephippia, to represent sexual reproduction. The relative proportion of ephippia of all chydorid species, i.e. total chydorid ephippia (TCE) indicates the relative proportion of sexual reproduction during the open-water season. This thesis is part of the EPHIPPIUM-project which aims to develop ephippium analysis towards a quantitative climate reconstruction tool. To be able to develop a valid climate model, the influence of the environmental stressors other than climate on contemporary sexual reproduction and its reflection in sediment assemblages must be clarified so they can be eliminated from the model. During contemporary monitoring a few sexual individuals were observed during summer, apparently forced to sexual reproduction by non-climatic local environmental factors, such as crowding or invertebrate predation. Monitoring also revealed that the autumnal chydorid sexual reproduction period was consistent between the different lakes and climate-related factors appeared to act as the main inducers and regulators of autumnal sexual reproduction. However, during autumn, chydorid species and populations among the lakes exhibited a wide variation in the intensity, induction time, and length of autumnal sexual reproduction. These variations apparently act as mechanisms for local adaptations due to the genetic variability provided by sexual reproduction that enhance the ecological flexibility of chydorid species, allowing them to inhabit a wide range of environments. A large variation was also detected in the abundance of parthenogenetic and gamogenetic individuals during the open-water season among the lakes. On the basis of surface sediment samples, the general level of the TCE is ca. 3-4% in southern Finland, reflecting an average proportion of sexual reproduction in this specific climate. The variation in the TCE was much lower than could be expected on the basis of the monitoring results. This suggests that some of the variation detected by monitoring may derive from differences between sampling sites and years smoothed out in the sediment samples, providing an average of the entire lake area and several years. The TCE is always connected to various ecological interactions in lake ecosystems and therefore is always lake-specific. Hypothetically, deterioration of climate conditions can be detected in the TCE as an increase in ephippia of all chydorid species, since a shortening open-water season is reflected in the relative proportions of the two reproduction modes. Such an increase was clearly detected for the time period of the Little Ice Age in a sediment core. The paleolimnological results also indicated that TCE can suddenly increase due to ephippia of one or two species, which suggests that at least some chydorids can somehow increase the production of resting eggs under local environmental stress. Thus, some environmental factors may act as species-specific environmental stressors. The actual mechanism of the increased sexual reproduction seen in sediments has been unknown but the present study suggests that the mechanism is probably the increased intensity of gamogenesis, i.e. that a larger proportion of individuals in autumnal populations reproduce sexually, which results in a larger proportion of ephippia in sediments and a higher TCE. The results of this thesis demonstrate the utility of ephippium analysis as a paleoclimatological method which may also detect paleolimnological changes by identifying species-specific environmental stressors. For a quantitative TCE-based climate reconstruction model, the natural variation in the TCE of surface sediments in different climates must be clarified with more extensive studies. In addition, it is important to recognize the lakes where the TCE is not only a reflection of the length of the open-water season, but is also non-climatically forced. The results of ephippium analysis should always be interpreted in a lake-specific manner and in the context of other paleoecological proxies.