UHP metamorphic rocks of the Eastern Rhodope Massif, NE Greece: new constraints from petrology, geochemistry and zircon ages

In this PhD thesis, a multidisciplinary study has been carried out on metagranitoids and paragneisses from the Eastern Rhodope Massif, northern Greece, to decipher the pre-Alpine magmatic and geodynamic evolution of the Rhodope Massif and to correlate the eastern part with the western/central parts of the orogen. The Rhodope Massif, which occupies the major part of NE Greece and S Bulgaria, represents the easternmost part of the Internal Hellenides. It is regarded as a nappe stack of high-grade units, which is classically subdivided into an upper unit and a lower unit, separated by a SSE-NNW trending thrust plane, the Nestos thrust. Recent research in the central Greek Rhodope Massif revealed that the two units correspond to two distinct terranes of different age, the Permo-Carboniferous Thracia Terrane, which was overthrusted by the Late Jurassic/Early Cretaceous Rhodope Terrane. These terranes are separated by the Nestos suture, a composite zone comprising metapelites, metabasites, metagranitoids and marbles, which record high-pressure and even ultrahigh-pressure metamorphism in places. Similar characteristic rock associations were investigated during this study along several well-constrained cross sections in vincity to the Ada, Sidiro and Kimi villages in the Greek Eastern Rhodope Massif. Field evidence revealed that the contact zone of the two terranes in the Eastern Rhodope Massif is characterized by a mélange of metapelites, migmatitic amphibolites/eclogites, strongly sheared orthogneisses and marbles. The systematical occurrence of this characteristic rock association between the terranes implies that the Nestos suture is a continuous belt throughout the Greek Rhodope Massif. In this study, a new UHP locality could be established and for the first time in the Greek Rhodope, metamorphic microdiamonds were identified in situ in their host zircons using Laser-Raman spectroscopy. The presence of the diamonds as well as element distribution patterns of the zircons, obtained by TOF-SIMS, indicate metamorphic conditions of T > 1000 °C and P > 4 GPa. The high-pressure and ultrahigh-pressure rocks of the mélange zone are considered to have formed during the subduction of the Nestos Ocean in Jurassic times at ~150 Ma. Melting of metapelitic rocks at UHP conditions facilitated the exhumation to lower crustal levels. To identify major crust forming events, basement granitoids were dated by LA-SF-ICPMS and SHRIMP-II U-Pb analyses of zircons. The geochronological results revealed that the Eastern Rhodope Massif consists of two crustal units, a structurally lower Permo-Carboniferous unit corresponding to the Thracia Terrane and a structurally upper Late Jurassic/Early Cretaceous unit corresponding to the Rhodope Terrane, like it was documented for the Central Rhodope Massif. Inherited zircons in the orthogneisses from the Thracia Terrane of the Eastern Rhodope Massif indicate the presence of a pre-existing Neoproterozoic and Ordovician-Silurian basement in this region. Triassic magmatism is witnessed by the zircons of few orthogneisses from the easternmost Rhodope Massif and is interpreted to be related to rifting processes. Whole-rock major and trace element analyses indicate that the metagranitoids from both terranes originated in a subduction-related magmatic-arc environment. The Sr-Nd isotope data for both terranes of the Eastern and Central Rhodope Massif suggest a mixed crust-mantle source with variable contributions of older crustal material as already indicated by the presence of inherited zircons. Geochemical and isotopic similarity of the basement of the Thracia Terrane and the Pelagonian Zone implies that the Thracia Terrane is a fragment of a formerly unique Permo-Carboniferous basement, separated by rifting and opening of the Meliata-Maliac ocean system in Triassic times. A branch of the Meliata-Maliac ocean system, the Nestos Ocean, subducted northwards in Late Jurassic times leading to the formation of the Late Jurassic/Early Cretaceous Rhodope magmatic arc on remnants of the Thracia Terrane as suggested by inherited Permo-Carboniferous zircons. The ~150 Ma zircon ages of the orthogneisses from the Rhodope Terrane indicate that subduction-related magmatism and HP/UHP metamorphism occurred during the same subduction phase. Subduction ceased due to the closure of the Nestos Ocean in the Late Jurassic/Early Cretaceous. The post-Jurassic evolution of the Rhodope Massif is characterized by the exhumation of the Rhodope core complex in the course of extensional tectonics associated with late granite intrusions in Eocene to Miocene times.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Nature and origin of ultramafic lamprophyres and carbonatites from the borders of the Labrador Sea

Die Ränder des Labrador Meeres wurden während des späten Neoproterozoikums intensiv von karbonatreichen silikatischen Schmelzen durchsetzt. Diese Schmelzen bildeted sich bei Drucken zwischen ca. 4-6 GPa (ca. 120-180 km Tiefe) an der Basis der kontinentalen Mantel-Lithosphäre. Diese Magmengenerierung steht in zeitlichem und räumlichem Zusammenhang mit kontinentalen Extensionsprozessen, welche zu beiden Seiten des sich öffnenden Iapetus-Ozeans auftraten.

Ophiolites of the eastern Vardar Zone, N. Greece

