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.

Variation of stable silicon isotopes: Analytical developments and applications in Precambrian geochemistry

Die vorliegende Dissertation behandelt die Gesamtgesteinsanalyse stabiler Siliziumisotope mit Hilfe einer „Multi Collector-ICP-MS“. Die Analysen fanden in Kooperation mit dem „Royal Museum for Central Africa“ in Belgien statt. Einer der Schwerpunkte des ersten Kapitels ist die erstmalige Analyse des δ30Si –Wertes an einem konventionellen Nu PlasmaTM „Multi-Collector ICP-MS“ Instrument, durch die Eliminierung der den 30Si “peak” überlagernden 14N16O Interferenz. Die Analyse von δ30Si wurde durch technische Modifikationen der Anlage erreicht, welche eine höherer Massenauflösung ermöglichten. Die sorgsame Charakterisierung eines adäquaten Referenzmaterials ist unabdingbar für die Abschätzung der Genauigkeit einer Messung. Die Bestimmung der „U.S. Geological Survey“ Referenzmaterialien bildet den zweiten Schwerpunkt dieses Kapitales. Die Analyse zweier hawaiianischer Standards (BHVO-1 and BHVO-2), belegt die präzise und genaue δ30Si Bestimmung und bietet Vergleichsdaten als Qualitätskontrolle für andere Labore. Das zweite Kapitel befasst sich mit kombinierter Silizium-/Sauerstoffisotope zur Untersuchung der Entstehung der Silizifizierung vulkanischer Gesteine des „Barberton Greenstone Belt“, Südafrika. Im Gegensatz zu heute, war die Silizifizierung der Oberflächennahen Schichten, einschließlich der „Chert“ Bildung, weitverbreitete Prozesse am präkambrischen Ozeanboden. Diese Horizonte sind Zeugen einer extremen Siliziummobilisierung in der Frühzeit der Erde. Dieses Kapitel behandelt die Analyse von Silizium- und Sauerstoffisotopen an drei unterschiedlichen Gesteinsprofilen mit unterschiedlich stark silizifizierten Basalten und überlagernden geschichteten „Cherts“ der 3.54, 3.45 und 3.33 Mill. Jr. alten Theespruit, Kromberg und Hooggenoeg Formationen. Siliziumisotope, Sauerstoffisotope und die SiO2-Gehalte demonstrieren in allen drei Gesteinsprofilen eine positive Korrelation mit dem Silizifizierungsgrad, jedoch mit unterschiedlichen Steigungen der δ30Si-δ18O-Verhältnisse. Meerwasser wird als Quelle des Siliziums für den Silizifizierungsprozess betrachtet. Berechnungen haben gezeigt, dass eine klassische Wasser-Gestein Wechselwirkung die Siliziumisotopenvariation nicht beeinflussen kann, da die Konzentration von Si im Meerwasser zu gering ist (49 ppm). Die Daten stimmen mit einer Zwei-Endglieder-Komponentenmischung überein, mit Basalt und „Chert“ als jeweilige Endglieder. Unsere gegenwärtigen Daten an den „Cherts“ bestätigen einen Anstieg der Isotopenzusammensetzung über der Zeit. Mögliche Faktoren, die für unterschiedliche Steigungen der δ30Si-δ18O Verhältnisse verantwortlich sein könnten sind Veränderungen in der Meerwasserisotopie, der Wassertemperatur oder sekundäre Alterationseffekte. Das letzte Kapitel beinhaltet potentielle Variationen in der Quellregion archaischer Granitoide: die Si-Isotopen Perspektive. Natriumhaltige Tonalit-Trondhjemit-Granodiorit (TTG) Intrusiva repräsentieren große Anteile der archaischen Kruste. Im Gegensatz dazu ist die heutige Kruste kaliumhaltiger (GMS-Gruppe: Granit-Monzonite-Syenite). Prozesse, die zu dem Wechsel von natriumhaltiger zu kaliumhaltiger Kruste führten sind die Thematik diesen Kapitels. Siliziumisotopenmessungen wurden hier kombiniert mit Haupt- und Spurenelementanalysen an unterschiedlichen Generationen der 3.55 bis 3.10 Mill. Yr. alten TTG und GMS Intrusiva aus dem Arbeitsgebiet. Die δ30Si-Werte in den unterschiedlichen Plutonit Generationen zeigen einen leichten Anstieg der Isotopie mit der Zeit, wobei natriumhaltige Intrusiva die niedrigste Si-Isotopenzusammensetzung aufweisen. Der leichte Anstieg in der Siliziumisotopenzusammensetzung über die Zeit könnte auf unterschiedliche Temperaturbedingungen in der Quellregion der Granitoide hinweisen. Die Entstehung von Na-reichen, leichten d30Si Granitoiden würde demnach bei höheren Temperaturen erfolgen. Die Ähnlichkeit der δ30Si-Werte in archaischen K-reichen Plutoniten und phanerozoischen K-reichen Plutoniten wird ebenfalls deutlich.

