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.

Magmatic evolution of topaz-bearing granite stocks within the Wiborg rapakivi granite batholith

This study brings new insights into the magmatic evolution of natural F-enriched peraluminous granitic systems. The Artjärvi, Sääskjärvi and Kymi granite stocks within the 1.64 Ga Wiborg rapakivi granite batholith have been investigated by petrographic, geochemical, experimental and melt inclusion methods. These stocks represent late-stage leucocratic and weakly peraluminous intrusive phases typical of rapakivi granites worldwide. The Artjärvi and Sääskjärvi stocks are multiphase intrusions in which the most evolved phase is topaz granite. The Kymi stock contains topaz throughout and has a well-developed zoned structure, from the rim to the center: stockscheider pegmatite equigranular topaz granite porphyritic topaz granite. Geochemically the topaz granites are enriched in F, Li, Be, Ga, Rb, Sn and Nb and depleted in Mg, Fe, Ti, Ba, Sr, Zr and Eu. The anomalous geochemistry and mineralogy of the topaz granites are essentially magmatic in origin; postmagmatic reactions have only slightly modified the compositions. The Kymi equigranular topaz granite shows the most evolved character, and the topaz granites at Artjärvi and Sääskjärvi resemble the less evolved porphyritic topaz granite of the Kymi stock. Stockscheiders are found at the roof contacts of the Artjärvi and Kymi stocks. The stockscheider at Artjärvi is composed of biotite-rich schlieren and pegmatite layers parallel to the contact. The schlieren layering is considered to have formed by velocity-gradient sorting mechanism parallel to the flow, which led to the accumulation of mafic minerals along the upper contact of the topaz granite. Cooling and contraction of the topaz granite formed fractures parallel to the roof contact and residual pegmatite magmas were injected along the fractures and formed the pegmatite layers. The zoned structure of the Kymi stock is the result of intrusion of highly evolved residual melt from deeper parts of the magma chamber along the fractured contact between the porphyritic granite crystal mush and country rock. The equigranular topaz granite and marginal pegmatite (stockscheider) crystallized from this evolved melt. Phase relations of the Kymi equigranular topaz granite have been investigated utilizing crystallization experiments at 100 to 500 MPa as a function of water activity and F content. Fluorite and topaz can crystallize as liquidus phases in F-rich peraluminous systems, but the F content of the melt should exceed 2.5 - 3.0 wt % to facilitate crystallization of topaz. In peraluminous F-bearing melts containing more than 1 wt % F, topaz and muscovite are expected to be the first F-bearing phases to crystallize at high pressure, whereas fluorite and topaz should crystallize first at low pressure. Overall, the saturation of fluorite and topaz follows the reaction: CaAl2Si2O8 (plagioclase) + 2[AlF3]melt = CaF2 (fluorite) + 2Al2SiO4F2 (topaz). The obtained partition coefficient for F between biotite and glass D(F)Bt/glass is 1.89 to 0.80 (average 1.29) and can be used as an empirical fluormeter to determine the F content of coexisting melts. In order to study the magmatic evolution of the Kymi stock, crystallized melt inclusions in quartz and topaz grains in the porphyritic and the equigranular topaz granites and the marginal pegmatite were rehomogenized and analyzed. The homogenization conditions for the melt inclusions from the granites were 700 °C, 300 MPa, and 24 h, and for melt inclusions from the pegmatite, 700 °C, 100 MPa, and 24/96 h. The majority of the melt inclusions is chemically similar to the bulk rocks (excluding H2O content), but a few melt inclusions in the equigranular granite show clearly higher F and low K2O contents (on average 11.6 wt % F, 0.65 wt % K2O). The melt inclusion compositions indicate coexistence of two melt fractions, a prevailing peraluminous and a very volatile-rich, possibly peralkaline. Combined petrological, experimental and melt inclusion studies of the Kymi equigranular topaz granite indicate that plagioclase was the liquidus phase at nearly water-saturated (fluid-saturated) conditions and that the F content of the melt was at least 2 wt %. The early crystallization of biotite and the presence of muscovite in crystallization experiments at 200 MPa contrasts with the late-stage crystallization of biotite and the absence of muscovite in the equigranular granite, indicating that crystallization pressure may have been lower than 200 MPa for the granite.

