Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

The history of subaquatic volcanism recorded in the sediments of Lake Kivu

Subaquatic volcanic activity has been ongoing in Lake Kivu since the early Holocene and has a dynamic effect on the biological productivity in the surface water, and the preservation of carbonate in the deep anoxic water. Groundwater discharge into the lake’s deepwater propels the upward advection of the water column that ultimately supplies nutrients to the surface water for biological production. The amount of nutrients supplied from the deepwater can be increased suddenly by (1) a cold meteorological event that drives deep seasonal mixing resulting in increased nutrients from below and oxygen from above, or (2) subaquatic volcanic activity that induces a buoyant hydrothermal plume, which entrains nutrients from the deepwater and results in anoxia or suboxic conditions in the surface water. Previous sedimentological studies in Lake Kivu have hypothesized that regional climatic changes are responsible for sudden changes in the preservation of carbonates in the Main Basin. Here we reveal that sublacustrine volcanic events most likely induce the abrupt changes to the geochemistry in the sediment in Lake Kivu. An unprecedented look into the sediment stratigraphy and geochemistry from high-resolution seismic-reflection, and 15N-isotope analyses was conducted in the Main Basin. The results reveal that buoyant hydrothermal plumes caused by subaquatic volcanic activity are a possible trigger for increased biological productivity and organic matter preservation, and that ongoing hydrothermal activity increases the alkalinity in the deepwater, leading to carbonate preservation. The onset of carbonate preservation since the 1970s that is currently observed in the sediment could indicate that hydrothermal discharge has recently increased in the lake.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

Diffusion of HTO, Br-, I-, Cs+, Sr-85(2+) and Co-60(2+) in a clay formation: Results and modelling from an in situ experiment in Opalinus Clay

The migration of radioactive and chemical contaminants in clay materials and argillaceous host rocks is characterised by diffusion and retention processes. Valuable information on such processes can be gained by combining diffusion studies at laboratory scale with field migration tests. In this work, the outcome of a multi-tracer in situ migration test performed in the Opalinus Clay formation in the Mont Terri underground rock laboratory (Switzerland) is presented. Thus, 1.16 x 10(5) Bq/L of HTO, 3.96 x 10(3) Bq/L of Sr-85, 6.29 x 10(2) Bq/L of Co-60, 2.01 x 10(-3) mol/L Cs, 9.10 x 10(-4) mol/L I and 1.04 x 10(-3) mol/L Br were injected into the borehole. The decrease of the radioisotope concentrations in the borehole was monitored using in situ gamma-spectrometry. The other tracers were analyzed with state-of-the-art laboratory procedures after sampling of small water aliquots from the reservoir. The diffusion experiment was carried out over a period of one year after which the interval section was overcored and analyzed. Based on the experimental data from the tracer evolution in the borehole and the tracer profiles in the rock, the diffusion of tracers was modelled with the numerical code CRUNCH. The results obtained for HTO (H-3), I- and Br- confirm previous lab and in situ diffusion data. Anionic fluxes into the formation were smaller compared to HTO because of anion exclusion effects. The migration of the cations Sr-85(2+), Cs+ and Co-60(2+) was found to be governed by both diffusion and sorption processes. For Sr-85(2+), the slightly higher diffusivity relative to HTO and the low sorption value are consistent with laboratory diffusion measurements on small-scale samples. In the case of Cs+, the numerically deduced high diffusivity and the Freundlich-type sorption behaviour is also supported by ongoing laboratory data. For Co, no laboratory diffusion data were yet available for comparison; however, the modelled data suggests that Co-60(2+) sorption was weaker than would be expected from available batch sorption data. Overall, the results demonstrate the feasibility of the experimental setup for obtaining high-quality diffusion data for conservative and sorbing tracers. (C) 2007 Elsevier Ltd. All rights reserved.

A combined vacuum crushing and sieving (CVCS) system designed to determine noble gas paleotemperatures from stalagmite samples

