Suomen maankuoren kivilajikoostumuksen tulkinta seismisistä aineistoista

Aim of this study is to investigate composition of the crust in Finland using seismic wide-angle velocity models and laboratory measurements on P- and S-wave velocities of different rock types. The velocities adopted from wide-angle velocity models were compared with laboratory velocities of different rock types corrected for the crustal PT conditions in the study area. The wide-angle velocity models indicate that the P-wave velocity does not only increase step-wise at boundaries of major crustal layers, but there is also gradual increase of velocity within the layers. On the other hand, the laboratory measurements of velocities indicate that no single rock type is able to provide the gradual downward increasing trends. Thus, there must be gradual vertical changes in rock composition. The downward increase of velocities indicates that the composition of the crust becomes gradually more mafic with increasing depth. Even though single rock types cannot simulate the wide-angle model velocities, it can be done with a mixture of rock types. There are a large number of rock type mixtures giving the correct P-wave velocities. Therefore, the inverse solution of rock types and their proportions from velocities is a non-unique problem if only P-wave velocities is available. Amount of the possible rock type mixtures can be limitted using S-wave velocities, reflection seismic results and other geological and geophysical results of the study area. Crustal model FINMIX-2 is presented in this study and it suggest that the crustal velocity profiles can be simulated with rock type mixtures, where the upper crust consists of felsic gneisses and granitic-granodioritic rocks with a minor contribution of quartzite, amphibolite and diabase. In the middle crust the amphibolite proportion increases. The lower crust consists of tonalitic gneiss, mafic garnet granulite, hornblendite, pyroxenite and minor mafic eclogite. This composition model is in agreement with deep crustal kimberlite-hosted xenolith data in eastern Finland and reflectivity of the FIRE (Finnish Reflection Experiment). According to FINMIX-2 model the Moho is deeper and the crustal composition is a more mafic than an average global continental model would suggest. Composition models of southern Finland are quite similar than FINMIX-2 model. However, there are minor differencies between the models, which indicates areal differences of composition. Models of northern Finland shows that the crustal thickness is smaller than southern Finland and composition of the upper crust is different. Density profiles calculated from the lithological models suggest that there is practically no density contrast at Moho in areas of the high-velocity lower crust. This implies that crustal thickness in the central Fennoscandian Shield may have been controlled by the densities of the lower crustal and upper mantle rocks.

Archaean crustal evolution in eastern Finland : New geochronological and geochemical constraints from the Kuhmo terrain and the Nurmes belt

This thesis summarises the results of four original papers concerning U-Pb geochronology and geochemical evolution of Archaean rocks from the Kuhmo terrain and the Nurmes belt, eastern Finland. The study area belongs to a typical Archaean granite-greenstone terrain, composed of metavolcanic and metasedimentary rocks in generally N-S trending greenstone belts as well as a granitoid-gneiss complex with intervening gneissic and migmatised supracrustal and plutonic rocks. U-Pb data on migmatite mesosomes indicate that the crust surrounding the Tipasjärvi-Kuhmo-Suomussalmi greenstone belt is of varying age. The oldest protolith detected for a migmatite mesosome from the granitoid-gneiss complex is 2.94 Ga, whereas the other dated migmatites protoliths have ages of 2.84 2.79 Ga. The latter protoliths are syngenetic with the majority of volcanic rocks in the adjacent Tipasjärvi-Kuhmo-Suomussalmi greenstone belt. This suggests that the genesis of some of the volcanic rocks within the greenstone belt and surrounding migmatite protoliths could be linked. Metamorphic zircon overgrowths with ages of 2.84 2.81 Ga were also obtained. The non-migmatised plutonic rocks in the Kuhmo terrain and in the Nurmes belt record secular geochemical evolution, typical of Archaean cratons. The studied tonalitic rocks have ages of 2.83 2.75 Ga and they have geochemical characteristics similar to low-Al and high-Al TTD (tonalite-trondhjemite-dacite). The granodiorites, diorites, and gabbros with high Mg/Fe and LILE-enriched characteristics were mostly emplaced between 2.74 2.70 Ga and they exhibit geochemical characteristics typical of Archaean sanukitoid suites. The latest identified plutonic episode took place at 2.70 2.68 Ga, when compositionally heterogeneous leucocratic granitoid rocks, with a variable crustal component, were emplaced. U-Pb data on migmatite leucosomes suggest that leucosome generation may have been coeval with this latest plutonic event. On the basis of available U-Pb and Sm-Nd isotopic data it appears that the plutonic rocks of the Kuhmo terrain and the Nurmes belt do not contain any significant input from Palaeoarchaean sources. A characteristic feature of the Nurmes belt is the presence of migmatised paragneisses, locally preserving primary edimentary structures, with sporadic amphibolite intercalations. U-Pb studies on zircons indicate that the precursors of the Nurmes paragneisses were graywackes that were deposited between 2.71 Ga and 2.69 Ga and that they had a prominent 2.75 2.70 Ga source. Nd isotopic and whole-rock geochemical data for the intercalated amphibolites imply MORB sources. U-Pb data on zircons from the plutonic rocks and paragneisses reveal that metamorphic zircon growth took place at 2.72 2.63 Ga. This was the last tectonothermal event related to cratonisation of the Archaean crust of eastern Finland.

