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     

Diffusion-controlled growth of albite and pyroxene reaction rims

The growth rates of albite and pyroxene (enstatite + diopside + spinel) reaction rims were measured at 1000°C and ˜700 MPa and found to be parabolic indicating diffusion-controlled growth. The parabolic rate constants for the pyroxene (+ spinel) rims in samples with 0.5 wt% H₂O added or initially vacuum dried at 25°C and 250°C are 1.68 ± 0.09, 0.54 ± 0.05 and 0.25 ± 0.06 μm²/h, respectively. The values for albite rim growth in samples initially dried at 60°C and with 0.1 wt% H₂O added are 0.25 ± 0.04 and 0.33 ± 0.03 μm²/h, respectively. The latter values were used to derive the product of the grain boundary diffusion coefficient D′_A, where A = SiO₂, NaAlO₂, or NaAlSi₋₁, and the grain boundary thickness δ in albite. The calculated D′_SIO2δ in the albite aggregate for the situations of two different water contents are about 9.9 × 10⁻²³ and 1.4 × 10⁻²² m³ s⁻¹, respectively. Both the rate constants and the calculated D′_Aδ demonstrate that the effect of water content on the grain boundary diffusion rate in monomineralic albite and polymineralic pyroxene (+ spinel) aggregates is small, consistent with recent studies of monomineralic enstatite and forsterite rims. Received: 1 July 1995 / Accepted: 1 August 1996     

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     

Mineralogy and geochemistry of El Dorado epithermal gold deposit, El Sauce district, central-northern Chile

Summary The El Dorado Au-Cu deposit is located in an extensive intra-caldera zone of hydrothermal alteration affecting Upper Cretaceous andesites of the Los Elquinos Formation at La Serena (≈ 29°47′S Lat., 70°43′W Long., Chile). Quartz-sulfide veins of economic potential are hosted by N25W and N20E fault structures associated with quartz-illite alteration (+supergene kaolinite). The main ore minerals in the deposit are pyrite, chalcopyrite ± fahlore (As/(As + Sb): 0.06−0.98), with electrum, sphalerite, galena, bournonite-seligmanite (As/(As + Sb): 0.21−0.31), marcasite, pyrrhotite being accessory phases. Electrum, with an Ag content between 32 and 37 at.%, occurs interstitial to pyrite aggregates or along pyrite fractures. Pyrite commonly exhibits chemical zonation with some zones up to 1.96 at.% As. Electron probe microanalyses of pyrite indicate that As-rich zones do not exhibit detectable Au values. Fluid inclusion microthermometry shows homogenization temperatures between 130 and 352 °C and salinities between 1.6 and 6.9 wt.% NaCl eq. Isotope data for quartz, ankerite and phyllosilicates and estimated temperatures show that δ¹⁸O and δD for the hydrothermal fluids were between 3 and 10‰ and between −95 and −75‰, respectively. These results suggest the mineralizing fluids were a mixture of meteoric and magmatic waters. An epithermal intermediate-sulfidation model is proposed for the formation of the El Dorado deposit. Author’s present address: J. Carrillo-Rosúa, Dpto. de Didáctica de las Ciencias Experimentales, Universidad de Granada, Campus de Cartuja, 18071, Granada, Spain     

The normal grain growth behaviour of nominally pure calcitic aggregates

The normal grain growth behaviour of four different, but all nominally pure, calcite powders (99%+ analytic grade calcite, 99.7% chalk, 99.97% crushed Iceland Spar, 99.95%+ chelometric grade calcite) has been investigated as a function of temperature (550, 600, 650, 700 °C) and confining pressure (100, 190 MPa) under both “dry” and hydrostatic (P _fluid = P _total) conditions. The initial particle size of both the analytic grade and chelometric grade calcite was about 5 μm, and that of the chalk was about 3 μm, while the experiments on the Iceland Spar were conducted on powders of three different initial particle sizes (3.4, 7.5, 38.5 μm). On each material, at each pressure/temperature condition 6 to 15 experiments, equally spaced in log time from 15 minutes to 50 days, were conducted. Under dry conditions all four materials recrystallized to aggregates which contained less than 2% porosity and which had a grain size of between 4 and 20 μm (depending on the initial particle size). Subsequently the aggregates coarsened by normal grain growth, with the kinetics of the growth process being controlled by the rate at which the grain boundaries could drag the residual pores with them as they migrated. Under nominally identical conditions both the mechanism and rates of pore drag differed greatly for the different materials, implying that this process is highly sensitive to trace solute impurity concentrations. This sensitivity renders the task of providing a systematic account of dry calcite grain growth kinetics highly problematic. Under hydrostatic conditions all the powders followed the same normal grain growth kinetics in which the growth process was rate-controlled by diffusion through the pore fluid on the grain boundaries. An activation enthalpy of 162.6 kJ mol⁻¹ and an activation volume of 34.35 cm³ mol⁻¹ was obtained for this process. Received: 23 May 1996 / Accepted: 8 July 1997     

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     

