Diffusion-controlled growth of wollastonite rims between quartz and calcite: comparison between nature and experiment

Growth rates of wollastonite reaction rims between quartz and calcite were experimentally determined at 0.1 and 1 GPa and temperatures from 850 to 1200 °C. Rim growth follows a parabolic rate law indicating that this reaction is diffusion‐controlled. From the rate constants, the D′δ‐values of the rate‐limiting species were derived, i.e. the product of grain boundary diffusion coefficient D′ and the effective grain boundary width, δ. In dry runs at 0.1 GPa, wollastonite grew exclusively on quartz surfaces. From volume considerations it is inferred that (D′_CaOδ)/(D′_SiO2δ)≥1.33, and that SiO₂ diffusion controls rim growth. D′_SiO2δ increases from about 10^?25 to 10^?23 m³ s^?1 as temperature increases from 850 to 1000 °C, yielding an apparent activation energy of 330±36 kJ mol^?1. In runs at 1 GPa, performed in a piston‐cylinder apparatus, there were always small amounts of water present. Here, wollastonite rims always overgrew calcite. Rims around calcite grains in quartz matrix are porous and their growth rates are controlled by a complex diffusion‐advection mechanism. Rim growth on matrix calcite around quartz grains is controlled by grain boundary diffusion, but it is not clear whether CaO or SiO₂ diffusion is rate‐limiting. D′δ increases from about 10^?21 to 10^?20 m³ s^?1 as temperature increases from 1100 to 1200 °C. D′_SiO2δ or D′_CaOδ in rims on calcite is c. 10 times larger than D′_SiO2δ in dry rims at the same temperature. Growth structures of the experimentally produced rims are very similar to contact‐metamorphic wollastonite rims between metachert bands and limestone in the Bufa del Diente aureole, Mexico, whereby noninfiltrated metacherts correspond to dry and brine‐infiltrated metacherts to water‐bearing experiments. However, the observed diffusivities were 4 to 5 orders of magnitude larger during contact‐metamorphism as compared to our experimental results.     

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     

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     

Formation of Wollastonite by Chemically Reactive Fluid Flow During Contact Metamorphism, Mt. Morrison Pendant, Sierra Nevada, California, USA

Quartz–calcite sandstones experienced the reaction calcite+ quartz = wollastonite + CO₂ during prograde contact metamorphismat P = 1500 bars and T = 560°C. Rocks were in equilibriumduring reaction with a CO₂–H₂O fluid with XCO₂ = 0·14.The transition from calcite-bearing, wollastonite-free to wollastonite-bearing,calcite-free rocks across the wollastonite isograd is only severalmillimeters wide. The wollastonite-forming reaction was drivenby infiltration of quartz–calcite sandstone by chemicallyreactive H₂O-rich fluids, and the distribution of wollastonitedirectly images the flow paths of reactive fluids during metamorphism.The mapped distribution of wollastonite and modeling of an O-isotopeprofile across a lithologic contact indicate that the principaldirection of flow was layer-parallel, directed upward, withany cross-layer component of flow <0·1% of the layer-parallelcomponent. Fluid flow was channeled at a scale of 1–100m by pre-metamorphic dikes, thrust and strike-slip faults, foldhinges, bedding, and stratigraphic contacts. Limits on the amountof fluid, based on minimum and maximum estimates for the displacementof the wollastonite reaction front from the fluid source, are(0·7–1·9) x 10⁵ cm³ fluid/cm² rock. Thesharpness of the wollastonite isograd, the consistency of mineralthermobarometry, the uniform measured ¹⁸O–¹⁶O fractionationsbetween quartz and calcite, and model calculations all arguefor a close approach to local mineral–fluid equilibriumduring the wollastonite-forming reaction. KEY WORDS: contact metamorphism, fluid flow, wollastonite, oxygen isotopes, reaction front     

