On the potential for CO2 mineral storage in continental flood basalts – PHREEQC batch- and 1D diffusion–reaction simulations

Continental flood basalts (CFB) are considered as potential CO₂ storage sites because of their high reactivity and abundant divalent metal ions that can potentially trap carbon for geological timescales. Moreover, laterally extensive CFB are found in many place in the world within reasonable distances from major CO₂ point emission sources. Based on the mineral and glass composition of the Columbia River Basalt (CRB) we estimated the potential of CFB to store CO₂ in secondary carbonates. We simulated the system using kinetic dependent dissolution of primary basalt-minerals (pyroxene, feldspar and glass) and the local equilibrium assumption for secondary phases (weathering products). The simulations were divided into closed-system batch simulations at a constant CO₂ pressure of 100?bar with sensitivity studies of temperature and reactive surface area, an evaluation of the reactivity of H₂O in scCO₂, and finally 1D reactive diffusion simulations giving reactivity at CO₂ pressures varying from 0 to 100?bar. Although the uncertainty in reactive surface area and corresponding reaction rates are large, we have estimated the potential for CO₂ mineral storage and identified factors that control the maximum extent of carbonation. The simulations showed that formation of carbonates from basalt at 40?C may be limited to the formation of siderite and possibly FeMg carbonates. Calcium was largely consumed by zeolite and oxide instead of forming carbonates. At higher temperatures (60 ?C 100?C), magnesite is suggested to form together with siderite and ankerite. The maximum potential of CO₂ stored as solid carbonates, if CO₂ is supplied to the reactions unlimited, is shown to depend on the availability of pore space as the hydration and carbonation reactions increase the solid volume and clog the pore space. For systems such as in the scCO₂ phase with limited amount of water, the total carbonation potential is limited by the amount of water present for hydration of basalt.     

Caprock interaction with CO2: A laboratory study of reactivity of shale with supercritical CO2 and brine

Crushed rock from two caprock samples, a carbonate-rich shale and a clay-rich shale, were reacted with a mixture of brine and supercritical CO₂ (CO₂–brine) in a laboratory batch reactor, at different temperature and pressure conditions. The samples were cored from a proposed underground CO₂ storage site near the town of Longyearbyen in Svalbard. The reacting fluid was a mixture of 1 M NaCl solution and CO₂ (110 bar) and the water/rock ratio was 20:1. Carbon dioxide was injected into the reactors after the solution had been bubbled with N₂, in order to mimic O₂-depleted natural storage conditions. A control reaction was also run on the clay-rich shale sample, where the crushed rock was reacted with brine (CO₂-free brine) at the same experimental conditions. A total of 8 batch reaction experiments were run at temperatures ranging from 80 to 250 °C and total pressures of 110 bar (∼40 bar for the control experiment). The experiments lasted 1–5 weeks.Fluid analysis showed that the aqueous concentration of major elements (i.e. Ca, Mg, Fe, K, Al) and SiO₂ increased in all experiments. Release rates of Fe and SiO₂ were more pronounced in solutions reacted with CO₂–brine as compared to those reacted with CO₂-free brine. For samples reacted with the CO₂–brine, lower temperature reactions (80 °C) released much more Fe and SiO₂ than higher temperature reactions (150–250 °C). Analysis by SEM and XRD of reacted solids also revealed changes in mineralogical compositions. The carbonate-rich shale was more reactive at 250 °C, as revealed by the dissolution of plagioclase and clay minerals (illite and chlorite), dissolution and re-precipitation of carbonates, and the formation of smectite. Carbon dioxide was also permanently sequestered as calcite in the same sample. The clay-rich shale reacted with CO₂–brine did not show major mineralogical alteration. However, a significant amount of analcime was formed in the clay-rich shale reacted with CO₂-free brine; while no trace of analcime was observed in either of the samples reacted with CO₂–brine.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

Terrestrial and freshwater carbonates in Hoxnian interglacial deposits, UK: micromorphology, stable isotopic composition and palaeoenvironmental significance