Miocene basin development and volcanism along a strike-slip to flat-slab subduction transition: Stratigraphy, geochemistry, and geochronology of the central Wrangell volcanic belt, Yakutat-North America collision zone

New geochronologic, geochemical, sedimentologic, and compositional data from the central Wrangell volcanic belt (WVB) document basin development and volcanism linked to subduction of overthickened oceanic crust to the northern Pacific plate margin. The Frederika Formation and overlying Wrangell Lavas comprise >3 km of sedimentary and volcanic strata exposed in the Wrangell Mountains of south-central Alaska (United States). Measured stratigraphic sections and lithofacies analyses document lithofacies associations that reflect deposition in alluvial-fluvial-lacustrine environments routinely influenced by volcanic eruptions. Expansion of intrabasinal volcanic centers prompted progradation of vent-proximal volcanic aprons across basinal environments. Coal deposits, lacustrine strata, and vertical juxtaposition of basinal to proximal lithofacies indicate active basin subsidence that is attributable to heat flow associated with intrabasinal volcanic centers and extension along intrabasinal normal faults. The orientation of intrabasinal normal faults is consistent with transtensional deformation along the Totschunda-Fairweather fault system. Paleocurrents, compositional provenance, and detrital geochronologic ages link sediment accumulation to erosion of active intrabasinal volcanoes and to a lesser extent Mesozoic igneous sources. Geochemical compositions of interbedded lavas are dominantly calc-alkaline, range from basaltic andesite to rhyolite in composition, and share geochemical characteristics with Pliocene-Quaternary phases of the western WVB linked to subduction-related magmatism. The U/Pb ages of tuffs and Ar-40/Ar-39 ages of lavas indicate that basin development and volcanism commenced by 12.5-11.0 Ma and persisted until at least ca. 5.3 Ma. Eastern sections yield older ages (12.5-9.3 Ma) than western sections (9.6-8.3 Ma). Samples from two western sections yield even younger ages of 5.3 Ma. Integration of new and published stratigraphic, geochronologic, and geochemical data from the entire WVB permits a comprehensive interpretation of basin development and volcanism within a regional tectonic context. We propose a model in which diachronous volcanism and transtensional basin development reflect progressive insertion of a thickened oceanic crustal slab of the Yakutat microplate into the arcuate continental margin of southern Alaska coeval with reported changes in plate motions. Oblique northwestward subduction of a thickened oceanic crustal slab during Oligocene to Middle Miocene time produced transtensional basins and volcanism along the eastern edge of the slab along the Duke River fault in Canada and subduction-related volcanism along the northern edge of the slab near the Yukon-Alaska border. Volcanism and basin development migrated progressively northwestward into eastern Alaska during Middle Miocene through Holocene time, concomitant with a northwestward shift in plate convergence direction and subduction collision of progressively thicker crust against the syntaxial plate margin.

Provenance discrimination of Lower Cretaceous synrift sandstones (eastern Iberian Chain Spain) Constraints from detrital modes heavy minerals and geochemistry