Iron Diffusion in Single-Crystal Diopside

 Iron tracer diffusion experiments in diopside have been performed using natural and synthetic single crystals of diopside, and stable iron tracers enriched in ⁵⁴Fe, at temperatures in the range 950–1100 °C, total pressure 1 atm, for times up to 29 days. Iron isotope diffusion profiles were determined with an ion microprobe. For experiments performed at log pO₂ = −13, in directions parallel to the c axis and the b axis of two natural, low iron (Fe ∼ 1.8 at %) diopsides, the data obey a single Arrhenius relationship of the form D = 6.22_−5.9 ^+49.6×10⁻¹⁵ exp(−161.5 ± 35.0 kJ mol⁻¹/RT) m² s⁻¹. A single datum for iron diffusion in iron-free, single-crystal diopside at 1050 °C, is approximately 1 order of magnitude slower than in the natural crystals. The pO₂ dependence of iron diffusion in natural crystals at 1050 °C (power exponent = 0.229 ± 0.036) indicates a vacancy mechanism; this is consistent with the results of unpublished atomistic simulation studies. There is no evidence of anisotropy for iron diffusion in diopside. Received: 16 March 1999 / Accepted: 10 April 2000     

Contrasting fluids and reservoirs in the contiguous Marcona and Mina Justa iron oxide–Cu (–Ag–Au) deposits,south-central Perú

The Marcona–Mina Justa deposit cluster, hosted by Lower Paleozoic metaclastic rocks and Middle Jurassic shallow marine andesites, incorporates the most important known magnetite mineralization in the Andes at Marcona (1.9 Gt at 55.4% Fe and 0.12% Cu) and one of the few major iron oxide–copper–gold (IOCG) deposits with economic Cu grades (346.6 Mt at 0.71% Cu, 3.8 g/t Ag and 0.03 g/t Au) at Mina Justa. The Middle Jurassic Marcona deposit is centred in Ica Department, Perú, and the Lower Cretaceous Mina Justa Cu (Ag, Au) prospect is located 3–4 km to the northeast. New fluid inclusion studies, including laser ablation time-of-flight inductively coupled plasma mass spectrometry (LA-TOF-ICPMS) analysis, integrated with sulphur, oxygen, hydrogen and carbon isotope analyses of minerals with well-defined paragenetic relationships, clarify the nature and origin of the hydrothermal fluid responsible for these contiguous but genetically contrasted deposits. At Marcona, early, sulphide-free stage M-III magnetite–biotite–calcic amphibole assemblages are inferred to have crystallized from a 700–800°C Fe oxide melt with a δ¹⁸O value from +5.2‰ to +7.7‰. Stage M-IV magnetite–phlogopite–calcic amphibole–sulphide assemblages were subsequently precipitated from 430–600°C aqueous fluids with dominantly magmatic isotopic compositions (δ³⁴S = +0.8‰ to +5.9‰; δ¹⁸O = +9.6‰ to +12.2‰; δD = −73‰ to −43‰; and δ¹³C = −3.3‰). Stages M-III and M-IV account for over 95% of the magnetite mineralization at Marcona. Subsequent non-economic, lower temperature sulphide–calcite–amphibole assemblages (stage M-V) were deposited from fluids with similar δ³⁴S (+1.8‰ to +5.0‰), δ¹⁸O (+10.1‰ to +12.5‰) and δ¹³C (−3.4‰), but higher δD values (average −8‰). Several groups of lower (<200°C, with a mode at 120°C) and higher temperature (>200°C) fluids can be recognized in the main polymetallic (Cu, Zn, Pb) sulphide stage M-V and may record the involvement of modified seawater. At Mina Justa, early magnetite–pyrite assemblages precipitated from a magmatic fluid (δ³⁴S = +0.8‰ to +3.9‰; δ¹⁸O = +9.5‰ to +11.5‰) at 540–600°C, whereas ensuing chalcopyrite–bornite–digenite–chalcocite–hematite–calcite mineralization was the product of non-magmatic, probably evaporite-sourced, brines with δ³⁴S ≥ +29‰, δ¹⁸O = 0.1‰ and δ¹³C = −8.3‰. Two groups of fluids were involved in the Cu mineralization stage: (1) Ca-rich, low-temperature (approx. 140°C) and high-salinity, plausibly a basinal brine and (2) Na (–K)-dominant with a low-temperature (approx. 140°C) and low-salinity probably meteoric water. LA-TOF-ICPMS analyses show that fluids at the magnetite–pyrite stage were Cu-barren, but that those associated with external fluids in later stages were enriched in Cu and Zn, suggesting such fluids could have been critical for the economic Cu mineralization in Andean IOCG deposits.     

