Geological setting of the Guelb Moghrein Fe oxide-Cu-Au-Co mineralization, Akjoujt area, Mauritania

The Guelb Moghrein Fe oxide-Cu-Au-Co deposit is located at the western boundary of the West African craton in NW Mauritania. The wall rocks to the mineralization represent a meta-volcanosedimentary succession typical of Archaean greenstone belts. Two types of meta-volcanic rocks are distinguished: (1) volcanoclastic rocks of rhyodacite-dacite composition (Sainte Barbe volcanic unit), which form the stratigraphic base; (2) tholeiitic andesites-basalts (Akjoujt meta-basalt unit). The trace element signature of both types is characteristic of a volcanic arc setting. A small meta-pelitic division belongs to the Sainte Barbe volcanic unit. A meta-carbonate body, which contains the mineralization, forms a tectonic lens in the Akjoujt meta-basalt unit. It can be defined by the high X(mg) (=36) of Fe-Mg carbonate, the REE pattern and the delta(13)C values of -18 to -17 parts per thousand as a marine precipitate similar to Archaean banded iron formation (BIF). Additionally, small slices of Fe-Mg clinoamphibole-chlorite schist in the meta-carbonate show characteristics of marine shale. This assemblage, therefore, does not represent an alteration product, but represents an iron formation unit deposited on a continental shelf, which probably belongs to the Lembeitih Formation. The hydrothermal mineralization at 2492 Ma was contemporaneous with regional D(2) thrusting of the Sainte Barbe volcanic unit and imbrications of the meta-carbonate in the upper greenschist facies. This resulted in the formation of an ore breccia in the meta-carbonate, which is enriched in Fe, Ni, Co, Cu, Bi, Mo, As and Au. Massive sulphide ore breccia contains up to 20 wt% Cu. The ore fluid was aqueous-carbonic in nature and either changed its composition from a Mg-rich oxidizing to an Fe-rich reducing fluid or the two fluid types mixed at the trap site. All lithologies at Guelb Moghrein were deformed by D(3) thrusting to the east in the lower greenschist facies. The mobility of REE in the retrogressed rocks explains the formation of a second generation of hydrothermal monazite, which was dated at c. 1742 Ma. Archaean rocks of the West African craton extend to the west to Guelb Moghrein. The active continental margin was deformed and mineralized in the Late Archaean-Early Proterozoic and again reactivated in the Mid-Proterozoic and Westphalian, showing that the western boundary of the craton was reactivated several times.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.

Origin of mineralizing fluids of the sediment-hosted Navachab gold mine, Namibia: Constraints from stable (O, H, C, S) isotopes

The Navachab gold mine in the Damara belt of central Namibia is characterized by a polymetallic Au-Bi-As-Cu-Ag ore assemblage, including pyrrhotite, chalcopyrite, sphalerite, arsenopyrite, bismuth, gold, bismuthinite, and bismuth tellurides. Gold is hosted by quartz sulfide veins and semimassive sulfide lenses that are developed in a near-vertical sequence of shelf-type metasedimentary rocks, including marble, calcsilicate rock, and biotite schist. The sequence has been intruded by abundant syntectonic lamprophyre, aplite, and pegmatite dikes, documenting widespread igneous activity coeval with mineralization. The majority of quartz from the veins has delta(18)O values of 14 to 15 per mil (V-SMOW). The total variations in delta(18)O values of the biotite schist and calcsilicate rock are relatively small (12-14 parts per thousand), whereas the marble records steep gradients in delta(18)O values (17-21 parts per thousand), the lowest values being recorded at the vein margins. Despite this, there is no correlation between delta(18)O and delta(13)C values and the carbonate content of the rocks, indicating that fluid-rock interaction alone cannot explain the isotopic gradients. In addition, the marble records increased delta(13)C values at the contact to the veins, possibly related to a change in the physicochemical conditions during fluid-rock interaction. Gold is interpreted to have precipitated in equilibrium with metamorphic find (delta(18)O 12-14 parts per thousand; delta D = -40 to -60 parts per thousand) at peak metamorphic conditions of ca. 550 degrees C and 2 kbars, consistent with isotopic fractionations between coexisting calcite, garnet, and clinopyroxene in the alteration halos. The most likely source of the mineralizing fluid was a midcrustal fluid in equilibrium with the Damaran metapelites that underwent prograde metamorphism at amphibolite- to granulite-facies grades. Although there is no isotopic evidence for the contribution of magmatic fluids, they may have been important in contributing to the overall hydraulic regime and high apparent geothermal gradients (ca. 80 degrees C/km(-1)) in the mine area.

