Geoelectrical characterization of Socompa lagoon area in the Andean Central Volcanic Zone from 3-D audiomagnetotelluric inversion

https://doi.org/10.1016/j.jvolgeores.2021.107246Get rights and content

Highlights

  • The Socompa Geothermal Prospect is located in the Andean Central Volcanic Zone.

  • An audiomagnetotelluric survey was carried out in the Socompa Geothermal Prospect.

  • A high conductive layer at 400–500 m and at least 200 m thick was determined by the 3D model.

  • The main scenarios for the high conductivity layer are: a saline aquifer; a Mio-Piocene? clay cap.

Abstract

The Socompa stratovolcano is located in the Andean Central Volcanic Zone, where active volcanism turns it into an area of interest for geothermal resources. In the surroundings of the volcano, the Socompa Geothermal Prospect was carried out looking for thermal manifestations related to an active geothermal system. The results of these preliminary studies are not conclusive.

Here we present a geoelectrical model from a 3-D inversion of audiomagnetotelluric data acquired in an area of 16 km2. It covers the topographic low where the Socompa lagoon and the thermal springs are located, in order to determine the reservoir of circulating geothermal waters. The most distinctive feature determined by the 3-D model is a high conductive layer (less than 10 Ωm) at a depth of about 400–500 m and at least 200 m thick. Due to the environment and the depth to which the conductive layer is found, there is more than one possibility to explain its presence. It could be due to a deep saline aquifer that would be hosting the circulating waters. But since the low topography could coincide with a cryptic caldera of Miocene - Pliocene(?) age, the high conductivity layer could also be explained by the presence of a clay cap developed during the period of activity. From the inversion of the audiomagnetotelluric data it was possible to carry out a geoelectrical characterization of the upper units present in the study area, as well as to delimit, in a first approximation, the depth to which the hydrothermal reservoir would be found.

Introduction

The Socompa is one of the largest stratovolcanoes of the Central Volcanic Zone (CVZ) of the Andes (see Fig. 1). The region is characterized by an active volcanic arc related to the subduction of the Nazca plate below the South American plate, in a high plateau environment, from semi-arid to hyper-arid weather conditions (Stern, 2004). Thermo-mechanical studies of the Andes state indicate that the region where the volcano is located shows the highest heat-flow and the shallowest depths of the Curie isotherm (Ibarra and Prezzi, 2019). The active volcanism makes the region attractive for the exploration of geothermal systems (e.g., Bissig et al., 2001; Tassi et al., 2009; McCoy-West et al., 2011; Lahsen et al., 2015; Munoz-Saez et al., 2018; Chiodi et al., 2019; Filipovich et al., 2020).

Previous studies in the Socompa area suggest favorable preliminary geological and hydrogeological features for a blind geothermal system characterized by heat source, related to the current magmatic system that feeds the Socompa volcano, reservoirs and seals within a possible calderic structure system, and water availability from different feeding zones and flow paths (Galliski et al., 1987; Seggiaro and Apaza, 2018; Lelli, 2018). The most recent studies were carried out within the framework of the Socompa Goethermal Prospect (SGP) which is located in the topographic low with mean altitude of 3500–3800 m.a.s.l. located between the southeastern flank of the volcano and the western boundary of the Caipe range, encompassing the Socompa lagoon. The main manifestations are the thermal springs in the Socompa lagoon and in Quebrada del Agua locality and pliocene rocks that are altered by the presence of circulating fluids on the Quebrada del Agua fault (see Fig. 1a) (Seggiaro and Apaza, 2018).