Timescales and mechanisms of fluid infiltration in a marble: an ion microprobe study

Using a recently developed ion microprobe technique, a detailed oxygen isotope map of calcite grains in a coarse-grained marble has been constructed, supported by trace element (Mn, Sr, Fe) analysis and cathodoluminescence (CL) imaging, in order to constrain scales of oxygen isotope equilibrium, timescales and mechanisms of metamorphic fluid infiltration, and fluid sources and pathways. Results are compared with a previous study of this sample (Wada 1988) carried out using a cryo-microtome technique and conventional oxygen isotope analysis. The marble, from the high temperature/low pressure Hida metamorphic belt in north-central Japan, underwent granulite facies followed by amphibolite facies metamorphic events, the latter associated with regional granite intrusion. The CL imaging indicates two types of calcite, a yellow luminescing (YLC) and a purple luminescing (PLC) variety. The YLC, which occupies grain boundaries, fractures, replacement patches, and most of the abundant deformation twin lamellae, post-dates the dominant PLC calcite and maps out fluid pathways. Systematic relationships were established between oxygen isotope and trace element composition, calcite type and texture, based on 74 ¹⁸O/¹⁶O and 17 trace element analyses with 20–30 μ m spatial resolution. The YLC is enriched in Mn and Fe, and depleted in ¹⁸O and Sr compared to PLC, and is much more ¹⁸O depleted than is indicated from conventional analyses. Results are interpreted to indicate infiltration of ¹⁸O-depleted (metamorphic or magmatic) fluid (initial δ¹⁸O = 9‰–10.5‰) along grain boundaries, fractures and deformation twin lamellae, depleting calcite grains in Sr and enriching them in Mn and Fe. The sample is characterised by gross isotopic and elemental disequilibrium, with important implications for the application of chromatographic theory to constrain fluid fluxes in metacarbonate rocks. Areas of PLC unaffected by “short-circuiting” fluid pathways contain oxygen diffusion profiles of ∼10‰/∼200 μm in grain boundary regions or adjacent to fractures/patches. When correction is made for estimated grain boundary/fracture and profile orientation in 3D, profiles are indistinguishable within error. Modelling of these profiles gives consistent estimates of Dt (where D is the diffusion coefficient and t is time) of ∼0.8 × 10⁻⁸ m², from which, using experimental data for oxygen diffusion in calcite, timescales of fluid transport along grain boundaries at amphibolite facies temperatures of ∼10³ to ∼10⁴ years are obtained. These short timescales, which are much shorter than plausible durations of metamorphism, imply that rock permeabilities may be transiently much higher during fluid flow than those calculated from time integrated fluid fluxes or predicted from laboratory measurements. The preservation of ¹⁸O/¹⁶O profiles requires either rapid cooling rates (∼100–600 °C/million years), or, more plausibly, loss of grain boundary fluid such that a dry cooling history followed the transient passage of fluid. The δ¹⁸O/trace element correlations are also consistent with volume diffusion-controlled transport in the PLC. Fluid transport and element exchange occurred by two inter-related mechanisms on short timescales and on different lengthscales – long-distance flow along cracks, grain boundaries and twin lamellae coupled to ∼200 μm-scale volume diffusion of oxygen. Received: 8 December 1997 / Accepted: 18 May 1998     