MIS 11 is often considered to be the best climatic analogue for the Holocene. Many studies have suggested, however, that it is a period of extreme climate warmth comparable in temperatures to the Middle and Late Pliocene. In Britain deposits of the Hoxnian interglacial are correlated to MIS 11 and multi-proxy techniques can be used to reconstruct the climate of this interglacial. Soil, groundwater and freshwater carbonates are common in Hoxnian deposits and the stable isotopic composition of these precipitates can be used to increase our understanding of MIS 11 environments in Britain. Carbonates from Marks Tey, Clacton, Swanscombe, Elveden and West Stow are studied, the stratigraphic context of which indicates that their formation is broadly synchronous (in the mid-Hoxnian, pollen zones Ho II to Ho III). The carbon isotopic composition of groundwater and pedogenic carbonates is typically depleted with respect to δ¹³C (ca −9 to −8‰ VPDB) reflecting uptake of plant respired CO₂ during water migration/recharge. The carbon isotopic composition of lacustrine carbonate is more enriched with respect to δ¹³C (ca 0-1‰VPDB) reflecting the equilibration of lake waters with atmospheric CO₂. The δ¹⁸O of groundwater and pedogenic carbonates is slightly more enriched than modern soil carbonates but not as enriched as soil carbonates formed under interglacials that were warmer than the Holocene (i.e. the Cromerian). The stable isotopic composition of Hoxnian carbonates does not, therefore, indicate that this interglacial was characterised by uniquely warm climates in the context of other Middle Pleistocene interglacials and the Holocene. This is contrary to many marine and littoral records from around the world but consistent with environmental records from Britain and Europe.     

Tufa and travertines of southern Italy: deep‐seated,fault‐related CO2 as the key control in precipitation

Continental carbonates of Quaternary age in southern Italy commonly exhibit the facies of calcareous tufa, often reported as related to shallow aquifers fed by meteoric waters and to organic processes. A close spatial relationship exists between the mappable tufa deposits and major Quaternary extensional faults. With respect to the Ca‐Mg‐HCO₃ composition of limestone aquifers’ springs, tufa‐depositing springs exhibit higher salinity and alkalinity, are slightly warmer, have lower pH and are enriched in SO₄ and CO₂. Their δ¹³C values are systematically positive and compatible with a deep‐seated carbon source. A clear input of soil‐derived organic carbon is indicated only for small, non‐mappable tufas deposited by perched springs. The dataset indicates that the large tufa deposits owe their origin to a supplementary source of CO₂ advected by degassing through active faults, as a necessary prerequisite for inducing a rise of total dissolved salts and alkalinity. Meteoric waters that have come from a shallow aquifer are able to precipitate only limited amount of carbonates.     

Mopping up leaking carbon: A natural analog at Wadi Namaleh,Jordan

Carbon capture and sequestration (CCS) is one of the important options available for partially stemming greenhouse gas emissions from large point sources. The possibility of leaking from deep storage needs to be addressed. The Wadi Namaleh area in southern Jordan provides an interesting case study of how excess CO₂ can be trapped in the form of carbonates in the near surface, even when the local geology is not obviously conducive for such a process.Carbonate veins are formed in surface alteration zones of rhyolite host rock in this arid region. The alteration zones are limited to areas where surface soil or colluvium are present. Oxygen, deuterium and carbon isotopes of the carbonates and near-surface ground water in the area suggest that the source of carbon is deep seated CO₂, and that the carbonate precipitated in local meteoric water under ambient temperature conditions. Analysis of strontium in the carbonate, fresh rhyolite and altered host shows that the source for calcium is aeolian. Trace elements show that metal and REE mobility are constrained to the alteration zone.Thus, interaction of H₂O, CO₂ and atmospheric wet and dry deposition lead to the formation of the clayey (montmorillonite) alteration zone. This zone acts to trap seeping CO₂ and water, and thus produces conditions of progressively more efficient trapping of carbon dioxide by means of a positive feedback mechanism. Replication of these conditions in other areas will minimize CO₂ leakage from man-made CCS sites.     

