Deep hydrothermal circulation in a granite intrusion beneath Larderello geothermal area (Italy): constraints from mineralogy, fluid inclusions and stable isotopes

https://doi.org/10.1016/S0377-0273(03)00150-1Get rights and content

Abstract

A shallow (∼2 km) granite body has been intersected by the MONT-4 well in new producing areas of the Larderello geothermal field. The granite at 2.2–2.5 km and below 3 km underwent pervasive propylitic alteration that is dominated by the assemblage chlorite+adularia+quartz+illite+calcite+albite±epidote±sphene. Fluid inclusion studies indicate two stages of alteration. The first one is characterised by high-temperature fluids that, at 2218 m, were trapped under boiling conditions at about 340°C and 15 MPa, whereas at 3520 m evidence of boiling is not revealed and the temperature of the fluid was calculated to be about 350°C. The second alteration stage shows a cooling of the system, particularly evident in the shallower portion of the intrusion (∼235°C at 2218 m and ∼315°C at 3520 m). The temperature of homogenisation for subsequent cooler inclusions approximates present-day down-hole temperatures, and may be considered as representative of the fluid actually circulating in the granite in recent times. Calculated δ18O values of 0±1‰, for fluids in equilibrium at 340°C (2218 m) and 350°C (3520 m) with feldspar, epidote and chlorite, and δD values of −51 to −61‰, for fluids in equilibrium with chlorite, together with directly measured δD values of −33 to −61‰ for fluid inclusions in quartz, indicate derivation from meteoric waters. The original meteoric signature of the fluids has been altered through water–rock interaction, but remains clearly distinguished from a primary magmatic signature, for which no isotopic evidence was found. Active fracturing and the high structural relief of granite emplacement promoted hydrothermal circulation in the granitic body drilled by MONT-4 well, suggesting that such characteristics may be a pre-requisite for such circulation at Larderello.

Introduction

The crust in Southern Tuscany (Italy) is characterised by high heat flow due to the Neogene magmatism of the Tuscan Magmatic Province and crustal thinning (Marinelli, 1963, Marinelli et al., 1993, Serri et al., 1993, Baldi et al., 1995, Mongelli et al., 1998). Since the late Miocene, the development of widespread hydrothermal activity connected with this magmatism is indicated by the presence of hydrothermal ore deposits and active geothermal systems (Marinelli, 1963, Marinelli, 1969, Marinelli, 1983, Cataldi et al., 1978, Tanelli, 1983, Gianelli et al., 1988, Gianelli et al., 1997, Lattanzi, 1999). The vapour-dominated Larderello geothermal field is the most striking feature of this active hydrothermal activity, and has been extensively exploited for several decades for electrical energy production by ENEL (Italian Electricity Board). Two main productive horizons are exploited: (1) a shallow reservoir, located in Mesozoic carbonate and anhydritic units, and (2) a deep reservoir in permeable levels of the Palaeozoic basement which has been metamorphosed regionally and by contact with the Neogene intrusions. Geological and geophysical data indicate the presence of these shallow intrusive bodies beneath the Larderello geothermal field, belonging to a composite batholith distributed at different depths within the crust (Gianelli et al., 1997). Radiometric ages of 3.8–2.5 Ma were measured in granites and granite dykes from deep Larderello geothermal wells (Batini et al., 1985, Villa et al., 1997). As a consequence, these shallow intrusions cannot be considered responsible for the actual geothermal anomaly. However, geophysical data support the presence of a still partially molten magmatic body below Larderello (Gianelli et al., 1997). This is considered to be responsible for the very high positive thermal anomaly of the area (geothermal gradient up to >300°C·km−1; Baldi et al., 1995). Aplitic dykes and two-mica granites have been penetrated by deep (3–4.5 km) geothermal wells, and associated high-grade contact metamorphic aureoles are common (Gianelli et al., 1997, Carella et al., 2000). Several studies revealed a complex evolution of the fluids circulated in the subterrane of Larderello (Cavarretta et al., 1982, Belkin et al., 1985, Valori et al., 1992, Cathelineau et al., 1994, Ruggieri et al., 1999, Ruggieri and Gianelli, 1999), from an early stage, dominated by magmatic and thermometamorphic fluids associated with the emplacement of granitic intrusions, to the present-day fluids, prevalently of meteoric origin.