Petrography, geochemical whole-rock composition, and chemical analyses of tourmaline were performed in order to determine the source areas of Lower Cretaceous Mora, El Castellar, and uppermost Camarillas Formation sandstones from the Iberian Chain, Spain. Sandstones were deposited in intraplate subbasins, which are bound by plutonic and volcanic rocks of Permian, Triassic, and Jurassic age, Paleozoic metamorphic rocks, and Triassic sedimentary rocks. Modal analyses together with petrographic and cathodoluminescence observations allowed us to define three quartz-feldspathic petrofacies and recognize diagenetic processes that modified the original framework composition. Results from average restored petrofacies are: Mora petrofacies = P/F >1 and Q(r)70 F(r)22 R(r)9; El Castellar petrofacies = P/F >1 and Q(r)57 F(r)25 R(r)18; and Camarillas petrofacies = P/F ∼ zero and Q(r)64 F(r)28 R(r)7 (P—plagioclase; F—feldspar; Q—quartz; R—rock fragments; r—restored composition). Trace-element and rare earth element abundances of whole-rock analyses discriminate well between the three petrofacies based on: (1) the Rb concentration, which is indicative of the K content and reflects the amount of K-feldspar modal abundance, and (2) the relative modal abundance of heavy minerals (tourmaline, zircon, titanite, and apatite), which is reproduced by the elements hosted in the observed heavy mineral assemblage (i.e., B and Li for tourmaline; Zr, Hf, and Ta for zircon; Ti, Ta, Nb, and their rare earth elements for titanite; and P, Y, and their rare earth elements for apatite). Tourmaline chemical composition for the three petrofacies ranges from Fe-tourmaline of granitic to Mg-tourmaline of metamorphic origin. The three defined petrofacies suggest a mixed provenance from plutonic and metamorphic source rocks. However, a progressively major influence of granitic source rocks was detected from the lowermost Mora petrofacies toward the uppermost Camarillas petrofacies. This provenance trend is consistent with the uplift and erosion of the Iberian Massif, which coincided with the development of the latest Berriasian synrift regional unconformity and affected all of the Iberian intraplate basins. The uplifting stage of Iberian Massif pluton caused a significant dilution of Paleozoic metamorphic source areas, which were dominant during the sedimentation of the lowermost Mora and El Castellar petrofacies. The association of petrographic data with whole-rock geochemical compositions and tourmaline chemical analysis has proved to be useful for determining source area characteristics, their predominance, and the evolution of source rock types during the deposition of quartz-feldspathic sandstones in intraplate basins. This approach ensures that provenance interpretation is consistent with the geological context.

Geochemistry and Stable Isotopes of Surface Water and Groundwater in the Continental Pit in Butte, Montana, USA

The Continental porphyry Cu‐Mo mine, located 2 km east of the famous Berkeley Pit lake of Butte, Montana, contains two small lakes that vary in size depending on mining activity. In contrast to the acidic Berkeley Pit lake, the Continental Pit waters have near-neutral pH and relatively low metal concentrations. The main reason is geological: whereas the Berkeley Pit mined highly‐altered granite rich in pyrite with no neutralizing potential, the Continental Pit is mining weakly‐altered granite with lower pyrite concentrations and up to 1‐2% hydrothermal calcite. The purpose of this study was to gather and interpret information that bears on the chemistry of surface water and groundwater in the active Continental Pit. Pre‐existing chemistry data from sampling of the Continental Pit were compiled from the Montana Bureau of Mines and Geology and Montana Department of Environmental Quality records. In addition, in March of 2013, new water samples were collected from the mine’s main dewatering well, the Sarsfield well, and a nearby acidic seep (Pavilion Seep) and analyzed for trace metals and several stable isotopes, including dD and d18O of water, d13C of dissolved inorganic carbon, and d34S of dissolved sulfate. In December 2013, several soil samples were collected from the shore of the frozen pit lake and surrounding area. The soil samples were analyzed using X‐ray diffraction to determine mineral content. Based on Visual Minteq modeling, water in the Continental Pit lake is near equilibrium with a number of carbonate, sulfate, and molybdate minerals, including calcite, dolomite, rhodochrosite (MnCO3), brochantite (CuSO4·3Cu(OH)2), malachite (Cu2CO3(OH)2), hydrozincite (Zn5(CO3)2(OH)6), gypsum, and powellite (CaMoO4). The fact that these minerals are close to equilibrium suggests that they are present on the weathered mine walls and/or in the sediment of the surface water ponds. X‐Ray Diffraction (XRD) analysis of the pond “beach” sample failed to show any discrete metal‐bearing phases. One of the soil samples collected higher in the mine, near an area of active weathering of chalcocite‐rich ore, contained over 50% chalcanthite (CuSO4·5H2O). This water‐soluble copper salt is easily dissolved in water, and is probably a major source of copper to the pond and underlying groundwater system. However, concentrations of copper in the latter are probably controlled by other, less‐soluble minerals, such as brochantite or malachite. Although the acidity of the Pavilion Seep is high (~ 11 meq/L), the flow is much less than the Sarsfield Well at the current time. Thus, the pH, major and minor element chemistry in the Continental Pit lakes are buffered by calcite and other carbonate minerals. For the Continental Pit waters to become acidic, the influx of acidic seepage (e.g., Pavilion Seep) would need to increase substantially over its present volume.