The low-pressure stability of phase A,Mg7Si2O8(OH)6

 Phase A, Mg₇Si₂O₈(OH)₆, is a dense hydrous magnesium silicate whose importance as a host of H₂O in the Earth’s mantle is a subject of debate. We have investigated the low-pressure stability of phase A in experiments on the reaction phase A=brucite+forsterite. Experiments were conducted in piston-cylinder and multi-anvil apparatus, using mixtures of synthetic phase A, brucite and forsterite. The reaction was bracketed between 2.60 and 2.75 GPa at 500° C, between 3.25 and 3.48 GPa at 600° C and between 3.75 and 3.95 GPa at 650^° C. These pressures are much lower than observed in the synthesis experiments of Yamamoto and Akimoto (1977). At 750^° C the stability field of brucite + chondrodite was entered. The enthalpy of formation and entropy of phase A at 1 bar (10⁵ Pa), 298 K, were derived from the experimental brackets on the reaction phase A=brucite+forsterite using a modified version of the thermodynamic dataset THERMOCALC of Holland and Powell (1990), which includes a new equation of state of H₂O derived from the molecular dynamics simulations of Brodholt and Wood (1993). The data for phase A are: ΔH ^o _f =−7126±8 kJ mol^-1, S ^o=351 J K^-1 mol^-1. Incorporating these data into THERMOCALC allows the positions of other reactions involving phase A to be calculated, for example the reaction phase A + enstatite=forsterite+vapour, which limits the stability of phase A in equilibrium with enstatite. The calculated position of this reaction (753° C at 7 GPa to 937° C at 10 GPa) is in excellent agreement with the experimental brackets of Luth (1995) between 7 and 10 GPa, supporting the choice of equation of state of H₂O used in THERMOCALC. Comparison of our results with calculated P-T paths of subducting slabs (Peacock et al. 1994) suggests that, in the system MgO–SiO₂–H₂O, phase A could crystallise in compositions with Mg/Si>2 at pressures as low as 3 GPa. In less Mg rich compositions phase A could crystallise at pressures above approximately 6 GPa. Received: 3 July 1995/Accepted: 14 December 1995     

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     

Volcanogenic-related origin of sulfide-rich quartz veins: evidence from O and S isotopes at the Géant Dormant gold mine, Abitibi belt, Canada

The Géant Dormant gold mine is a sulfide-rich quartz vein gold deposit hosted by a volcano-sedimentary sequence and an associated felsic endogenous dome and dikes. The auriferous quartz-sulfide veins were preceded by two synvolcanic gold-bearing mineralizing events: early sulfidic seafloor-related and later disseminated pyrite in the felsic dome. This deposit differs from classical Archean auriferous quartz vein deposits by the low carbonate and high sulfide contents of the veins and by their formation prior to ductile penetrative deformation. The δ¹⁸O values of quartz associated with seafloor-related auriferous sulfides average 11.9 ± 0.6‰ (n = 3). The seafloor hydrothermal fluids had a δ¹⁸O value of 3.2‰ calculated at 250 °C. The oxygen isotope composition of quartz and chlorite from veins average 12.5 ± 0.3‰ (n = 20) and 5.9 ± 1.1‰ (n = 4) respectively. Assuming oxygen isotope equilibrium between quartz and chlorite, the veins formed at a temperature of ∼275 °C, which is consistent with the calculated temperature of 269 ± 10 °C from chlorite chemistry. The gold-bearing fluids had a δ¹⁸O value of 4.7‰ calculated at 275 °C. The δ³⁴S values of sulfides from the three gold events range from 0.6 to 2.8‰ (n = 32) and are close to magmatic values. Sulfur isotope geothermometry constrains the sulfide precipitation in the gold-bearing veins at a temperature of ∼350 °C. The similarity of the isotope data, the calculated δ¹⁸O of the mineralizing fluids and the likely seawater fluid source suggest that the three mineralizing events are genetically related to a volcanogenic hydrothermal system. The high value of the auriferous fluids (δ¹⁸O = 4.7‰) is attributed to a significant magmatic fluid contribution to the evolved seawater-dominated convective hydrothermal system. The two-stage filling of veins at increasing temperature from quartz-chlorite (275 °C) to sulfides (350 °C) may reflect the progressive maturation of volcanogenic hydrothermal systems. These results, together with field and geochemical data, suggest that formation of gold-rich volcanogenic systems require specific conditions that comprise a magmatic fluid contribution and gold from arc-related felsic rocks, coeval with the mineralizing events. This study shows that some auriferous quartz-vein orebodies in Archean terranes are formed in volcanogenic rather than mesothermal systems. Received: 12 December 1998 / Accepted: 5 July 1999     

