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.

Long-term fluid circulation in extensional faults in the central Catalan Coastal Ranges: P-T constraints from neoformed chlorite and K-white mica

The neoformation of chlorite and K-white mica in fault rocks from two main faults of the central Catalan Coastal Ranges, the Vallès and the Hospital faults, has allowed us to constrain the P–T conditions during fault evolution using thermodynamic modeling. Crystallization of M1 and M2 muscovite and microcline occured as result of deuteric alteration during the exhumation of the pluton (290 °C > T > 370 °C) in the Permian. After that, three tectonic events have been distinguished. The first tectonic event, attributed to the Mesozoic rifting, is characterized by precipitation of M3 and M4 phengite together with chlorite and calcite C1 at temperatures between 190 and 310 °C. The second tectonic event attributed to the Paleogene compression has only been identified in the Hospital fault with precipitation of low-temperature calcite C2. The shortcut produced during inversion of the Vallès fault was probably the responsible for the lack of neoformed minerals within this fault. Finally, the third tectonic event, which is related to the Neogene extension, is characterized in the Vallès fault by a new generation of chlorite, associated with calcite C4 and laumontite, formed at temperatures between 125 and 190 °C in the absence of K-white mica. Differently, the Hospital fault is characterized by the precipitation of calcite C3 during the syn-rift stage at temperatures around 150 °C and by low-temperature fluids precipitating calcites C5, C6 and PC1 during the post-rift stage. During the two extensional events (Mesozoic and Neogene), faults acted as conduits for hot fluids producing anomalous high geothermal gradients (50 °C/km minimum).

P-T-t estimation of syn-kinematic strain in low-grade quartz-feldspar bearing rocks using thermodynamic modeling and 40Ar/39Ar dating techniques: example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)

Pressure–Temperature–time (P–T–t) estimates of the syn-kinematic strain at the peak-pressure conditions reached during shallow underthrusting of the Briançonnais Zone in the Alpine subduction zone was made by thermodynamic modelling and 40Ar/39Ar dating in the Plan-de-Phasy unit (SE of the Pelvoux Massif, Western Alps). The dated phengite minerals crystallized syn-kinematically in a shear zone indicating top-to-the-N motion. By combining X-ray mapping with multi-equilibrium calculations, we estimate the phengite crystallization conditions at 270 ± 50 °C and 8.1 ± 2 kbar at an age of 45.9 ± 1.1 Ma. Combining this P–T–t estimate with data from the literature allows us to constrain the timing and geometry of Alpine continental subduction. We propose that the Briançonnais units were scalped on top of the slab during ongoing continental subduction and exhumed continuously until collision.

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.

Experimental dissolution of molybdenum-sulphides at low oxygen concentrations: A first-order approximation of late Archean atmospheric conditions

The abundance of atmospheric oxygen and its evolution through Earth's history is a highly debated topic. The earliest change of the Mo concentration and isotope composition of marine sediments are interpreted to be linked to the onset of the accumulation of free O2 in Earth's atmosphere. The O2 concentration needed to dissolve significant amounts of Mo in water is not yet quantified, however. We present laboratory experiments on pulverized and surface-cleaned molybdenite (MoS2) and a hydrothermal breccia enriched in Mo-bearing sulphides using a glove box setup. Duration of an experiment was 14 days, and first signs of oxidation and subsequent dissolution of Mo compounds start to occur above an atmospheric oxygen concentration of 72 ± 20 ppmv (i.e., 2.6 to 4.6 × 10−4 present atmospheric level (PAL)). This experimentally determined value coincides with published model calculations supporting atmospheric O2 concentrations between 1 × 10−5 to 3 × 10−4 PAL prior to the Great Oxidation Event and sets an upper limit to the molecular oxygen needed to trigger Mo accumulation and Mo isotope variations recorded in sediments. In combination with the published Mo isotope composition of the rock record, this result implies an atmospheric oxygen concentration prior to 2.76 Ga of below 72 ± 20 ppmv.

Long-term summer temperature variations in the Pyrenees from detrended stable carbon isotopes

Substantial effort has recently been put into the development of climate reconstructions from tree-ring stable carbon isotopes, though the interpretation of long-term trends retained in such timeseries remains challenging. Here we use detrended δ13C measurements in Pinus uncinata tree-rings, from the Spanish Pyrenees, to reconstruct decadal variations in summer temperature back to the 13th century. The June-August temperature signal of this reconstruction is attributed using decadally as well as annually resolved, 20th century δ13C data. Results indicate that late 20th century warming has not been unique within the context of the past 750 years. Our reconstruction contains greater am-plitude than previous reconstructions derived from traditional tree-ring density data, and describes particularly cool conditions during the late 19th century. Some of these differences, including early warm periods in the 14th and 17th centuries, have been retained via δ13C timeseries detrending - a novel approach in tree-ring stable isotope chronology development. The overall reduced variance in earlier studies points to an underestimation of pre-instrumental summer temperature variability de-rived from traditional tree-ring parameters.