In the last few years, research into the exploitation of high enthalpy fluids at Larderello has moved the exploration toward the southeastern part of the field (Barelli et al., 2000). In the Montieri zone (Fig. 1) a deep well (MONT-4) encountered a granitic body at relatively shallow depth (∼2000 m below ground level) and traversed through to the depth of 3721 m. This granite shows clear evidence of hydrothermal alteration throughout the well. Our work is focussed on the characterisation of the fluids involved in the hydrothermal alteration of this granite, through mineralogical, fluid inclusion and stable isotope (δ18O and δD) investigations. We show that the granitic intrusion was first affected by the circulation of meteoric-sourced fluids at around 340°C, and a more recent event – perhaps still on-going – associated with meteoric-sourced fluids at temperatures ranging between 235°C (2218 m) and 315°C (3520 m).

Section snippets

Geological setting and drill-hole lithostratigraphy

The MONT-4 well is located on 100 m of shale belonging to the Ligurian flysch complex. The borehole then passes through a 1000-m sequence of tectonic slices, formed by Triassic dolomitic limestone and anhydrite of the Tuscan Nappe, which is underlain by 900 m of Palaeozoic metamorphic basement consisting of chlorite–graphite phyllites. Thermometamorphic biotite begins to nucleate at about 1300 m, whilst hornfels rocks form below 1900 m (Barelli et al., 2000). The granite intrusion is cut from a

Analytical methods

Minerals were analysed by SEM–EDS using a Philips XL30 micro-analytical system coupled with EDAX DX4 at 20-kV accelerating energy, 0.1-nA beam current, 100-s live time (Earth Science Department, University of Pisa). Before each analytical session, calibration and standardisation using international mineral standards were performed. Analyses are normalised to 100 due to EDAX software used.

Microthermometries were performed on double polished thin sections, using a Linkam THMS 600 heating/freezing

Hydrothermal mineral distribution

