The role of the grain boundary at chemical and isotopic fronts in marble during contact metamorphism

Carbon and oxygen isotopic profiles around a low pressure metasomatic wollastonite reaction front in a marble of the Hida metamorphic terrain, central Japan, display typical metamorphic fluid-enhanced isotopic zonations. Isotopic profiles obtained from detailed microscale analyses perpendicular to the chemical reaction front in calcite marble show that diffusion-enhanced isotopic exchange may control these profiles. Carbon and oxygen isotopic behaviour in grain boundaries is remarkably different. Oxygen isotopic troughs (¹⁸O depleted rims) around the calcite-grain boundaries are widely observed in this contact aureole, demonstrating that diffusion of oxygen in calcite grain boundary dominates over lattice diffusion in calcite. In contrast, no difference is observed in carbon isotopic profiles obtained from grain cores and rims. There is thus no specific role of the grain boundary for diffusion of carbonic species in the metamorphic fluid during transportation. Carbon chemical species such as CO₂ and CO₃ ions in metamorphic fluid migrate mainly through lattice diffusion. The carbon and oxygen isotope profiles may be modelled by diffusion into a semi-infinite medium. Empirically lattice diffusion of oxygen isotopes is almost six times faster than that of carbon isotopes, and oxygen grain-boundary diffusion is ten times faster than oxygen lattice diffusion. Oxygen isotopic results around the wollastonite vein indicate that migration of the metamorphic fluid into calcite marble was small and was parallel to the aquifer. From the stability of wollastonite and the attainment of oxygen isotopic equilibrium, we suggest that diffusion of oxygen occurred through an aqueous fluid phase. The timescale of formation of the oxygen isotopic profile around the wollastonite vein is calculated to be about 0.76 × 10⁶ years using the experimentally determined diffusion constant. Received: 14 January 1997 / Accepted: 23 April 1998     

Oxygen isotope exchange and disequilibrium between calcite and tremolite in the absence and presence of an experimental C–O–H fluid

Oxygen isotope exchange between minerals during metamorphism can occur in either the presence or the absence of aqueous fluids. Oxygen isotope partitioning among minerals and fluid is governed by both chemical and isotopic equilibria during these processes, which progress by intragranular and intergranular diffusion as well as by surface reactions. We have carried out isotope exchange experiments in two- and three-phase systems, respectively, between calcite and tremolite at high temperatures and pressures. The two-phase system experiments were conducted without fluid either at 1 GPa and 680 °C for 7 days or at 500 MPa and 560 °C for 20 days. Extrapolated equilibrium fractionations between calcite and tremolite are significantly lower than existing empirical estimates and experimental determinations in the presence of small amounts of fluid, but closely match calculated fractionations by means of the increment method for framework oxygen in tremolite. The small fractionations measured in the direct calcite–tremolite exchange experiments are interpreted by different rates of oxygen isotope exchange between hydroxyl oxygen, framework oxygen and calcite during the solid–solid reactions where significant recrystallization occurs. The three-phase system experiments were accomplished in the presence of a large amount of fluid (CO₂+H₂O) at 500 MPa and 560 °C under conditions of phase equilibrium for 5, 10, 20, 40, 80, 120, 160, and 200 days. The results show that oxygen isotope exchange between minerals and fluid proceeds in two stages: first, through a mechanism of dissolution-recrystallization and very rapidly; second, through a mechanism of diffusion and very slowly. Synthetic calcite shows a greater rate of isotopic exchange with fluid than natural calcite in the first stage. The rate of oxygen diffusion in calcite is approximately equal to or slightly greater than that in tremolite in the second stage. A calculation using available diffusion coefficients for calcite suggests that grain boundary diffusion, rather than volume diffusion, has been the dominant mechanism of oxygen transport between the fluid and the mineral grains in the later stage.Editorial responsibility: T.L. Grove     

Fast diffusion along mobile grain boundaries in calcite