Hydrogen diffusion in spinel grain boundaries and consequences for chemical homogenization in hydrous peridotite

Hydrogen can be stored in the structure of nominally anhydrous minerals as point defects, and these impurities substantially modify many physical properties of Earth’s mantle minerals. However, mantle rocks are composed of mineral grains separated by grain boundaries and interphase grains boundaries. Therefore, as a potential hydrogen reservoir, grain boundaries should be given proper attention. Here, I report an experimental investigation into hydrogen diffusion through grain boundaries in polycrystalline aggregates. Sintering and diffusion experiments were performed using a gas-medium high-pressure vessel at under pressure of 300 MPa and over a temperature range of 900–1,250°C. The diffusion assembly consisted of a polycrystalline cylinder of aluminous spinel + olivine crystals with a talc cylinder as the main hydrogen source. A Ni capsule was used to buffer the oxygen fugacity at Ni–NiO. Experimental durations varied from 5 min to 5 h. The presence of hydrogen in the crystals was measured by Fourier-transform infrared spectroscopy. The calculation of the diffusion coefficients was based on the estimation of the characteristic distance. The absence or presence of hydrogen recorded by the ‘hydrogen sensor’ olivines embedded in the aggregate allows the estimation of bounds on this characteristic distance. Results presented here suggest that hydrogen effective diffusion coefficients are only one order of magnitude faster (~10⁻⁹ m²s⁻¹ at 1,000°C) than in an olivine single crystal along the [100] axis. Resulting diffusion coefficients for hydrogen in grain boundary are four orders of magnitude faster than in a single crystal, but this diffusivity is not fast enough to affect hydrogen mobility in mantle rocks with grain sizes greater than ~1 mm. Thus, very limited chemical homogenization would occur using grain boundaries diffusion in mantle hydrous peridotite for incompatible and volatile element, such as hydrogen.     

The partitioning of copper between silicate melts and two-phase aqueous fluids: An experimental investigation at 1 kbar, 800° C and 0.5 kbar, 850° C

 Experiments were performed in the three phase system high-silica rhyolite melt+low-salinity aqueous vapor+hydrosaline brine, to investigate the partitioning equilibria for copper in magmatic-hydrothermal systems at 800° C and 1 kbar, and 850° C and 0.5 kbar. D^aqm/mlt _Cu and apparent equilibrium constants, K^aqm/mlt _Cu,Na, between the aqueous mixture (aqm=quenched vapor+brine) and the silicate melt (mlt) are calculated. D^aqm/mlt _Cu increases with increasing aqueous chloride concentration and is a function of pressure. K^aqm/mlt _Cu,Na=215(±73) at 1 kbar and 800° C and K^aqm/mlt _Cu,Na=11(±6) at 0.5 kbar and 850°C. Decreasing pressure from 1 to 0.5 kbar lowers K^aqm/mlt _Cu,Na by a factor of approximately 20. Data revealed no difference in K^aqm/mlt _Cu,Na or D^aqm/mlt _Cu as a function of the melt aluminium saturation index. Within the 2-phase field the K^aqm/mlt _Cu,Na show no variation with total aqueous chloride, indicating that copper-sodium exchange between the vapor, brine and silicate melt is independent of the mass proportion of vapor and brine. Model copper-sodium apparent equilibrium constants for the hydrosaline brine and the silicate melt revealed a negative dependence on pressure. Model apparent equilibrium constants for copper-sodium exchange between the brine and vapor were close to unity at 1 kbar and 800° C. Received: 27 June 1994/Accepted: 30 March 1995     

Oxygen diffusion in anhydrous sanidine feldspar

Oxygen exchange experiments have been performed between single crystals of sanidine feldspar and oxygen gas enriched in ¹⁸O, at temperatures in the range 869–1053 °C, total pressure 1 atmosphere, for times up to 28 days. Oxygen isotope diffusion profiles in a direction perpendicular to (001) were determined with an ion microprobe. The experimental data obey a single Arrhenius relationship of the form D = 8.4 × 10⁻¹¹ exp. (−245 ± 15 kJ mol⁻¹/RT) m²s⁻¹. The results indicate that oxygen diffusion in anhydrous sanidine feldspar is marginally slower than oxygen diffusion in anhydrous anorthite. Comparison with published atomistic simulation studies suggests that oxygen transport in feldspar is by an interstitial mechanism. Received: 17 October 1997 / Accepted: 6 July 1998     