This paper presents a novel extraction device for water and noble gases from speleothem samples for noble gas paleotemperature determination. The “combined vacuum crushing and sieving (CVCS) system” was designed to reduce the atmospheric noble gas contents from air inclusions in speleothem samples by up to 2 orders of magnitude without adsorbing atmospheric noble gases onto the freshly produced grain surfaces, a process that had often hampered noble gas temperature (NGT) determination in the past. We also present the results from first performance tests of the CVCS system processing stalagmite samples grown at a known temperature. This temperature is reliably reproduced by the NGTs derived from Ar, Kr, and Xe extracted from the samples. The CVCS system is, therefore, suitable for routine determinations of accurate NGTs. In combination with stalagmite dating, these NGTs will allow reconstructing past regional temperature evolutions, and also support the interpretation of the often complex stable isotope records preserved in the stalagmites' calcite.

Taxon-specific δ13C analysis of chitinous invertebrate remains in sediments from Strandsjön, Sweden

Taxon-specific stable carbon isotope (δ13C) analysis of chitinous remains of invertebrates can provide valuable information about the carbon sources used by invertebrates living in specific habitats of lake ecosystems (for example, sediments, water column, or aquatic vegetation). This is complementary to δ13C of sedimentary organic matter (SOM), which provides an integrated signal of organic matter produced in a lake and its catchment, and of diagenetic processes within sediments. In a sediment record from Strandsjön (Sweden) covering the past circa 140 years, we analyzed SOM geochemistry (δ13C, C:Natomic, organic carbon content) and δ13C of chitinous invertebrate remains in order to examine whether taxon-specific δ13C records could be developed for different invertebrate groups and whether these analyses provide insights into past changes of organic carbon sources for lacustrine invertebrates available in benthic and planktonic compartments of the lake. Invertebrate taxa included benthic chironomids (Chironomus, Chironomini excluding Chironomus, Tanytarsini, and Tanypodinae), filter-feeders on suspended particulate organic matter (Daphnia, Plumatella and Cristatella mucedo), and Rhabdocoela. δ13C of chironomid remains indicated periodic availability of 13C-depleted carbon sources in the benthic environment of the lake as δ13C values of the different chironomid taxa fluctuated simultaneously between -34.7 and -30.5‰ (VPDB). Daphnia and Bryozoa showed parallel changes in their δ13C values which did not coincide with variations in δ13C of chironomids, though, and a 2-3‰ decrease since circa AD 1960. The decrease in δ13C of Daphnia and Bryozoa could indicate a decrease in phytoplankton δ13C as a result of lower lake productivity, which is in accordance with historical information about the lake that suggests a shift to less eutrophic conditions after AD 1960. In contrast, Rhabdocoela cocoons were characterized by relatively high δ13C values (-30.4 to -28.2‰) that did not show a strong temporal trend, which could be related to the predatory feeding mode and wide prey spectrum of this organism group. The taxon-specific δ13C analyses of invertebrate remains indicated that different carbon sources were available for the benthic chironomid larvae than for the filter-feeding Daphnia and bryozoans. Our results therefore demonstrate that taxon-specific analysis of δ13C of organic invertebrate remains can provide complementary information to measurements on bulk SOM and that δ13C of invertebrate remains may allow the reconstruction of past changes in carbon sources and their δ13C in different habitats of lake ecosystems.

The Yanque prospect (Peru): from polymetallic Zn-Pb mineralization to a nonsulfide deposit