Experimental measurements of grain boundary diffusion are usually conducted on static boundaries, despite the fact that grain boundaries deep in the Earth are frequently mobile. In order to explore the possible effect of boundary mobility on grain boundary diffusion rates we have measured the uptake of ⁴⁴Ca from a layer of ⁴⁴Ca-enriched calcite powder during the static recrystallization of a single crystal of calcite at 900°C. A region about 500 μm wide adjacent to the powder layer is heterogeneously enriched in ⁴⁴Ca, and complex zoning patterns, including sharp steps in composition and continuous increases and decreases in ⁴⁴Ca content, are developed. In metamorphic rocks, these would normally be interpreted in terms of changes in pressure or temperature, Rayleigh fractionation, or episodic fluid infiltration. These explanations cannot apply to our experiments, and instead the zoning patterns are interpreted as being due to variations in grain boundary migration rate. We have applied an analytical model which allows the product of grain boundary diffusion coefficient and grain boundary width (D _GB δ) to be calculated from the grain boundary migration rate and the compositional gradient away from the powder layer. The value of D _GB δ in the mobile grain boundaries is at least five orders of magnitude greater than the published value for static boundaries under the same conditions. In order to allow the scale of chemical equilibrium (and hence textural evolution) to be predicted under both experimental and geological conditions, we present quantitative diffusion-regime maps for static and mobile boundaries in calcite, using both published values and our new values for grain boundary diffusion in mobile boundaries. Enhanced diffusion in mobile boundaries has wide implications for the high temperature rheology of Earth materials, for geochronology, and for interpretations of the length- and time-scales of chemical mass-transport. Moreover, zones of anomalously high electrical conductivity in the crust and mantle could be regions undergoing recrystallization such as active shear zones, rather than regions of anomalous mineralogy, water- or melt-content as is generally suggested.     

Calcite solubility and speciation in supercritical NaCl-HCl aqueous fluids

The solubility of calcite in NaCl-H₂O and in HCl-H₂O fluids was measured using an extraction-quench hydrothermal apparatus. Experiments were conducted at 2 kbar, between 400° C and 600° C. Measurements in NaCl-H₂O were conducted in two ways: 1) at constant pressure and NaCl concentration, as a function of temperature; and 2) at constant pressure and temperature, as a function of NaCl concentration. In both the NaCl-H₂O and the HCl-H₂O systems, the solubility of calcite increases with increasing chlorine concentrations. For example, the log calcium molality in equilibrium with calcite increases from –3.75 at 2 kbar and 500° C, in pure H₂O to –3.10 at 2 kbar and 500° C at log NaCl molality=–1.67. At fixed pressure and NaCl molality, the solubility of calcite is almost constant from 400° C to 550° C, but increases somewhat at higher temperatures. The results can be used to determine the dominant calcium species in the experimental solutions as a function of NaCl concentration and to obtain values for the second dissociation constant of CaCl_2(aq). At 2 kbar, 400° C, 500° C, and 600° C, we calculate values for the log of the dissociation constant of CaCl⁺ of –2.1, –3.2, and –4.3, respectively. The 400° C and 500° C values are consistent with those obtained by Frantz and Marshall (1982) using electrical conductance techniques. However, our 600° C value is 0.8 log units higher than that reported by Frantz and Marshall. The calcite solubilities in the NaCl-H₂O and HCl-H₂O systems are inconsistent with the solubilities of calcite in pure H₂O reported by Walther and Long (1986). They are, however, consistent with the measurements of calcite solubilities in pure H₂O presented in this study. These results allow for the calculation of the solubilities of calcium silicates and carbonates in fluids that contain CO₂ and NaCl.     

Compositional zoning of calcite in a high-pressure metamorphic calc-schist: clues to heterogeneous grain-scale fluid distribution during exhumation

