Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature‐induced grain coarsening. Despite this large‐scale trend, a considerable grain size scatter occurs on the outcrop‐scale indicating local influence of second‐order effects such as thermal perturbations, fluid flow and second‐phase particles. Second‐phase particles, whose sizes range from nano‐ to the micron‐scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron‐scale second‐phase particles, or to samples containing a large number of nano‐scale particles. The obtained data set suggests that the second phases induce a temperature‐controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D = C₂ e^Q*/RT(d_p/f_p)^m*, where C₂ is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio d_p/f_p, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C₂ and Q*. Besides variations in the average grain size, the presence of second‐phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian‐type distributions of the pure samples. In contrast, fluid‐enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable δ¹⁸O and δ¹³C isotope ratios in fluid‐affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid‐induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1. To summarize, considerable grain size variations of up to one order of magnitude can locally result from second‐order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.     

Petrogenesis and Mineralization of Two‐Stage A‐Type Granites in Jiuyishan,South China: Constraints from Whole‐rock Geochemistry,Mineral Composition and Zircon U‐Pb‐Hf Isotopes

The Jiuyishan complex massif, located in the northern section of the Nanling region, is a combination of five plutons, namely, the Xuehuading, Jinjiling, Pangxiemu, Shaziling and Xishan plutons. Whole-rock geochemistry, mineral electron microprobe analysis, zircon U-Pb dating and Hf isotope analysis were carried out for the Jinjiling and Pangxiemu plutons. The zircon U-Pb dating yields weighted mean ages of 152.9±0.9 Ma for the Jinjiling pluton and 151.7±1.5 Ma for the Pangxiemu pluton, with a narrow gap between them. The Jinjiling and Pangxiemu plutons both have geochemical characteristics of high SiO2, Al2O3, Na2O, K2O and low TiO2, MgO, CaO, P2O5 contents, with intense depletionS in Sr, Ba, Ti, Eu and enrichments?in Ga, FeOT and HFSE, and these characteristics reflect an A-type affinity. From the Jinjiling to the Pangxiemu plutons, the mineral composition of mica changes from lepidomelane to zinnwaldite, with increases in F, Li2O and Rb2O contents. The mineral composition of zircon changes from low Zr/Hf to high Zr/Hf, with increasing HfO2, P2O5 and UO2+ThO2+Y2O3 contents. The mineral compositions of feldspar indicate that the Pangxiemu pluton contains more alkali feldspar than the Jinjiling pluton. The whole-rock geochemistry and mineral compositions reveal a higher degree of differentiation for the Pangxiemu pluton. The nearly uniform εHf(t) indicates the same source region for the two plutons: both were derived from partial melting of the lower crust, with small contributions of mantle materials. In addition, higher F, lower Nb/Ta and Zr/Hf ratios in the Pangxiemu Pluton suggest a closer relationship with the rare metal mineralization than for the Jinjiling pluton.     

White mica 40Ar/39Ar age spectra and the timing of multiple episodes of high‐P metamorphic mineral growth in the Cycladic eclogite–blueschist belt,Syros, Aegean Sea,Greece

New geochronology from Syros in the Cycladic eclogite–blueschist belt, Aegean Sea, Greece, shows that ⁴⁰Ar/³⁹Ar geochronology consistently dates microstructural events in metamorphic rocks. We demonstrate that the age spectra depend on microstructure in a predictable and systematic way. Ages can be inferred by applying the method of asymptotes and limits to data from the step ‐ heating experiments. The results are consistent with previously published estimates for the timing of a sequence of distinct and discrete episodes of high ‐ P metamorphic mineral growth observed regionally across this belt. Arrhenius plots from these experiments imply that phengitic white mica is highly retentive of argon, and therefore (if these data can be extrapolated to the natural environment) the ages can be interpreted as recording the timing of episodic deformation and metamorphism. Porphyroblastic growth begins: (i) for omphacite – jadeite–eclogite facies parageneses at c. 53 Ma; and (ii) for garnet – glaucophane facies parageneses at c. 47 Ma. The Kini Shear Zone started as an extensional post ‐ epidote–albite‐transitional–blueschist facies shear zone that had completed operation by c. 31 Ma. The scatter in ages is due to the effect of deformation, recrystallization and multiple growth events in shear zones that continued operating for 3 – 6 million years from the start of each episode.     

Reactive flow of mixed CO2–H2O fluid and progress of calc-silicate reactions in contact metamorphic aureoles: insights from two-dimensional numerical modelling