Metamorphic reactions between fluid inclusions and mineral hosts. I. Progress of the reaction calcite+quartz=wollastonite+CO2 in natural wollastonite-hosted fluid inclusions

 At the Bufa del Diente contact-metamorphic aureole, brine infiltration through metachert layers embedded in limestones produced thick wollastonite rims, according to Cc+Qz=Wo+CO₂. Fluid inclusions trapped in recrystallized quartz hosts include: (1) high salinity four phase inclusions [Th(V-L)=460–573° C; Td(salts)=350–400° C; (Na+K)_Cleq=64–73 wt%; X _{CO 2}≤0.02]; (2) low density vapour-rich CO₂-bearing inclusions [Th(L-V)≈500±100° C; X _{CO 2}=0.22–0.44; X _NaCl≤0.01], corresponding to densities of 0.27± 0.05 gcm⁻³. Petrographical observations, phase compositions and densities show that the two fluids were simultaneously trapped in the solvus of the H₂O-CO₂-salts system at 500–600° C and 700±200 bars. The low density fluid was generated during brine infiltration at the solvus via the wollastonite producing reaction. Identical fluid types were also trapped as inclusion populations in wollastonite hosts 3 cm adjacent to quartz crystals. At room temperature, both fluid types additionally contain one quartz and one calcite crystal, generated by the back-reaction Wo+CO₂=Cc+Qz of the host with the CO₂-proportion of the fluid during retrogression. All of the CO₂ was removed from the fluid. On heating in the microstage, the reaction progress of the prograde reaction was estimated via volume loss of the calcites. In vapour-rich fluids, 50% progress is reached at 490–530° C; 80% at 520–560° C; and 100% at 540–590° C, the latter representing the trapping temperatures of the original fluid at the two fluid solvus. The progress is volume controlled. With knowledge of compositions and densities from unmodified inclusions in quartz and using the equation of state of Duan et al. (1995) for H₂O-CO₂-NaCl, along with f _{CO 2}-values extracted from it, the reaction progress curve was recalculated in the P-T-X-space. The calculated progress curve passes through the two fluid solvus up to 380° C/210 bars, continues in the one fluid field and meets the solvus again at trapping conditions. The P-T slope is steep, most of the reaction occurs above 450° C and there is high correspondence between calculated and measured reaction progress. We emphasize that with the exception of quartz, back-reactions between inclusion fluids and mineral hosts is a common process. For almost any prograde metamorphic mineral that was formed by a devolatilization reaction and that trapped the equilibrium fluid or any peak metamorphic fluid as an inclusion, a fluid-host back-reaction exists which must occur somewhere along the retrograde path. Such retrograde reactions may cause drastic changes in density and composition of the fluid. In most cases, however, evidence of the evolving mineral assemblages is not given for they might form submicroscopical layers at the inclusion walls. Received: 15 March 1995 / Accepted: 1 June 1995     

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     

Fluid flow in marbles at Jervois, central Australia: oxygen isotope disequilibrium and zoning produced by decoupling of mineralogical and isotopic resetting

The Jervois region of the Arunta Inlier, central Australia, contains para- and orthogneisses that underwent low-pressure amphibolite facies metamorphism (P = 200–300 MPa, T = 520–600 °C). Marble layers cut by metre-wide quartz + garnet ± epidote veins comprise calcite, quartz, epidote, clinopyroxene, grandite garnet, and locally wollastonite. The marbles also contain locally discordant decimetre-thick garnet and epidote skarn layers. The mineral assemblages imply that the rocks were infiltrated by water-rich fluids (XCO₂ = 0.1–0.3) at ∼600 °C. The fluids were probably derived from the quartz-garnet vein systems that represent conduits for fluids exsolved from crystallizing pegmatites emplaced close to the metamorphic peak. At one locality, the marble has calcite (Cc) δ¹⁸O values of 9–18‰ and garnet (Gnt) δ¹⁸O values of 10–14‰. The δ¹⁸O(Gnt) values are only poorly correlated with δ¹⁸O(Cc), and the δ¹⁸O values of some garnet cores are higher than the rims. The isotopic disequilibrium indicates that garnet grew before the δ¹⁸O values of the rock were reset. The marbles contain  ≤15% garnet and, for water-rich fluids, garnet-forming reactions are predicted to propagate faster than O-isotopes are reset. The Sm-Nd and Pb-Pb ages of garnets imply that fluid flow occurred at 1750–1720 Ma. There are no significant age differences between garnet cores and rims, suggesting that fluid flow was relatively rapid. Texturally late epidote has δ¹⁸O values of 1.5–6.2‰ implying δ¹⁸O(H₂O) values of 2–7‰. Waters with such low-δ¹⁸O values are probably at least partly meteoric in origin, and the epidote may be recording the late influx of meteoric water into a cooling hydrothermal system. Received: 29 April 1996 / Accepted: 12 March 1997     

Wollastonite--Scapolite Assemblages as Indicators of Granulite Pressure-Temperature-Fluid Histories: The Rauer Group, East Antarctica

