Fluid flow and mass transport at the Valentine wollastonite deposit, Adirondack Mountains, New York State

The Valentine wollastonite skarn in the north-west Adirondack Mountains, New York, is a seven million ton deposit which resulted from channellized infiltration of H₂O-rich, silica-bearing fluids. The wollastonite formed by reaction of these fluids with non-siliceous calcite marble. The skarn formed at the contact of the syenitic Diana Complex and was subsequently overprinted by Grenville-age granulite facies metamorphism and retrograde hydrothermal alteration during uplift. Calcite marbles adjacent to the deposit have generally high δ¹⁸O values (c. 21‰), typical of Grenville marbles which have not exchanged extensively with externally derived fluids. Carbon isotopic fractiona-tions between coexisting calcite and graphite in the marbles indicate equilibration at 675d? C, consistent with the conditions of regional metamorphism. Oxygen isotopic ratios from wollastonite skarn are lower than in the marbles and show a 14‰ variation (-1‰ to 13‰). Some isotopic heterogeneity is preserved from skarn formation, and some represents localized exchange with low-δ¹⁸O retrograde fluids. Detailed millimetre- to centimetre-scale isotopic profiles taken across skarn/marble contacts reveal steep δ¹⁸O gradients in the skarn, with values increasing towards the marble. The gradients reflect isotopic evolution of the fluid as it reacted with high δ¹⁸O calcite to form wollastonite. Calcite in the marble preserves high δ¹⁸O values to within <5 mm of the skarn contact. The preservation of high δ¹⁸O values in marbles at skarn contacts and the disequilibrium fractionation between wollastonite skarn and calcite marble across these contacts indicate that the marbles were not infiltrated with significant quantities of the fluid. Thus, the marbles were relatively impermeable during both the skarn formation and retrograde alteration. Skarn formation may have been episodic and fluid flow was either chaotic or dominantly parallel to lithological contacts. Although these steep isotope gradients resemble fluid infiltration fronts, they actually represent the sides of the major flow system. Because chromatographic infiltration models of mass transport require the assumption of pervasive fluid flow through a permeable rock, such models are not applicable to this hydrothermal system and, by extension, to many other metamorphic systems where low-permeability rocks restrict fluid migration pathways. Minimum time-integrated fluid fluxes have been calculated at the Valentine deposit using oxygen isotopic mass balance, reaction progress of fluid buffering reactions, and silica mass balance. All three approaches show that large volumes of fluid were necessary to produce the skarn, but silica mass balance calculations yield the largest minimum flux and are hence the most realistic.     

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     

Phase equilibrium and stable isotope constraints on the formation of metasomatic garnet-vesuvianite veins (SW Adamello,N Italy)

 Metasomatic garnet-vesuvianite veins occur within the contact metamorphic marble sequence of the Lower Triassic Prezzo formation in a narrow, 1–5 m wide zone along an intrusive marble-granodiorite contact at the southwestern border of the Tertiary Adamello batholith. The metasomatic mineral assemblage is comprised of garnet, vesuvianite, clinopyroxene, wollastonite, and pyrrhotite, which were precipitated from the vein-forming fluid in a preexisting calcite matrix at conditions of about 2800 bars and 630° C. The veins are enriched in silicon, aluminum, iron, magnesium, titanium and depleted in calcium with respect to the unaltered contact metamorphic marble. Graphite, which is present in the unaltered Prezzo Marble is absent in the veins. Irregularly shaped mineralogically distinct zones with different degrees of silicification can be distinguished within the veins. The isotopic compositions of calcite (cc) in the unaltered marble are about δ¹⁸O (SMOW; Standard mean Ocean Water)=21.0‰ and δ¹³C(PDB; Peedee belemnite)=0.0‰. They are reset to significantly lower values within the veins, where δ¹⁸O_cc is 15.0 to 16.0‰ and δ¹³C_cc is −4.5 to −3.5‰. The isotopic front coincides with an abrupt change in the microscopic texture of matrix carbonate which occurs at the sharp boundary between graphite-bearing and graphite-free material. Within the veins the oxygen isotope fractionation between calcite and garnet (gar) varies systematically with distance from highly silicified zones. The variations in Δ¹⁸O_cc-gar are as large as 2‰, on a millimeter scale, indicating garnet-calcite isotopic disequilibrium. Vein formation was due to the infiltration of a water rich fluid of magmatic provenance into the carbonate country rock along fractures. Removal of graphite from the wall rock by dissolution through the metasomatic fluid induced recrystallization of matrix calcite. Permeability was enhanced during calcite recrystallization facilitating material transport into the wall rock and metasomatic alteration. Vein garnet was precipitated in isotopic equilibrium with the metasomatic fluid. The isotopic composition of preexisting calcite was initially out of equilibrium with the vein-forming fluid and it was shifted towards equilibrium by surface-reaction controlled calcite-fluid isotopic exchange during calcite recrystallization. Due to the short lifetime of the metasomatic system, calcite-fluid isotopic equilibrium was generally not attained. Within the veins, oxygen and carbon transport was fast relative to mineral-fluid exchange of their isotopes and the geometry of the isotopic pattern is largely controlled by the kinetics of mineral-fluid exchange. Received: 16 June 1994/Accepted: 20 May 1995     