Weathering geochemistry and Sr-Nd fingerprints of equatorial upper Nile and Congo muds

This study investigates processes of sediment generation in equatorial central Africa. An original, complete and integrated mineralogical-geochemical database on silt-sized sediments derived from different parent rocks (basalt, granite, gneiss, metapsammite, sandstone) along the East African Rift from 5°S in Tanzania to 5°N in Sudan is presented and used to assess the incidence of diverse factors controlling sediment composition (source-rock lithology, geomorphology, hydraulic sorting, grain size, recycling), with particular emphasis on chemical weathering.

Geochronology and isotope geochemistry of Eocene dykes intruding the Ladakh Batholith

In order to determine the extent and timing of dyke formation in the Ladakh Batholith we examined about 30 mostly andesitic dykes intruding the Ladakh batholith in a ca. 50 km wide area to the west of Leh (NW India). The dykes in the east of the area trend E-NE and those in the west trend N-NW. The difference in orientation is also evident in the petrography and isotopic signatures. The eastern dykes contain corroded quartz xenocrysts and show negative ε0(Nd) and positive ε0(Sr) values, where as the western dykes do not contain quartz xenocrysts and exhibit positive ε0(Nd) and near-zero ε0(Sr) values. The variability in Sr-Nd isotopes (ε0(Nd) = 3.6 to −9.6, ε0(Sr) = 0.4 to 143) and the quartz xenocrysts can best be explained by (differing degrees of) crustal assimilation of the parent magma of the dykes. Separated minerals from five dykes were dated by 40Ar-39Ar incremental heating: amphibole ages range between 50 and 54 Ma, and one biotite dated both by Rb-Sr and by 40Ar-39Ar gave an age of 45 Ma. One dated pseudotachylyte sample attests to brittle faulting at ca. 54 Ma. The combination of structural field evidence with petrographic, isotopic and geochronological analyses demonstrates that the dykes did not form from a single, progressively differentiating magma chamber, despite having formed in the same tectonic setting around the same time, and that processes such as crustal assimilation and magma mixing/mingling also played a significant role in magma petrogenesis.

Redox front geochemistry and weathering: theory with application to the Osamu Utsumi uranium mine, Poços de Caldas, Brazil

Numerical calculations describing weathering of the Poços de Caldas alkaline complex (Minas Gerais, Brazil) by infiltrating groundwater are carried out for time spans up to two million years in the absence of pyrite, and up to 500,000 years with pyrite present. Deposition of uranium resulting from infiltration of oxygenated, uranium bearing groundwater through the hydrothermally altered phonolitic host rock at the Osamu Utsumi uranium mine is also included in the latter calculation. The calculations are based on the quasi-stationary state approximation to mass conservation equations for pure advective transport. This approximation enables the prediction of solute concentrations, mineral abundances and porosity as functions of time and distance over geologic time spans. Mineral reactions are described by kinetic rate laws for both precipitation and dissolution. Homogeneous equilibrium is assumed to be maintained within the aqueous phase. No other constraints are imposed on the calculations other than the initial composition of the unaltered host rock and the composition of the inlet fluid, taken as rainwater modified by percolation through a soil zone. The results are in qualitative agreement with field observations at the Osamu Utsumi uranium mine. They predict a lateritic cover followed by a highly porous saprolitic zone, a zone of oxidized rock with pyrite replaced by iron-hydroxide, a sharp redox front at which uranium is deposited, and the reduced unweathered host rock. Uranium is deposited in a narrow zone located on the reduced side of the redox front in association with pyrite, in agreement with field observations. The calculations predict the formation of a broad dissolution front of primary kaolinite that penetrates deep into the host rock accompanied by the precipitation of secondary illite. Secondary kaolinite occurs in a saprolitic zone near the surface and in the vicinity of the redox front. Gibbsite forms a bi-modal distribution consisting of a maximum near the surface followed by a thin tongue extending downward into the weathered profile in agreement with field observations. The results are found to be insensitive to the kinetic rate constants used to describe mineral reactions.