Electromagnetic methods are widely used for exploration and characterization of environments associated with geothermal systems (Arnason et al., 2000; Harinarayana et al., 2006; Farquharson and Craven, 2009; Amatyakul et al., 2015; Blake et al., 2016; Barcelona et al., 2019). In particular, the audiomagnetotelluric technique has been used to explore hydrogeological resources in geothermal environments due to its ability to detect water-bearing rocks in the subsurface (Blake et al., 2016; Siniscalchi et al., 2019). Several factors strongly influence the conductivity of the rock such as temperature, porosity, permeability, poral content and the presence of clay minerals (Spichak and Manzella, 2009). The conductivity of the fluid in the pore space is a factor that influences the resistivity of the rock. In turn, fluid conductivity generally increases with salt concentration, taking extreme values of 0.01 Ωm for hypersaline fluids with 25 wt% salt concentration (Bedrosian, 2007). On the other hand, in rocks with clay minerals, the mechanism that causes the high conductivity is the electrical double layer, that is formed at the interface of the clay mineral surface and the water. This additional conductivity contribution is called surface conductivity (Ussher et al., 2000). Besides, the effect of surface conductivity is generally more important when porosity or permeability is low (Spichak and Manzella, 2009). In turn, temperature is the main factor influencing mineralogical alteration.

This paper presents an AMT model of the first few hundred meters of the region between the Socompa lagoon and Quebrada del Agua. The main objectives are the geoelectrical characterization of the geological formations present in the topographic low, as well as the structures that facilitate the circulation of hydrothermal waters in order to determine the reservoir of such waters at depth.

Section snippets

Geological setting

The generalized stratigraphic column of the region is composed of a pre-Cenozoic crystalline basement over which different Cenozoic units are supported, ranging from the Paleogene to the Quaternary. The regional pre-Cenozoic basement cannot be observed locally in the SGP. It is complex and includes Lower Paleozoic granitic and metamorphic rocks, covered by Devonian quartzities, Permian acidic volcanic rocks, and Upper Cretaceous sands and clays (Gardeweg et al., 1993; Reutter and Munier, 2006).

Method

