Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.

Softening the Lower Crust: Modes of Syn-Transport Transposition Around and Adjacent to a Deep Crustal Granulite Nappe, Parry Sound Domain, Grenville Province, Ontario, Canada

The Parry Sound domain is a granulite nappe-stack transported cratonward during reactivation of the ductile lower and middle crust in the late convergence of the Mesoproterozoic Grenville orogeny. Field observations suggest the following with respect to the ductile sheath: (1) Formation of a carapace of transposed amphibolite facies gneiss derived from and enveloping the western extremity of the Parry Sound domain and separating it from high-strain gneiss of adjacent allochthons. This ductile sheath formed dynamically around the moving granulite nappe through the development of systems of progressively linked shear zones. (2) Transposition initiated by hydration (amphibolization) of granulite facies gneiss by introduction of fluid along cracks accompanying pegmatite emplacement. Shear zones nucleated along pegmatite margins and subsequently linked and rotated. The source of the pegmatites was most likely subjacent migmatitic and pegmatite-rich units or units over which Parry Sound domain was transported. Comparison of gneisses of the ductile sheath with high-strain layered gneiss of adjacent allochthons show the mode of transposition of penetratively layered gneiss depended on whether or not the gneiss protoliths were amphibolite or granulite facies tectonites before initiation of transposition, resulting in, e.g., folding before shearing, no folding before shearing, respectively. Meter-scale truncation along high-strain gradients at the margins of both types of transposition-related shear zones observed within and marginal to Parry Sound domain mimic features at kilometer scales, implying that apparent truncation by transposition originating in a manner similar to the ductile sheath may be a common feature of deep crustal ductile reworking. Citation: Culshaw, N., C. Gerbi, and J. Marsh (2010), Softening the lower crust: Modes of syn-transport transposition around and adjacent to a deep crustal granulite nappe, Parry Sound domain, Grenville Province, Ontario, Canada, Tectonics, 29, TC5013, doi:10.1029/2009TC002537.

Isotopic and chemical compositions of lower crust rocks and minerals from ODP Hole 176-735B

The composition of gabbroic rocks from the drill core of Hole 735B (ODP Leg 176) at the 11 Ma Atlantis II bank close to the slow spreading Southwest Indian Ridge (SWIR) has been analyzed for major and trace elements and Sr, Nd and Pb isotopic composition. The samples are thought to represent much of the mineralogical and geochemical variation in a vertical 1-km section (500-1500 m below the sea floor) of the lower ocean crust. Primitive troctolitic gabbros, olivine gabbros and gabbros that have Mg#=84-70, Ca#>61 and low Na# (Na/(Na+Al)) (8-17) are intruded by patches or veins of more evolved FeTi-oxide rich gabbroic and dioritic rocks with Mg# to 20, Ca# to 32, Na#=14-23, TiO2<7 wt.% and FeOtotal<18 wt.%. All rocks are acdcumulates, and incompatible element concentrations are low, e.g. Pb=0.1-0.7 ppm and Ulower values, 0.70276. Separated clinopyroxene has elevated 87Sr/86Sr up to 0.7035, while plagioclase generally has close to whole rock Sr. Leaching reduced 87Sr/86Sr in clinopyroxene and in two (out of nine) cases leached separates and whole rock display isotopic equilibrium. Relatively minor hydrothermal seawater alteration is thought to have increased 87Sr/86Sr in the rocks, while a secondary high temperature percolation of a mantle-derived agent is thought to be the cause for the trend towards radiogenic Pb. This material had intermediate 87Sr/86Sr and may have originated from non-MORB off axis mantle. The main primary igneous isotopic variation of the gabbros is suggested to have been derived from the MORB-mantle and is defined mainly by leached samples from both ODP Leg 176 and Leg 118 and can be explained by two-component mixing of an end-member with composition like Central Indian Ridge basalts and an end-member with composition unlike any MORB. The latter is characterized by very unradiogenic Pb, in particular 207Pb/204Pb, and may have an origin with affinity to old depleted mantle (DM). The isotopic composition of the magmas parental to the FeTi-oxide rich rocks cannot be distinguished from the magmas parental to the primitive gabbros and an intimate relationship is indicated. The small-scale inhomogeneity indicated for the SWIR MORB-mantle at the Atlantis II Fracture Zone was probably inherited by the lower crustal rocks due to small-scale melting and monogenetic magma chambers at this slow spreading ridge.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.

Mineral and oxygen isotopic compositions in a profile through the altered lower oceanic crust, ODP Hole 176-735B