Snow characteristics in Dronning Maud Land, Antarctica

Snow cover is very sensitive to climate change and has a large feedback effect on the climate system due to the high albedo. Snow covers almost all surfaces in Antarctica and small changes in snow properties can mean large changes in absorbed radiation. In the ongoing discussion of climatic change, the mass balance of Antarctica has received increasing focus during recent decades, since its reaction to global warming strongly influences sea-level change. The aim of the present work was to examine the spatial and temporal variations in the physical and chemical characteristics of surface snow and annual accumulation rates in western Dronning Maud Land, Antarctica. The data were collected along a 350-km-long transect from the coast to the plateau during the years 1999-2004 as a part of the Finnish Antarctic Research Programme (FINNARP). The research focused on the most recent annual accumulation in the coastal area. The results show that the distance from the sea, and the moisture source, was the most predominant factor controlling the variations in both physical (conductivity, grain size, oxygen isotope ratio and accumulation) and chemical snow properties. The sea-salt and sulphur-containing components predominated in the coastal region. The local influences of nunataks and topographic highs were also visible on snow. The variations in all measured properties were wide within single sites mostly due to redistribution by winds and sastrugi topography, which reveals the importance of the spatially representative measurements. The mean accumulations occurred on the ice shelf, in the coastal region and on the plateau: 312 ± 28, 215 ± 43 and 92 ± 25 mm w.e., respectively. Depth hoar layers were usually found under the thin ice crust and were associated with a low dielectric constant and high concentrations of nitrate. Taking into account the vast size of the Antarctic ice sheet and its geographic characteristics, it is important to extend investigation of the distribution of surface snow properties and accumulation to provide well-documented data.

Improving GPS time series for geodynamic studies

Accurate and stable time series of geodetic parameters can be used to help in understanding the dynamic Earth and its response to global change. The Global Positioning System, GPS, has proven to be invaluable in modern geodynamic studies. In Fennoscandia the first GPS networks were set up in 1993. These networks form the basis of the national reference frames in the area, but they also provide long and important time series for crustal deformation studies. These time series can be used, for example, to better constrain the ice history of the last ice age and the Earth s structure, via existing glacial isostatic adjustment models. To improve the accuracy and stability of the GPS time series, the possible nuisance parameters and error sources need to be minimized. We have analysed GPS time series to study two phenomena. First, we study the refraction in the neutral atmosphere of the GPS signal, and, second, we study the surface loading of the crust by environmental factors, namely the non-tidal Baltic Sea, atmospheric load and varying continental water reservoirs. We studied the atmospheric effects on the GPS time series by comparing the standard method to slant delays derived from a regional numerical weather model. We have presented a method for correcting the atmospheric delays at the observational level. The results show that both standard atmosphere modelling and the atmospheric delays derived from a numerical weather model by ray-tracing provide a stable solution. The advantage of the latter is that the number of unknowns used in the computation decreases and thus, the computation may become faster and more robust. The computation can also be done with any processing software that allows the atmospheric correction to be turned off. The crustal deformation due to loading was computed by convolving Green s functions with surface load data, that is to say, global hydrology models, global numerical weather models and a local model for the Baltic Sea. The result was that the loading factors can be seen in the GPS coordinate time series. Reducing the computed deformation from the vertical time series of GPS coordinates reduces the scatter of the time series; however, the long term trends are not influenced. We show that global hydrology models and the local sea surface can explain up to 30% of the GPS time series variation. On the other hand atmospheric loading admittance in the GPS time series is low, and different hydrological surface load models could not be validated in the present study. In order to be used for GPS corrections in the future, both atmospheric loading and hydrological models need further analysis and improvements.

Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Seismic Tomography and Earthquake Mechanism Beneath the Central Fennoscandian Shield