The origin of mineral waters in Kiseljak near Sarajevo,Bosnia and Herzegovina

The mineral water deposits in Kiseljak are located in the central Dinarids, Bosnia and Herzegovina, in the southwestern edge of Sarajevo–Zenica basin that was formed in the zone of Busova?a fault. Busova?a fault reaches deep into the Earth’s crust and is characterised by the presence of mineral and thermomineral water enriched with CO₂ and CO₂ springs (mofetes) in the direction of Ilid?a–Kiseljak–Busova?a. Deposits are constructed of layers of Palaeozoic to Cretaceous age. Primary aquifer of mineral waters is Permian clastites and evaporites and secondary Anisian carbonates. Mineral water and CO₂ are of different origin. The water is of atmospheric origin. Due to slow circulation, water descends in the primary aquifer where it becomes enriched with CO₂ and minerals. Due to high pressure in the primary aquifer mineral water ascends along Busova?a fault, mounts into the secondary aquifer and rises at spring sources. Water is a mixture of two or more waters of different mineralization. Mixing of water occurs in the zone of secondary aquifer even at greater depths without the influence of contemporary climatic factors. Intensive water mixing is indicated by the high ratio of Ca/Sr, Na/Cl and Ca/SO₄ and the mixing diagram. CO₂ is thermometamorphic, arising from the catalytic activity of SiO₂ on carbonates in the deeper layers of the Earth’s crust, where quartz porphyry broke through Palaeozoic formations.     

Open to closed system transition traced through the TDIC isotopic signature at the aquifer recharge stage, implications for groundwater C dating

¹⁴C dating models are limited when considering recent groundwater for which the carbon isotopic signature of the total dissolved inorganic carbon (TDIC) is mainly acquired in the unsaturated zone. Reducing the uncertainties of dating thus implies a better identification of the processes controlling the carbon isotopic composition of the TDIC during groundwater recharge. Geochemical interactions between gas, water and carbonates in the unsaturated zone were investigated for two aquifers (the carbonate-free Fontainebleau sands and carbonate-bearing Astian sands, France) in order to identify the respective roles of CO₂ and carbonates on the carbon isotopic signatures of the TDIC; this analysis is usually approached using open or closed system terms. Under fully open system conditions, the seasonality of the ¹³C values in the soil CO₂ can lead to important uncertainties regarding the so-called “initial ¹⁴C activity” used in ¹⁴C correction models. In a carbonate-bearing unsaturated zone such as in the Astian aquifer, we show that an approach based on fully open or closed system conditions is not appropriate. Although the chemical saturation between water and calcite occurs rapidly within the first metre of the unsaturated zone, the carbon isotopic contents (δ¹³C) of the CO₂ and the TDIC evolve downward, impacted by the dissolution-precipitation of the carbonates. In this study, we propose a numerical approach to describe this evolution. The δ¹³C and the A¹⁴C (radiocarbon activity) of the TDIC at the base of the carbonate-bearing unsaturated zone depends on (i) the δ¹³C and the A¹⁴C of the TDIC in the soil determined by the soil CO_2, (ii) the water’s residence time in the unsaturated zone and (iii) the carbonate precipitation-dissolution fluxes. In this type of situation, the carbonate δ¹³C-A¹⁴C evolutions indicate the presence of secondary calcite and permit the calculation of its accretion flux, equal to . More generally, for other sites under temperate climate and with similar properties to the Astian sands site, this approach allows for a reliable determination of the carbon isotopic composition at the base of the unsaturated zone as the indispensable “input function” data of the carbon cycle into the aquifer.     

Reactivity of micas and cap-rock in wet supercritical CO2 with SO2 and O2 at CO2 storage conditions

