Geochemistry, petrogenesis and tectonic setting of ophiolites and mafic-ultramafic complexes in the Northeastern Aegean Region

In this study two ophiolites and a mafic-ultramafic complexes of the northeastern Aegean Sea, Greece, have been investigated to re-evaluate their petrogenetic evolution and tectonic setting. These complexes are: the mafic-ultramafic complex of Lesvos Island and the ophiolites of Samothraki Island and the Evros area. In order to examine these complexes in detail whole-rock major- and trace-elements as well as Sr and Nd isotopes, and minerals were analysed and U-Pb SHRIMP ages on zircons were determined. The mafic-ultramafic complex of Lesvos Island consists of mantle peridotite thrusted over a tectonic mélange containing metasediments, metabasalts and a few metagabbros. This succession had previously been interpreted as an ophiolite of Late Jurassic age. The new field and geochemical data allow a reinterpretation of this complex as representing an incipient continental rift setting that led to the subsequent formation of the Meliata-Maliac-Vardar branches of Neotethys in Upper Permian times (253 ± 6 Ma) and the term “Lesvos ophiolite” should be abandoned. With proceeding subduction and closure of the Maliac Ocean in Late Jurassic times (155 Ma) the Lesvos mafic-ultramafic complex was obducted. Zircon ages of 777, 539 and 338 Ma from a gabbro strongly suggest inheritance from the intruded basement and correspond to ages of distinct terranes recently recognized in the Hellenides (e.g. Florina terrane). Geochemical similar complexes which contain rift associations with Permo-Triassic ages can be found elsewhere in Greece and Turkey, namely the Teke Dere Thrust Sheet below the Lycian Nappes (SW Turkey), the Pindos subophiolitic mélange (W Greece), the Volcanosedimentary Complex on Central Evia Island (Greece) and the Karakaya Complex (NW Turkey). This infers that the rift-related rocks from Lesvos belong to an important Permo-Triassic rifting episode in the eastern Mediterranean. The ‘in-situ’ ophiolite of Samothraki Island comprises gabbros, sparse dykes and basalt flows as well as pillows cut by late dolerite dykes and had conventionally been interpreted as having formed in an ensialic back-arc basin. The results of this study revealed that none of the basalts and dolerites resemble mid-ocean ridge or back-arc basin basalts thus suggesting that the Samothraki ophiolite cannot represent mature back-arc basin crust. The age of the complex is regarded to be 160 ± 5 Ma (i.e. Oxfordian; early Upper Jurassic), which precludes any correlation with the Lesvos mafic-ultramafic complex further south (253 ± 6 Ma; Upper Permian). Restoration of the block configuration in NE Greece, before extensional collapse of the Hellenic hinterland and exhumation of the Rhodope Metamorphic Core Complex (mid-Eocene to mid-Miocene), results in a continuous ophiolite belt from Guevgueli in the NW to Samothraki in the SE, thus assigning the latter to the Innermost Hellenic Ophiolite Belt. In view of the data of this study, the Samothraki ophiolite represents a rift propagation of the Sithonia ophiolite spreading ridge into the Chortiatis calc-alkaline arc. The ophiolite of the Evros area consists of a plutonic sequence comprising cumulate and non-cumulate gabbros with plagiogranite veins, and an extrusive sequence of basalt dykes, massive and pillow lavas as well as pyroclastic rocks. Furthermore, in the Rhodope Massif tectonic lenses of harzburgites and dunites can be found. All rocks are spatially separated. The analytical results of this study revealed an intra-oceanic island arc setting for the Evros ophiolitic rocks. During late Middle Jurassic times (169 ± 2 Ma) an intra-oceanic arc has developed above a northwards directed intra-oceanic subduction zone of the Vardar Ocean in front of the Rhodope Massif. The boninitic, island arc tholeiitic and calc-alkaline rocks reflect the evolution of the Evros island arc. The obduction of the ophiolitic rocks onto the Rhodope basement margin took place during closure of the Vardar ocean basins. The harzburgites and dunites of the Rhodope Massif are strongly depleted and resemble harzburgites from recent oceanic island arcs. After melt extraction they underwent enrichment processes by percolating melts and fluids from the subducted slab. The relationship of the peridotites and the Evros ophiolite is still ambiguous, but the stratigraphic positions of the peridotites and the ophiolitic rocks indicate separated origin. The harzburgites and dunites most probably represent remnants of the mantle wedge of the island arc of the Rhodope terrane formed above subducted slab of the Nestos Ocean in late Middle Jurassic times. During collision of the Thracia terrane with the Rhodope terrane thrusting of the Rhodope terrane onto the Thracia terrane took place, whereas the harzburgites and dunites were pushed between the two terranes now cropping out on top of the Thracia terrane of the Rhodope Massif.

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.

Petrogenetic evolution of Uralian-Alaskan-type mafic-ultramafic complexes in the southern and middle Urals, Russia