The Yanque nonsulfide Pb-Zn deposit (inferred resources 12.5 Mt @ 3.7% Pb and @ 3.5% Zn) is located in the Andahuaylas-Yauri ore province (Cuzco, southern Peru). The deposit occurs within a base metal mineralized district, centered on the medium-sized Dolores porphyry copper. A thorough geological, mineralogical and geochemical study has carried out in order to define: the relationships between the Dolores Cu-porphyry ore and the Yanque Zn-Pb polymetallic mineralization, and the characteristics of the economic nonsulfide concentrations. Both sedimentary and igneous rocks constitute the backbone of the Yanque-Dolores area. The sedimentary lithologies belong to the Soraya, Mara and Ferrobamba Fms. (upper Jurassic-middle Cretaceous). The Yanque orebody is hosted by the Mara Fm., which prevailingly consists of a siliciclastic sedimentary breccia. The original sulfide mineralization consisted of galena, pyrite and sphalerite. The host rock has been affected by a strong hydrothermal alteration, characterized by prevailing sericite/illite, as in the typical porphyry-related phyllic-argillic alteration stage, and by minor kaolinite, dolomite and quartz. Minor element geochemistry, characterized by Sb, As, Mn, Ag and locally also by Cu, points to magmatic-hydrothermal related mineralizing fluids. The Pb isotopic compositions from Dolores and Yanque sulfides are similar, and are typical of the Tertiary magmatically-derived ores in this part of Peru. The hydrothermally altered rocks at Yanque have the same Pb isotopic compositions as the sulfides, thus confirming the hypothesis that the Yanque primary Zn-Pb mineralization may have been produced by hydrothermal circulation related to the emplacement of the Dolores Cu-porphyry, as it is the case of other porphyry Cu systems associated with polymetallic mineralization elsewhere. However, no simple genetic model for the mineralization involving just one fluid circulation episode is able to explain the data. The Yanque economic nonsulfide ore association consists of sauconite, hemimorphite, smithsonite and cerussite, which result from the weathering and alteration of the original sulfide mineralization. Zinc is allocated mainly in sauconite (Zn-smectite), rather than in carbonates: a factor strictly related to the prevailing siliciclastic character of the host rock. Distinctive features of the Yanque orebody are the comparable ore grades for both Pb and Zn (3.5% Zn and 3.7% Pb), and the inverse supergene chemical zoning. In fact, contrary to other supergene ores of this type, zinc prevails in the top zone of the Yanque deposit, whereas lead content increases with depth. Considering the different mobility of the two metals in solution, it may be assumed that most of the primary zinc that was the source for the Yanque nonsulfides was originally located far from the position occupied by the galena mineralization, whose remnants have been observed on site. Zinc sulfides may have been originally contained in the eroded rock volumes that surrounded the actual deposit: the zinc-rich solutions have possibly migrated through the siliciclastic Mara Fm. and precipitated the nonsulfide minerals by porosity filling and replacement processes. In this sense, the Yanque secondary Zn-Pb deposit could be considered as a special type of “Exotic” mineralization.

Provenance versus weathering control on the composition of tropical river mud (southern Africa)

This study presents an integrated mineralogical-geochemical data base on fine-grained sediments transported by all major rivers of southern Africa, including the Zambezi, Okavango, Limpopo, Olifants, Orange and Kunene. Clay mineralogy, bulk geochemistry, Sr and Nd isotopic signatures of river mud, considered as proxy of suspended load, are used to investigate the influence of source-rock lithology and weathering intensity on the composition of clay and silt produced in subequatorial to subtropical latitudes. Depletion in mobile alkali and alkaline-earth metals, minor in arid Namibia, is strong in the Okavango, Kwando and Upper Zambezi catchments, where recycling is also extensive. Element removal is most significant for Na, and to a lesser extent for Sr. Depletion in K, Ca and other elements, negligible in Namibia, is moderate elsewhere. The most widespread clay minerals are smectite, dominant in muds derived from Karoo or Etendeka flood basalts, or illite and chlorite, dominant in muds derived from metasedimentary rocks of the Damara Orogen or Zimbabwe Craton. Kaolinite represents 30-40% of clay minerals only in Okavango and Upper Zambezi sediments sourced in humid subequatorial Angola and Zambia. After subtracting the effects of recycling and of local accumulation of authigenic carbonates in soils, the regional distribution of clay minerals and chemical indices consistently reflect weathering intensity primarily controlled by climate. Bulk geochemistry identifies most clearly volcaniclastic sediments and mafic sources in general, but cannot discriminate the other sources of detritus in detail. Instead, Sr and Nd isotopic fingerprints are insensitive to weathering, and thus mirror faithfully the tectonic structure of the southern African continent. Isotopic tools thus represent a much firmer basis than bulk geochemistry or clay mineralogy in the provenance study of mudrocks.