Seal or cap-rock integrity is a safety issue during geological carbon dioxide capture and storage (CCS). Industrial impurities such as SO₂, O₂, and NOx, may be present in CO₂ streams from coal combustion sources. SO₂ and O₂ have been shown recently to influence rock reactivity when dissolved in formation water. Buoyant water-saturated supercritical CO₂ fluid may also come into contact with the base of cap-rock after CO₂ injection. Supercritical fluid-rock reactions have the potential to result in corrosion of reactive minerals in rock, with impurity gases additionally present there is the potential for enhanced reactivity but also favourable mineral precipitation.The first observation of mineral dissolution and precipitation on phyllosilicates and CO₂ storage cap-rock (siliciclastic reservoir) core during water-saturated supercritical CO₂ reactions with industrial impurities SO₂ and O₂ at simulated reservoir conditions is presented. Phyllosilicates (biotite, phlogopite and muscovite) were reacted in contact with a water-saturated supercritical CO₂ containing SO₂, or SO₂ and O₂, and were also immersed in the gas-saturated bulk water. Secondary precipitated sulfate minerals were formed on mineral surfaces concentrated at sheet edges. SO₂ dissolution and oxidation resulted in solution pH decreasing to 0.74 through sulfuric acid formation. Phyllosilicate dissolution released elements to solution with ∼50% Fe mobilized. Geochemical modelling was in good agreement with experimental water chemistry. New minerals nontronite (smectite), hematite, jarosite and goethite were saturated in models. A cap-rock core siltstone sample from the Surat Basin, Australia, was also reacted in water-saturated supercritical CO₂ containing SO₂ or in pure supercritical CO₂. In the presence of SO₂, siderite and ankerite were corroded, and Fe-chlorite altered by the leaching of mainly Fe and Al. Corrosion of micas in the cap-rock was however not observed as the pH was buffered by carbonate dissolution. Ca-sulfate, and Fe-bearing precipitates were observed post SO₂-CO₂ reaction, mainly centered on surface cracks and an illite rich illite-smectite precipitate quantified. Water saturated impure supercritical CO₂ was observed to have reactivity to rock-forming biotite, muscovite and phlogopite mineral separates. In the cap-rock core however carbonates and chlorite were the main reacting minerals showing the importance of assessing actual whole core.     

On the interaction of pure and impure supercritical CO2 with rock forming minerals in saline aquifers: An experimental geochemical approach

The aim of this experimental study was to evaluate and compare the geochemical impact of pure and impure CO₂ on rock forming minerals of possible CO₂ storage reservoirs. This geochemical approach takes into account the incomplete purification of industrial captured CO₂ and the related effects during injection, and provides relevant data for long-term storage simulations of this specific greenhouse gas. Batch experiments were conducted to investigate the interactions of supercritical CO₂, brine and rock-forming mineral concentrates (albite, microcline, kaolinite, biotite, muscovite, calcite, dolomite and anhydrite) using a newly developed experimental setup. After up to 42 day (1000 h) experiments using pure and impure supercritical CO₂ the dissolution and solution characteristics were examined by XRD, XRF, SEM and EDS for the solid, and ICP–MS and IC for the fluid reactants, respectively. Experiments with mixtures of supercritical CO₂ (99.5 vol.%) and SO₂ or NO₂ impurities (0.5 vol.%) suggest the formation of H₂SO₄ and HNO₃, reflected in pH values between 1 and 4 for experiments with silicates and anhydrite and between 5 and 6 for experiments with carbonates. These acids should be responsible for the general larger amount of cations dissolved from the mineral phases compared to experiments using pure CO₂. For pure CO₂ a pH of around 4 was obtained using silicates and anhydrite, and 7–8 for carbonates. Dissolution of carbonates was observed after both pure and impure CO₂ experiments. Anhydrite was corroded by approximately 50 wt.% and gypsum precipitated during experiments with supercritical CO₂ + NO₂. Silicates do not exhibit visible alterations during all experiments but released an increasing amount of cations in the reaction fluid during experiments with impure CO₂. Nonetheless, precipitated secondary carbonates could not be identified.     

CO2-water-basalt interaction. Numerical simulation of low temperature CO2 sequestration into basalts