The PhD thesis at hand consists of three parts and describes the petrogenetic evolution of Uralian-Alaskan-type mafic ultramafic complexes in the Ural Mountains, Russia. Uralian-Alaskan-type mafic-ultramafic complexes are recognized as a distinct class of intrusions. Characteristic petrologic features are the concentric zonation of a central dunite body grading outward into wehrlite, clinopyroxenite and gabbro, the absence of orthopyroxene and frequently occurring platinum group element (PGE) mineralization. In addition, the presence of ferric iron-rich spinel discriminates Uralian-Alaskan-type complexes from most other mafic ultramafic rock assemblages. The studied Uralian-Alaskan-type complexes (Nizhnii Tagil, Kytlym and Svetley Bor) belong to the southern part of a 900 km long, N–S-trending chain of similar intrusions between the Main Uralian Fault to the west and the Serov-Mauk Fault to the east. The first chapter of this thesis studies the evolution of the ultramafic rocks tracing the compositional variations of rock forming and accessory minerals. The comparison of the chemical composition of olivine, clinopyroxene and chromian spinel from the Urals with data from other localities indicates that they are unique intrusions having a characteristic spinel and clinopyroxene chemistry. Laser ablation-ICPMS (LA-ICPMS ) analyses of trace element concentrations in clinopyroxene are used to calculate the composition of their parental melt which is characterized by enriched LREE (0.5-5.2 prim. mantle) and other highly incompatible elements (U, Th, Ba, Rb) relative to the HREE (0.25-2.0 prim. mantle). A subduction-related geotectonic setting is indicated by a positive anomaly for Sr and negative anomalies for Ti, Zr and Hf. The mineral compositions monitor the evolution of the parental magmas and decipher differences between the studied complexes. In addition, the observed variation in LREE/HREE (for example La/Lu = 2-24) can be best explained with the model of an episodically replenished and erupted open magma chamber system with the extensive fractionation of olivine, spinel and clinopyroxene. The data also show that ankaramites in a subduction-related geotectonic setting could represent parental magmas of Uralian-Alaskan-type complexes. The second chapter of the thesis discusses the chemical variation of major and trace elements in rock-forming minerals of the mafic rocks. Electron microprobe and LA-ICPMS analyses are used to quantitatively describe the petrogenetic relationship between the different gabbroic lithologies and their genetic link to the ultramafic rocks. The composition of clinopyroxene identifies the presence of melts with different trace element abundances on the scale of a thin section and suggests the presence of open system crustal magma chambers. Even on a regional scale the large variation of trace element concentrations and ratios in clinopyroxene (e.g. La/Lu = 3-55) is best explained by the interaction of at least two fundamentally different magma types at various stages of fractionation. This requires the existence of a complex magma chamber system fed with multiple pulses of magmas from at least two different coeval sources in a subduction-related environment. One source produces silica saturated Island arc tholeiitic melts. The second source produces silica undersaturated, ultra-calcic, alkaline melts. Taken these data collectively, the mixing of the two different parental magmas is the dominant petrogenetic process explaining the observed chemical variations. The results further imply that this is an intrinsic feature of Uralian-Alaskan-type complexes and probably of many similar mafic-ultramafic complexes world-wide. In the third chapter of this thesis the major element composition of homogeneous and exsolved spinel is used as a petrogenetic indicator. Homogeneous chromian spinel in dunites and wehrlites monitors the fractionation during the early stages of the magma chamber and the onset of clinopyroxene fractionation as well as the reaction of spinel with interstitial liquid. Exsolved spinel is present in mafic and ultramafic rocks from all three studied complexes. Its composition lies along a solvus curve which defines an equilibrium temperature of 600°C, given that spinel coexists with olivine. This temperature is considered to be close to the temperature of the host rocks into which the studied Uralian-Alaskan-type complexes intruded. The similarity of the exsolution temperatures in the different complexes over a distance of several hundred kilometres implies a regional tectonic event that terminated the exsolution process. This event is potentially associated with the final exhumation of the Uralian-Alaskan-type complexes along the Main Uralian Fault and the Serov-Mauk Fault in the Uralian fold belt.

Aeolian geomorphodynamics in endorheic basins of the Mongolian Gobi Desert