Calcite in former aragonite–dolomite-bearing calc-schists from the ultrahigh-pressure metamorphic (UHPM) oceanic complex at Lago di Cignana, Valtournanche, Italy, preserved different kinds of zoning patterns at calcite grain and phase boundaries. These patterns are interpreted in terms of lattice diffusion and interfacial mass transport linked with a heterogeneous distribution of fluid and its response to a changing state of stress. The succession of events that occurred during exhumation is as follows: As the rocks entered the calcite stability field at T=530–550 °C, P ca. 1.2 GPa, aragonite occurring in the matrix and as inclusions in poikilitic garnet was completely transformed to calcite. Combined evidence from microstructures and digital element distribution maps (Mn-, Mg-, Fe- and Ca–Kα radiation intensity patterns) indicates that transformation rates have been much higher than rates of compositional equilibration of calcite (involving resorption of dolomite and grain boundary transport of Mg, Fe and Ca). This rendered the phase transformation an isochemical process. During subsequent cooling to T ca. 490 °C (where lattice diffusion effectively closed), grains of matrix calcite have developed diffusion-zoned rims, a few hundred micrometres thick, with Mg and Fe increasing and Ca decreasing towards the phase boundary. Composition profiles across concentrically zoned, large grains in geometrically simple surroundings can be successfully modelled with an error function describing diffusion into a semi-infinite medium from a source of constant composition. The diffusion rims in matrix calcite are continuous with quartz, phengite, paragonite and dolomite in the matrix. This points to an effective mass transport on phase boundaries over a distance of several hundred micrometres, if matrix dolomite has supplied the Mg and Fe needed for incorporation in calcite. In contrast, diffusion rims are lacking at calcite–calcite and most calcite–garnet boundaries, implying that only very minor mass transport has occurred on these interfaces over the same T–t interval. From available grain boundary diffusion data and experimentally determined fluid–solid grain boundary structures, inferred large differences in transport rates can be best explained by the discontinuous distribution of aqueous fluid along grain/phase boundaries. Observed patterns of diffusion zoning indicate that fluid was distributed not only along grain-edge channels, but spread out along most calcite–white mica and calcite–quartz two-grain junctions. On the other hand, the inferred non-wetting of calcite grain boundaries in carbonate-rich domains is compatible with fluid–calcite–calcite dihedral angles >60° determined by Holness and Graham (1995) for a wide range of fluid compositions under the P–T conditions of interest. Whereas differential stress has been very low at the stage of diffusion zoning (T > 490 °C), it increased as the rocks were cooling below 440 °C (at 0.3–0.5 GPa). Dislocation creep and the concomitant increase of strain energy in matrix calcite induced migration recrystallisation of high-angle grain boundaries. For that stage, the compositional microstructure of recrystallised calcite grain boundary domains indicates significant mass transport along calcite two-grain junctions, which at the established low temperatures is likely to have been accomplished by ionic diffusion within a hydrous grain boundary fluid film (“dynamic wetting” of migrating grain boundaries). Received: 10 January 2000 / Accepted: 10 April 2000     

Volume and grain boundary diffusion of calcium in natural and hot-pressed calcite aggregates

 Calcium self-diffusion rates in natural calcite single crystals were experimentally determined at 700 to 900° C and 0.1 MPa in a stream of CO₂. Diffusion coefficients (D) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The Arrhenius parameters yield an activation energy (Q)=382±37 kJ/mol and pre-exponential factor (D₀)=0.13 m²/s, and there is no measurable anisotropy. Calcium grain boundary diffusion rates were experimentally determined in natural (Solnhofen) limestone and hot-pressed calcite aggregates at 650° to 850° C and 0.1 to 100 MPa pressure. The Solnhofen limestone was first pre-annealed for 24 h at 700° C and 100 MPa confining pressure under anhydrous conditions to produce an equilibrium microstructure for the diffusion experiments. Values for the product of the grain boundary diffusion coefficient (D′) and the effective grain boundary diffusion width (δ) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The results show that there is no measurable difference between D′δ values obtained for pre-annealed Solnhofen samples at 0.1 and 100 MPa or between hot-pressed calcite aggregates and pre-annealed Solnhofen samples. The temperature dependence for calcium grain boundary diffusion in Solnhofen samples annealed at 0.1 MPa is described by the Arrhenius parameters D^′ ₀δ=1.5×10⁻⁹ m³/s and Q=267±47 kJ/mol. Comparison of the results of this study with previously published data show that calcium is the slowest volume diffusing species in calcite. The calcium diffusivities measured in this study place constraints on several geological processes that involve diffusive mass transfer including diffusion-accommodated mechanisms in the deformation of calcite rocks. Received: 19 December 1994/Accepted: 30 June 1995     

Reaction-induced fluid flow in synthetic quartz-bearing marbles