Previous models of hydrodynamics in contact metamorphic aureoles assumed flow of aqueous fluids, whereas CO₂ and other species are also common fluid components in contact metamorphic aureoles. We investigated flow of mixed CO₂–H₂O fluid and kinetically controlled progress of calc‐silicate reactions using a two‐dimensional, finite‐element model constrained by the geological relations in the Notch Peak aureole, Utah. Results show that CO₂ strongly affects fluid‐flow patterns in contact aureoles. Infiltration of magmatic water into a homogeneous aureole containing CO₂–H₂O sedimentary fluid facilitates upward, thermally driven flow in the inner aureole and causes downward flow of the relatively dense CO₂‐poor fluid in the outer aureole. Metamorphic CO₂‐rich fluid tends to promote upward flow in the inner aureole and the progress of devolatilization reactions causes local fluid expulsion at reacting fronts. We also tracked the temporal evolution of P‐T‐X_CO2conditions of calc‐silicate reactions. The progress of low‐ to medium‐grade (phlogopite‐ to diopside‐forming) reactions is mainly driven by heat as the CO₂ concentration and fluid pressure and temperature increase simultaneously. In contrast, the progress of the high‐grade wollastonite‐forming reaction is mainly driven by infiltration of chemically out‐of‐equilibrium, CO₂‐poor fluid during late‐stage heating and early cooling of the inner aureole and thus it is significantly enhanced when magmatic water is involved. CO₂‐rich fluid dominates in the inner aureole during early heating, whereas CO₂‐poor fluid prevails at or after peak temperature is reached. Low‐grade metamorphic rocks are predicted to record the presence of CO₂‐rich fluid, and high‐grade rocks reflect the presence of CO₂‐poor fluid, consistent with geological observations in many calc‐silicate aureoles. The distribution of mineral assemblages predicted by our model matches those observed in the Notch Peak aureole.     

The effect of pre‐tectonic reaction and annealing extent on behaviour during subsequent deformation: insights from paired shear zones in the lower crust of Fiordland,New Zealand

Granulite facies pargasite orthogneiss is partially to completely reacted to garnet granulite either side of narrow (<20 mm) felsic dykes, in Fiordland, New Zealand, forming ~10–80 mm wide garnet reaction zones. The metamorphic reaction changed the abundance of minerals, and their shape and grain size distribution. The extent of reaction and annealing (temperature‐related coarsening and nucleation) is greatest close to the dykes, whereas further away the reaction is incomplete. As a consequence, grain size and the abundance of garnet decreases away from the felsic dykes over a few centimetres. The aspect ratios of clusters of S1 pyroxene and pargasite in the orthogneiss, which are variably reacted to post‐S1 garnet, decrease from high in the host, to near equidimensional close to the dyke. Post‐reaction deformation localized in the fine‐grained partially reacted areas. This produced a pattern of ‘paired’ shear zones located at the outer parts of the garnet reaction zone. Our study shows that grain size sensitive deformation occurs where the grain size is sufficiently reduced by metamorphic reaction. The weakening of the rock due to the change in grain size distribution outweighs the addition of nominally stronger garnet to the assemblage.     

Bed‐scale spatial patterns in dolomite abundance: Part II. Effect of varied fluid chemistry,flow rate,precursor mineralogy,temperature, textural heterogeneity,nucleation density and bed geometry

Reaction‐transport modelling shows that a lateral, metre‐scale pattern in dolomite abundance can form during replacive dolomitization of calcite and aragonite precursors by Mississippian, Triassic and modern seawaters advecting at >20 cm year^?1 in low temperature (≤60°C) hydrogeological systems. The modelled pattern develops best in beds >1 m thick with a relatively uniform nucleate density and a low variance in reactive surface area (for example, grainstones). These conditions suggest that a lateral pattern in dolomite abundance should be expected in dolomites formed by any shallow hydrodynamic system that advects dolomitizing fluid horizontally, including the down‐gradient portion of a reflux system, the oceanward reaches of geothermal (Kohout) convection systems, and the seaward reaches of a seawater entrainment zone below a freshwater aquifer. However, even in those systems, a pattern will not form in all dolomites. Pattern will be muted in thinner (≤ ca 50 cm thick) beds and non‐emergent where the precursor had a high variance in reactive surface area (for example, a skeletal wackestone) and/or large variation in nucleate density. Pattern also did not form at temperatures above ca 60°C, which implies that pattern should not be expected in dolomites formed at intermediate to deep burial depths. As the patterns are horizontal, they also should not be expected where dolomitizing fluids moved vertically (for example, reflux immediately below a brine source). Pattern metrics (short‐range correlation length, and wavelength and amplitude of the longer cyclic component) vary with flow rate, fluid chemistry, bed thickness, porosity, grain size, temperature and the flow complexity (one‐dimensional, two‐dimensional or three‐dimensional) within the bed. However, no modelled pattern produced metrics larger than those documented on dolomite outcrops. The results thus constrain the length scales of porosity and permeability variance to include in petrophysical models of dolomite reservoirs, as well as the geological scenarios in which to consider metre‐scale lateral petrophysical variability.