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.

Oceanic crust and island arc formation in Central Asia during late Neoproterozoic times - evidence from petrological and geochemical studies

Die vorliegende Arbeit behandelt die Entwicklung des 570 Ma alten, neoproterozoischen Agardagh - Tes-Chem Ophioliths (ATCO) in Zentralasien. Dieser Ophiolith liegt südwestlich des Baikalsees (50.5° N, 95° E) und wurde im frühen Stadium der Akkretion des Zentralasiatischen Mobilgürtels auf den nordwestlichen Rand des Tuvinisch-Mongolischen Mikrokontinentes aufgeschoben. Bei dem Zentralasiatische Mobilgürtel handelt es sich um einen riesigen Akkretions-Subduktionskomplex, der heute das größte zusammenhängende Orogen der Erde darstellt. Im Rahmen dieser Arbeit wurden eine Reihe plutonischer und vulkanischer Gesteine, sowie verschiedene Mantelgesteine des ATCO mittels mikroanalytischer und geochemischer Verfahren untersucht (Elektronenstrahlmikrosonde, Ionenstrahlmikrosonde, Spurenelement- und Isotopengeochemie). Die Auswertung dieser Daten ermöglichte die Entwicklung eines geodynamisch-petrologischen Modells zur Entstehung des ATCO. Die vulkanischen Gesteine lassen sich aufgrund ihrer Spurenelement- und Isotopenzusammensetzung in inselbogenbezogene und back-arc Becken bezogene Gesteine (IA-Gesteine und BAB-Gesteine) unterscheiden. Darüber hinaus gibt es eine weitere, nicht eindeutig zuzuordnende Gruppe, die hauptsächlich mafische Gänge umfasst. Der grösste Teil der untersuchen Vulkanite gehört zur Gruppe der IA-Gesteine. Es handelt sich um Al-reiche Basalte und basaltische Andesite, welche aus einem evolvierten Stammmagma mit Mg# 0.60, Cr ~ 180 µg/g und Ni ~ 95 µg/g hauptsächlich durch Klinopyroxenfraktionierung entstanden sind. Das Stammmagma selbst entstand durch Fraktionierung von ca. 12 % Olivin und geringen Anteilen von Cr-Spinell aus einer primären, aus dem Mantel abgeleiteten Schmelze. Die IA-Gesteine haben hohe Konzentrationen an inkompatiblen Spurenelementen (leichte-(L)- Seltenerdelement-(SEE)-Konzentrationen etwa 100-fach chondritisch, chondrit-normierte (La/Yb)c von 14.6 - 5.1), negative Nb-Anomalien (Nb/La = 0.37 - 0.62) und niedrige Zr/Nb Verhältnisse (7 - 14) relativ zu den BAB-Gesteinen. Initiale eNd Werte liegen bei etwa +5.5, initiale Bleiisotopenverhältnisse sind: 206Pb/204Pb = 17.39 - 18.45, 207Pb/204Pb = 15.49 - 15.61, 208Pb/204Pb = 37.06 - 38.05. Die Anreicherung lithophiler inkompatibler Spurenelemente (LILE) in dieser Gruppe ist signifikant (Ba/La = 11 - 130) und zeigt den Einfluss subduzierter Komponenten an. Die BAB-Gesteine repräsentieren Schmelzen, die sehr wahrscheinlich aus der gleichen Mantelquelle wie die IA-Gesteine stammen, aber durch höhere Aufschmelzgrade (8 - 15 %) und ohne den Einfluss subduzierter Komponenten entstanden sind. Sie haben niedrigere Konzentrationen an inkompatiblen Spurenelementen, flache SEE-Muster ((La/Yb)c = 0.6 - 2.4) und höhere initiale eNd Werte zwischen +7.8 und +8.5. Nb Anomalien existieren nicht und Zr/Nb Verhältnisse sind hoch (21 - 48). Um die geochemische Entwicklung der vulkanischen Gesteine des ATCO zu erklären, sind mindestens drei Komponenten erforderlich: (1) eine angereicherte, ozeaninselbasalt-ähnliche Komponente mit hoher Nb Konzentration über ~ 30 µg/g, einem niedrigen Zr/Nb Verhältnis (ca. 6.5), einem niedrigen initialen eNd Wert (um 0), aber mit radiogenen 206Pb/204Pb-, 207Pb/204Pb- und 208Pb/204Pb-Verhältnissen; (2) eine N-MORB ähnliche back-arc Becken Komponente mit flachem SEE-Muster und einem hohen initialen eNd Wert von mindestens +8.5, und (3) eine Inselbogen-Komponente aus einer verarmten Mantelquelle, welche durch die abtauchende Platte geochemisch modifiziert wurde. Die geochemische Entstehung der ATCO Vulkanite lässt sich dann am besten durch eine Kombination aus Quellenkontamination, fraktionierte Kristallisation und Magmenmischung erklären. Geodynamisch gesehen entstand der ATCO sehr wahrscheinlich in einem intraozeanischen Inselbogen - back-arc System. Bei den untersuchten Plutoniten handelt es sich um ultramafische Kumulate (Wehrlite und Pyroxenite) sowie um gabbroische Plutonite (Olivin-Gabbros bis Diorite). Die geochemischen Charakteristika der mafischen Plutonite sind deutlich unterschiedlich zu denen der vulkanischen Gesteine, weshalb sie sehr wahrscheinlich ein späteres Entwicklungsstadium des ATCO repräsentieren. Die Spurenelement-Konzentrationen in den Klinopyroxenen der ultramafischen Kumulate sind extrem niedrig, mit etwa 0.1- bis 1-fach chondritischen SEE-Konzentrationen und mit deutlich LSEE-verarmten Mustern ((La/Yb)c = 0.27 - 0.52). Berechnete Gleichgewichtsschmelzen der ultramafischen Kumulate zeigen grosse Ähnlichkeit zu primären boninitischen Schmelzen. Die primären Magmen waren daher boninitischer Zusammensetzung und entstanden in dem durch vorausgegangene Schmelzprozesse stark verarmten Mantelkeil über einer Subduktionszone. Niedrige Spurenelement-Konzentrationen zeigen einen geringen Einfluss der abtauchenden Platte an. Die Spurenelement-Konzentrationen der Gabbros sind ebenfalls niedrig, mit etwa 0.5 - 10-fach chondritischen SEE-Konzentrationen und mit variablen SEE-Mustern ((La/Yb)c = 0.25 - 2.6). Analog zu den Vulkaniten der IA-Gruppe haben alle Gabbros eine negative Nb-Anomalie mit Nb/La = 0.01 - 0.31. Die initialen eNd Werte der Gabbros variieren zwischen +4.8 und +7.1, mit einem Mittelwert von +5.9, und sind damit identisch mit denen der IA-Vulkanite. Bei den untersuchten Mantelgesteinen handelt es sich um teilweise serpentinisierte Dunite und Harzburgite, die alle durch hohe Mg/Si- und niedrige Al/Si-Verhältnisse gekennzeichnet sind. Dies zeigt einen refraktären Charakter an und steht in guter Übereinstimmung mit den hohen Cr-Zahlen (Cr#) der Spinelle (bis zu Cr# = 0.83), auf deren Basis der Aufschmelzgrad der residuellen Mantelgesteine berechnet wurde. Dieser beträgt etwa 25 %. Die geochemische Zusammensetzung und die petrologischen Daten der Ultramafite und Gabbros lassen sich am besten erklären, wenn man für die Entstehung dieser Gesteine einen zweistufigen Prozess annimmt. In einer ersten Stufe entstanden die ultramafischen Kumulate unter hohem Druck in einer Magmenkammer an der Krustenbasis, hauptsächlich durch Klinopyroxen-Fraktionierung. Bei dieser Magmenkammer handelte es sich um ein offenes System, dem von unten laufend neue Schmelze zugeführt wurde, und aus dem im oberen Bereich evolviertere Schmelzen geringerer Dichte entwichen. Diese evolvierten Schmelzen stiegen in flachere krustale Bereiche auf und bildeten dort meist isolierte Intrusionskörper. Diese Intrusionskörper erstarrten ohne Magmen-Nachschub, weshalb petrographisch sehr unterschiedliche Gesteine entstehen konnten. Eine geochemische Modifikation der abkühlenden Schmelzen erfolgte allerdings durch die Assimilation von Nebengestein. Da innerhalb der Gabbros keine signifikante Variation der initalen eNd Werte existiert, handelte es sich bei dem assimilierten Material hauptsächlich um vulkanische Gesteine des ATCO und nicht um ältere, möglicherweise kontinentale Kruste.