Features of the geomagnetic field during the Holocene and Proterozoic

The geomagnetic field is one of the most fundamental geophysical properties of the Earth and has significantly contributed to our understanding of the internal structure of the Earth and its evolution. Paleomagnetic and paleointensity data have been crucial in shaping concepts like continental drift, magnetic reversals, as well as estimating the time when the Earth's core and associated geodynamo processes begun. The work of this dissertation is based on reliable Proterozoic and Holocene geomagnetic field intensity data obtained from rocks and archeological artifacts. New archeomagnetic field intensity results are presented for Finland, Estonia, Bulgaria, Italy and Switzerland. The data were obtained using sophisticated laboratory setups as well as various reliability checks and corrections. Inter-laboratory comparisons between three laboratories (Helsinki, Sofia and Liverpool) were performed in order to check the reliability of different paleointensity methods. The new intensity results fill up considerable gaps in the master curves for each region investigated. In order to interpret the paleointensity data of the Holocene period, a novel and user-friendly database (GEOMAGIA50) was constructed. This provided a new tool to independently test the reliability of various techniques and materials used in paleointensity determinations. The results show that archeological artifacts, if well fired, are the most suitable materials. Also lavas yield reliable paleointensity results, although they appear more scattered. This study also shows that reliable estimates are obtained using the Thellier methodology (and its modifications) with reliability checks. Global paleointensity curves during Paleozoic and Proterozoic have several time gaps with few or no intensity data. To define the global intensity behavior of the Earth's magnetic field during these times new rock types (meteorite impact rocks) were investigated. Two case histories are presented. The Ilyinets (Ukraine) impact melt rocks yielded a reliable paleointensity value at 440 Ma (Silurian), whereas the results from Jänisjärvi impact melts (Russian Karelia, ca. 700 Ma) might be biased towards high intensity values because of non-ideal magnetic mineralogy. The features of the geomagnetic field at 1.1 Ga are not well defined due to problems related to reversal asymmetries observed in Keweenawan data of the Lake Superior region. In this work new paleomagnetic, paleosecular variation and paleointensity results are reported from coeval diabases from Central Arizona and help understanding the asymmetry. The results confirm the earlier preliminary observations that the asymmetry is larger in Arizona than in Lake Superior area. Two of the mechanisms proposed to explain the asymmetry remain plausible: the plate motion and the non-dipole influence.

Physical properties of a boreal clay soil under differently managed perennial vegetation