In the present thesis, the geochemistry, petrology and geochronology of ophiolite complexes from central northern Greece were studied in detail in order to gain insights on the petrogenetic pathways and geodynamic processes that lead to their formation and evolution. The major- and trace-element content of minerals and whole rocks from all four ophiolite complexes was determined using high-precision analytical equipment. These results were then coupled with Nd and Sr isotopic measurements. In order to precisely place the evolution of these ophiolites in time, U-Pb geochronology on zircons was conducted using a SHRIMP-II. The data obtained suggest that the ophiolites studied invariably show typical characteristics of subduction-zone magmatism (e.g. negative Nb anomalies, Th enrichment). In N-MORB-normalised multielement profiles the high field-strength elements display patterns that vary from depleted to N-MORB-like. Chondrite-normalised rare-earth element (REE) profiles show flat heavy-REE patterns suggesting a shallow regime of source melting for all the ophiolites, well within the stability field of spinel lherzolite. The majority of the samples have light-REE depleted patterns. 87Sr/86Sr isotopic ratios range from 0.703184 to 0.715853 and are in cases influenced by alteration. The εNd values are positive (the majority of the mafic samples is typically 7.1-3.1) but lower than N-MORB and depleted mantle. With the exception of the Thessaloniki ophiolite that has uniform island-arc tholeiitic chemical characteristics, the rest of the ophiolites show dual chemistry consisting of rocks with minor subduction-zone characteristics that resemble chemically back-arc basin basalts (BABB) and rocks with more pronounced subduction-zone characteristics. Tectonomagmatic discrimination schemes classify the samples as island-arc tholeiites and back-arc basin basalts or N-MORB. Melting modelling carried out to evaluate source properties and degree of melting verifies the dual nature of the ophiolites. The samples that resemble back-arc basin basalts require very small degrees of melting (<10%) of fertile sources, whereas the rest of the samples require higher degrees (25-15%) of melting. As deduced from the present geochemical and petrological investigation, the ophiolites from Guevguely, Oraeokastro, Thessaloniki, and Chalkidiki represent relics of supra-subduction zone crust that formed in succeeding stages of island-arc rifting and back-arc spreading as well as in a fore arc setting. The geochronological results have provided precise determination of the timing of formation of these complexes. The age of the Guevguely ophiolite has been determined as 167±1.2 Ma, that of Thessaloniki as 169±1.4 Ma, that of Kassandra as 167±2.2 Ma and that of Sithonia as 160±1.2 Ma.

Deciphering petrological signatures of reactive melt stagnation and cooling in the oceanic mantle underneath ultraslow-spreading ridges

The global mid-ocean ridge system creates oceanic crust and lithosphere that covers more than two-thirds of the Earth. Basalts are volumetrically the most important rock type sampled at mid-ocean ridges. For this reason, our present understanding of upper mantle dynamics and the chemical evolution of the earth is strongly influenced by the study of mid-ocean ridge basalts (MORB). However, MORB are aggregates of polybarically generated small melt increments that can undergo a variety of physical and chemical processes during their ascent and consequently affect their derivative geochemical composition. Therefore, MORB do not represent “direct” windows to the underlying upper mantle. Abyssal peridotites, upper mantle rocks recovered from the ocean floor, are the residual complement to MORB melting and provide essential information on melt extraction from the upper mantle. In this study, abyssal peridotites are examined to address these overarching questions posed by previous studies of MORB: How are basaltic melts formed in the mantle, how are they extracted from the mantle and what physical and chemical processes control mantle melting? The number of studies on abyssal peridotites is small compared to those on basalts, in part because seafloor exposures of abyssal peridotites are relatively rare. For this reason, abyssal peridotite characteristics need to be considered in the context of subaerially exposed peridotites associated with ophiolites, orogenic peridotite bodies and basalt-hosted xenoliths. However, orogenic peridotite bodies are mainly associated with passive continental margins, most ophiolites are formed in supra-subduction zone settings, and peridotite xenoliths are often contaminated by their host magma. Therefore, studies of abyssal peridotites are essential to understanding the primary characteristics of the oceanic upper mantle free from the influence of continental rifting, subduction and tectonic emplacement processes. Nevertheless, numerous processes such as melt stagnation and cooling-induced, inter-mineral exchange can affect residual abyssal peridotite compositions after the cessation of melting. The aim of this study is to address these post-melting modifications of abyssal peridotites from a petrological-geochemical perspective. The samples in this study were dredged along the axis of the ultraslow-spreading Gakkel Ridge in the Arctic Ocean within the “Sparsely Magmatic Zone”, a 100 km ridge section where only mantle rocks are exposed. During two expeditions (ARK XVII-2 in 2001 and ARK XX-2 in 2004), exceptionally fresh peridotites were recovered. The boulders and cobbles collected cover a range of mantle rock compositions, with most characterized as plagioclase-free spinel peridotites or plagioclase- spinel peridotites. This thesis investigates melt stagnation and cooling processes in the upper mantle and is divided into two parts. The first part focuses on processes in the stability field of spinel peridotites (>10 kb) such as melt refertilization and cooling related trace element exchange, while the second part investigates processes in the stability field of plagioclase peridotites (< 10 kb) such as reactive melt migration and melt stagnation. The dissertation chapters are organized to follow the theoretical ascent of a mantle parcel upwelling beneath the location where the samples were collected.