Transport of Pb and Sr in leaky aquifers of the Bufa del Diente contact metamorphic aureole, north-east Mexico

The Sr and Pb isotopes from the 31.6 ± 0.3 Ma (2σ) old Diente del Bufa alkali syenite, northeastern Mexico, and marbles of its contact aureole were used to trace the sources and the mobility of these metals during hydrothermal activity. Chert layers form aquifers within the marbles. The marbles represent aquitards. During fluid-wallrock reaction, the chert layers developed wollastonite rims. Early wollastonite rims have Sr and Pb isotopic compositions similar to those of their immediate host marbles, which indicates that the isotopic composition of Sr and Pb is initially buffered by the marble. Later wollastonite and other replacement minerals rimming the aquifer have Sr and Pb isotopic compositions that carry with time increasingly larger contributions from the high-salinity magmatic brine. The Sr and Pb contributions from the alkali syenite can be traced isotopically for more than 90 m away from the contact of the intrusion. In contrast, Sr and Pb originating from the alkali syenite are traceable within the marbles only for 3 to 5 cm from the aquifer-marble boundary. This distance is comparable to the spatial distribution of isotopic alterations of C and O implying that Sr and Pb were transported into the marbles through a fluid phase. The isotopic variation of Sr, Pb, C, and O across the aquifer-marble profiles reflects infiltration as a transport mechanism rather than diffusion. Because Sr and Pb are minor components in both the infiltrating fluid and the rock and because their concentrations are strongly affected by the distribution coefficients among the solid phases present, there is little correlation between the isotopic compositions of the trace elements Sr and Pb and those of C and O, which are major components in fluid and rock. Very thin meta-argillite rinds at the outer margin of the aquifer represent residual material after the dissolution of calcite. They are distinctly enriched in Rb, Sr, and U. The Rb and Sr are to some extent residual from the original limestone mineralogy, whereas U is dominantly derived from the magmatic fluid and leaked from the aquifer with the escaping immiscible CO₂-rich H₂O-CO₂ fluid that was produced by decarbonation. The ²³⁸U/²⁰⁴Pb values ranging from 100 to 250 and distinctly lowered Th/U in the meta-argillite rims (1) demonstrate that U was transported with the magmatic fluid along the aquifer and (2) imply that during unmixing of the highly saline magmatic fluid U fractionated into the CO₂-rich H₂O-CO₂ fluid from which it precipitated selectively in the meta-argillite band across the aquifer. Radioautographs demonstrate that the upper meta-argillite rim has 20 to 40 times more U than the lower rim, which implies that 20 to 40 times more CO₂-rich H₂O-CO₂ fluid has left through the upper aquifer contact. Received: 30 September 1997 / Accepted: 15 December 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     

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.     

Prograde and retrograde fluid-rock interaction in calc-silicates northwest of the Idaho batholith: stable isotopic evidence