The aim of this thesis was to study the seismic tomography structure of the earth s crust together with earthquake distribution and mechanism beneath the central Fennoscandian Shield, mainly in southern and central Finland. The earthquake foci and some fault plane solutions are correlated with 3-D images of the velocity tomography. The results are discussed in relation to the stress field of the Shield and with other geophysical, e.g. geomagnetic, gravimetric, tectonic, and anisotropy studies of the Shield. The earthquake data of the Fennoscandian Shield has been extracted from the Nordic earthquake parameter data base which was founded at the time of inception of the earthquake catalogue for northern Europe. Eight earlier earthquake source mechanisms are included in a pilot study on creating a novel technique for calculating an earthquake fault plane solution. Altogether, eleven source mechanisms of shallow, weak earthquakes are related in the 3-D tomography model to trace stresses of the crust in southern and central Finland. The earthquakes in the eastern part of the Fennoscandian Shield represent low-active, intraplate seismicity. Earthquake mechanisms with NW-SE oriented horizontal compression confirm that the dominant stress field originates from the ridge-push force in the North Atlantic Ocean. Earthquakes accumulate in coastal areas, in intersections of tectonic lineaments, in main fault zones or are bordered by fault lines. The majority of Fennoscandian earthquakes concentrate on the south-western Shield in southern Norway and Sweden. Onwards, epicentres spread via the ridge of the Shield along the west-coast of the Gulf of Bothnia northwards along the Tornio River - Finnmark fault system to the Barents Sea, and branch out north-eastwards via the Kuusamo region to the White Sea Kola Peninsula faults. The local seismic tomographic method was applied to find the terrane distribution within the central parts of the Shield the Svecofennian Orogen. From 300 local explosions a total of 19765 crustal Pg- and Sg-wave arrival times were inverted to create independent 3-D Vp and Vs tomographic models, from which the Vp/Vs ratio was calculated. The 3-D structure of the crust is presented as a P-wave and for the first time as an S-wave velocity model, and also as a Vp/Vs-ratio model of the SVEKALAPKO area that covers 700x800 km2 in southern and central Finland. Also, some P-wave Moho-reflection data was interpolated to image the relief of the crust-mantle boundary (i.e. Moho). In the tomography model, the seismic velocities vary smoothly. The lateral variations are larger for Vp (dVp =0.7 km/s) than for Vs (dVs =0.4 km/s). The Vp/Vs ratio varies spatially more distinctly than P- and S-wave velocities, usually from 1.70 to 1.74 in the upper crust and from 1.72 to 1.78 in the lower crust. Schist belts and their continuations at depth are associated with lower velocities and lower Vp/Vs ratios than in the granitoid areas. The tomography modelling suggests that the Svecofennian Orogen was accreted from crustal blocks ranging in size from 100x100 km2 to 200x200 km2 in cross-sectional area. The intervening sedimentary belts have ca. 0.2 km/s lower P- and S-wave velocities and ca. 0.04 lower Vp/Vs ratios. Thus, the tomographic model supports the concept that the thick Svecofennian crust was accreted from several crustal terranes, some hidden, and that the crust was later modified by intra- and underplating. In conclusion, as a novel approach the earthquake focal mechanism and focal depth distribution is discussed in relation to the 3-D tomography model. The schist belts and the transformation zones between the high- and low-velocity anomaly blocks are characterized by deeper earthquakes than the granitoid areas where shallow events dominate. Although only a few focal mechanisms were solved for southern Finland, there is a trend towards strike-slip and oblique strike-slip movements inside schist areas. The normal dip-slip type earthquakes are typical in the seismically active Kuusamo district in the NE edge of the SVEKALAPKO area, where the Archean crust is ca. 15-20 km thinner than the Proterozoic Svecofennian crust. Two near vertical dip-slip mechanism earthquakes occurred in the NE-SW junction between the Central Finland Granitoid Complex and the Vyborg rapakivi batholith, where high Vp/Vs-ratio deep-set intrusion splits the southern Finland schist belt into two parts in the tomography model.