The reaction kinetics and fluid expulsion during the decarbonation reaction of calcite+quartz=wollastonite+CO₂ in water-absent conditions were experimentally investigated using a Paterson-type gas apparatus. Starting materials consisted of synthetic calcite/quartz rock powders with variable fractions of quartz (10, 20, and 30 wt%) and grain sizes of 10 µm (calcite) and 10 and 30 µm (quartz). Prior to reaction, samples were HIPed at 700 °C and 300 MPa confining pressure and varying pore pressures. Initial porosity was low at 2.7–6.3%, depending on pore pressure during HIP and the amount and grain size of quartz particles. Samples were annealed at reaction temperatures of 900 and 950 °C at 150 and 300 MPa confining pressures, well within the wollastonite stability field. Run durations were between 10 min and 20 h. SEM micrographs of quenched samples show growth of wollastonite rims on quartz grains and CO₂-filled pores between rims and calcite grains and along calcite grain boundaries. Measured widths of wollastonite rims vs. time indicate a parabolic growth law. The reaction is diffusion-controlled and reaction progress and CO₂ production are continuous. Porosity increases rapidly at initial stages of the reaction and attains about 10–12% after a few hours. Permeability at high reaction temperatures is below the detection limit of 10^–21 m² and not affected by increased porosity. This makes persistent pore connectivity improbable, in agreement with observed fluid inclusion trails in form of unconnected pores in SEM micrographs. Release of CO₂ from the sample was measured in a downstream reservoir. The most striking observation is that fluid release is not continuous but occurs episodic and in pulses. Ongoing continuous reaction produces increase in pore pressure, which is, once having attained a critical value (P_crit), spontaneously released. Connectivity of the pore space is short-lived and transient. The resulting cycle includes pore pressure build-up, formation of a local crack network, pore pressure release and crack closure. Using existing models for plastic stretching and decrepitation of pores along with critical stress intensity factors for the calcite matrix and measured pore widths, it results that P_crit is about 20 MPa. Patterns of fluid flow based on mineralogical and stable isotope evidence are commonly predicted using the simplifying assumption of a continuous and constant porosity and permeability during decarbonation of the rock. However, simple flow models, which assume constant pore pressure, constant fluid filled porosity, and constant permeability may not commonly apply. Properties are often transient and it is most likely that fluid flow in a specific reacting rock volume is a short-lived episodic process.Editorial responsibility: J. Hoefs     

P-T-X effects on equilibrium carbonate-H2O-CO2-NaCl dihedral angles: constraints on carbonate permeability and the role of deformation during fluid infiltration

Fluid-solid-solid dihedral angles in the NaCl-H₂O-CO₂-calcite-dolomite-magnesite system have been determined at pressures ranging from 0.5 to 7 kbar and temperatures from 450°C to 750°C. At 1 kbar and 650°C, both dolomite and magnesite exhibit a dihedral angle minimum for intermediate H₂O-CO₂ fluids similar to that previously determined by the present authors for calcite, but the depth of the minimum is smaller, being above the critical value of 60° for both dolomite and magnesite for all fluid compositions. Calcite-calcite-brine dihedral angles at 650°C have been determined in the pressure range 1–5 kbar. Angles decrease with increasing salt content of the fluid, tending towards a constant value of about 65° for strong brines at pressures above 2 kbar. There is a general increase of angle with increasing pressure which is most marked for strong brines. A positive correlation of angle with pressure is also observed in calcite-H₂O-CO₂ fluids, the position of the minimum moving towards higher angles and towards H₂O-rich fluids with increasing pressure. The permeability window previously observed by the present authors at 1 kbar and intermediate fluid compositions closes at about 1.5 kbar. The results demonstrate that the permeability of carbonates to grain edge fluid flow is only possible at low pressures and for fluids of restricted H₂O-CO₂-NaCl compositions. However, geochemical evidence from metamorphic terrains suggests that pervasive infiltration does occur under conditions where impermeability is predicted. From examination of published studies of infiltrated carbonates we conclude that deformation plays a critical role in enhancing carbonate permeability. Possible mechanisms for this include shear-enhanced dilatancy (micro-cracking), fluid inclusion drag by deformation-controlled grain boundary migration, and dynamically maintained transient grain boundary fluid films.     

The effect of bulk composition on the solidus of carbonated eclogite from partial melting experiments at 3 GPa