The interaction between CO₂-rich waters and basaltic glass was studied using reaction path modeling in order to get insight into the water-rock reaction process including secondary mineral composition, water chemistry and mass transfer as a function of CO₂ concentration and reaction progress (ξ). The calculations were carried out at 25-90 °C and pCO₂ to 30 bars and the results were compared to recent experimental observations and natural systems. A thermodynamic dataset was compiled from 25 to 300 °C in order to simulate mineral saturations relevant to basalt alteration in CO₂-rich environment including revised key aqueous species for mineral dissolution reactions and apparent Gibbs energies for clay and carbonate solid solutions observed to form in nature. The dissolution of basaltic glass in CO₂-rich waters was found to be incongruent with the overall water composition and secondary mineral formation depending on reaction progress and pH. Under mildly acid conditions in CO₂ enriched waters (pH <6.5), SiO₂ and simple Al-Si minerals, Ca-Mg-Fe smectites and Ca-Mg-Fe carbonates predominated. Iron, Al and Si were immobile whereas the Mg and Ca mobility depended on the mass of carbonate formed and water pH. Upon quantitative CO₂ mineralization, the pH increased to >8 resulting in Ca-Mg-Fe smectite, zeolites and calcite formation, reducing the mobility of most dissolved elements. The dominant factor determining the reaction path of basalt alteration and the associated element mobility was the pH of the water. In turn, the pH value was determined by the concentration of CO₂ and extent of reaction. The composition of the carbonates depended on the mobility of Ca, Mg and Fe. At pH <6.5, Fe was in the ferrous oxidation state resulting in the formation of Fe-rich carbonates with the incorporation of Ca and Mg. At pH >8, the mobility of Fe and Mg was limited due to the formation of clays whereas Ca was incorporated into calcite, zeolites and clays. Competing reactions between clays (Ca-Fe smectites) and carbonates at low pH, and zeolites and clays (Mg-Fe smectites) and carbonates at high pH, controlled the availability of Ca, Mg and Fe, playing a key role for low temperature CO₂ mineralization and sequestration into basalts. Several problems of the present model point to the need of improvement in future work. The determinant factors linking time to low temperature reaction path modeling may not only be controlled by the primary dissolving phase, which presents challenges concerning non-stoichiometric dissolution, the leached layer model and reactive surface area, but may include secondary mineral precipitation kinetics as rate limiting step for specific reactions such as retrieved from the present reaction path study.     

Dissolution kinetics of selected natural minerals relevant to potential CO2-injection sites − Part 1: A review

This publication provides a literature review on experimental studies of dissolution kinetics of mainly carbonates and feldspar group minerals, i.e. most common minerals at potential CO₂-injection and/or storage sites. Geochemical interaction processes between injected CO₂ and coexisting phases, namely reservoir and cap rock minerals and formation fluids close to the CO₂-injection site can be simulated by flow-through or mixed flow reactors, while processes far from the injection site and long-term processes after termination actual CO₂-injection can be mimicked by batch reactors. At sufficient small stirring rates or fluid flow rates as well as low solute concentrations flow-through reactors are also able to simulate processes far from the injection site. The experimental parameter temperature not only intensifies the dissolution process, the dominant dissolution mechanisms are also influenced by temperature. The dissolution mechanisms change from incongruent and surface controlled mechanisms at lower temperatures to congruent and transport controlled mechanisms at higher temperatures. The CO₂ partial pressure has only a second order influence on dissolution behavior compared to the influence of pH-value and ionic strength of the CO₂-bearing brine. Minerals exposed to CO₂-bearing brines at elevated temperatures and pressures are subject of alteration, leading to severe changes of reactive surfaces and potential precipitation of secondary minerals.Computational simulations of mineral reactions at potential CO₂ storage sites have therefore to include not only the time-resolved changes of dissolution behavior and hence kinetics of mineral dissolution, but also the influence of secondary minerals on the interaction of the minerals with CO₂-enriched brines.     

Fate of carbonates within oceanic plates subducted to the lower mantle, and a possible mechanism of diamond formation