Physical Properties of Deep Drill Cores : Implications for Meteorite Impact Effects and Crustal Structures

Physical properties provide valuable information about the nature and behavior of rocks and minerals. The changes in rock physical properties generate petrophysical contrasts between various lithologies, for example, between shocked and unshocked rocks in meteorite impact structures or between various lithologies in the crust. These contrasts may cause distinct geophysical anomalies, which are often diagnostic to their primary cause (impact, tectonism, etc). This information is vital to understand the fundamental Earth processes, such as impact cratering and associated crustal deformations. However, most of the present day knowledge of changes in rock physical properties is limited due to a lack of petrophysical data of subsurface samples, especially for meteorite impact structures, since they are often buried under post-impact lithologies or eroded. In order to explore the uppermost crust, deep drillings are required. This dissertation is based on the deep drill core data from three impact structures: (i) the Bosumtwi impact structure (diameter 10.5 km, 1.07 Ma age; Ghana), (ii) the Chesapeake Bay impact structure (85 km, 35 Ma; Virginia, U.S.A.), and (iii) the Chicxulub impact structure (180 km, 65 Ma; Mexico). These drill cores have yielded all basic lithologies associated with impact craters such as post-impact lithologies, impact rocks including suevites and breccias, as well as fractured and unfractured target rocks. The fourth study case of this dissertation deals with the data of the Paleoproterozoic Outokumpu area (Finland), as a non-impact crustal case, where a deep drilling through an economically important ophiolite complex was carried out. The focus in all four cases was to combine results of basic petrophysical studies of relevant rocks of these crustal structures in order to identify and characterize various lithologies by their physical properties and, in this way, to provide new input data for geophysical modellings. Furthermore, the rock magnetic and paleomagnetic properties of three impact structures, combined with basic petrophysics, were used to acquire insight into the impact generated changes in rocks and their magnetic minerals, in order to better understand the influence of impact. The obtained petrophysical data outline the various lithologies and divide rocks into four domains. Based on target lithology the physical properties of the unshocked target rocks are controlled by mineral composition or fabric, particularly porosity in sedimentary rocks, while sediments result from diverse sedimentation and diagenesis processes. The impact rocks, such as breccias and suevites, strongly reflect the impact formation mechanism and are distinguishable from the other lithologies by their density, porosity and magnetic properties. The numerous shock features resulting from melting, brecciation and fracturing of the target rocks, can be seen in the changes of physical properties. These features include an increase in porosity and subsequent decrease in density in impact derived units, either an increase or a decrease in magnetic properties (depending on a specific case), as well as large heterogeneity in physical properties. In few cases a slight gradual downward decrease in porosity, as a shock-induced fracturing, was observed. Coupled with rock magnetic studies, the impact generated changes in magnetic fraction the shock-induced magnetic grain size reduction, hydrothermal- or melting-related magnetic mineral alteration, shock demagnetization and shock- or temperature-related remagnetization can be seen. The Outokumpu drill core shows varying velocities throughout the drill core depending on the microcracking and sample conditions. This is similar to observations by Kern et al., (2009), who also reported the velocity dependence on anisotropy. The physical properties are also used to explain the distinct crustal reflectors as observed in seismic reflection studies in the Outokumpu area. According to the seismic velocity data, the interfaces between the diopside-tremolite skarn layer and either serpentinite, mica schist or black schist are causing the strong seismic reflectivities.

Geochemistry and petrology of the ferropicrite dikes and associated rocks of Vestfjella, western Dronning Maud Land, Antarctica