Carbon and oxygen isotopic analyses of silicate and carbonate minerals indicate that isotopic compositions in metasediments of the Wallace Formation (Belt Supergroup) exposed northwest of the Idaho batholith have been affected by both prograde and retrograde fluid-rock interaction. Silicates retain isotopic fractionations that reflect equilibration at peak metamorphic temperatures. In contrast, calcite oxygen isotopic compositions range from δ¹⁸O(Calcite)=+2.3 to +18.6‰ SMOW (standard mean oceanic water) and indicate that some calcites have exchanged with low-δ¹⁸O meteorichydrothermal fluids. Values of Δ¹⁸O (Quartz-Calcite) as large as +15.5 clearly indicate that the isotopic depletion of these calcites postdates the peak of regional metamorphism. Carbon isotopic compositions of ¹⁸O-depleted calcites are not significantly shifted relative to δ¹³C values in undepleted calcites, suggesting that the retrograde fluid was carbon-poor. Petrographically, retrograde fluid-rock interaction is associated with the occurrence of fine-grained, highly-luminescent calcite overgrowths on less-luminescent, metamorphic calcites, slight to moderate argillic alteration, and pseudomorphing of scapolite porphyroblasts by fine-grained albite. Retrograde isotopic depletions may be related to shallow meteoric-hydrothermal fluid flow developed around the Idaho batholith after intrusion and rapid uplift of the terrane. Peak metamorphic isotopic compositions in the Wallace Formation reflect mineralogically heterogeneous protolith compositions and isotopic fractionation due to devolatilization and/or infiltration. Variability in oxygen isotopic compositions on the order of 4–6‰ within the same rock type can be attributed to the combined effects of inherited isotopic compositions and isotopic shifts resulting from prograde devolatilization. Isotopic and compositional heterogeneity on the scale of mm to m precludes generalization of isotopic gradients on a regional scale. The isotopic data presented here, and metamorphic fluid compositions determined in previous studies, are best reconciled with heterogeneous bulk compositions, dominantly channelized prograde and retrograde fluid flow, and locally low fluid-rock ratios.     

Strontium and oxygen isotope profiles across marble-silicate contacts,Lizzies Basin,East Humboldt Range,Nevada: constraints on metamorphic permeability contrasts and fluid flow

 Sr isotope profiles across marble-silicate rock contacts are used in conjunction with previously published oxygen isotope profiles to constrain fluid movement, porosity and permeability contrasts in migmatitic rocks from Lizzies Basin in the East Humboldt Range, Nevada. The ¹⁸O/¹⁶O systematics in the high-grade sequence have been interpreted to reflect infiltration of ∼2×10² m³/m² of a relatively low ¹⁸O hydrous fluid through the sequence, but with preservation of δ¹⁸O anomalies in thin marble bands due to a 30-fold lower porosity in the marble compared with silicate rocks (Wickham and Peters 1992). The Sr isotope profiles confirm that tracer exchange between marble and silicate rock was primarily by diffusion, and in one case, indicate that porosities differed by less than a factor of four in the ∼10 cm boundary layer which exhibits diffusive modification of ⁸⁷Sr/⁸⁶Sr ratios. This contrasts with modelling of the oxygen isotope profiles which imply porosity contrasts >10 for one marble band and >50 for a second marble band. Either strontium and oxygen isotope diffusion reflect different events (possible if fluid Sr contents varied with time) or porosity varied substantially with the silicate rocks. Oxygen isotope profiles in the deeper part of the metamorphic section in which δ¹⁸O values of silicate rocks have been homogenised and lowered, indicate similar diffusion distances (and thus porosity-time evolution) to oxygen isotopic profiles higher in the section. Comparison of strontium and oxygen isotope diffusion distances constrains fluid Sr contents to between ∼50 and ∼500 ppm deep in the section, but less than ∼10 ppm higher in the section. The difference is related to release of relatively saline, Sr-rich fluids, by the abundant leucogranites and associated skarns deep in the section (cf. Peters and Wickham 1995). Received: 9 December 1994/Accepted: 13 April 1995     