We report on high-pressure and high-temperature experiments involving carbonates and silicates at 30–80 GPa and 1,600–3,200 K, corresponding to depths within the Earth of approximately 800–2,200 km. The experiments are intended to represent the decomposition process of carbonates contained within oceanic plates subducted into the lower mantle. In basaltic composition, CaCO₃ (calcite and aragonite), the major carbonate phase in marine sediments, is altered into MgCO₃ (magnesite) via reactions with Mg-bearing silicates under conditions that are 200–300°C colder than the mantle geotherm. With increasing temperature and pressure, the magnesite decomposes into an assemblage of CO₂ + perovskite via reactions with SiO₂. Magnesite is not the only host phase for subducted carbon—solid CO₂ also carries carbon in the lower mantle. Furthermore, CO₂ itself breaks down to diamond and oxygen under geotherm conditions over 70 GPa, which might imply a possible mechanism for diamond formation in the lower mantle.     

Carbon and oxygen stable isotopes in the middle-upper miocene and lower pliocene carbonates of the Eastern Paratethys (Kerch-Taman Region): Palaeoenvironments and post-sedimentation changes

C and O isotope composition of Middle-Upper Miocene and Lower Pliocene carbonates from Kerch-Taman Region (Eastern Paratethys) have been studied in order to reconstruct palaeoenvironmental variability and post-sedimentation changes. The δ¹³C and δ¹⁸О values of the Upper Sarmatian to Lower Pliocene organogenic carbonates reflect the desalinization of paleobasins, global Late Miocene Cooling, and increase in seasonal temperature fluctuations. Isotopic composition of the Middle Sarmatian organogenic carbonates was strongly influenced by evaporation processes, high bioproductivity, and local submarine methane emissions. Warm climate and low bioproductivity together with unstable hydrological regime during the Late Chokrakian and the Karaganian times influenced the isotope composition of primary carbonates. Calcite shell of Spiratella sp. (δ¹³C =–0.4‰ and δ¹⁸О =–0.4‰) from Tarkhanian sediments was formed in warm marine environment. Dolomitization prevails over other secondary mineralization in the studied carbonate rocks. Two groups of secondary dolomites that are characterized by negative and positive δ¹³C values have been recognized. Lowe δ¹³C values (up to–31.4‰) in dolomites indicate the influence of both dissolved inorganic carbon (DIC) from oxidized organic matter (С_org) and methane. Dolomites with positive δ¹³C values (7.0 and 7.8‰) associat with migration of CO₂- and CH₄-containing saline groundwater.     

Microscale carbon isotope variability in ALH84001 carbonates and a discussion of possible formation environments

The carbonates in martian meteorite ALH84001 preserve a record of aqueous processes on Mars at 3.9 Ga, and have been suggested to contain signatures of ancient martian life. The conditions of the carbonate formation environment are critical for understanding possible evidence for life on Mars, the history of water on Mars, and the evolution of the martian atmosphere. Despite numerous studies of petrographic relationships, microscale oxygen isotope compositions, microscale chemical compositions, and other minerals associated with the carbonates, formation models remain relatively unconstrained. Microscale carbon isotope analyses of ALH84001 carbonates reveal variable δ¹³C values ranging from +27 to +64‰. The isotopic compositions are correlated with chemical composition and extent of crystallization such that the Mg-poor, early-formed carbonates are relatively ¹³C depleted and the Mg-rich, later forming carbonates, are ¹³C enriched. These data are inconsistent with many of the previously proposed environments for carbonate formation, and a new set of hypotheses are proposed. Specifically, two new models that account for the data involve low temperature (<100°C) aqueous processes: (1) the carbonates formed during mixing of two fluids derived from separate chemical and isotopic reservoirs; or (2) the carbonates formed from high pH fluids that are exposed to a CO₂-rich atmosphere and precipitate carbonate, similar to high pH springs on Earth.     

Secondary siderite-oxide-sulphide and carbonate-andalusite assemblages in cordierite granulites from Sri Lanka: post-granulite facies fluid evolution during uplift