This study provides insights into the composition and origin of ferropicrite dikes (FeOtot = 13 17 wt. %; MgO = 13 19 wt. %) and associated meimechite, picrite, picrobasalt, and basalt dikes found at Vestfjella, western Dronning Maud Land, Antarctica. The dikes crosscut Jurassic Karoo continental flood basalts (CFB) that were emplaced during the early stages of the breakup of the Gondwana supercontinent ~180 Ma ago. Selected samples (31 overall from at least eleven dikes) were analyzed for their mineral chemical, major element, trace element, and Sr, Nd, Pb, and Os isotopic compositions. The studied samples can be divided into two geochemically distinct types: (1) The depleted type (24 samples from at least nine dikes) is relatively depleted in the most incompatible elements and exhibits isotopic characteristics (e.g., initial εNd of +4.8 to +8.3 and initial 187Os/188Os of 0.1256 0.1277 at 180 Ma) similar to those of mid-ocean ridge basalts (MORB); (2) The enriched type (7 samples from at least two dikes) exhibits relatively enriched incompatible element and isotopic characteristics (e.g., initial εNd of +1.8 to +3.6 and initial 187Os/188Os of 0.1401 0.1425 at 180 Ma) similar to those of oceanic island basalts. Both magma types have escaped significant contamination by the continental crust. The depleted type is related to the main phase of Karoo magmatism and originated as highly magnesian (MgO up to 25 wt. %) partial melts at high temperatures (mantle potential temperature >1600 °C) and pressures (~5 6 GPa) from a sublithospheric, water-bearing, depleted peridotite mantle source. The enriched type sampled pyroxene-bearing heterogeneities that can be traced down to either recycled oceanic crust or melt-metasomatized portions of the sublithospheric or lithospheric mantle. The source of the depleted type represents a sublithospheric end-member source for many Karoo lavas and has subsequently been sampled by the MORBs of the Indian Ocean. These observations, together with the purported high temperatures, indicate that the Karoo CFBs were formed in an extensive melting episode caused mainly by internal heating of the upper mantle beneath the Gondwana supercontinent. My research supports the view that ferropicritic melts can be generated in several ways: the relative Fe-enrichment of mantle partial melts is most readily achieved by (1) relatively low degree of partial melting, (2) high pressure of partial melting, and (3) melting of enriched source components (e.g., pyroxenite and metasomatized peridotite). Ferropicritic whole-rock compositions could also result from accumulation, secondary alteration, and fractional crystallization, however, and caution is required when addressing the parental magma composition.

Svekofennialaisen orogenian ylityönnön terminen mallinnus Itä-Suomen arkeeisella kratonilla

This study reviews the thermophysical properties of rocks and the lithosphere and describes a one dimensional thermal numeric model of hypothetical 20 km thick overthrust plate obtruded on to the Archeaen craton in the Svecofennian orogeny (1.92-1.77 Ga). The objective is to find out if the overthrust plate and its radiogenic heat sources were able to produce the thermal effects observed on the current erosion level of the Archaean craton. Heat transfer in lithosphere is assumed conductive, and advective heat transfer due to melting and melt transfer is supposed negligible. The study area is located in the Eastern Finland, approximately current Kainuu and Northern Karelia regions, east from the most active orogenic belt (Raahe-Ladoga zone), so that orogenic magmatism can be neglected. Physical parameters and boundary conditions for the model are from different earlier published sources: deep seismic profiles (rock variation in depth), laboratory measurements (heat production and conductivity of rocks), field measurements (heat flow densities), and pT(t) estimations from the Finnish precambrian to estimate the size and thickness of the sheet. Comparison of the modelling results to previous K-Ar datings and other pTt estimations show, that the effect of the overthrust sheet has been adequate (max. T 450°C at 4 kbar) to produce the K-Ar resetting ages measured from the Archaean bedrock at current erosion level. No other kind of thermal activation in lithosphere is required. Results show possibly very minor partial melting in upper middle crust underneath the overthrust sheet.

The prehistory of Suomussalmi, eastern Finland; the first billion years as revealed by isotopes and the composition of granitoid suites