The pre-alpine evolution of the basement of the Pelagonian Zone and the Vardar Zone, Greece

The Pelagonian Zone and the Vardar Zone in Greece represent the western part of the Hellenide hinterland (Internal Hellenides). While the Pelagonian Zone comprises predominantly crystalline basement and sedimentary cover rocks, the Vardar Zone has long been regarded as an ophiolite-decorated suture zone separating the Pelagonian Zone from the Serbo-Macedonian Massif to the east. Felsic basement rocks from both areas, with the main focus put on the Pelagonian Zone, were dated in order to identify the major crust-forming episodes and to improve the understanding of the evolutionary history of the region. The interpretation of the single-zircon geochronology results was aided by geochemical investigations. The majority of the basement rocks from the Pelagonian Zone yielded Permo-Carboniferous intrusion ages around 300 Ma, underlining the importance of this crust-forming event for the Internal Hellenides of Greece. Geochemically these basement rocks are classified as subduction-related granitoids, which formed in an active continental margin setting. An important result was the identification of a Precambrian crustal unit within the crystalline basement of the Pelagonian Zone. Orthogneisses from the NW Pelagonian Zone yielded Neoproterozoic ages of c. 700 Ma and are so far the oldest known rocks in Greece. These basement rocks, which are also similar to active margin granitoids, were interpreted as remnants of a terrane, the Florina Terrane, which can be correlated to a Pan-African or Cadomian arc. Since the gneisses contain inherited zircons of Middle to Late Proterozoic ages, the original location of the Florina Terrane was probably at the northwestern margin of Gondwana. In the Vardar Zone an important phase of Upper Jurassic felsic magmatism is documented by igneous formation ages ranging from 155 to 164 Ma. The chemical and isotopic composition of these rocks is also in accord with their formation in a volcanic-arc setting at an active continental margin. Older continental material incorporated in the Vardar Zone is documented by 319-Ma-old gneisses and by inherited zircons of mainly Middle Palaeozoic ages. The prevalence of subduction-related igneous rocks indicates that arc formation and accretion orogeny were the most important processes during the evolution of this part of the Internal Hellenides. The geochronological results demonstrate that most of the Pelagonian Zone and the Vardar Zone crystalline basement formed during distinct pre-Alpine episodes at c. 700, 300 and 160 Ma with a predominance of the Permo-Carboniferous magmatic phase.

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.