Fluid-rock history of granulite facies humite-marbles from Ambasamudram, southern India.

An extensive humite‐bearing marble horizon within a supracrustal sequence at Ambasamudram, southern India, was studied using petrological and stable isotopic techniques to define its metamorphic history and fluid characteristics. At peak metamorphic temperatures of 775±73°C, based on calcite‐graphite carbon isotope thermometry, the mineral assemblages suggest layer‐by‐layer control of fluid compositions. Clinohumite + calcite‐bearing assemblages suggest X_CO2 < 0.4 (at 700°C and 5 kbar), calcite + forsterite + K‐feldspar‐bearing assemblages suggest X_CO2>0.9 (at 790°C); and local wollastonite + scapolite + grossular‐bearing zones formed at X_CO2 of c. 0.3. Retrograde reaction textures such as scapolite + quartz symplectites after feldspar and calcite and replacement of dolomite + diopside or tremolite+dolomite after calcite+forsterite or calcite+clinohumite are indicative of retrogression under high X_CO2 conditions. Calcite preserves late Proterozoic carbon and oxygen isotopic signatures and the marble lacks evidence for extensive retrograde fluid infiltration, while during prograde metamorphism the possible infiltration of aqueous fluids did not produce significant isotopic resetting. Isotopic zonation of calcite and graphite grains was likely produced by localized CO₂ fluid infiltration during retrogression. Contrary to the widespread occurrence of humite‐marbles related to retrograde aqueous fluid infiltration, the Ambasamudram humite‐marbles record a prograde‐to‐peak metamorphic humite formation and retrogression under conditions of low X_H2O.     

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     

The formation of columnar fiber texture in wollastonite rims by induced stress and implications for diffusion-controlled corona growth

 The evolution of columnar fiber texture was studied in wollastonite reaction rims synthesized by the reaction calcite + quartz=wollastonite + CO₂. Experiments were performed at 850 to 950 °C at 100 MPa in dry CO₂ and were evaluated by scanning and transmission electron microscopy. Rim growth rates are interpreted as controlled by the diffusion of the SiO₂ component through the rims from the quartz–wollastonite to the wollastonite–calcite interface. The temperature dependence of rim growth rates yields an apparent activation energy of 314 ± 53 kJ mol⁻¹. The columnar fibrous wollastonite crystallizes at the quartz–wollastonite interface and comprises the largest parts of the rims. Ultimately, at the growth front strain contrast centers are present in the quartz. The strained volume extends about 200 nm into the quartz grains. We suggest that this might signify deformation of the quartz lattice due to wollastonite crystallization. Wollastonite fiber thickness was measured from TEM images along traverses that represent intermediate positions of the growth front during the experiments. The average thickness is in the 100–200 nm range. Fiber thickness increases with increasing growth temperature. At a given temperature, the thickness of the fibers at the growth front slightly decreases with time, i.e., the number of fiber tips per unit area in the growth front increases. The decrease of the fiber thickness is well fitted by a parabolic rate law. The generation of the columnar fiber texture is interpreted as an effect of induced stresses at the growth front, resulting from the volume increase due to the local reaction. This volume increase forces SiO₂ to diffuse along the growth front to the grain boundaries between the wollastonite fibers. These serve as fast diffusion pathways through the rims. The fiber thickness monitors the diffusion distances in the growth front and thus the height of the induced stress gradients. Since interface reactions are usually associated with volume changes, growth rates of reaction rims and zones in coronas are not only controlled by the diffusive mobility of the components but also by the volume restraints on the interface reactions. Received: 19 July 2002 / Accepted: 14 February 2003     

Elemental dispersion and stable isotope fractionation during reactive fluid-flow and fluid immiscibility in the Bufa del Diente aureole,NE-Mexico: evidence from radiographies and Li,B, Sr,Nd, and Pb isotope systematics