The physical properties of surface soil horizons, essentially pore size, shape, continuity and affinity for water, regulate water entry into the soil. These properties are prone to changes caused by natural forces and human activity. The hydraulic properties of the surface soil greatly impact the generation of surface runoff and accompanied erosion, the major concern of agricultural water protection. The general target of this thesis was to improve our understanding of the structural and hydraulic properties of boreal clay soils. Physical properties of a clayey surface soil (0 - 10 cm, clay content 51%), with a micaceous/illitic mineralogy subjected to three different management practices of perennial vegetation, were studied. The study sites were vegetated buffer zones located side by side in SW Finland: 1) natural vegetation with no management, 2) harvested once a year, and 3) grazed by cattle. The soil structure, hydraulic properties, shrinkage properties and soil water repellency were determined at all sites. Two distinct flow domains were evident. The surface soil was characterized by subangular blocky, angular blocky and platy aggregates. Hence, large, partially accommodated, irregular elongated pores dominated the macropore domain at all sites. The intra-aggregate pore system was mostly comprised of pores smaller than 30 μm, which are responsible for water storage. Macropores at the grazed site, compacted by hoof pressure, were horizontally oriented and pore connectivity was poorest, which decreased water and air flux compared with other sites. Drying of the soil greatly altered its structure. The decrease in soil volume between wet and dry soil was 7 - 10%, most of which occurred in the moisture range of field conditions. Structural changes, including irreversible collapse of interaggregate pores, began at matric potentials around -6 kPa indicating, instability of soil structure against increasing hydraulic stress. Water saturation and several freezethaw cycles between autumn and spring likely weakened the soil structure. Soil water repellency was observed at all sites at the time of sampling and when soil was dryer than about 40 vol.%. (matric potential < -6 kPa). Therefore, water repellency contributes to water flow over a wide moisture range. Water repellency was also observed in soils with low organic carbon content (< 2%), which suggests that this phenomenon is common in agricultural soils of Finland due to their relatively high organic carbon content. Aggregate-related pedofeatures of dense infillings described as clay intrusions were found at all sites. The formation of these intrusions was attributed to clay dispersion and/or translocation during spring thaw and drying of the suspension in situ. These processes generate very new aggregates whose physical properties are most probably different from those of the bulk soil aggregates. Formation of the clay infillings suggested that prolonged wetness in autumn and spring impairs soil structure due to clay dispersion, while on the other hand it contributes to the pedogenesis of the soil. The results emphasize the dynamic nature of the physical properties of clay soils, essentially driven by their moisture state. In a dry soil, fast preferential flow is favoured by abundant macropores including shrinkage cracks and is further enhanced by water repellency. Increase in soil moisture reduces water repellency, and swelling of accommodated pores lowers the saturated hydraulic conductivity. Moisture- and temperature-related processes significantly alter soil structure over a time span of 1 yr. Thus, the pore characteristics as well as the hydraulic properties of soil are time-dependent.

Suomussalmen varhaisproterotsooisten diabaasien mineralogia tektonisessa tulkinnassa

The object of this research is to study the mineralogy of the diabase dykes in Suomussalmi and the relevance of the mineralogy to tectonic events, specifically large block movements in the Archaean crust. Sharp tectonic lines separate two anomalies in the dyke swarms, shown on a geomagnetic map as positive anomalies. In one of these areas, the Toravaara anomaly, the diabases seem to contain pyroxenes as a main component. Outside the Toravaara anomaly hornblende is the main ferromagnesian mineral in diabases. The aim of this paper is to research the differences in the diabases inside and outside the anomalies and interpret the processes that formed the anomalies. The data for this sudy consist of field observations, 120 thin sections, 334 electron microprobe analyses, 19 whole-rock chemical analyses, a U-Pb age analysis and geomagnetic low-altitude aerial survey maps. The methods are interpretation of field observations, chemical analyses, microprobe analyses of single minerals and radiometric age determination, microscopic studies of the thin sections, geothermometers and geobarometers. On the basis of field observations and petrographic studies the diabases in the area are divided into pyroxene diabases, hornblende diabases and the Lohisärkkä porphyritic dyke swarm. Hornblende diabases are found in the entire study area, while the pyroxene diabases concentrate on the area of the Toravaara geomagnetic anomaly. The Lohisärkkä swarm transects the whole area as a thin line from east to west. The diabases are fairly homogenous both chemically and by mineral composition. The few exceptions are part of rarer older swarms or are significantly altered. The Lohisärkkä dyke swarm was dated as 2,21 Ga old, significantly older than the most common 1,98 Ga swarm in the area. The geothermometers applied showed that the diabases on the Toravaara anomaly were stabilized at a much higher temperature than the dykes outside the anomaly. The geobarometers showed the pyroxenes to have crystallized at varying depths. The research showed the Toravaara anomaly to have formed by a vertical block movement, and the fault on its west side to have a total lateral transfer of only a few kilometers. The formation of the second anomaly was also interpreted to be tectonic in nature. In addition, the results of the geothermobarometry uncovered necessary conditions for the study of diabase emplacement depth: the minerals for the study must be chosen by minimum crystallization depth, and a geobarometer capable of determining the magmatic temperature must be used. In addition, it would be more suitable to conduct this kind of study in an area where the dykes are more exposed.