Cut-first branch-and-price-second for the capacitated arc-routing problem

This paper presents the first full-fledged branch-and-price (bap) algorithm for the capacitated arc-routing problem (CARP). Prior exact solution techniques either rely on cutting planes or the transformation of the CARP into a node-routing problem. The drawbacks are either models with inherent symmetry, dense underlying networks, or a formulation where edge flows in a potential solution do not allow the reconstruction of unique CARP tours. The proposed algorithm circumvents all these drawbacks by taking the beneficial ingredients from existing CARP methods and combining them in a new way. The first step is the solution of the one-index formulation of the CARP in order to produce strong cuts and an excellent lower bound. It is known that this bound is typically stronger than relaxations of a pure set-partitioning CARP model.rnSuch a set-partitioning master program results from a Dantzig-Wolfe decomposition. In the second phase, the master program is initialized with the strong cuts, CARP tours are iteratively generated by a pricing procedure, and branching is required to produce integer solutions. This is a cut-first bap-second algorithm and its main function is, in fact, the splitting of edge flows into unique CARP tours.

On arc-routing problems

Das Basisproblem von Arc-Routing Problemen mit mehreren Fahrzeugen ist das Capacitated Arc-Routing Problem (CARP). Praktische Anwendungen des CARP sind z.B. in den Bereichen Müllabfuhr und Briefzustellung zu finden. Das Ziel ist es, einen kostenminimalen Tourenplan zu berechnen, bei dem alle erforderlichen Kanten bedient werden und gleichzeitig die Fahrzeugkapazität eingehalten wird. In der vorliegenden Arbeit wird ein Cut-First Branch-and-Price Second Verfahren entwickelt. In der ersten Phase werden Schnittebenen generiert, die dem Master Problem in der zweiten Phase hinzugefügt werden. Das Subproblem ist ein kürzeste Wege Problem mit Ressourcen und wird gelöst um neue Spalten für das Master Problem zu liefern. Ganzzahlige CARP Lösungen werden durch ein neues hierarchisches Branching-Schema garantiert. Umfassende Rechenstudien zeigen die Effektivität dieses Algorithmus. Kombinierte Standort- und Arc-Routing Probleme ermöglichen eine realistischere Modellierung von Zustellvarianten bei der Briefzustellung. In dieser Arbeit werden jeweils zwei mathematische Modelle für Park and Loop und Park and Loop with Curbline vorgestellt. Die Modelle für das jeweilige Problem unterscheiden sich darin, wie zulässige Transfer Routen modelliert werden. Während der erste Modelltyp Subtour-Eliminationsbedingungen verwendet, werden bei dem zweiten Modelltyp Flussvariablen und Flusserhaltungsbedingungen eingesetzt. Die Rechenstudie zeigt, dass ein MIP-Solver den zweiten Modelltyp oft in kürzerer Rechenzeit lösen kann oder bei Erreichen des Zeitlimits bessere Zielfunktionswerte liefert.