To explore the effect of bulk composition on the solidus of carbonated eclogite, we determined near-solidus phase relations at 3 GPa for four different nominally anhydrous, carbonated eclogites. Starting materials (SLEC1, SLEC2, SLEC3, and SLEC4) were prepared by adding variable proportions and compositions of carbonate to a natural eclogite xenolith (66039B) from Salt Lake crater, Hawaii. Near-solidus partial melts for all bulk compositions are Fe–Na calcio-dolomitic and coexist with garnet + clinopyroxene + ilmenite ± calcio-dolomitic solid solution. The solidus for SLEC1 (Ca#=100 × molar Ca/(Ca + Mg + Fe_T)=32, 1.63 wt% Na₂O, and 5 wt% CO₂) is bracketed between 1,050°C and 1,075°C (Dasgupta et al. in Earth Planet Sci Lett 227:73–85, 2004), whereas initial melting for SLEC3 (Ca# 41, 1.4 wt% Na₂O, and 4.4 wt% CO₂) is between 1,175°C and 1,200°C. The solidus for SLEC2 (Ca# 33, 1.75 wt% Na₂O, and 15 wt% CO₂) is estimated to be near 1,100°C and the solidus for SLEC3 (Ca# 37, 1.47 wt% Na₂O, and 2.2 wt% CO₂) is between 1,100°C and 1,125°C. Solidus temperatures increase with increasing Ca# of the bulk, owing to the strong influence of the calcite–magnesite binary solidus-minimum on the solidus of carbonate bearing eclogite. Bulk compositions that produce near-solidus crystalline carbonate closer in composition to the minimum along the CaCO₃-MgCO₃ join have lower solidus temperatures. Variations in total CO₂ have significant effect on the solidus if CO₂ is added as CaCO₃, but not if CO₂ is added as a complex mixture that maintains the cationic ratios of the bulk-rock. Thus, as partial melting experiments necessarily have more CO₂ than that likely to be found in natural carbonated eclogites, care must be taken to assure that the compositional shifts associated with excess CO₂ do not unduly influence melting behavior. Near-solidus dolomite and calcite solid solutions have higher Ca/(Ca + Mg) than bulk eclogite compositions, owing to Ca–Mg exchange equilibrium between carbonates and silicates. Carbonates in natural mantle eclogite, which have low bulk CO₂ concentration, will have Ca/Mg buffered by reactions with silicates. Consequently, experiments with high bulk CO₂ may not mimic natural carbonated eclogite phase equilibria unless care is taken to ensure that CO₂ enrichment does not result in inappropriate equilibrium carbonate compositions. Compositions of eclogite-derived carbonate melt span the range of natural carbonatites from oceanic and continental settings. Ca#s of carbonatitic partial melts of eclogite vary significantly and overlap those of partial melts of carbonated lherzolite, however, for a constant Ca-content, Mg# of carbonatites derived from eclogitic sources are likely to be lower than the Mg# of those generated from peridotite.     

Grain growth kinetics of dolomite,magnesite and calcite: a comparative study

The rates of grain growth of stoichiometric dolomite [CaMg(CO₃)₂] and magnesite (MgCO₃) have been measured at temperatures T of 700–800°C at a confining pressure P _c of 300 MPa, and compared with growth rates of calcite (CaCO₃). Dry, fine-grained aggregates of the three carbonates were synthesized from high purity powders by hot isostatic pressing (HIP); initial mean grain sizes of HIP-synthesized carbonates were 1.4, 1.1, and 17 μm, respectively, for CaMg(CO₃)₂, MgCO₃, and CaCO₃, with porosities of 2, 28, and 0.04% by volume. Grain sizes of all carbonates coarsened during subsequent isostatic annealing, with mean values reaching 3.9, 5.1, and 27 μm for CaMg(CO₃)₂, MgCO₃, and CaCO₃, respectively, in 1 week. Grain growth of dolomite is much slower than the growth rates of magnesite or calcite; assuming normal grain growth and n = 3 for all three carbonates, the rate constant K for dolomite (≃5 × 10⁻⁵ μm³/s) at T = 800°C is less than that for magnesite by a factor of ~30 and less than that for calcite by three orders of magnitude. Variations in carbonate grain growth may be affected by differences in cation composition and densities of pores at grain boundaries that decrease grain boundary mobility. However, rates of coarsening correlate best with the extent of solid solution; K is the largest for calcite with extensive Mg substitution for Ca, while K is the smallest for dolomite with negligible solid solution. Secondary phases may nucleate at advancing dolomite grain boundaries, with implications for deformation processes, rheology, and reaction kinetics of carbonates.     

Contact metamorphism and metasomatism at a dolerite-limestone contact in the Gebel Yelleq area, Northern Sinai, Egypt

Summary At the northeastern flank of Gebel Yelleq, northern Sinai, pure limestones of Upper Cretaceous age were subjected to a thermal overprint, caused by a c. 80m thick Tertiary olivine dolerite sill. Metasomatic supply of Si, Al, Fe, Mg and Ti was greater to the c. 7m wide upper than to the c. 25m wide lower thermal aureole. The greater width of the lower aureole is possibly due to a longer duration of the thermal overprint at this contact. Mineral assemblages in both aureoles are (from the contact outward):(i) clinopyroxene + garnet ± wollastonite + calcite(ii) garnet ± wollastonite + calcite;(iii) wollastonite + calcite.In places, late stage xenoblasts of apophyllite and witherite overgrow these assemblages. Garnets are grandites to melanites with Grs_56–86Adr_14–42Sch_0–2Sps_0–0.2Prp₀ in the lower, and Grs_29–94Adr_5–64Sch_0–12Sps_0–0.2Prp_0–1.7 in the upper aureole. Close to the upper contact, clinopyroxene is virtually pure diopside with X _Mg = Mg/(Mg + Fe²⁺) = 0.97–1.0, whereas clinopyroxenes farther away from the upper contact and in the lower aureole have X _Mg-values of 0.49 and 0.53, respectively.The minimum temperatures reached during contact metamorphism in the upper and lower aureole are defined by the lower stability limit of wollastonite. The temperatures are inferred with a calculated T-X(CO₂) projection in the system CMASCH and are estimated at c. 290 °C and 380 °C for X(CO₂) values of 0.05 and 0.25, respectively. A pressure of roughly 100 bar is estimated for the lower dolerite-limestone contact. As indicated by one-dimensional thermal modelling, a maximum temperature of 695 °C was attained at this contact, assuming a magma temperature of 1150 °C. Further modelling results indicate (i) wollastonite, which occurs first 13 m away from the lower contact, formed at a maximum temperature of c. 575 °C, (ii) there, wollastonite formation lasted for approximately 170 years and, (iii) at the outer rim of the lower aureole, the maximum temperature reached was 480 °C, and temperatures sufficient for wollastonite formation lasted for about 140 years.     

Oxygen isotope fractionation between calcite and tremolite: an experimental study

Oxygen isotope partitioning between calcite and tremolite was experimentally calibrated in the presence of small amounts of a supercritical CO₂–H₂O fluid at temperatures from 520 to 680° C and pressures from 3 to 10 kbar. The experiments were carried out within the stability field of the calcite-tremolite assemblage based on phase equilibrium relationships in the system CaO–MgO–SiO₂–CO₂–H₂O, so that decomposition of calcite and tremolite was avoided under the experimental conditions. Appropriate proportions of carbon dioxide to water were used to meet this requirement. Large weight ratios of mineral to fluid were employed in order to make the isotopic exchange between calcite and tremolite in the presence of a fluid close to that without fluid. The data processing method for isotopic exchange in a three-phase system has been applied to extrapolate partial equilibrium data to equilibrium values. The determined fractionation factors between calcite (Cc) and tremolite (Tr) are expressed as:10³1n _Cc-Tr=3.80 × 10⁶/T ²-1.67By combining the present data with the experimental calibrations of Clayton et al. (1989) on the calcite-quartz system, we obtain the fractionation for the quartztremolite system: 10³1n _Qz-Tr=4.18 × 10⁶/T ²-1.67Our experimental calibrations are in good agreement with the theoretical calculations of Hoffbauer et al. (1994) and the empirical estimates of Bottinga and Javoy (1975) based on isotopic data from naturall assemblages. At 700 C good agreement also exists between our experimental data and theoretical values calculated by Zheng (1993b). With decreasing temperature, however, an increasing difference between these data appears.Retrograde isotopic reequilibration by oxygen diffusion may be common for amphibole relative to diopside in metamorphic rocks. However, isotopic equilibrium in amphibole can be preserved in cases of rapid cooling.     

Calculated greenschist facies mineral equilibria in the system CaO−FeO−MgO−Al2O3−SiO2−CO2−H2O

The kinetic problems associated with the experimental determination of reactions among complex solidsolution phases at low temperatures have hindered our understanding of the phase relations in greenschist facies rocks. In the absence of reliable experimental data, we have used the new, expanded internally-consistent thermodynamic dataset of Holland and Powell (1990), to present calculated phase equilibria for the system CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ (CaFMASCH) with quartz in eccess, in the range 400°–500°C at low to intermediate pressures, involving the minerals amphibole, chlorite, anorthite, clinozoisite, dolomite, chloritoid, garnet, margarite, andalusite, and calcite. By solving independent sets of non-linear equations formed from equilibrium relationships, we calculate not only the loci of reactions in pressuretemperature-x(CO₂) space, but also the compositions of coexisting minerals in terms of the substitutions, FeMg_-1 and (Fe,Mg)SiAl_-1Al_-1. Invariant, univariant and divariant equilibria are calculated and discussed in relation to naturally-occurring greenschist facies metabasic and siliceous dolomitic mineral assemblages. We thus avoid the use of activity-corrected curves so commonly presented in the literature as a substitute for genuine univariant phase diagram boundaries.     

Regional fluid and metal mobility in the Dalradian metamorphic belt,Southern Grampian Highlands,Scotland

A prominent set of veins was formed during post-metamorphic deformation of the Caledonian Dalradian metamorphic belt. These veins are concentrated in dilational zones in fold hinges, but apophyses follow schistosity and fold axial surface fractures. The veins are most common in the cores of regional structures, especially the Dalradian Downbend and consist of quartz, calcite, chlorite and metallic sulphides and oxides. Metals, including gold, have been concentrated in the veins. The fluid which formed the veins was low salinity (1–5 wt% NaCl and KCl) CO₂-bearing (3–16 wt% CO₂) water of metamorphic origin. The fluid varies slightly in composition within and between samples, but is essentially uniform in composition over several hundred km². Vein formation occurred at about 350±50 °C and 200–300 MPa pressure. Further quartz mineralization occurred in some dilational zones at lower temperatures (160–180 °C). This later mineralization was accompanied by CO₂ immiscibility. Dilution and oxidation of the metamorphic fluid occurred due to mixing with meteoric water as the rocks passed through the brittle-ductile transition. A similar metamorphic fluid is thought to have been responsible for gold mineralization in the nearby Tyndrum Fault at a later stage in the Dalradian uplift.     

Fluid inclusion study of Xihuashan tungsten deposit in the southern Jiangxi province,China

The Xihuashan tungsten deposit is closely related to a small highly evolved granitic intrusion. The fluid phases associated with the wolframite-bearing quartz veins have been investigated using microthermometry and the Raman microprobe; they are highly variable in density and composition. The earlier fluids are low-density and low-salinity CO₂-bearing aqueous solutions circulating at temperatures up to 420 °C, and low-salinity (2–3 equiv. wt% NaCl) aqueous solutions without traces of CO₂ circulating at high temperatures 280°–400 °C) involved in a specific hydrothermal fracturing event; limited unmixing occurs at 380 °C and 200–100 bar in response to a sudden pressure drop. The second types of fluids related to deposition of idiomorphic drusy quartz are typical CO₂-bearing aqueous solutions with low salinity (2.5 equiv. wt% NaCl) homogenizing at low to moderate temperatures (180°–340 °C). The late fluids characterize the sulfide deposition stage; they are aqueous fluids with variable salinities homogenizing in the liquid phase between 100° and 275 °C. The Xihuashan hydrothermal evolution resulted from a discontinuous sequence of specific events occurring between 420° and 150 °C and during a continuous hydrothermal evolution of the system during cooling. The role played by the CO₂-rich fluids in the transport and deposition of tungsten in the hydrothermal environment is discussed.     