The effect of dissolved magnesium on creep of calcite II: transition from diffusion creep to dislocation creep

We extended a previous study on the influence of Mg solute impurity on diffusion creep in calcite to include deformation under a broader range of stress conditions and over a wider range of Mg contents. Synthetic marbles were produced by hot isostatic pressing (HIP) mixtures of calcite and dolomite powders for different intervals (2–30 h) at 850°C and 300 MPa confining pressure. The HIP treatment resulted in high-magnesian calcite aggregates with Mg content ranging from 0.5 to 17 mol%. Both back-scattered electron images and chemical analysis suggested that the dolomite phase was completely dissolved, and that Mg distribution was homogeneous throughout the samples at the scale of about two micrometers. The grain size after HIP varied from 8 to 31 μm, increased with time at temperature, and decreased with increasing Mg content (>3.0 mol%). Grain size and time were consistent with a normal grain growth equation, with exponents from 2.4 to 4.7, for samples containing 0.5–17.0 mol% Mg, respectively. We deformed samples after HIP at the same confining pressure with differential stresses between 20 and 200 MPa using either constant strain rate or stepping intervals of loading at constant stresses in a Paterson gas-medium deformation apparatus. The deformation tests took place at between 700 and 800°C and at strain rates between 10⁻⁶ and 10⁻³ s⁻¹. After deformation to strains of about 25%, a bimodal distribution of large protoblasts and small recrystallized neoblasts coexisted in some samples loaded at higher stresses. The deformation data indicated a transition in mechanism from diffusion creep to dislocation creep. At stresses below 40 MPa, the strength was directly proportional to grain size and decreased with increasing Mg content due to the reductions in grain size. At about 40 MPa, the sensitivity of log strain rate to log stress, (n), became greater than 1 and eventually exceeded 3 for stresses above 80 MPa. At a given strain rate and temperature, the stress at which that transition occurred was larger for samples with higher Mg content and smaller grain size. At given strain rates, constant temperature, and fixed grain size, the strength of calcite in the dislocation creep regime increased with solute content, while the strength in the diffusion creep regime was independent of Mg content. The results suggest that chemical composition will be an important element to consider when solid substitution can occur during natural deformation.     

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.     

Effect of impurities on grain growth in synthetic calcite aggregates

 We investigated grain growth of calcite aggregates fabricated from crushed natural single crystals with different impurity content. The total trace-element concentration of the starting powders varied from about 170 ppm to more than 930 ppm with Mn as the major component. Samples were produced by hot-isostatic pressing of the different powders at 300 MPa confining pressure at 600 °C for 2 h. The starting material for the anneals was dry and had a uniform microstructure with an average grain size of about 5 μm and a porosity of <2.1%. Three disks with Mn concentrations of 10, 350, and 670 ppm, respectively, were annealed in the same run at a confining pressure of 300 MPa, and temperatures between 700 and 900 °C for up to 20 h. Grain growth was fastest in samples with the highest Mn concentrations. A multivariable fit to the data yields grain-growth exponents of 2.0 ± 0.3 for samples with 10 ppm Mn and 2.3 ± 0.2 for those with 670 ppm Mn. The activation energies for grain growth vary from 99 ± 12 kJ mol⁻¹ to 147 ± 14 kJ mol⁻¹ for the respective calcite compositions. Received: 22 August 2000 / Accepted: 12 March 2001     

A fluid inclusion and stable isotope study of 200 Ma of fluid evolution in the Galway Granite, Connemara, Ireland