Water-use efficiency and transpiration across European forests during the Anthropocene

The Earth’s carbon and hydrologic cycles are intimately coupled by gas exchange through plant stomata1, 2, 3. However, uncertainties in the magnitude4, 5, 6 and consequences7, 8 of the physiological responses9, 10 of plants to elevated CO2 in natural environments hinders modelling of terrestrial water cycling and carbon storage11. Here we use annually resolved long-term δ13C tree-ring measurements across a European forest network to reconstruct the physiologically driven response of intercellular CO2 (Ci) caused by atmospheric CO2 (Ca) trends. When removing meteorological signals from the δ13C measurements, we find that trees across Europe regulated gas exchange so that for one ppmv atmospheric CO2 increase, Ci increased by ~0.76 ppmv, most consistent with moderate control towards a constant Ci/Ca ratio. This response corresponds to twentieth-century intrinsic water-use efficiency (iWUE) increases of 14 ± 10 and 22 ± 6% at broadleaf and coniferous sites, respectively. An ensemble of process-based global vegetation models shows similar CO2 effects on iWUE trends. Yet, when operating these models with climate drivers reintroduced, despite decreased stomatal opening, 5% increases in European forest transpiration are calculated over the twentieth century. This counterintuitive result arises from lengthened growing seasons, enhanced evaporative demand in a warming climate, and increased leaf area, which together oppose effects of CO2-induced stomatal closure. Our study questions changes to the hydrological cycle, such as reductions in transpiration and air humidity, hypothesized to result from plant responses to anthropogenic emissions.

Lithium isotope constraints on crust–mantle interactions and surface processes on Mars

Lithium abundances and isotope compositions are reported for a suite of martian meteorites that span the range of petrological and geochemical types recognized to date for Mars. Samples include twenty-one bulk-rock enriched, intermediate and depleted shergottites, six nakhlites, two chassignites, the orthopyroxenite Allan Hills (ALH) 84001 and the polymict breccia Northwest Africa (NWA) 7034. Shergottites unaffected by terrestrial weathering exhibit a range in δ7Li from 2.1 to 6.2‰, similar to that reported for pristine terrestrial peridotites and unaltered mid-ocean ridge and ocean island basalts. Two chassignites have δ7Li values (4.0‰) intermediate to the shergottite range, and combined, these meteorites provide the most robust current constraints on δ7Li of the martian mantle. The polymict breccia NWA 7034 has the lowest δ7Li (−0.2‰) of all terrestrially unaltered martian meteorites measured to date and may represent an isotopically light surface end-member. The new data for NWA 7034 imply that martian crustal surface materials had both a lighter Li isotope composition and elevated Li abundance compared with their associated mantle. These findings are supported by Li data for olivine-phyric shergotitte NWA 1068, a black glass phase isolated from the Tissint meteorite fall, and some nakhlites, which all show evidence for assimilation of a low-δ7Li crustal component. The range in δ7Li for nakhlites (1.8 to 5.2‰), and co-variations with chlorine abundance, suggests crustal contamination by Cl-rich brines. The differences in Li isotope composition and abundance between the martian mantle and estimated crust are not as large as the fractionations observed for terrestrial continental crust and mantle, suggesting a difference in the styles of alteration and weathering between water-dominated processes on Earth versus possibly Cl–S-rich brines on Mars. Using high-MgO shergottites (>15 wt.% MgO) it is possible to estimate the δ7Li of Bulk Silicate Mars (BSM) to be 4.2 ± 0.9‰ (2σ). This value is at the higher end of estimates for the Bulk Silicate Earth (BSE; 3.5 ± 1.0‰, 2σ), but overlaps within uncertainty.

The Basement of the Deccan Traps and Its Madagascar Connection: Constraints from Xenoliths