Interdiffusion of H2O and CO2 in metamorphic fluids at ∼490 to 690°C and 1 GPa

Interdiffusion coefficients have been determined for H₂O-CO₂ mixtures by quantifying the flux of CO₂ between two fluid-filled chambers in a specially designed piston-cylinder cell. The two chambers, which are maintained at 1.0 GPa and at temperatures differing by ∼100°C, each contain the X_CO2-buffering assemblage calcite + quartz + wollastonite, in H₂O. The positive dependence of X_CO2 on temperature results in a down-temperature, steady-state flux of CO₂ through a capillary tube that connects the two chambers. This flux drives the wollastonite = calcite + quartz equilibrium to the right in the cooler chamber, producing a measurable amount of calcite that is directly related to CO₂-H₂O interdiffusion rates. Diffusivities calculated from seven experiments range from 1.0 × 10⁻⁸ to 6.1 × 10⁻⁸ m²/s for mean capillary temperatures between ∼490 and 690°C. The data set can be approximated by an Arrhenius-type relation:

     

A petrographic and fluid inclusion comparison of silver-poor and silver-rich zinc-lead ores,N.E. Washington State

Many of the zinc-lead deposits of NE Washington State are poorly known examples of Mississippi Valley Type (MVT) mineralization. This study compares inclusion fluids from the Josephine Breccia ores with the later cross-cutting sulfide-bearing quartz veins. The breccia ores are cemented mainly by open space fillings of dolomite, sphalerite, quartz, galena, jasperoid and calcite. Replacement is of minor importance. Ore and gangue deposition occurred over the range 150–250 °C with most of the temperatures less than 200 °C. The aqueous brines typically contain 17–23 equivalent weight percent NaCl with often substantial amounts of Ca and/or Mg chlorides. Homogenization temperatures do not delineate any cooling or paragenetic sequence. The cross-cutting vein quartz contains CO₂-rich inclusions with overall densities usually less than 0.7 g/cc and homogenization temperatures from 250–325 °C. Sulfur isotope analyses yield two populations with the quartz vein ores being lighter (<13 permil CDT) than the average for the conformable ores. The later veins are not remobilized MVT sulfides but represent a separate, high-silver period of mineralization.     