The distribution of hydrothermal minerals is reported in Fig. 3. Microanalytical analyses on hydrothermal minerals were carried out at 2218 m (core), 2850 m (cuttings) and 3520 m (cuttings). Pervasive hydrothermal alteration was observed in the intervals 2200–2500 m and below 3000 m, whereas between 2500 and 3000 m the rocks are less altered. The secondary paragenesis is dominated by the assemblage chlorite+adularia+quartz+illite+calcite+albite±epidote±sphene (Fig. 3; Table 1, Table 2, Table 3

Fluid inclusions

Fluid inclusions in hydrothermal and magmatic quartz were examined at two different depths (2218 and 3520 m). The size of the analysed fluid inclusions are in the range 5–30 μm. Microscopic observation at room temperature reveals three main types of fluid inclusions.

Type 1: two phase liquid-rich fluid inclusions (L+V). Two sub-populations can be distinguished on the basis of optical observation: Types 1a and 1b. Type 1a is the most widespread population; the fluid inclusions are primary in

Stable isotopes

To better constrain the origin of the fluids involved in the hydrothermal alteration of the granite, primary and hydrothermal silicate minerals were analysed for δ18O and δD, and δD was analysed also in decrepitated Type 1a fluid inclusions hosted in quartz (Table 7). The samples come from three depths (2218, 2740, 3520 m) of the MONT-4 well. The high δ18O values of magmatic quartz (between 10.4‰ and 12.7‰; Table 7) are compatible with the 18O-enriched magmas of the Tuscan Province (+11 to +15;

Discussion

The granitic intrusion in the MONT-4 well underwent propylitic alteration at 240–350°C. At 2218 m the coexistence of liquid-rich (Type 1a) and vapour-rich (Type 2) fluid inclusions and their similar range of Th (Table 5) suggest that boiling may have occurred during deposition of hydrothermal alteration mineral assemblage. The mode of Th of Type 2 fluid inclusions is slightly, but significantly higher than that of the coexisting Type 1a fluid inclusions. This evidence, coupled with the relative

Conclusions

This study shows that the shallow-emplaced granite, cut by the MONT-4 geothermal well in the southeastern part of the Larderello geothermal field, is pervasively propylitically altered by hydrothermal fluids. Two main stages were observed and are clearly distinguishable in the shallower part of the intrusion. The first stage is characterised by a temperature of about 340–350°C, with clear evidence for boiling at 2218 m. The second stage reveals a cooling of the system, particularly evident in

Acknowledgements

We thank ENEL Green Power SpA (formerly ENEL) for access to the samples, for encouraging this work, and for permission to publish the data. The authors are grateful to Robert Criss and an anonymous reviewer for their comments, suggestions and criticisms that helped to improve the quality of the manuscript. This research was supported by a research grant to A.S. in the framework of the project MURST 2000 (Relationship between continental crust structure and geothermal resources in southern

References (62)

  • P Baldi et al.

    Geothermal anomalies and structural features of southern Tuscany (Italy)

    Proc. World Geotherm. Congr., Florence

    (1995)
  • Barelli, A., Bertini, G., Buonasorte, G., Cappetti, G., Fiordelisi, A., 2000. Recent deep exploration results at the...
  • P.B Barton et al.

    Refinement of the evaluation of the role of CO2 in modifying estimates of pressure of epithermal mineralization

    Econ. Geol.

    (1993)
  • F Batini et al.

    Deep structure, age and evolution of the Larderello–Travale geothermal field

    Geotherm. Res. Counc. Trans.

    (1985)
  • H.E Belkin et al.

    Fluid inclusions in minerals from the geothermal fields of Tuscany, Italy

    Geothermics

    (1985)
  • D.K Bird et al.

    Calc–silicate mineralization in active geothermal systems

    Econ. Geol.

    (1984)
  • B.P Bishop et al.

    Variation in sericite compositions from fracture zones within the Coso Hot Springs geothermal system

    Geochim. Cosmochim. Acta

    (1987)
  • R.J Bodnar

    Revised equationand table for determining the freezing point depression of H2O–NaCl solutions

    Geochim. Cosmochim. Acta

    (1993)
  • Bodnar, R.J., Vityk, M.O., 1994. Interpretation of microthermometric data for H2O–NaCl fluid inclusions. In: De Vivo,...
  • Bodnar, R.J., Reynolds, T.J., Kuehn, C.A., 1985. Fluid-inclusion systematics in epithermal systems. In: Berger, B.R.,...
  • P.R.L Browne et al.

    The Ohaki-Broadlands hydrothermal area, New Zealand: Mineralogy and related geochemistry

    Am. J. Sci.

    (1970)
  • M Carella et al.

    Metamorphic consequences of Neogene thermal anomaly in the northern Apennines (Radicondoli–Travale area, Larderello geothermal field – Italy)

    Geodin. Acta

    (2000)
  • R Cataldi et al.

    Assessment of geothermal potential of central and Southern Tuscany

    Geothermics

    (1978)
  • M Cathelineau

    Cation site occupancy in chlorites and illites as a function of temperature

    Clays Clay Miner.

    (1988)
  • M Cathelineau et al.

    A chlorite solid solution geothermometer: The Los Azufres (Mexico) geothermal system

    Contrib. Miner. Petrol.

    (1985)
  • M Cathelineau et al.

    Evidence for Li-rich brines and early magmatic fluid–rock interaction in the Larderello geothermal system

    Geochim. Cosmochim. Acta

    (1994)
  • G Cavarretta et al.

    Formation of authigenic minerals and their use as indicators of physico–chemical parameters of the fluid in the Larderello–Travale geothermal field

    Econ. Geol.

    (1982)
  • D.R Cole et al.

    Oxygen isotope fractionation between chlorite and water from 170 to 350°C: A preliminary assessment based on partial exchange and fluid/rock experiments

    Geochim. Cosmochim. Acta

    (1998)
  • Crawford, M.L., 1981. Phase equilibria in aqueous fluid inclusions. In: Hollister, L.S., Crawford, M.L. (Eds.), Short...
  • R.E Criss et al.

    An 18O/16O and D/H study of Tertiary hydrothermal systems in the southern half of the Idaho batholith

    Geol. Soc. Am. Bull.

    (1983)
  • F D’Amore et al.

    Isotopic evidence for a magmatic contribution to fluids of the geothermal systems of Larderello, Italy, and the Geyser, California

    Geothermics

    (1994)
  • A Dallegno et al.

    Pyrite deposits of the Gavorrano area, Grosseto

    Atti Soc. Tosc. Sci. Nat.

    (1979)
  • P De Caritat et al.

    Chlorite geothermometry: A review

    Clays Clay Miner.

    (1993)
  • T Donnelly et al.

    Hydrogen isotope analysis of natural abundance and deuterium-enriched waters by reduction over chromium on-line to a dynamic dual inlet isotope-ratio mass spectrometer

    Rapid Commun. Mass Spectrom.

    (2001)
  • Field, C.W., Fifarek, R.H., 1985. Light stable-isotope systematics in the epithermal environment. In: Berger, B.R.,...
  • P Fulignati et al.

    Fluid evolution and minerogenetic implications in the Botro ai Marmi monzogranitic intrusion (Campiglia Marittima, Italy)

    Terra Nostra

    (1999)
  • G Gianelli et al.

    Structural setting of the Larderello–Travale geothermal region

    Mem. Soc. Geol. Ital.

    (1978)
  • G Gianelli et al.

    Geological model of a young volcano–plutonic system: The geothermal region of Monte Amiata (Tuscany, Italy)

    Geothermics

    (1988)
  • G Gianelli et al.

    Crustal models of the geothermal areas of southern Tuscany (Italy)

    Tectonophysics

    (1997)
  • C.M Graham et al.

    Experimental study of hydrogen-isotope exchange between aluminous chlorite and water and of hydrogen diffusion in chlorite

    Am. Miner.

    (1987)
  • F.M Haynes et al.

    Synthetic fluid inclusions in natural quartz. IV. Chemical analyses of fluid inclusions by SEM/EDA. Evaluation of method

    Geochim. Cosmochim. Acta

    (1988)
  • Cited by (32)

    • Episodic hydrothermal alteration recorded by microscale oxygen isotope analysis of white mica in the Larderello-Travale Geothermal Field, Italy

      2020, Chemical Geology
      Citation Excerpt :

      Numerous mineralogical, stable isotope (δ18O, δD) and fluid inclusion studies have established and improved the current reservoir model (Batini et al., 1985, D'Amore et al., 1987; Cavarretta and Puxeddu, 1990; Valori et al., 1992; Cathelineau et al., 1994; Petrucci et al., 1993, 1994; Gianelli et al., 1997; Ruggieri and Gianelli, 1999; Ruggieri et al., 1999; Carella et al., 2000; Gianelli and Ruggieri, 2002; Boyce et al., 2003; Dallai et al., 2005). In general, two major stages of hydrothermal fluid alteration are distinguished: (1) an early high-temperature (HT) phase, coeval with the emplacement of the granitic intrusions and the contact metamorphic overprint, with fluids of mostly local magmatic provenance (δ18O = 11–15‰), and (2) a late to present day low-temperature (LT) hydrothermal stage, dominated by mixing between infiltrating meteoric fluids and deep high-δ18O fluids which could be residual magmatic fluids and/or a meteoric fluids that undergone oxygen isotope exchange with the host rocks (Petrucci et al., 1993, 1994; D'Amore and Bolognesi, 1994; Gianelli and Ruggieri, 2002; Boyce et al., 2003; Dallai et al., 2005). Stable isotope data (δ18O, δD, δ13C) from the LTGF are available for a wide range of lithologies, minerals (primary and authigenic), fluid inclusions, dry steam used for energy production, and fluids of variable sources (condensate, reinjection, spring water) (Celati et al., 1973; Panichi et al., 1974; D'Amore et al., 1987; Petrucci et al., 1993, 1994; D'Amore and Bolognesi, 1994; Ruggieri and Gianelli, 1999; Gianelli and Ruggieri, 2002; Boyce et al., 2003; Dallai et al., 2005; Farina et al., 2018).

    • Reflection seismic sections across the Geothermal Province of Tuscany from reprocessing CROP profiles

      2015, Geothermics
      Citation Excerpt :

      It may correspond to the horizon H shown in the profile CROP-18A (Fig. 3) and in the interpretation of 2D (Fiordelisi et al., 2005) and 3D profiles (Casini et al., 2010b) in the geothermal area of Travale, that is not far off the north-west termination of CROP-18B. The horizon may correspond to a thermo-metamorphic level (skarn and hornfels layers) drilled at around 1300 m by the well Montieri-4 (Barelli et al., 2000; Boyce et al., 2003). The presence of tectonic slices can be easily inferred over this marker.

    View all citing articles on Scopus
    View full text