The 31.6±0.3 Ma old Bufa del Diente alkali-syenite (NE Mexico) intruded a sequence of Cretaceous limestones with intercalated sub-horizontal chert layers. The cherts acted as aquifers that facilitated transport of brines and pegmatitic melts within the shallow-level (<1 kbar) contact-metamorphic aureole. Fluid-driven reactions between chert and marble wallrock, and the influx of late melts and various fluids gave rise to distinct chemical and isotopic signatures within the aquifer and across the zones of infiltration and fluid-driven reaction. Aqueous brines of magmatic origin produced thick wollastonite mantles around the chert layers. Wollastonite formation occurred at the expense of limestone and chert and generated CO₂. This CO₂-induced fluid unmixing into an aqueous brine and a low-density CO₂-rich fluid, which was lost to the overlying marble where it oxidized organic matter and caused ¹³C and ¹⁸O shifts in a zone some 5–10 cm wide. After wollastonite formation, the chert aquifers were locally intruded by pegmatite veins carrying alkali feldspar, quartz, aegirine-augite, eudialyte, zircon, and apatite. Aqueous fluids that exsolved during crystallization of the pegmatite veins escaped along late cross-fractures and migrated along the inner and outer borders of the wollastonite margins. Chemical dispersion patterns of U, Al, Na + K, P, S, Fe, and REE across the chert-to-marble boundary and its metasomatic rims are shown by autoradiography and neutron-induced radiography. Scavenging of cations at mineralogical contacts and cation transport into the marbles occurred only on the mm to cm scale. Isotopic data for Pb and Sr across a simple metachert-marble boundary and for Pb, Sr, Nd, B, and Li across a metachert-pegmatite-marble sequence demonstrate the following: (1) The Pb and Sr isotopic signature of early fluids was buffered by the carbonate wallrock. Only late fluids, shielded from wallrock interaction by a wollastonite mantle, variably preserved a memory of their initial magmatic signature. (2) Since the Nd isotope signature of marble and chert is bound to calcite and clay minerals, systematic shifts to unradiogenic Nd in marble reflect loss of carbonate-bound Nd as the wollastonite margin is approached. Nd in the wollastonite margin is dominated by Nd originally bound to clay minerals. The later emplacement of the pegmatite, which carried the Nd isotope signature of its alkali-syenite source, had little effect on the Nd isotopic composition of the wollastonite rim. (3) Although the Li and B isotopic compositions reflect the alkali-syenite source, they are also affected by isotopic fractionation and partitioning between melt, fluid, and solids.Editorial responsibility J. Hoefs     

Mechanisms of oxygen isotopic exchange and isotopic evolution of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O-depleted periclase zone marbles in the Alta aureole,Utah: insights from ion microprobe analysis of calcite