We report here that some of the pelitic rocks from the Wanni and Highland Complexes of Sri Lanka reacted with CO₂-rich fluids to produce a wide range of unusual secondary carbonate-silicate-oxide-sulphide assemblages. These enable the depth, temperature and fluid compositions of CO₂ reactions to be calculated more rigorously than is generally possible for the patches of arrested charnockite that have been described from Sri Lanka. Magnesite-andalusite-quartz has partially replaced primary cordierite, and siderite-rutile replaced ilmenite. Paragenetic sequences involving primary pyrrhotite, ilmenite and magnetite and secondary pyrite-siderite-rutile-magnetite-(hematite) demonstrate the control which carbonate equilibria have upon evolving fluid compositions during cooling. Direct evidence for the role of graphite as a source of CO₂ is found in the Highland Complex where primary graphite partially reacted with silicates to form secondary siderite assemblages. It is proposed that following peak metamorphism, continued uplift along a clockwise P-T-t path was accompanied by a series of devolatilization reactions involving breakdown of graphite and the continuous production of secondary CO₂-rich fluids. The limited extent of disseminated secondary carbonate reflects the small amount of graphite inferred to have been present in the source rocks. These rocks demonstrate that CO₂-rich fluids, as found in disseminated fluid inclusions, need not form during peak granulite metamorphism but may be an inevitable consequence of continued uplift along a clockwise P-T-t path. The arrested charnockite which overprinted some of the hornblende-bearing felsic-intermediate composition rocks in Sri Lanka most likely formed by the same process. Received: 4 May 1994 / Accepted: 25 October 1996     

Isotope and chemical composition of gases from mud volcanoes in the Taman Peninsula and problem of their genesis

