Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

Climatic and anthropogenic influence on the stable isotope record from bulk carbonates and ostracodes in Lake Neuchatel, Switzerland, during the last two millennia

Lake Neuchatel is a medium sized, hard-water lake, lacking varved sediments, situated in the western Swiss Lowlands at the foot of the Jura Mountains. Stable isotope data (delta(18)O and delta(13)C) from both bulk carbonate and ostracode calcite in an 81 cm long, radiocarbon-dated sediment core represent the last 1500 years of Lake Neuchatel's environmental history. Comparison between this isotopic and other palaeolimnologic data (mineralogical, geochemical, palynological, etc.) helps to differentiate between anthropogenic and natural factors most recently affecting the lake. An increase in lacustrine productivity (450-650AD ca), inferred from the positive trend in delta(13)C values of bulk carbonate, is related to medieval forest clearances and the associated nutrient budget changes. A negative trend in both the bulk carbonate and ostracode calcite delta(18)O values between approximately 1300 and 1500AD, is tentatively interpreted as due to a cooling in mean air temperature at the transition from the Medieval Warm Period to the Little Ice Age. Negative trends in bulk carbonate delta(18)O and delta(13)C values through the uppermost sediments, which have no equivalent in ostracode calcite isotopic values, are concomitant with the recent onset of eutrophication in the lake. Isotopic disequilibrium during calcite precipitation, probably due to kinetic factors in periods of high productivity is postulated as the mechanism to explain the associated negative isotopic trends, although the effect of a shift of the calcite precipitation towards the warmer months cannot be excluded.

Apparent stable-isotope heterogeneities in gangue carbonates of the Mississippi Valley-type Zn-Pb deposits of San-Vicente, Central Peru

The aim of the present communication is to emphasize that some variations of the measured delta(13)C and delta(18)O values are apparent, and due to analytical interferences caused by the presence of sulfur and organosulfur compounds in the analyzed carbonates. This is particularly relevant for isotopic studies on carbonate-hosted mineral deposits, where the nearly ubiquitous association of the host carbonates with organic matter and sulfides can certainly affect the metallogenetic interpretations. In this work two methods were used to overcome the disturbing effects of sulfides and organic matter: (1) sample pretreatment following the method proposed by Charef and Sheppard (1984), combining the oxidation of organic matter with sodium hypochlorite and trapping of the sulfur species with silver phosphate; and (2) laser-based microprobe extraction. Apparent isotopic variations in sparry dolomite from a single hand sample of zebra ore from the MVT Zn-Pb deposit, San Vicente, central Peru, are as large as 6 parts per thousand delta(13)C and 4 parts per thousand delta(18)O. These variations are reduced to several tenths of a per mil when the samples are pretreated. A careful examination of the effects of treatment with NaOCl and/or Ag3PO4 in relation to the concentration of sulfide inclusions indicates that the main disturbing effects for delta(13)C values are the presence of sulfur species and organic matter, whereas the delta(18)O values are mainly affected by the presence of sulfides. Fine- and medium-grained replacement carbonates from MVT and other sediment-hosted base metal deposits are potentially the most affected during isotope analysis, due to the common presence of organic matter and sulfides. Using in situ laser microprobe techniques, it is possible to determine isotopic variations at a sub-millimeter scale. Our results show that laser extraction analysis allows a more precise sampling of the carbonate minerals, and minimizes contamination of the sample with sulfides and to some extent with intergrown organic matter. However, there is an isotopic shift associated with the laser extraction technique, of the order of 0.5-1 parts per thousand for delta(13)C and delta(18)O values.

Modélisation par automate cellulaire des phénomènes diagénétiques des plateformes carbonatées : calibration et paramétrage à partir de deux cas d'études : l'Urgonien du Vercors (Crétacé inférieur, SE France) et les Calcaires Gris du Mont Compomolon (Lias, NE Italie)

Une fois déposé, un sédiment est affecté au cours de son enfouissement par un ensemble de processus, regroupé sous le terme diagenèse, le transformant parfois légèrement ou bien suffisamment pour le rendre méconnaissable. Ces modifications ont des conséquences sur les propriétés pétrophysiques qui peuvent être positives ou négatives, c'est-à-dire les améliorer ou bien les détériorer. Une voie alternative de représentation numérique des processus, affranchie de l'utilisation des réactions physico-chimiques, a été adoptée et développée en mimant le déplacement du ou des fluides diagénétiques. Cette méthode s'appuie sur le principe d'un automate cellulaire et permet de simplifier les phénomènes sans sacrifier le résultat et permet de représenter les phénomènes diagénétiques à une échelle fine. Les paramètres sont essentiellement numériques ou mathématiques et nécessitent d'être mieux compris et renseignés à partir de données réelles issues d'études d'affleurements et du travail analytique effectué. La représentation des phénomènes de dolomitisation de faible profondeur suivie d'une phase de dédolomitisation a été dans un premier temps effectuée. Le secteur concerne une portion de la série carbonatée de l'Urgonien (Barrémien-Aptien), localisée dans le massif du Vercors en France. Ce travail a été réalisé à l'échelle de la section afin de reproduire les géométries complexes associées aux phénomènes diagénétiques et de respecter les proportions mesurées en dolomite. De plus, la dolomitisation a été simulée selon trois modèles d'écoulement. En effet, la dédolomitisation étant omniprésente, plusieurs hypothèses sur le mécanisme de dolomitisation ont été énoncées et testées. Plusieurs phases de dolomitisation per ascensum ont été également simulées sur des séries du Lias appartenant aux formations du groupe des Calcaire Gris, localisées au nord-est de l'Italie. Ces fluides diagénétiques empruntent le réseau de fracturation comme vecteur et affectent préférentiellement les lithologies les plus micritisées. Cette étude a permis de mettre en évidence la propagation des phénomènes à l'échelle de l'affleurement. - Once deposited, sediment is affected by diagenetic processes during their burial history. These diagenetic processes are able to affect the petrophysical properties of the sedimentary rocks and also improve as such their reservoir capacity. The modelling of diagenetic processes in carbonate reservoirs is still a challenge as far as neither stochastic nor physicochemical simulations can correctly reproduce the complexity of features and the reservoir heterogeneity generated by these processes. An alternative way to reach this objective deals with process-like methods, which simplify the algorithms while preserving all geological concepts in the modelling process. The aim of the methodology is to conceive a consistent and realistic 3D model of diagenetic overprints on initial facies resulting in petrophysical properties at a reservoir scale. The principle of the method used here is related to a lattice gas automata used to mimic diagenetic fluid flows and to reproduce the diagenetic effects through the evolution of mineralogical composition and petrophysical properties. This method developed in a research group is well adapted to handle dolomite reservoirs through the propagation of dolomitising fluids and has been applied on two case studies. The first study concerns a mid-Cretaceous rudist and granular platform of carbonate succession (Urgonian Fm., Les Gorges du Nan, Vercors, SE France), in which several main diagenetic stages have been identified. The modelling in 2D is focused on dolomitisation followed by a dédolomitisation stage. For the second study, data collected from outcrops on the Venetian platform (Lias, Mont Compomolon NE Italy), in which several diagenetic stages have been identified. The main one is related to per ascensum dolomitisation along fractures. In both examples, the evolution of the effects of the mimetic diagenetic fluid on mineralogical composition can be followed through space and numerical time and help to understand the heterogeneity in reservoir properties. Carbonates, dolomitisation, dédolomitisation, process-like modelling, lattice gas automata, random walk, memory effect.