Oxygen isotope ratios have been measured by ion microprobe and millimeter-scale dental drill along detailed sampling traverses across the boundary between periclase-bearing (δ¹⁸O = 11.8‰) and periclase-free (δ¹⁸O = 17.2‰) marble layers in the periclase (Per) zone of the Alta Stock aureole, Utah. These data define a steep, coherent gradient in δ¹⁸O that is displaced a short distance (~4 cm) into the periclase-free (Cal + Fo) layer. SEM and ion microprobe analyses show two isotopically and texturally distinct types of calcite at the grain scale. Clear (well polished) calcite grains are isotopically homogeneous (within analytical uncertainty; ±0.27‰, 2SD). More poorly polished (pitted), texturally retrograde ‘turbid’-looking calcite has lower and more variable δ¹⁸O values, and replaces clear calcite along fractures, cleavage traces or grain boundaries. Despite significant lowering of the δ¹⁸O values in calcite throughout both layers during prograde metamorphism, ion microprobe analyses indicate that individual clear calcite grains are now isotopically homogeneous across the entire gradient in δ¹⁸O. Diffusion calculations indicate that conservative time scales required for isotopic homogenization of calcite grains by volume diffusion, 30,000–62,000 years at 575–600°C, exceed significantly the timescale (~1,250 years) estimated for the prograde development of the δ¹⁸O gradient at the boundary between these two marble layers. The ion microprobe data and these diffusion calculations suggest instead that surface reaction mechanisms accompanying recrystallization are responsible for the observed oxygen isotope homogeneity of these calcite grains. Thus, the ion microprobe data are consistent with the formation of calcite in oxygen isotope exchange equilibrium with infiltrating fluid during prograde reaction and recrystallization. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The correlation of reaction and isotope fronts and the mechanism of metamorphic fluid flow

 Infiltration of a metabasite sill from Islay, Scotland by an H₂O-CO₂ fluid caused (1) modification of δ¹⁸O and (2) carbonation at the sill margins. Maps of δ¹⁸O and reaction progress were constructed from a 20 × 47.7 metre sample grid across the sill. The grid consisted of 300 samples, spaced at m, dm and cm intervals, many of which were analysed for both δ¹⁸O and reaction progress. The δ¹⁸O was determined by laser fluorination of whole rock silicate powders and reaction progress was determined by rapid field-based measurement of % calcite (“fizz-o-meter”, Skelton et al. 1995). Reaction and isotope fronts outlined tube-like features that emanate from the sill margin and discrete nodes that, although detached from the sill margin in two dimensions, are thought to represent sections through similar tubes in three dimensions. We envisage that these protrusions are the fossil record of metamorphic “fluid pathways” whereby fluid permeated the sill. Isotope and reaction fronts are found to correlate spatially as predicted by a modified form of the chromatographic equation which describes this envisaged geometry, that is where isotopic and reactive transport in the fluid phase are facilitated by advection along specific fluid pathways and transverse diffusion in the surrounding rock. These fluid pathways consist of bundles of anastomosing grain boundary channels or micro-cracks, which are thought to propagate through transient cyclic infiltration, reaction, porosity enhancement and fracturing. This mechanism is self-perpetuating and accentuates random perturbations at the sill margin to form the observed tubes. We argue that this is the earliest stage of the infiltration process which has affected metabasites of the SW Scottish Highlands and that subsequent shear deformation of the reacted rims of these pathways, has caused their re-orientation and juxtaposition to form the reacted sill margins described by Skelton et al. (1995). Received: 17 February 1998 / Accepted: 6 December 1999     

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     

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     

Isotope geochemistry and fluid inclusion study of skarns from Vesuvius

Summary We present new mineral chemistry, fluid inclusion, stable carbon and oxygen, as well as Pb, Sr, and Nd isotope data of Ca-Mg-silicate-rich ejecta (skarns) and associated cognate and xenolithic nodules from the Mt. Somma-Vesuvius volcanic complex, Italy. The typically zoned skarn ejecta consist mainly of diopsidic and hedenbergitic, sometimes “fassaitic” clinopyroxene, Mg-rich and Ti-poor phlogopite, F-bearing vesuvianite, wollastonite, gehlenite, meionite, forsterite, clinohumite, anorthite and Mg-poor calcite with accessory apatite, spinell, magnetite, perovskite, baddeleyite, and various REE-, U-, Th-, Zr- and Ti-rich minerals. Four major types of fluid inclusions were observed in wollastonite, vesuvianite, gehlenite, clinopyroxene and calcite: a) primary silicate melt inclusions (T_HOM = 1000–1050 °C), b) CO₂ ± H₂S-rich fluid inclusions (T_HOM = 20–31.3 °C into the vapor phase), c) multiphase aqueous brine inclusions (T_HOM = 720–820 °C) with mainly sylvite and halite daughter minerals, and d) complex chloride-carbonate-sulfate-fluoride-silicate-bearing saline-melt inclusions (T_HOM = 870–890 °C). The last inclusion type shows evidence for immiscibility between several fluids (silicate melt – aqueous chloride-rich liquid – carbonate/sulfate melt?) during heating and cooling below 870 °C. There is no evidence for fluid circulation below 700 °C and participation of externally derived meteoric fluids in skarn formation. Skarns have considerably variable ²⁰⁶Pb/²⁰⁴Pb (19.047–19.202), ²⁰⁷Pb/²⁰⁴Pb (15.655–15.670), and ²⁰⁸Pb/²⁰⁴Pb (38.915–39.069) and relatively low ¹⁴³Nd/¹⁴⁴Nd (0.51211–0.51244) ratios. The carbon and oxygen isotope compositions of skarn calcites (δ¹³C_V-PDB = −5.4 to −1.1‰; δ18O_V-SMOW = 11.7 to 16.4‰) indicate formation from a ¹⁸O- and ¹³C-enriched fluid. The isotope composition of skarns and the presence of silicate melt inclusion-bearing wollastonite nodules suggests assimilation of carbonate wall rocks by the alkaline magma at moderate depths (< 5 km) and consequent exsolution of CO₂-rich vapor and complex saline melts from the contaminated magma that reacted with the carbonate rocks to form skarns. Received March 1, 2000; revised version accepted November 2, 2000     