The four papers summarized in this thesis deal with the Archean and earliest Paleoproterozoic granitoid suites observed in the Suomussalmi district, eastern Finland. Geologically, the area belongs to the Kianta Complex of the Western Karelian Terrane in the Karelian Province of the Fennoscandian shield. The inherited zircons up to 3440 Ma old together with Sm Nd and Pb Pb data confirm the existence of previously anticipated Paleoarchean protocrust in Suomussalmi. The general timeline of granitoid magmatism is similar to that of the surrounding areas. TTG magmatism occurred in three distinct phases: ca 2.95 Ga, 2.83 2.78 Ga and 2.76 2.74 Ga. In Suomussalmi the TTGs sensu stricto (K2O/Na2O less than 0.5) belong to the low-HREE type and are interpreted as partial melts of garnet amphibolites, which did not significantly interact with mantle peridotites. Transitional TTGs (K2O/Na2O more than 0.5), present in Suomussalmi and absent from surrounding areas, display higher LILE concentrations, but otherwise closely resemble the TTGs sensu stricto and indicate that recycling of felsic crust commenced in Suomussalmi 200 Ma earlier than in surrounding areas. The youngest TTG phase was coeval with the intrusion of the Likamännikkö quartz alkali feldspar syenite (2741 ± 2 Ma) complex. The complex contains angular fragments of ultrabasic rock, which display considerable compositional heterogeneity and are interpreted as cumulates containing clinopyroxene (generally altered to actinolite), apatite, allanite, epidote, and albite. The quartz alkali feldspar syenite cannot be regarded as alkaline sensu stricto, despite clear alkaline affinities. Within Likamännikkö there are also calcite carbonatite patches, which display mantle-like O- and C-isotope values, as well as trace element characteristics consistent with a magmatic origin, and could thus be among the oldest known carbonatites in the world. Sanukitoid (2.73 2.71 Ga) and quartz diorite suites (2.70 Ga) overlap within error margins and display compositional similarities, but can be differentiated from each other on the basis of higher Ba, K2O and LREE contents of the sanukitoids. The Likamännikkö complex, sanukitoids and quartz diorites are interpreted as originating from the metasomatized mantle and mark the diversification of the granitoid clan after 200 Ma of evolution dominated by the TTG suite. Widespread migmatization and the intrusion of anatectic leucogranitoids as dykes and intrusions of varying size took place at 2.70 2.69 Ga, following collisional thickening of the crust. The leucogranitoids and leucosomes of migmatized TTGs are compositionally alike and characterized by high silica contents and a leucocratic appearance. Due to compositional overlap, definitive discrimination between leucogranitoids and transitional TTGs requires isotope datings and/or knowledge of field relationships. Leucogranitoids represent partial melts of the local TTGs, both the sensu stricto and transitional types, mostly derived under water fluxed conditions, with possible fluid sources being late sanukitoids and quartz diorites as well as dehydrating lower crust. The Paleoproterozoic 2.44 2.39 Ga A-type granitoids of the Kianta Complex emplaced in an extensional environment are linked to the coeval and more widespread mafic intrusions and dykes observed over most of the Archean nucleus of the Fennoscandian shield. The A-type intrusions in the Suomussalmi area are interpreted as partial melts of the Archean lower crust and display differences in composition and magnetite content, which indicate differences in the composition and oxidation state of the source.

Neoarchean sanukitoid series in the Karelian Province, Finland

Sanukitoid series intrusions can be found throughout the Archean Karelian Province of the Fennoscandian shield. All sanukitoids share the same controversial elemental characteristics: they have high content of incompatible elements such as K, Ba, and Sr as well as high content of the compatible elements Mg, Cr, and Ni, and high Mg#. This composition is explained by an enriched mantle wedge origin in a Neoarchean subduction setting. This study concentrates on sanukitoid intrusions and tonalite-trondhjemite-granodiorite series (TTGs) from Finnish part of the Karelian Province. The collected rock samples have been studied in the field and under microscope as well as for their whole-rock (including isotopes) and mineral compositions. The new data together with previously published analyses help us to better understand the petrogenesis, tectonic setting and reworking of the Archean rock units. TTGs from the Karelian Province form a voluminous series of granitoids and reworked migmatites. This study divides TTG series into two subgroups based on their elemental composition: low-HREE (heavy rare earth element) TTGs and high-HREE TTGs indicating pressure differences in their source. Sanukitoid series is a minor, divergent group of intrusions. These intrusions are variable sized, and the texture varies from even-grained to K-feldspar porphyritic. The elemental composition differentiates sanukitoids from more voluminous TTG groups, the SiO2 in sanukitoids varies to include series of gabbro, diorite, and granodiorite. U Pb age determinations from sanukitoid series show temporally limited emplacement between ~ 2745 2715 Ma after the main crust forming period in the area. Hafnium, neodymium, common lead, and oxygene isotopes indicate well homogenized characteristics. Recycled crust has made a variable, yet minor, contribution to sanukitoids, as evidenced by oxygene isotopes and inherited zircon cores. A proposed tectonic setting for the formation of the sanukitoid series is slab breakoff of oceanic lithosphere in subduction setting, with sanukitoids deriving from an enriched mantle wedge. The proposed setting explains some of the peculiar features of sanukitoids, such as their temporally limited occurrence and controversial elemental composition. Sanukitoids would occur after cessation of the regional growth of Archean crust, and they could be derived from mantle wedge previously enriched by melts and fluids from oceanic crust and sediments. A subsequent event during the Paleoproterozoic Svecofennian orogeny at ~1.9 Ga affected the appearance and microstructures of the rocks as well as caused redistribution of lead between minerals and whole rock. However, the deformation was not able to obliterate the original geochemical characteristics of these sanukitoids.