Fluid inclusions in granite quartz and three generations of veins indicate that three fluids have affected the Caledonian Galway Granite. These fluids were examined by petrography, microthermometry, chlorite thermometry, fluid chemistry and stable isotope studies. The earliest fluid was a H₂O-CO₂-NaCl fluid of moderate salinity (4–10 wt% NaCl eq.) that deposited late-magmatic molybdenite mineralised quartz veins (V₁) and formed the earliest secondary inclusions in granite quartz. This fluid is more abundant in the west of the batholith, corresponding to a decrease in emplacement depth. Within veins, and to the east, this fluid was trapped homogeneously, but in granite quartz in the west it unmixed at 305–390 °C and 0.7–1.8 kbar. Homogeneous quartz δ¹⁸O across the batholith (9.5 ± 0.4‰n = 12) suggests V₁ precipitation at high temperatures (perhaps 600 °C) and pressures (1–3 kbar) from magmatic fluids. Microthermometric data for V₁ indicate lower temperatures, suggesting inclusion volumes re-equilibrated during cooling. The second fluid was a H₂O-NaCl-KCl, low-moderate salinity (0–10 wt% NaCl eq.), moderate temperature (270–340 °C), high δD (−18 ± 2‰), low δ¹⁸O (0.5–2.0‰) fluid of meteoric origin. This fluid penetrated the batholith via quartz veins (V₂) which infill faults active during post-consolidation uplift of the batholith. It forms the most common inclusion type in granite quartz throughout the batholith and is responsible for widespread retrograde alteration involving chloritization of biotite and hornblende, sericitization and saussuritization of plagioclase, and reddening of K-feldspar. The salinity was generated by fluid-rock interactions within the granite. Within granite quartz this fluid was trapped at 0.5–2.3 kbar, having become overpressured. This fluid probably infiltrated the Granite in a meteoric-convection system during cooling after intrusion, but a later age cannot be ruled out. The final fluid to enter the Granite and its host rocks was a H₂O-NaCl-CaCl₂-KCl fluid with variable salinity (8–28 wt% NaCl eq.), temperature (125–205 °C), δD (−17 to −45‰), δ¹⁸O (−3 to + 1.2‰), δ¹³C_CO2 (−19 to 0‰) and δ³⁴S_sulphate (13–23‰) that deposited veins containing quartz, fluorite, calcite, barite, galena, chalcopyrite sphalerite and pyrite (V₃). Correlations of salinity, temperature, δD and δ¹⁸O are interpreted as the result of mixing of two fluid end-members, one a high-δD (−17 to −8‰), moderate-δ¹⁸O (1.2–2.5‰), high-δ¹³C_CO2 (> −4‰), low-δ³⁴S_sulphate (13‰), high-temperature (205–230 °C), moderate-salinity (8–12 wt% NaCl eq.) fluid, the other a low-δD (−61 to −45‰), low-δ¹⁸O (−5.4 to −3‰), low-δ¹³C (<−10‰), high-δ³⁴S_sulphate (20–23‰) low-temperature (80–125 °C), high-salinity (21–28 wt% NaCl eq.) fluid. Geochronological evidence suggests V₃ veins are late Triassic; the high-δD end-member is interpreted as a contemporaneous surface fluid, probably mixed meteoric water and evaporated seawater and/or dissolved evaporites, whereas the low-δD end-member is interpreted as a basinal brine derived from the adjacent Carboniferous sequence. This study demonstrates that the Galway Granite was a locus for repeated fluid events for a variety of reasons; from expulsion of magmatic fluids during the final stages of crystallisation, through a meteoric convection system, probably driven by waning magmatic heat, to much later mineralisation, concentrated in its vicinity due to thermal, tectonic and compositional properties of granite batholiths which encourage mineralisation long after magmatic heat has abated. Received: 3 April 1996 / Accepted: 5 May 1997     

Ion microprobe analysis of oxygen isotope ratios in quartz from Skye granite: healed micro-cracks,fluid flow,and hydrothermal exchange

 Quartz grains in hydrothermally altered granites from the Isle of Skye are highly heterogeneous and not equilibrated in oxygen isotope ratio at the 20 μm scale. Ion microprobe analysis of one grain shows a gradient of 13‰ over 400 μm and a greater range in δ ¹⁸O than all quartz previously analyzed on the Isle of Skye. Other crystals from the same outcrop are homogeneous. Digitized cathodoluminescence images reveal patterns of magmatic zoning and brittle fracturing not otherwise detectable. The ion probe analysis correlates low δ ¹⁸O values on a micro-scale to one set of healed cracks. Thus, quartz exchanges oxygen isotopes primarily by solution and reprecipitation along fractures, in contrast to more reactive feldspar that appears to exchange from the grain boundary inward. Macroscopic models of isotope exchange are not realistic for these rocks; the flow of hydrothermal fluids was heterogeneous, anisotropic and crack controlled. Received: 23 October 1995/Accepted: 9 April 1996