Limited fluid-rock interaction at marble-gneiss contacts during Cretaceous granulite-facies metamorphism,Seward Peninsula,Alaska

Stable-isotope profiles show that flat-lying marble units acted as impermeable barriers to upward fluid flow in transitional amphibolite-granulite grade rocks of the Kigluaik Mountains, Seward Peninsula, Alaska. The degree of permeability is related to the composition of the marble. The margin of a thick pure dolomite marble chemically reacted with underlying metasyenite (aH₂O=0.2) to form a 2 cm boundary layer of calcite + forsterite by introduction of SiO₂. No fluid penetrated past this reaction front, although the high temperature of metamorphism (800°C) allowed transport of carbon and oxygen isotopes for an additional 2 cm by diffusion through the solid dolomite. A second marble with a higher silica content underwent more decarbonation, which enhanced porosity and lead to a greater extent of isotope transport (2–3 m) in contact with quartzo-feld-spathic gneiss below. An estimate of total fluid flux across the bottom of this marble layer based on the shape of the isotope profile is 1 cm³/cm² directed down, out of the marble. At two other marble-gneiss contacts steep isotopic gradients coincide with lithologic contacts, indicating very little cross-lithology fluid flow. The extent of diffusional transport of isotopes in the marbles is limited and interpreted as indicating the transient presence of a pore fluid, generated by thermally driven devolatilization reactions. No wholesale pervasive advection of C-O-H fluid occurred across the thick, continuous, marble units near the exposed base of the Kigluaik Group section during the entire regional metamorphic cycle. Activities of pore-fluid species were controlled by internal processes. Movement of volatiles and stable-isotopes between contrasting rock-types was dominantly diffusive. Channelized fluid pathways through the marble units developed during uplift and cooling but were not present during peak metamorphism. Heating of the section occurred by conduction, probably from an underlying magma source, and not by advection of a C-O-H fluid.     

Wollastonite at the Sterling Hill Fe–Zn–Mn ore body, Ogdensburg, New Jersey

Summary Wollastonite occurs abundantly at the Sterling Hill Fe–Zn–Mn ore deposit, Ogdensburg, New Jersey, one of the few occurrences of wollastonite in regionally metamorphosed rocks; it is absent from the surrounding Franklin marble. Wollastonite occurs in two distinct bands along the inner margins of the synclinal ore deposit. Minerals associated with wollastonite include calcite, grossular-andradite, diopsidic pyroxene, alkali feldspar, and rarely vesuvianite, quartz or bustamite. Assuming the generally accepted values of 750°C at 5kbar at Sterling Hill during metamorphism in the Grenville Orogeny, thermodynamic modeling of reactions involving garnet and wollastonite suggest XCO₂ 0.35 in the wollastonite-bearing rocks. Infiltrating metamorphic fluid rich in H₂O was necessary for the formation of wollastonite; at XCO₂ of 0.35, the calculated minimum volumetric water:rock ratio is 0.51. The source of the water is believed to be the dehydration of water-rich phases in adjacent ores or mafic rocks. The chemical compositions, textures, stratigraphy, and calculated metamorphic conditions show that wollastonite formed from calcite and quartz at the peak of the Grenville Orogeny.Present address: Maryland State Highway AdministrationReceived August 18, 2002; revised version accepted February 5, 2003     

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