Local systems, global impacts : Using life cycle assessment to analyse the potential and constraints of industrial symbioses

Human activities extract and displace different substances and materials from the earth s crust, thus causing various environmental problems, such as climate change, acidification and eutrophication. As problems have become more complicated, more holistic measures that consider the origins and sources of pollutants have been called for. Industrial ecology is a field of science that forms a comprehensive framework for studying the interactions between the modern technological society and the environment. Industrial ecology considers humans and their technologies to be part of the natural environment, not separate from it. Industrial operations form natural systems that must also function as such within the constraints set by the biosphere. Industrial symbiosis (IS) is a central concept of industrial ecology. Industrial symbiosis studies look at the physical flows of materials and energy in local industrial systems. In an ideal IS, waste material and energy are exchanged by the actors of the system, thereby reducing the consumption of virgin material and energy inputs and the generation of waste and emissions. Companies are seen as part of the chains of suppliers and consumers that resemble those of natural ecosystems. The aim of this study was to analyse the environmental performance of an industrial symbiosis based on pulp and paper production, taking into account life cycle impacts as well. Life Cycle Assessment (LCA) is a tool for quantitatively and systematically evaluating the environmental aspects of a product, technology or service throughout its whole life cycle. Moreover, the Natural Step Sustainability Principles formed a conceptual framework for assessing the environmental performance of the case study symbiosis (Paper I). The environmental performance of the case study symbiosis was compared to four counterfactual reference scenarios in which the actors of the symbiosis operated on their own. The research methods used were process-based life cycle assessment (LCA) (Papers II and III) and hybrid LCA, which combines both process and input-output LCA (Paper IV). The results showed that the environmental impacts caused by the extraction and processing of the materials and the energy used by the symbiosis were considerable. If only the direct emissions and resource use of the symbiosis had been considered, less than half of the total environmental impacts of the system would have been taken into account. When the results were compared with the counterfactual reference scenarios, the net environmental impacts of the symbiosis were smaller than those of the reference scenarios. The reduction in environmental impacts was mainly due to changes in the way energy was produced. However, the results are sensitive to the way the reference scenarios are defined. LCA is a useful tool for assessing the overall environmental performance of industrial symbioses. It is recommended that in addition to the direct effects, the upstream impacts should be taken into account as well when assessing the environmental performance of industrial symbioses. Industrial symbiosis should be seen as part of the process of improving the environmental performance of a system. In some cases, it may be more efficient, from an environmental point of view, to focus on supply chain management instead.  

Inverse problem for the wave equation : partial data and novel boundary sources

An inverse problem for the wave equation is a mathematical formulation of the problem to convert measurements of sound waves to information about the wave speed governing the propagation of the waves. This doctoral thesis extends the theory on the inverse problems for the wave equation in cases with partial measurement data and also considers detection of discontinuous interfaces in the wave speed. A possible application of the theory is obstetric sonography in which ultrasound measurements are transformed into an image of the fetus in its mother's uterus. The wave speed inside the body can not be directly observed but sound waves can be produced outside the body and their echoes from the body can be recorded. The present work contains five research articles. In the first and the fifth articles we show that it is possible to determine the wave speed uniquely by using far apart sound sources and receivers. This extends a previously known result which requires the sound waves to be produced and recorded in the same place. Our result is motivated by a possible application to reflection seismology which seeks to create an image of the Earth s crust from recording of echoes stimulated for example by explosions. For this purpose, the receivers can not typically lie near the powerful sound sources. In the second article we present a sound source that allows us to recover many essential features of the wave speed from the echo produced by the source. Moreover, these features are known to determine the wave speed under certain geometric assumptions. Previously known results permitted the same features to be recovered only by sequential measurement of echoes produced by multiple different sources. The reduced number of measurements could increase the number possible applications of acoustic probing. In the third and fourth articles we develop an acoustic probing method to locate discontinuous interfaces in the wave speed. These interfaces typically correspond to interfaces between different materials and their locations are of interest in many applications. There are many previous approaches to this problem but none of them exploits sound sources varying freely in time. Our use of more variable sources could allow more robust implementation of the probing.

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.