Hydrothermal alteration of Tertiary igneous rocks from the Isle of Skye,northwest Scotland

Hydrothermal alteration of Tertiary gabbros from Skye involved the reaction of igneous olivine, augite, hypersthene, plagioclase, magnetite, and ilmenite with aqueous fluid primarily to combinations of talc, chlorite, montmorillonite, calcic amphibole, biotite, and secondary magnetite. Lesser amounts of calcite, epidote, quartz, sphene, prehnite, and garnet also developed. During mineralogical alteration of gabbro there was a net addition to rock of K, Na, Sr, and H₂O and a net loss of Mg. Gabbro was oxidized early in the hydrothermal event and later reduced. Iron and silicon were probably initially lost and later added. There is no evidence for significant change in the Al or Ca content of the gabbros. Hydrothermal alteration of Skye gabbro involved not only large-scale migration of ¹⁸O, ¹⁶O, D and H but also of K, Na, Sr, Mg, and probably Fe and Si.Mineral thermometry indicates that pyroxenes in the gabbros crystallized at 1000° C–1150° C and were very resistent chemically as well as isotopically to later hydrothermal alteration. Hypothetical equilibrium between primary and secondary mafic silicates suggests that mineralogical alteration of gabbro occurred at 450°–550° C. The lack of correlation between mineralogical and isotopic alteration of gabbro requires that much isotopic alteration occurred at temepratures above those at which the secondary minerals developed, 550°–1000° C. The chemical alteration of gabbro is correlated with its mineralogical alteration and therefore occurred at 450°–550° C.Measured progress of the mineral-fluid reactions was used to estimate the amount of H₂O fluid that infiltrated the gabbro as primary olivine was converted to talc+magnetite at 525°–550° C. Calculated fluid-rock ratios are in the range 0.2–6 (volume basis) and are smaller than values estimated from isotopic data (fluid/rock 1–10, volume basis). Both isotopic and petrologic data point to pervasive flow of fluid through crystalline rock at elevated temperatures of 500°–1000° C. Isotopic fluid-rock ratios are larger than petrologic fluid-rock ratios because isotopic alteration of cooling gabbro began earlier and at higher temperatures than did the mineralogical alteration.     

