Morphological ripening of fluid inclusions and coupled zone-refining in quartz crystals revealed by cathodoluminescence imaging: Implications for CL-petrography, fluid inclusion analysis and trace-element geothermometry

Cathodoluminescence (CL) studies have previously shown that some secondary fluid inclusions in luminescent quartz are surrounded by dark, non-luminescent patches, resulting from fracture-sealing by late, trace-element-poor quartz. This finding has led to the tacit generalization that all dark CL patches indicate influx of low temperature, late-stage fluids. In this study we have examined natural and synthetic hydrothermal quartz crystals using CL imaging supplemented by in-situ elemental analysis. The results lead us to propose that all natural, liquid-water-bearing inclusions in quartz, whether trapped on former crystal growth surfaces (i.e., of primary origin) or in healed fractures (i.e., of pseudosecondary or secondary origin), are surrounded by three-dimensional, non-luminescent patches. Cross-cutting relations show that the patches form after entrapment of the fluid inclusions and therefore they are not diagnostic of the timing of fluid entrapment. Instead, the dark patches reveal the mechanism by which fluid inclusions spontaneously approach morphological equilibrium and purify their host quartz over geological time. Fluid inclusions that contain solvent water perpetually dissolve and reprecipitate their walls, gradually adopting low-energy euhedral and equant shapes. Defects in the host quartz constitute solubility gradients that drive physical migration of the inclusions over distances of tens of μm (commonly) up to several mm (rarely). Inclusions thus sequester from their walls any trace elements (e.g., Li, Al, Na, Ti) present in excess of equilibrium concentrations, thereby chemically purifying their host crystals in a process analogous to industrial zone refining. Non-luminescent patches of quartz are left in their wake. Fluid inclusions that contain no liquid water as solvent (e.g., inclusions of low-density H2O vapor or other non-aqueous volatiles) do not undergo this process and therefore do not migrate, do not modify their shapes with time, and are not associated with dark-CL zone-refined patches. This new understanding has implications for the interpretation of solids within fluid inclusions (e.g., Ti- and Al-minerals) and for the elemental analysis of hydrothermal and metamorphic quartz and its fluid inclusions by microbeam methods such as LA-ICPMS and SIMS. As Ti is a common trace element in quartz, its sequestration by fluid inclusions and its depletion in zone-refined patches impacts on applications of the Ti-in-quartz geothermometer.

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Noble gas, CFC and other geochemical evidence for the age and origin of the Bath thermal waters, UK

The city of Bath is a World Heritage site and its thermal waters, the Roman Baths and new spa development rely on undisturbed flow of the springs (45 °C). The current investigations provide an improved understanding of the residence times and flow regime as basis for the source protection. Trace gas indicators including the noble gases (helium, neon, argon, krypton and xenon) and chlorofluorocarbons (CFCs), together with a more comprehensive examination of chemical and stable isotope tracers are used to characterise the sources of the thermal water and any modern components. It is shown conclusively by the use of 39Ar that the bulk of the thermal water has been in circulation within the Carboniferous Limestone for at least 1000 years. Other stable isotope and noble gas measurements confirm previous findings and strongly suggest recharge within the Holocene time period (i.e. the last 12 kyr). Measurements of dissolved 85Kr and chlorofluorocarbons constrain previous indications from tritium that a small proportion (<5%) of the thermal water originates from modern leakage into the spring pipe passing through Mesozoic valley fill underlying Bath. This introduces small amounts of O2 into the system, resulting in the Fe precipitation seen in the King’s Spring. Silica geothermometry indicates that the water is likely to have reached a maximum temperature of between 69–99 °C, indicating a most probable maximum circulation depth of ∼3 km, which is in line with recent geological models. The rise to the surface of the water is sufficiently indirect that a temperature loss of >20 °C is incurred. There is overwhelming evidence that the water has evolved within the Carboniferous Limestone formation, although the chemistry alone cannot pinpoint the geometry of the recharge area or circulation route. For a likely residence time of 1–12 kyr, volumetric calculations imply a large storage volume and circulation pathway if typical porosities of the limestone at depth are used, indicating that much of the Bath-Bristol basin must be involved in the water storage.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

Tectono-metamorphic evolution of a nappe stack: A case study of the Swiss Alps

Fold-and-thrust belts are prominent structures that occur at the front of compressional orogens. To unravel the tectonic and metamorphic evolution of such complexes, kinematic investigations, quantitative microstructural analysis and geothermometry (calcite–graphite, calcite–dolomite) were performed on carbonate mylonites from thrust faults of the Helvetic nappe stack in Central Switzerland. Paleo-isotherms of peak temperature conditions and cooling stages (fission track) of the nappe pile were reconstructed in a vertical section and linked with the microstructural and kinematic evolution. Mylonitic microstructures suggest that under metamorphic conditions close to peak temperature, strain was highly localized within thrust faults where deformation temperatures spatially continuously increased in both directions, from N to S within each nappe and from top–down in the nappe stack, covering a temperature range of 180–380 °C. Due to the higher metamorphic conditions, thrusting of the lowermost nappe, the Doldenhorn nappe, was accompanied by a much more pronounced nappe internal ductile deformation of carbonaceous rock types than was the case for the overlying Wildhorn- and Gellihorn nappes. Ongoing thrusting brought the Doldenhorn nappe closer to the surface. The associated cooling resulted in a freezing in of the paleo-isotherms of peak metamorphic conditions. Contemporaneous shearing localized in the basal thrust, initially still in the ductile deformation regime and finally as brittle faulting and cataclasis inducing ultimately an inverse metamorphic zonation. With ongoing exhumation and the formation of the Helvetic antiformal nappe stack, a bending of large-scale tectonic structures (thrusts, folds), peak temperature isotherms and cooling isotherms occurred. While this local bending can directly be attributed to active deformation underneath the section investigated up to times of 2–3 ma, a more homogeneous uplift of the entire region is suggested for the very late and still active exhumation stage.