Paleogeographic reconstructions of India and Madagascar before their late Cretaceous rifting juxtapose the Antongil Block of Madagascar against the Deccan Traps of India, indicating that the Western Dharwar Craton extends below the Deccan lavas. Some recent studies have suggested that the South Maharashtra Shear Zone along the northern Konkan coast of India limits the northern extent of the Western Dharwar Craton, implying that the craton does not extend below the Deccan Traps, raising a question mark on paleogeographic reconstructions of India and Madagascar. The continuity of the Western Dharwar Craton north of the South Maharashtra Shear Zone below the Deccan Traps—or its lack thereof—is critical for validating tectonic models correlating Madagascar with India. In this study, zircons in tonalitic basement xenoliths hosted in Deccan Trap dykes were dated in situ, using the U-Pb isotope system. The data furnish U-Pb ages that define three populations at 2527 ± 6, 2456 ± 6, and 2379 ± 9 Ma. The 2527 ± 6 Ma ages correspond to the igneous crystallization of the tonalites, whereas the 2456 ± 6 and 2379 ± 9 Ma ages date metamorphic overprints. The results help to establish for the first time that the basement is a part of the Neoarchean granitoid suite of the Western Dharwar Craton, which extends northward up to at least Talvade in central and Kihim beach in the western Deccan. By implication, the South Maharashtra Shear Zone cannot be the northern limit of the Western Dharwar Craton. The granitoids are correlated with the Neoarchean felsic intrusions (2.57–2.49) of the Masaola suite in the Antongil Block of Madagascar, supporting the existence of a Neoarchean Greater Dharwar Craton comprising the Western Dharwar Craton and the Antongil-Masora Block.

Gas composition of the main volatile elements in protoplanetary discs and its implication for planet formation

Context. Direct observations of gaseous exoplanets reveal that their gas envelope has a higher C/O ratio than that of the host star (e.g., Wasp 12-b). This has been explained by considering that the gas phase of the disc could be inhomogeneous, exceeding the stellar C/O ratio in regions where these planets formed; but few studies have considered the drift of the gas and planet migration. Aims. We aim to derive the gas composition in planets through planet formation to evaluate if the formation of giant planets with an enriched C/O ratio is possible. The study focusses on the effects of different processes on the C/O ratio, such as the disc evolution, the drift of gas, and planet migration. Methods. We used our previous models for computing the chemical composition, together with a planet formation model, to which we added the composition and drift of the gas phase of the disc, which is composed of the main volatile species H2O, CO, CO2, NH3, N2, CH3OH, CH4, and H2S, H2 and He. The study focusses on the region where ice lines are present and influence the C/O ratio of the planets. Results. Modelling shows that the condensation of volatile species as a function of radial distance allows for C/O enrichment in specific parts of the protoplanetary disc of up to four times the solar value. This leads to the formation of planets that can be enriched in C/O in their envelope up to three times the solar value. Planet migration, gas phase evolution and disc irradiation enables the evolution of the initial C/O ratio that decreases in the outer part of the disc and increases in the inner part of the disc. The total C/O ratio of the planets is governed by the contribution of ices accreted, suggesting that high C/O ratios measured in planetary atmospheres are indicative of a lack of exchange of material between the core of a planet and its envelope or an observational bias. It also suggests that the observed C/O ratio is not representative of the total C/O ratio of the planet.

Daily and seasonal thermal stresses in tilings: a field survey combined with numeric modeling

This study investigates thermally induced tensile stresses in ceramic tilings. Daily and seasonal thermal cycles, as well as, rare but extreme events, such as a hail-storm striking a heated terrace tiling, were studied in the field and by numerical modeling investigations. The field surveys delivered temperature– time diagrams and temperature profiles across tiling systems. These data were taken as input parameters for modeling the stress distribution in the tiling system in order to detect potential sites for material failure. Dependent on the thermal scenario (e.g., slow heating of the entire structure during morning and afternoon, or a rapid cooling of the tiles by a rain storm) the modeling indicates specific locations with high tensile stresses. Typically regions along the rim of the tiling field showed stresses, which can become critical with respect to the adhesion strength. Over the years, ongoing cycles of thermal expansion–contraction result in material fatigue promoting the propagation of cracks. However, the installation of flexible waterproofing membranes (applied between substrate and tile adhesive) represents an efficient technical innovation to reduce such crack propagation as confirmed by both numerical modeling results and microstructural studies on real systems.

Antarctic Zone nutrient conditions during the last two glacial cycles

In a sediment core from the Pacific sector of the Antarctic Zone (AZ) of the Southern Ocean, we report diatom-bound N isotope (δ15Ndb) records for total recoverable diatoms and two distinct diatom assemblages (pennate and centric rich). These data indicate tight coupling between the degree of nitrate consumption and Antarctic climate across the last two glacial cycles, with δ15Ndb (and thus the degree of nitrate consumption) increasing at each major Antarctic cooling event. Coupled with evidence from opal- and barium-based proxies for reduced export production during ice ages, the δ15Ndb increases point to ice age reductions in the supply of deep ocean-sourced nitrate to the AZ surface. The two diatom assemblages and species abundance data indicate that the δ15Ndb changes are not the result of changing species composition. The pennate and centric assemblage δ15Ndb records indicate similar changes but with a significant decline in their difference during peak ice ages. A tentative seasonality-based interpretation of the centric-to-pennate δ15Ndb difference suggests that late summer surface waters became nitrate free during the peak glacials.