The oceanic base of slope record of the Permian-Triassic crisis: view from Tethys (Oman)

The Oman Mountains provide some of the best sections of Permian and Triassic sediments from ocean sea floor to base-of-slope environments related to the distal South Tethyan margin. The central part of the range exposes the Buday'ah section of oceanic sediments in the so-called "Hawasina allochtons". The locality of Wadi Maqam in the north-western part of the Oman Mountains is among places where the thick Permian-Triassic base-of-slope sediments is exposed (Baud et al., 2001). Overlying 400 m of middle Permian limestones and dolomites, the upper Permian sediments consist of 50 m of ≈ 10 cm thick beds of cherts and dolomites rich in sponge spicules. The top of the Permian units is well bioturbated lime mudstone-wackestone, devoid of cherts and dated as late Changhsingian (Krystyn in Richoz et al., 2005). The boundary yellow shales are overlain by very thinly bedded, laminated microbial platy lime mudstone with H. parvus. The dramatic loss of the burrowing infauna indicates the appearance of oxygen-poor water. These Induan sediments are about 25 m thick and show at the top the first calcirudites, commonly clast-supported (edge-wise conglomerates), and are characterized by tabular clasts representing the sub- in situ reworking of the laminated, platy calcilutite. The very thick Smithian overlying litho-unit (up to 900 m) marks the onset on the base-of-slope of a deep-marine basin in which carbonate submarine fan deposits developed This very thick unit consists essentially of platy limestones, calcarenites and calcirudites. It comprises mainly grey-beige calcilutite, laminated and flaggy, interbedded with sparse beds of fine-grained calcarenite in cm beds. Channelized beds of intraformational calcirudite are also part of this succession which constitutes the greater part of the outcrop available. During the Spathian to Anisian, the sedimentation changes to terrigenous mudstone and siltstone that ended with Ladinian radiolarites.

The Permo-Triassic boundary in the Antalya Nappes (Western Taurides - Turkey)

The Antalya nappes (western Taurides-Turkey)*consist of several tectonic units which document the*Southern Neotethyan paleomargin from the Arabo-African*shallow shelf to the oceanic crust.*The Kerner Gorge Units (Upper Antalya nappes)*show a full stratigraphical succession from Ordovician*to Late Cretaceous. A carbonate platform regime*appeared during Late Permian times and existed up*to the early Middle Triassic. For detailed investigations*on the Permo-Triassic boundary, two lithostratigraphic*profiles have been selected: the Curuk*dag and the Kerner Gorge sections.*The main results presented in this paper are : 1) the Late Permian Pamucak Formation (Midian-*Dzhulfian) consists of calcareous algae-foraminiferae*bearing black limestones, locally rich in*brachiopods, crinoids and bryozoae : 2) this black limestone is overlain by an oolitic*grainstone. In the Curuk dag section, the oolitic horizon*is capped by a thin level of calcrete type; emersive*conditions are also inferred by a strong diagenetic*change within the oolitic deposit : 3) the first Early Triassic fossils, appearing*within or above the oolitic grainstone, are microforaminifera*and Pseudoclaraia wangi (late Griesbachian*in age). In the Curuk dag, a rich Early Triassic*microforaminifera association seems linked to a microbiallite*boundstone facies : 4) the overlying unfossiliferous lime mudstone,*the oolitical thick bedded grainstone, the variegated*marly limestone and the vermicular limestone facies*are present. not only in southern Turkey but seem*constant through ' the entire peri arabo-african platform.*Striking similarities appear between the studied*profiles and the Bellerophon-Werfen succession in the*Southern Alps..

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.