Geochemical variability in distal and proximal glass from the Youngest Toba Tuff eruption

The Youngest Toba Tuff (YTT, erupted ca. 74 ka ago) is a distinctive and widespread tephra marker across south and southeast Asia. The climatic, human and environmental consequences of the YTT eruption are widely debated. Although a considerable body of geochemical data is available for this unit, there has not been a systematic study of the variability of the ash geochemistry. Intrinsic (magmatic) and extrinsic (post-depositional) chemical variations bring fundamental information regarding the petrogenesis of the magma, the distribution of the tephra and the interaction between the ash and the receiving environment. Considering the importance of the geochemistry of the YTT for stratigraphic correlations and eruptive models, it is central to the YTT debate to quantify and interpret such variations. Here we collate all published geochemical data on the YTT glass, including analyses from 68 sites described in the literature and three new samples. Two principal sources of chemical variation are investigated: (i) compositional zonation of the magma reservoir, and (ii) post-depositional alteration. Post-depositional leaching is responsible for up to ca. 11% differences in Na2O/K2O and ca. 1% differences in SiO2/Al2O3 ratios in YTT glass from marine sites. Continental tephra are 2% higher in Na2O/K2O and 3% higher in SiO2/Al2O3 respect to the marine tephra. We interpret such post-depositional glass alteration as related to seawater induced alkali migration in marine environments, or to site-specific water pH. Crystal fractionation and consequential magmatic differentiation, which produced order-of-magnitude variations in trace element concentrations reported in the literature, also produced major element differences in the YTT glass. FeO/Al2O3 ratios vary by about 50 %, which is analytically significant. These variations represent magmatic fractionation involving Fe-bearing phases. We also compared major element concentrations in YTT and Oldest Toba Tuff (OTT) ash samples, to identify potential compositional differences that could constrain the stratigraphic identity of the Morgaon ash (Western India); no differences between the OTT and YTT samples were observed.

Anisotropy of magnetic susceptibility in natural olivine single crystals

Mantle flow dynamics can cause preferential alignment of olivine crystals that results in anisotropy of physical properties. To interpret anisotropy in mantle rocks, it is necessary to understand the anisotropy of olivine single crystals. We determined anisotropy of magnetic susceptibility (AMS) for natural olivine crystals. High-field AMS allows for the isolation of the anisotropy due to olivine alone. The orientations of the principal susceptibility axes are related to the olivine’s crystallographic structure as soon as it contains >3 wt % FeO. The maximum susceptibility is parallel to the c axis both at room temperature (RT) and at 77 K. The orientation of the minimum axis at RT depends on iron content; it is generally parallel to the a axis in crystals with 3–5 wt % FeO, and along b in samples with 6–10 wt % FeO. The AMS ellipsoid is prolate and the standard deviatoric susceptibility, k0, is on the order of 8*10210 m3/kg for the samples with <1wt % FeO, and ranges from 3.1*1029 m3/kg to 5.7*1029 m3/kg for samples with 3–10 wt % FeO. At 77 K, the minimum susceptibility is along b, independent of iron content. The shape of the AMS ellipsoid is prolate for samples with <5 wt % FeO, but can be prolate or oblate for higher iron content. The degree of anisotropy increases at 77 K with p0 7757.160.5. The results from this study will allow AMS fabrics to be used as a proxy for olivine texture in ultramafic rocks with high olivine content.