Die vorliegende Arbeit ist im Zuge des DFG Projektes Spätpleistozäne, holozäne und aktuelle Geomorphodynamik in abflusslosen Becken der Mongolischen Gobi´´ entstanden. Das Arbeitsgebiet befindet sich in der südlichen Mongolei im nördlichen Teil der Wüste Gobi. Neben einigen Teilen der Sahara (Heintzenberg, 2009), beispielsweise das Bodélé Becken des nördlichen Tschads (z.B. Washington et al., 2006a; Todd et al., 2006; Warren et al., 2007) wird Zentralasien als ein Hauptliefergebiet für Partikel in die globale Zirkulation der Atmosphäre gesehen (Goudie, 2009). Hauptaugenmerk liegt hierbei besonders auf den abflusslosen Becken und deren Sedimentablagerungen. Die, der Deflation ausgesetzten Flächen der Seebecken, sind hauptsächliche Quelle für Partikel die sich in Form von Staub respektive Sand ausbreiten. Im Hinblick auf geomorphologische Landschaftsentwicklung wurde der Zusammenhang von Beckensedimenten zu Hangdepositionen numerisch simuliert. Ein von Grunert and Lehmkuhl (2004) publiziertes Model, angelehnt an Ideen von Pye (1995) wird damit in Betracht gezogen. Die vorliegenden Untersuchungen modellieren Verbreitungsmechanismen auf regionaler Ebene ausgehend von einer größeren Anzahl an einzelnen punktuellen Standorten. Diese sind repräsentativ für die einzelnen geomorphologischen Systemglieder mit möglicherweise einer Beteiligung am Budget aeolischer Geomorphodynamik. Die Bodenbedeckung durch das charakteristische Steinpflaster der Gobi - Region, sowie unter anderem Korngrößenverteilungen der Oberflächensedimente wurden untersucht. Des Weiteren diente eine zehnjährige Zeitreihe (Jan 1998 bis Dez 2007) meteorologischer Daten als Grundlage zur Analyse der Bedingungen für äolische Geomorphodynamik. Die Daten stammen von 32 staatlichen mongolischen Wetterstationen aus der Region und Teile davon wurden für die Simulationen verwendet. Zusätzlich wurden atmosphärische Messungen zur Untersuchung der atmosphärischen Stabilität und ihrer tageszeitlichen Variabilität mit Mess-Drachenaufstiegen vorgenommen. Die Feldbefunde und auch die Ergebnisse der Laboruntersuchungen sowie der Datensatz meteorologischer Parameter dienten als Eingangsparameter für die Modellierungen. Emissionsraten der einzelnen Standorte und die Partikelverteilung im 3D Windfeld wurden modelliert um die Konvektivität der Beckensedimente und Hangdepositionen zu simulieren. Im Falle hoher mechanischer Turbulenz der bodennahen Luftschicht (mit einhergehender hoher Wind Reibungsgeschwindigkeit), wurde generell eine neutrale Stabilität festgestellt und die Simulationen von Partikelemission sowie deren Ausbreitung und Deposition unter neutraler Stabilitätsbedingung berechnet. Die Berechnung der Partikelemission wurde auf der Grundlage eines sehr vereinfachten missionsmodells in Anlehnung an bestehende Untersuchungen (Laurent et al., 2006; Darmenova et al., 2009; Shao and Dong, 2006; Alfaro, 2008) durchgeführt. Sowohl 3D Windfeldkalkulationen als auch unterschiedliche Ausbreitungsszenarien äolischer Sedimente wurden mit dem kommerziellen Programm LASAT® (Lagrange-Simulation von Aerosol-Transport) realisiert. Diesem liegt ein Langargischer Algorithmus zugrunde, mittels dessen die Verbreitung einzelner Partikel im Windfeld mit statistischer Wahrscheinlichkeit berechnet wird. Über Sedimentationsparameter kann damit ein Ausbreitungsmodell der Beckensedimente in Hinblick auf die Gebirgsfußflächen und -hänge generiert werden. Ein weiterer Teil der Untersuchungen beschäftigt sich mit der geochemischen Zusammensetzung der Oberflächensedimente. Diese Proxy sollte dazu dienen die simulierten Ausbreitungsrichtungen der Partikel aus unterschiedlichen Quellregionen nach zu verfolgen. Im Falle der Mongolischen Gobi zeigte sich eine weitestgehende Homogenität der Minerale und chemischen Elemente in den Sedimenten. Laser Bebohrungen einzelner Sandkörner zeigten nur sehr leichte Unterschiede in Abhängigkeit der Quellregionen. Die Spektren der Minerale und untersuchten Elemente deuten auf graitische Zusammensetzungen hin. Die, im Untersuchungsgebiet weit verbreiteten Alkali-Granite (Jahn et al., 2009) zeigten sich als hauptverantwortlich für die Sedimentproduktion im Untersuchungsgebiet. Neben diesen Mineral- und Elementbestimmungen wurde die Leichtmineralfraktion auf die Charakteristik des Quarzes hin untersucht. Dazu wurden Quarzgehalt, Kristallisation und das Elektronen-Spin-Resonanz Signal des E’1 - Centers in Sauerstoff Fehlstellungen des SiO2 Gitters bestimmt. Die Untersuchungen sind mit dem Methodenvorschlag von Sun et al. (2007) durchgeführt worden und sind prinzipiell gut geeignet um Herkunftsanalysenrndurchzuführen. Eine signifikante Zuordnung der einzelnen Quellgebiete ist jedoch auch in dieser Proxy nicht zu finden gewesen.