ODP Hole 735B located on the Southwest Indian Ridge at 57°E is an in situ sampled long, continuous section of lower oceanic crust. Oxygen isotope compositions of constituent minerals of Leg 176 gabbros have been measured by UV-laser oxygen isotope microprobe. Together with existing data from Leg 118, a complete oxygen isotope profile through the lower oceanic crust has been obtained. Most clinopyroxenes and olivines have normal mantle values of ~5.5 per mil and ~5.2 per mil, respectively, while plagioclases show slight d18O enrichment relative to its mantle value of 6.1per mil. Down-hole variations of Hole 735B gabbro indicate a downward decreasing d18O profile, with a kink at a depth of about 800 m below sea floor. Above this depth, gabbros are depleted in 18O relative to unaltered basalts, while below ~800 m they show nearly unmodified d18O values. Abundant seawater penetration appears to be limited to the upper part of the lower crust at ODP site 735 (~800 m into the gabbroic layer and ~2-2.5 km into the oceanic crust from the top of pillow basalts). Mass balance calculations show that the lower crust formed under this ultra-slow-spreading ridge has an average d18O value of 5.5 per mil. The whole crust at Site 735 has an overall 18O enrichment with d18O values of 6.0 per mil to 7.8 per mil, depending on the possible variation of the d18O values of the upper pillow basalts and sheeted dykes. The apparent difference in oxygen isotope compositions of ocean crusts formed with different spreading rates has important implications on the buffering of ocean water over geological time, as well as on the oxygen recycling between crust and mantle through subduction. The difference of seawater penetration between fast- and slow-spreading ridges could be related to their particular magmatic-tectonic history during the formation and aging of the crust. However, more analyses on continuous sections through oceanic and ophiolitic crust in different tectonic settings are required to derive any predictive models.

Peraluminous sapphirine-cordierite pods in Mg-rich orthopyroxene granulite from southern India: Implications for lower crustal processes

Sapphirine-cordierite intergrowths occur as pods within garnet-absent, high-Mg orthopyroxene-granulite xenoliths in the Kambam valley, Madurai Block, southern India. Whereas the cores of the pods are composed of sapphirine (X-Mg = 0.871-0.897) - cordierite (X-Mg = 0.892-0.931) intergrowth along with rutile, zircon and monazite, the rims are characterized by cordierite, apatite, plagioclase, K-feldspar, quartz and minor calcite. The surrounding matrix comprises orthopyroxene (maximum Al2O3 4.1 wt.%, X-Mg 0.848-0.850), plagioclase, biotite and quartz, similar to the assemblage in the surrounding charnockites. Sapphirine in the Kambam rocks is characterized by high Al contents with an end-member composition in the range of 7:9:3 and 3:5:1. The occurrence of peraluminous sapphirine in association with cordierite and in the absence of phases such as sillimanite and garnet is distinct from ultrahigh-temperature assemblages in other localities within the Madurai Block. The peraluminous composition of the pods suggests that these domains could represent cryptic pathways through which aluminous melts migrated. The reaction of such peraluminous melts with Mg-rich orthopyroxene in the host granulite at temperatures of 1025 degrees C and pressures around 8 kbar as computed from phase equilibria modeling followed by an isobaric cooling is inferred to have generated the sapphirine-cordierite pods. The unusual high-Mg orthopyroxene granulite suggests interaction of supracrustal rocks with mafic magmas, which probably acted as the heat source for the partial melting of lower crust and UHT metamorphism.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Mineraly of oceanic crust from the Kane Transform, Mid-Atlantic Ridge

Gabbroic cumulates drilled south of the Kane Transform Fault on the slow-spread Mid-Atlantic Ridge preserve up to three discrete magnetization components. Here we use absolute age constraints derived from the paleomagnetic data to develop a model for the magmatic construction of this section of the lower oceanic crust. By comparing the paleomagnetic data with mineral compositions, and based on thermal models of local reheating, we infer that magmas that began crystallizing in the upper mantle intruded into the lower oceanic crust and formed meter-scale sills. Some of these magmas were crystal-laden and the subsequent expulsion of interstitial liquid from them produced '"cumulus" sills. These small-scale magmatic injections took place over at least 210 000 years and at distances of ~3 km from the ridge axis and may have formed much of the lower crust. This model explains many of the complexities described in this area and can be used to help understand the general formation of oceanic crust at slow-spread ridges.

Stability of pargasite during ultrahigh-temperature metamorphism:A consequence of titanium and REE partitioning

Orthopyroxene-clinopyroxene-plagioclase needles and symplectite along the cleavage planes and grain boundaries of fluorine-bearing titanian-ferroan pargasite from the Highland Complex, Sri Lanka, are interpreted as evidence for dehydration melting at ultrahigh-temperature conditions. High Ti (up to 0.4 pfu) and F (XF up to 0.56) content in pargasite extends its stability to higher temperatures, and the composition indicates the dehydration melting reaction may take place at ultrahigh-temperatures (~950 °C) at a pressure around 10 kbar, close to peak metamorphic conditions. The increase of Ti content close to the grain boundaries and cleavage planes in pargasite indicates titanium partitioning from the melt during dehydration melting enhanced the stability of the mineral toward ultrahigh-temperature conditions. The REE content in the pargasite shows a similar behavior to that of titanium. The cores with no breakdown assemblage consist of low and flat REE concentrations with respect to the high and Eu-depleted rim. Clinopyroxene in symplectite and needle-shaped lamellae within the pargasite porphyroblasts have similar REE patterns with slightly low-concentrations relative to that of pargasite. In the breakdown assemblage, LREEs are partitioned mainly into plagioclase while the HREEs are partitioned into orthopyroxene. The REE enrichment in the pargasite rims signals their relative partitioning between pargasite rims and melt. Modeling of the partitioning of Ti and REEs associated with pargasite breakdown demonstrates that its stability is greatly enhanced at UHT conditions. This investigation implies that the stability of hydrous minerals such as amphibole can be extended to UHT conditions, and expands our knowledge of metamorphism in the lower crust.

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.

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.

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.