Mid-Proterozoic ( 1000 Ma) granulite facies calc-silicates fromthe Rauer Group, East Antarctica, contain grossular-wollastonite-scapolite-dinopyroxene( + quartz or calcite) assemblages which preserve symplectiteand corona textures typically involving the growth of secondarywollastonite. The textures include (1) wollastonite rims betweenquartz and calcite; (2) wollastonite-plagioclase rims and intergrowthsbetween quartz and scapolite; (3) wollastonite-scapolite-clinopyroxeneinter-growths replacing grossular; and (4) wollastonite-plagioclasesymplectites replacing grossular or earlier symplectites (3). Reactions between grossular, scapolite, wollastonite, calcite,quartz, anorthite, and vapour, have been modelled in the CaO-Al₂O₃SiO₂-H₂O-CO₂and more complex systems using the internally consistent data-setof Holland & Powell (1990). Reactions producing scapoliteand wollastonite consume vapour as temperature increases (i.e., carbonation), in agreement with the results of Moecher &Essene (1990). These calc-silicates can therefore behave asfluid sinks under high-grade conditions. Conversely, they maybe important fluid sources on cooling and contribute to theformation of post-metamorphic CO₂rich fluid inclusions in isobaricallycooled granulites. P-T-_CO2 diagrams calculated for typical phase compositions (e.g., garnet, scapolite) demonstrate that the observed texturesare a record of near-isothermal decompression at 800–850 C, consistent with P—rpath determinations based on otherrock types from the Rauer Group. For example, texture (2) resultsfrom crossing the reaction Scapolite + Quartz = Wollastonite + Plagioclase + V on decompression, at 6. 5–7 kb, 820 C, and a_CO2 of0–4–0–5. Furthermore, correlations betweenmodes of product phases (e. g., wollastonitexlinopyroxene) andreactant garnet composition preclude open-system behaviour inthe formation of these textures, consistent with post-peak vapour-absentreactions such as Grossular + Calcite + Quartz = Wollastonite + Scapolite occurring on decomposition at high temperatures (>800C). Reaction textures developed in calc-silicates from other granuliteterranes often involve the formation of grossular ( + quartz calcite) as rims on wollastonite-scapolite, or replacementof wollastonite by calcite-quartz. These textures have developedprincipally in response to cooling below 780–810 C andmay be signatures of near-isobaric cooling. Infiltration ofhydrous fluid is not a necessary condition for the productionof garnet coronas in wollastonite-scapolite granulites. ^*Present address: Department of Earth Sciences, University ofMelbourne, Parkville, Victoria 3052, Australia     