Variations in the carbon isotope composition in gases and waters of mud volcanoes in the Taman Peninsula are studied. The δ¹³C values in CH₄ and CO₂ vary from ?59.5 to ?44.0‰ (δ¹³C_av = ?52.4 ± 5.4‰) and from ?17.8 to +22.8‰ (δ¹³C_av = +6.9 ± 9.3‰), respectively. In waters from most mud volcanoes of the peninsula, this parameter ranges from +3.3 to +33.1‰, although locally lower values are also recorded (up to ?12‰. Fractionation of carbon isotopes in the CO₂-HCO₃ system corresponds to the isotope equilibrium under Earth’s surface temperatures. The growth of carbon dioxide concentration in the gaseous phase and increase in the HCO₃ ion concentration in their water phase is accompanied by the enrichment of the latter with the heavy ¹³C isotope. The δ¹³C_TDIC value in the water-soluble carbon depends on the occurrence time of water on the Earth’s surface (exchange with atmospheric CO₂, methane oxidation, precipitation of carbonates, and other processes), in addition to its primary composition. In this connection, fluctuations in δ¹³C_TDIC values in mud volcanoes with stagnant waters may amount to 10–20‰. In the clayey pulp, concentrations of carbonate matter recalculated to CaCO₃ varies from 1–4 to 36–50 wt %. The δ¹³C value in the latter ranges from ?3.6 to +8.4‰. Carbonate matter of the clayey pulp represents a mixture of sedimentogenic and authigenic carbonates. Therefore, it is usually unbalanced in terms of the carbon isotope composition with the water-soluble CO₂ forms.     

Evolution of gold-bearing veins in dykes of the Woods Point dyke swarm,Victoria

The Woods Point dyke swarm comprises hundreds of narrow, subparallel igneous dykes and dozens of pipe-shaped dyke bulges within strongly deformed early Palaeozoic turbidites of the Melbourne trough. Porpylitic alteration accompanied dyke emplacement and was followed by microfracturing induced by high fluid pressures, involving CO₂ of magmatic origin, as the dykes solidified. Further stress caused through-going faults having ladder and other patterns. Isotopic studies suggest that metamorphically or geothermally-derived solutions filled the faults and other fractures with quartz and carbonate and altered immediately adjacent dyke rock. However earlier-formed vein and wall rock carbonates retained their magmatic isotopic composition. Fluid inclusions indicate vein deposition began at approximately 400°C with salinities up to 9 weight percent NaCl. Nine sulfide minerals and gold were deposited in the veins after ankerite, sericite and albite, while quartz deposition continued through all stages. Sulfur isotopic determinations indicate the vein sulfur could not have been derived from adjacent sedimentary rocks, nor exclusively from the dykes. Metamorphic waters of marine origin is a viable source for sulfur. Saline and CO₂-rich alkaline solutions reacted with the dyke wall rocks and probably evolved chemically prior to deposition of gold. Vug carbonates deposited by meteoric water that leached vein carbonates mark the end of vein formation.Present Adress: 631 Station Street, North Carlton 3054, Victoria, Australia     

Fluid inclusion and stable isotope (O, H, C, and S) constraints on the genesis of the Serrinha gold deposit, Gurupi Belt, northern Brazil

The Serrinha gold deposit of the Gurupi Belt, northern Brazil, belongs to the class of orogenic gold deposits. The deposit is hosted in highly strained graphitic schist belonging to a Paleoproterozoic (∼2,160 Ma) metavolcano-sedimentary sequence. The ore-zones are up to 11 m thick, parallel to the regional NW–SE schistosity, and characterized by quartz-carbonate-sulfide veinlets and minor disseminations. Textural and structural data indicate that mineralization was syn- to late-tectonic and postmetamorphic. Fluid inclusion studies identified early CO₂ (CH₄-N₂) and CO₂ (CH₄-N₂)-H₂O-NaCl inclusions that show highly variable phase ratios, CO₂ homogenization, and total homogenization temperatures both to liquid and vapor, interpreted as the product of fluid immiscibility under fluctuating pressure conditions, more or less associated with postentrapment modifications. The ore-bearing fluid typically has 18–33mol% of CO₂, up to 4mol% of N₂, and less than 2mol% of CH₄ and displays moderate to high densities with salinity around 4.5wt% NaCl equiv. Mineralization occurred around 310 to 335°C and 1.3 to 3.0 kbar, based on fluid inclusion homogenization temperatures and oxygen isotope thermometry with estimated oxygen fugacity indicating relatively reduced conditions. Stable isotope data on quartz, carbonate, and fluid inclusions suggest that veins formed from fluids with δ¹⁸O_H2O and δD_H2O (310–335°C) values of +6.2 to +8.4‰ and −19 to −80‰, respectively, which might be metamorphic and/or magmatic and/or mantle-derived. The carbon isotope composition (δ¹³C) varies from −14.2 to −15.7‰ in carbonates; it is −17.6‰ in fluid inclusion CO₂ and −23.6‰ in graphite from the host rock. The δ³⁴S values of pyrite are −2.6 to −7.9‰. The strongly to moderately negative carbon isotope composition of the carbonates and inclusion fluid CO₂ reflects variable contribution of organic carbon to an originally heavier fluid (magmatic, metamorphic, or mantle-derived) at the site of deposition and sulfur isotopes indicate some oxidation of the originally reduced fluid. The deposition of gold is interpreted to have occurred mainly in response to phase separation and fluid-rock interactions such as CO₂ removal and desulfidation reactions that provoked variations in the fluid pH and redox conditions.     

The formation of graphite upon the interaction of subducted carbonates and sulfur with metal-bearing rocks of the lithospheric mantle

Experimental studies of the Fe⁰–(Mg, Ca)CO₃–S system were carried out during 18–20 h at 6.3 GPa, 900–1400°C. It is shown that the major processes resulting in the formation of free carbon include reduction of carbonates upon redox interaction with Fe⁰ (or Fe₃C), extraction of carbon from iron carbide upon interaction with a sulfur melt/fluid, and reduction of the carbonate melt by Fe–S and Fe?S–C melts. Reconstruction of the processes of graphite formation indicates that carbonates and iron carbide may be potential sources of carbon under the conditions of subduction, and participation of the sulfur melt/fluid may result in the formation of mantle sulfides.