Syn-metamorphic gold mineralisation,Invincible Vein,NW Otago Schist,New Zealand

The Invincible Vein fills a fault zone which strikes northeast and dips steeply southeast in the lower Rees Valley, NW Otago. The vein cuts north striking foliation in lower greenschist facies Otago Schist. Structures associated with the fault zone are both brittle and ductile, and the fault zone has had a complex history of post-mineralisation reactivation. Mineralised vein material filling parts of the fault zone consist of quartz, albite, muscovite, chlorite, calcite, pyrite, arsenopyrite and minor gold. These minerals have been strained and locally recrystallised during ductile deformation. Fluid inclusion homogenisation temperatures (140–175°C) and ice melting temperatures (0 to –1°C) indicate that the mineralising fluid was low salinity, low CO₂ water with a density between 0.88 and 0.93 g/cm₃. Arsenopyrite geothermometry implies a temperature of mineralisation of 370 ± 70°C. Mineralisation pressure lay between 2 and 5 kbar. Mineralisation pressure-temperature conditions and mineralogy are essentially the same as for metamorphism of the host schist. Vein calcite oxygen isotope ratios (+12 to +15 per mil) are similar to host schist values. Carbon isotope ratios of vein calcite (– 3 to –5 per mil) are distinctly different from ratios in host schist (–7 to –10 per mil). Elevated vein Cr contents, and isotopically depleted carbon data, are consistent with some degree of equilibration with metavolcanic rocks. It is inferred that metavolcanic rocks of the underlying Aspiring Terrane were a significant source for mineralising fluid and metals. Invincible mineralisation occurred in the latter stages of metamorphism, and is the earliest recognised gold-bearing vein system in the Otago Schist.     

Note on the genesis of brucite in contact metamorphism of dolomite

Consideration of available thermodynamic data and the published results of direct experiments relating to (1) formation. of periclase from dolomite and (2) hydration of periclase to brucite, permits the following conclusions to be drawn: (1) At very low partial pressures of CO₂ (perhaps of the order of 1 bar) and relatively high partial pressures of water (up to 2000 bars), dolomite can break down directly to brucite and calcite at temperatures above about 400° C, and below temperatures on the brucite dehydration curve. (2) The reaction dolomite calcite + periclase + CO₂ in contact metamorphism near granitic bodies is likely to occur only at low partial pressures of CO₂ (perhaps 10 or 20 bars); this can be achieved without direct formation of brucite, by maintaining a partial pressure of water of the order of 1000 bars or more. (3) At low CO₂ pressures dolomite may re-form in the cooling stages of metamorphism by reaction between calcite, brucite, and CO₂ at temperatures below about 400° C.