Rates of grain boundary diffusion through enstatite and forsterite reaction rims

 The growth rates of enstatite rims produced by reaction of Fo₉₂ and SiO₂ were determined at 250–1500 MPa and 900–1100°C for a wide range of water contents. Growth rates were also determined for forsterite rims between MgO and Mg₂Si₂O₆ and between MgO and SiO₂. Rim growth rates are parabolic indicating diffusion-controlled growth of the polycrystalline rims which are composed of ˜ 2 μm diameter grains. Rim growth rates were used to calculate the product of the grain boundary diffusion coefficient (D'_A) times the effective grain boundary thickness (δ) assuming in turn that MgO, SiO₂, and Mg₂Si₋₁ are the diffusing components (coupled diffusion of a cation and oxygen or interdiffusion of Mg and Si). The values for D'_MgOδ, D', and D' for enstatite at 1000°C and 700 MPa confining pressure with about 0.1 wt %  water are about five times larger than the corresponding D'_Aδ values for samples initially vacuum dried at 250°C. Most of the increase in D'_Aδ occurs with the first 0.1 wt %  water. The activation energy for diffusion through the enstatite rims (1100–950°C) is 162 ± 30 kJ/mole. The diffusion rate through enstatite rims is essentially unchanged for confining pressures from 210–1400 MPa, but the nucleation rate is greatly reduced at low confining pressure (for  ≤ 1.0 wt % water present) and limits the conditions at which rim growth can be measured. The corresponding values for D'_Aδ through forsterite rims are essentially identical for the two forsterite-producing reactions when 0.1 wt % water is added and similar to the D'_Aδ values for enstatite at the same conditions. The D'_Aδ values for forsterite are ˜ 28 times larger for samples starting with 0.1 wt %  water compared to samples that were first vacuum dried. Thus water enhances these grain boundary diffusion rates by a factor of 5–30 depending on the mineralogy, but the total range in D'_Aδ is only slightly more than an order of magnitude for as wide a range of water contents as expected for most crustal conditions. Received: 1 July 1995 / Accepted: 1 August 1996     

Stable isotope and fluid inclusion evidence for the origin of the Brandberg West area Sn–W vein deposits, NW Namibia

The Brandberg West region of NW Namibia is dominated by poly-deformed turbidites and carbonate rocks of the Neoproterozoic Damara Supergoup, which have been regionally metamorphosed to greenschist facies and thermally metamorphosed up to mid-amphibolite facies by Neoproterozoic granite plutons. The meta-sedimentary rocks host Damaran-age hydrothermal quartz vein-hosted Sn–W mineralization at Brandberg West and numerous nearby smaller deposits. Fluid inclusion microthermometric studies of the vein quartz suggests that the ore-forming fluids at the Brandberg West mine were CO₂-bearing aqueous fluids represented by the NaCl–CaCl₂–H₂O–CO₂ system with moderate salinity (mean=8.6 wt% NaCl_equivalent).Temperatures determined using oxygen isotope thermometry are 415–521°C (quartz–muscovite), 392–447°C (quartz–cassiterite), and 444–490°C (quartz–hematite). At Brandberg West, the oxygen isotope ratios of quartz veins and siliciclastic host rocks in the mineralized area are lower than those in the rocks and veins of the surrounding areas suggesting that pervasive fluid–rock interaction occurred during mineralization. The O- and H-isotope data of quartz–muscovite veins and fluid inclusions indicate that the ore fluids were dominantly of magmatic origin, implying that mineralization occurred above a shallow granite pluton. Simple mass balance calculations suggest water/rock ratios of 1.88 (closed system) and 1.01 (open system). The CO₂ component of the fluid inclusions had similar δ ¹³C to the carbonate rocks intercalated with the turbidites. It is most likely that mineralization at Brandberg West was caused by a combination of an impermeable marble barrier and interaction of the fluids with the marble. The minor deposits in the area have quartz veins with higher δ ¹⁸O values, which is consistent with these deposits being similar geological environments exposed at higher erosion levels.     

Reaction kinetics of dolomite rim growth

Reaction rims of dolomite (CaMg[CO₃]₂) were produced by solid-state reactions at the contacts of oriented calcite (CaCO₃) and magnesite (MgCO₃) single crystals at 400 MPa pressure, 750–850 °C temperature, and 3–146 h annealing time to determine the reaction kinetics. The dolomite reaction rims show two different microstructural domains. Elongated palisades of dolomite grew perpendicular into the MgCO₃ interface with length ranging from about 6 to 41 µm. At the same time, a 5–71 µm wide rim of equiaxed granular dolomite grew at the contact with CaCO₃. Platinum markers showed that the original interface is located at the boundary between the granular and palisade-forming dolomite. In addition to dolomite, a 12–80 µm thick magnesio-calcite layer formed between the dolomite reaction rims and the calcite single crystals. All reaction products show at least an axiotactic crystallographic relationship with respect to calcite reactant, while full topotaxy to calcite prevails within the granular dolomite and magnesio-calcite. Dolomite grains frequently exhibit growth twins characterized by a rotation of 180° around one of the $[11\bar{2}0]$ equivalent axis. From mass balance considerations, it is inferred that the reaction rim of dolomite grew by counter diffusion of MgO and CaO. Assuming an Arrhenius-type temperature dependence, activation energies for diffusion of CaO and MgO are E _a (CaO) = 192 ± 54 kJ/mol and E _a (MgO) = 198 ± 44 kJ/mol, respectively.     

Growth kinetics of enstatite reaction rims studied on nano-scale,Part I: Methodology,microscopic observations and the role of water

The kinetics of (Mg, Fe)SiO₃ pyroxene layer growth within silicate thin films with total thickness <1 μm was studied experimentally at 0.1 MPa total pressure, controlled fO₂ and temperatures from 1,000 to 1,300°C. The starting samples were produced by pulsed laser deposition. Layer thickness before and after the experiments and layer composition as well as microstructures, grain size and shape of the interfaces were determined by Rutherford back scattering and transmission electron microscopy assisted by focused ion beam milling. Due to the miniaturization of the starting samples and the use of high resolution analytical methods the experimentally accessible temperature range for rim growth experiments was extended by about 300°C towards lower temperatures. The thickness of the layers at a given temperature increases proprotional to the square root of time, indicating a diffusion-controlled growth mechanism. The temperature dependence of rim growth yields an apparent activation energy of 426 ± 34 kJ/mol. The small grain size in the orthopyroxene rims implies a significant contribution of grain boundary diffusion to the bulk diffusion properties of the polycrystalline rims. Based on microstructural observations diffusion scenarios are discussed for which the SiO₂ component behaves immobile relative to the MgO component. Volume diffusion data for Mg in orthopyroxene from the literature indicate that the measured diffusivity is probably controlled by the mobility of oxygen. The observed reaction rates are consistent with earlier results from dry high-temperature experiments on orthopyroxene rim growth. Compared to high pressure experiments at 1,000°C and low water fugacities, reaction rates are 3–4 orders of magnitude smaller. This observation is taken as direct evidence for a strong effect of small amounts of water on diffusion in silicate polycrystals. In particular SiO₂ changes from an immobile component at dry conditions to an extremely mobile component even at very low water fugacities.     

The behavior of Mg, Fe, and Ni during the replacement of olivine by orthopyroxene: experiments relevant to mantle metasomatism

Replacement of olivine by orthopyroxene is a frequently observed phenomenon in mantle metasomatism. In order to study element redistribution in SiO₂ metasomatism we synthesised orthopyroxene reaction rims at the contacts between forsterite-rich olivine and quartz. The orthopyroxene rims grew from the original quartz-olivine interface into both directions implying counterdiffusion of iron/magnesium and silicon. Following local equilibrium partitioning the X_Fe is lower in the orthopyroxene than in the reactant olivine at the olivine-orthopyroxene replacement front. The resulting local iron excess is compensated by formation of orthopyroxene with a higher X_Fe at the quartz-orthopyroxene interface, which is out of equilibrium with the reactant olivine. This is facilitated through short circuit diffusion along grain boundaries within the orthopyroxene rim. Due to the low capacity of orthopyroxene to accommodate Ni, this component is forced to diffuse back into the olivine producing a Ni enriched zone ahead of the replacement front. This leads to Ni contents in the orthopyroxene rim, which are higher than what is expected in equilibrium with the unaltered olivine. Taking quartz as a proxy for a silica rich fluid or liquid metasomatising agent, we conclude that the overall element fractionation between olivine and the silica rich phase may deviate from equilibrium partitioning so that the Fe and Ni concentrations in the orthopyroxene which is in contact with quartz are higher than in equilibrium with the reactant olivine. This indicates that kinetic fractionation is important for the chemical evolution of both the mantle rocks and the metasomatising agents.     

Relative timing of deformation and two-stage gold mineralization at the Hutti Mine, Dharwar Craton, India

Gold mineralization at Hutti is confined to a series of nine parallel, N–S to NNW–SSE trending, steeply dipping shear zones. The host rocks are amphibolites and meta-rhyolites metamorphosed at peak conditions of 660±40°C and 4±1 kbar. They are weakly foliated (S₁) and contain barren quartz extension veins. The auriferous shear zones (reefs) are typically characterized by four alteration assemblages and laminated quartz veins, which, in places, occupy the entire reef width of 2–10 m, and contain the bulk of gold mineralization. A <1.5 m wide distal chlorite-sericite (+biotite, calcite, plagioclase) alteration zone can be distinguished from a 3–5 m wide proximal biotite-plagioclase (+quartz, muscovite, calcite) alteration zone. Gold is both spatially and temporally associated with disseminated arsenopyrite and pyrite mineralization. An inner chlorite-K-feldspar (+quartz, calcite, scheelite, tourmaline, sphene, epidote, sericite) alteration halo, which rims the laminated quartz veins, is characterized by a pyrrhotite, chalcopyrite, sphalerite, ilmenite, rutile, and gold paragenesis. The distal chlorite-sericite and proximal biotite-plagioclase alteration assemblages are developed in microlithons of the S₂–S₃ crenulation cleavage and are replaced along S₃ by the inner chlorite-K-feldspar alteration, indicating a two-stage evolution for gold mineralization. Ductile D₂ shearing, alteration, and gold mineralization formed the reefs during retrograde evolution and fluid infiltration under upper greenschist to lower amphibolite facies conditions (560±60°C, 2±1 kbar). The reefs were reactivated in the D₃ dextral strike-slip to oblique-slip environment by fault-valve behavior at lower greenschist facies conditions (ca. 300–350°C), which formed the auriferous laminated quartz veins. Later D₄ crosscutting veins and D₅ faults overprint the gold mineralization. The alteration mineralogy and the structural control of the deposit clearly points to an orogenic style of gold mineralization, which took place either during isobaric cooling or at different levels of the Archean crust. From overlaps in the tectono-metamorphic history, it is concluded that gold mineralization occurred during two tectonic events, affecting the eastern Dharwar craton in south India between ca. 2550 – 2530 Ma: (1) The assemblage of various terranes of the eastern block, and (2) a tectono-magmatic event, which caused late- to posttectonic plutonism and a thermal perturbation. It differs, however, from the pre-peak metamorphic gold mineralization at Kolar and the single-stage mineralization at Ramagiri. Notably, greenschist facies gold mineralization occurred at Hutti 35–90 million years later than in the western Dharwar craton. Editorial handling: G. Beaudoin     

Petrology of calc-silicate rocks from Koduru,Andhra Pradesh,India

Mineral assemblages, rock and mineral chemistry, and mineral reactions, in calc-silicate rocks from Koduru area, Andhra Pradesh, India are discussed. Mineralogical and bulk chemical differences indicate 3 calc-silicate rock types — type I with K feldspar+calcite+wollastonite+quartz+scapolite+diopside_ss +andradite_ss+sphene, has relatively high rock oxidation ratios. Type II is a highly calcic variety with high rock MgFe ratios, and has K feldspar+calcite+wollastonite+quartz+scapolite + diopside_ss±grossularite_ss+sphene+zoisite. Type III has K feldspar +calcite+wollastonite+quartz+scapolite+diopside_ss +sphene+hornblende+magnetite, and has relatively low oxidation ratio and low MgFe ratio. The 3 calc-silicate rock types have originated as mixtures of limestone/dolomite/marl.Diopside was produced by a reaction involving Ca-amphibole +calcite+quartz, and reversed during retrogression. Andradite_ss in type I rocks was produced at the expense of hedenbergitic component of pyroxene in a continuous reaction as a consequence of increase in the oxygen content of the original sediment relative to type III. Calcite+quartz reacted to give wollastonite. During cooling an influx of water caused scapolite to alter to zoisite.     

Hydrocarbon- and ore-bearing basinal fluids: a possible link between gold mineralization and hydrocarbon accumulation in the Youjiang basin,South China

The Youjiang basin, which flanks the southwest edge of the Yangtze craton in South China, contains many Carlin-type gold deposits and abundant paleo-oil reservoirs. The gold deposits and paleo-oil reservoirs are restricted to the same tectonic units, commonly at the basinal margins and within the intrabasinal isolated platforms and/or bioherms. The gold deposits are hosted by Permian to Triassic carbonate and siliciclastic rocks that typically contain high contents of organic carbon. Paragenetic relationships indicate that most of the deposits exhibit an early stage of barren quartz ± pyrite (stage I), a main stage of auriferous quartz + arsenian pyrite + arsenopyrite + marcasite (stage II), and a late stage of quartz + calcite + realgar ± orpiment ± native arsenic ± stibnite ± cinnabar ± dolomite (stage III). Bitumen in the gold deposits is commonly present as a migrated hydrocarbon product in mineralized host rocks, particularly close to high grade ores, but is absent in barren sedimentary rocks. Bitumen dispersed in the mineralized rocks is closely associated and/or intergrown with the main stage jasperoidal quartz, arsenian pyrite, and arsenopyrite. Bitumen occurring in hydrothermal veins and veinlets is paragenetically associated with stages II and III mineral assemblages. These observations suggest an intimate relationship between bitumen precipitation and gold mineralization. In the paleo-petroleum reservoirs that typically occur in Permian reef limestones, bitumen is most commonly observed in open spaces, either alone or associated with calcite. Where bitumen occurs with calcite, it is typically concentrated along pore/vein centers as well as along the wall of pores and fractures, indicating approximately coeval precipitation. In the gold deposits, aqueous fluid inclusions are dominant in the early stage barren quartz veins (stage I), with a homogenization temperature range typically of 230°C to 270°C and a salinity range of 2.6 to 7.2 wt% NaCl eq. Fluid inclusions in the main and late-stage quartz and calcite are dominated by aqueous inclusions as well as hydrocarbon- and CO₂-rich inclusions. The presence of abundant hydrocarbon fluid inclusions in the gold deposits provides evidence that at least during main periods of the hydrothermal activity responsible for gold mineralization, the ore fluids consisted of an aqueous solution and an immiscible hydrocarbon phase. Aqueous inclusions in the main stage quartz associated with gold mineralization (stage II) typically have a homogenization temperature range of 200–230°C and a modal salinity around 5.3 wt% NaCl eq. Homogenization temperatures and salinities of aqueous inclusions in the late-stage drusy quartz and calcite (stage III) typically range from 120°C to 160°C and from 2.0 to 5.6 wt% NaCl eq., respectively. In the paleo-oil reservoirs, aqueous fluid inclusions with an average homogenization temperature of 80°C are dominant in early diagenetic calcite. Fluid inclusions in late diagenetic pore- and fissure-filling calcite associated with bitumen are dominated by liquid C₂H₆, vapor CH₄, CH₄–H₂O, and aqueous inclusions, with a typical homogenization temperature range of 90°C to 180°C and a salinity range of 2–8 wt% NaCl eq. It is suggested that the hydrocarbons may have been trapped at relatively low temperatures, while the formation of gold deposits could have occurred under a wider and higher range of temperatures. The timing of gold mineralization in the Youjiang basin is still in dispute and a wide range of ages has been reported for individual deposits. Among the limited isotopic data, the Rb–Sr date of 206 ± 12 Ma for Au-bearing hydrothermal sericite at Jinya as well as the Re–Os date of 193 ± 13 Ma on auriferous arsenian pyrite and ⁴⁰Ar/³⁹Ar date of 194.6 ± 2 Ma on vein-filling sericite at Lannigou may provide the most reliable age constraints on gold mineralization. This age range is comparable with the estimated petroleum charging age range of 238–185 Ma and the Sm–Nd date of 182 ± 21 Ma for the pore- and fissure-filling calcite associated with bitumen at the Shitouzhai paleo-oil reservoir, corresponding to the late Indosinian to early Yanshanian orogenies in South China. The close association of Carlin-type gold deposits and paleo-oil reservoirs, the paragenetic coexistence of bitumens with ore-stage minerals, the presence of abundant hydrocarbons in the ore fluids, and the temporal coincidence of gold mineralization and hydrocarbon accumulation all support a coeval model in which the gold originated, migrated, and precipitated along with the hydrocarbons in an immiscible, gold- and hydrocarbon-bearing, basinal fluid system.     

Cracking mechanisms during galena mineralization in a sandstone-hosted lead–zinc ore deposit: case study of the Jinding giant sulfide deposit, Yunnan, SW China

There are two types of lead–zinc ore bodies, i.e., sandstone-hosted ores (SHO) and limestone-hosted ores (LHO), in the Jinding giant sulfide deposit, Yunnan, SW China. Structural analysis suggests that thrust faults and dome structures are the major structural elements controlling lead–zinc mineralization. The two types of ore bodies are preserved in two thrust sheets in a three-layered structural profile in the framework of the Jinding dome structure. The SHO forms the cap of the dome and LHO bodies are concentrated beneath the SHO cap in the central part of the dome. Quartz, feldspar and calcite, and sphalerite, pyrite, and galena are the dominant mineral components in the sandstone-hosted lead–zinc ores. Quartz and feldspar occur as detrital clasts and are cemented by diagenetic calcite and epigenetic sulfides. The sulfide paragenetic sequence during SHO mineralization is from early pyrite to galena and late sphalerite. Galena occurs mostly in two types of cracks, i.e., crescent-style grain boundary cracks along quartz–pyrite, or rarely along pyrite–pyrite boundaries, and intragranular radial cracks in early pyrite grains surrounding quartz clasts. The radial cracks are more or less perpendicular to the quartz–pyrite grain boundaries and do not show any overall (whole rock) orientation pattern. Their distribution, morphological characteristics, and geometrical relationships with quartz and pyrite grains suggest the predominant role of grain-scale cracking. Thermal expansion cracking is one of the most important mechanisms for the generation of open spaces during galena mineralization. Cracking due to heating or cooling by infiltrating fluids resulted from upwelling fluid phases through fluid passes connecting the SHO and LHO bodies, provided significant spaces for crystallization of galena. The differences in coefficients of thermal expansion between pyrite and quartz led to a difference in volume changes between quartz grains and pyrite grains surrounding them and contributed to cracking of the pyrite grains when temperature changed. Combined thermal expansion and elastic mismatch due to heating and subsequent cooling resulted in the radial and crescent cracking in the pyrite grains and along the quartz–pyrite grain boundaries.     

Humite- and scapolite-bearing assemblages in marbles and calcsilicates of Dronning Maud Land, Antarctica: new data for Gondwana reconstructions

This study investigates marbles and calcsilicates in Central Dronning Maud Land (CDML), East Antarctica. The paleogeographic positioning of CDML as part of Gondwana is still unclear; however, rock types, mineral assemblages, textures and P–T conditions observed in this study are remarkably similar to the Kerala Khondalite Belt in India. The CDML marbles and calcsilicates experienced a Pan-African granulite facies metamorphism at c. 570 Ma and an amphibolite facies retrogression at c. 520 Ma. The highest grade assemblage in marbles is forsterite+spinel+calcite+dolomite, in calcsilicates the assemblages are diopside+spinel, diopside+garnet, scapolite+wollastonite+clinopyroxene±quartz, scapolite±anorthite±calcite+clinopyroxene+wollastonite. These assemblages constrain the peak metamorphic conditions to 830±20 °C, 6.8±0.5 kbar and X _CO2>0.46. During retrogression, highly fluoric humite-group minerals (humite, clinohumite, chondrodite) replaced forsterite, and garnet rims formed at the expense of scapolite during reactions with wollastonite, calcite or clinopyroxene but without involvement of anorthite. Metamorphic conditions were about 650 °C, 4.5±0.7 kbar, 0.2< X _CO2^fluid<0.36, and the co-existence of garnet, clinopyroxene, wollastonite and quartz constrains f_O2 to FMQ-1.5 log units. Mineral textures indicate a very limited influx of H₂O-rich fluid during amphibolite facies retrogression and point to significant variations of fluid composition in mm-sized areas of the rock. Gypsum was observed in two samples; it probably replaced metamorphic anhydrite which appears to have formed under amphibolite facies conditions. The observed extensive anorogenic magmatism (anorthosites, A-type granitoids) and the character of metamorphism between 610 and 510 Ma suggest that the crustal thermal structure was characterized by a long-lived (50–100 Ma) rise of the crustal geotherm probably caused by magmatic underplating.     

Growth of multilayered polycrystalline reaction rims in the MgO–SiO<Subscript>2</Subscript> system,part I: experiments

Growth of transport-controlled reaction layers between single crystals of periclase and quartz, and forsterite and quartz was investigated experimentally at 1.5 GPa, 1100°C to 1400°C, 5 min to 72 h under dry and melt-free conditions using a piston-cylinder apparatus. Starting assemblies consisting of Per | Qtz | Fo sandwiches produced polycrystalline double layers of forsterite and enstatite between periclase and quartz, and enstatite single layers between forsterite and quartz. The position of inert Pt-markers initially deposited at the interface of the reactants and inspection of mass balance confirmed that both layer-producing reactions are controlled by MgO diffusion, while SiO₂ is relatively immobile. BSE and TEM imaging revealed thicknesses from 0.6 μm to 14 μm for double layers and from 0 to 6.8 μm for single layers. Both single and double layers displayed non-parabolic growth together with pronounced grain coarsening. Textural evolution and growth rates for each reaction are directly comparable. Forsterite–enstatite double layers are always wider than enstatite single layers, and the growth of enstatite in the double layer is slower than that in the single layer. In double layers, the enstatite/forsterite layer thickness ratio significantly increases with temperature, reflecting different MgO mobilities as temperature varies. Thus, thickness ratios in multilayered reaction zones may contain a record of temperature, but also that of any physico-chemical parameter that modifies the mobilities of the chemical components between the various layers. This potential is largely unexplored in geologically relevant systems, which calls for further experimental studies of multilayered reaction zones.     

Prograde and retrograde fluid flow during contact metamorphism of siliceous carbonate rocks from the Ballachulish aureole,Scotland

 Siliceous dolomites and limestones contain abundant retrograde minerals produced by hydration-carbonation reactions as the aureole cooled. Marbles that contained periclase at the peak of metamorphism bear secondary brucite, dolomite, and serpentine; forsterite-dolomite marbles have retrograde tremolite and serpentine; wollastonite limestones contain secondary calcite and quartz; and wollastonite-free limestones have retrograde tremolite. Secondary tremolite never appears in marbles where brucite has replaced periclase or in wollastonite-bearing limestones. A model for infiltration of siliceous carbonates by CO₂-H₂O fluid that assumes (a) vertical upwardly-directed flow, (b) fluid flux proportional to cooling rate, and (c) flow and reaction under conditions of local equilibrium between peak temperatures and ≈400 °C, reproduces the modes of altered carbonate rocks, observed reaction textures, and the incompatibility between tremolite and brucite and between tremolite and wollastonite. Except for samples from a dolomite xenolith, retrograde time-integrated flux recorded by reaction progress is on the order of 1000 mol fluid/cm² rock. Local focusing of flow near the contact is indicated by samples from the xenolith that record values an order of magnitude greater. Formation of periclase, forsterite, and wollastonite at the peak of metamorphism also required infiltration with prograde time-integrated flux approximately 100–1000 mol/cm². The comparatively small values of prograde and retrograde time-integrated flux are consistent with lack of stable isotope alteration of the carbonates and with the success of conductive thermal models in reproducing peak metamorphic temperatures recorded by mineral equilibria. Although isobaric univariant assemblages are ubiquitous in the carbonates, most formed during retrograde metamorphism. Isobaric univariant assemblages observed in metacarbonates from contact aureoles may not record physical conditions at the peak of metamorphism as is commonly assumed. Received: 19 September 1995 / Accepted: 14 March 1996