Audiomagnetotellurics is a passive, frequency-dependent method that utilizes the natural variations of electric and magnetic fields as a source for electromagnetic induction on the Earth (Simpson and Bahr, 2005). The observed total field contains information about the subsurface resistivity structure. The penetration depth achieved is greater while the sounding period and the resistivity of the subsurface increase. The theory, developed in the wave zone (i.e. where the electromagnetic field in

Softwares, grid and data

The 3-D inversion was carried out with the ModEM code (Egbert and Kelbert, 2012; Kelbert et al., 2014). The code uses a non-linear conjugate gradient (NLCG) method as a strategy of minimization. The penalty functional depends on observed and predicted data, and prior model and model covariance. The input files required for the code were constructed using the 3DGrid software provided by Náser Meqbel (personal communication).

The total model length is 16x16x10 Km in x, y and z directions,

Discussion

For a description and correlation of the 3-D model with the geology formations and structures, horizontal slices of the final model shown in Fig. 7 were used. The high electrical resistivity observed at 100 m deep is considered the base of the Quaternary sedimentary cover. Below, between 150 and 250 m depth the lower electrical resistivity determined by the model could be due to the presence of a shallow aquifer. The topographic low where the Socompa lagoon is located corresponds to the hanging

Conclusion

An AMT study was carried out in the area where the Socompa Geothermal Prospect is located. The goal was the geoelectric characterization of the first hundreds meters (1000 m) to try to determine the depth at which the reservoir of the circulating thermal waters would be found. The 3 D model determines a shallow resistivity layer correlated to Quaternary sedimentary cover, a shallow low-resistivity layer associated to a shallow aquifer correlated with some unit of the upper member of CVSQA, and

Authorship conformation form

Guevara L., Favetto A. and Pomposiello C. designed the AMT campaign and analyzed, interpreted and inverted the data. Guevara L. and Apaza F.D. participated in the registration of the data in the field. Apaza and Seggiaro carried out the local geological survey. All authors contributed to the interpretation and approved the final manuscript.

Declaration of Competing Interest

None.

Acknowledgments

We thank RAMSA and SEGEMAR for allowing us to use the data, Gabriel Giordanengo (INGEIS) and Francisco Baigorria (SEGEMAR) for their invaluable technical assistance during the surveys. We also thanks A. Kelbert, and G. Egbert for providing their inversion algorithm and Náser Meqbel for providing his Grid3D software. The TanDEM-X DEMs were provided by the German Aerospace Center (DLR) through proposal DEM GEOL1342 thanks to Pablo Grosse.

References (48)

  • J. Stimac et al.

    Intrusion-related geothermal systems

  • G. Wadge et al.

    The Socompa collapse and avalanche event

    J. Volcanol. Geotherm. Res.

    (1995)
  • K. Arnason et al.

    The resistivity structure of hightemperature geothermal systems in Iceland

  • H. Barcelona et al.

    Audiomagnetotelluric survey at the Bañitos-Gollete geothermal area, main Andes Cordillera of San Juan, Argentina

    Geol. Acta

    (2019)
  • P.A. Bedrosian

    MT+, integrating magnetotellurics to determine Earth structure, physical state, and processes

    Surv. Geophys.

    (2007)
  • T. Bissig et al.

    The cenozoic history of volcanism and hydrothermal alteration in the central andean flat-slab region: new 40Ar-39Ar constraints from the El Indio–Pascua Au (-Ag, Cu) belt, 29 20–30 30 S

    Int. Geol. Rev.

    (2001)
  • J.R. Booker

    The magnetotelluric phase tensor: a critical review

    Surv. Geophys.

    (2014)
  • A. Burd et al.

    Three-dimensional electrical conductivity in the mantle beneath the Payún Matrú volcanic field in the Andean backarc of Argentina near 36.5 S: evidence for decapitation of a mantle plume by resurgent upper mantle shear during slab steepening

    Geophys. J. Int.

    (2014)
  • T.G. Caldwell et al.

    The magnetotelluric phase tensor

    Geophys. J. Int.

    (2004)
  • A. Conde Serra et al.

    Modelo Conceptual Geotermico Preliminar Del Volcan Socompa, Departamento de los Andes, Provincia de Salta, Argentina. Technical Report

    (2020)
  • B.V.W. De Vries et al.

    A gravitational spreading origin for the Socompa debris avalanche

    Journal of Volcanology and Geothermal Re-search

    (2001)
  • P.F. Dobson

    A review of exploration methods for discovering hidden geothermal systems

    GRC Trans.

    (2016)
  • G. Egbert et al.

    Computational recipes for electromagnetic inverse problems

    Geophys. J. Int.

    (2012)
  • M.E. Farías et al.

    The discovery of stromatolites developing at 3570 m above sea level in a high-altitude volcanic lake Socompa, Argentinean Andes

    PloS one

    (2013)
  • Cited by (6)

    • Resource assessment of the Cerro Blanco geothermal system

      2023, Journal of South American Earth Sciences
    • Electrical resistivity structure in the Tocomar geothermal system obtained from 3-D inversion of audio-magnetotelluric data (Central Puna, NW Argentina)

      2022, Geothermics
      Citation Excerpt :

      3-D inversion has emerged as a promising technique to model and image geothermal reservoirs in a single self-consistent manner at presumably optimal accuracy and resolution. Recent advances in 3-D magnetotelluric (MT) inversion codes (Mackie et al., 1994; Siripunvaraporn et al., 2005; Egbert and Kelbert, 2012; Kelbert et al., 2014) and the use of high performance parallel computing allowed to perform 3-D inversions (Heise et al., 2008; Arango et al., 2009; Árnason et al., 2010; Díaz et al., 2012; Bertrand et al., 2012; Piña-Varas et al., 2014; Rosenkjaer et al., 2015; Samrock et al., 2015; Comeau et al., 2015, Díaz et al., 2015; Abera et al., 2018; Guevara et al., 2020, 2021; Slezak et al., 2021). The improvement of MT technology and 3-D inversion facilitates the construction of more realistic subsurface resistivity models of structurally complex areas, such as the Tocomar Geothermal System.

    View full text