The San Diego Area, A New Geothermal Prospect On The Colombian Andes
The San Diego Area, A New Geothermal Prospect On The Colombian Andes
The San Diego Area, A New Geothermal Prospect On The Colombian Andes
The San Diego Area, a New Geothermal Prospect on the Colombian Andes
Maria L. Monsalve1, Ivan D. Ortiz and Claudia M. Alfaro
Servicio Geolgico Colombiano. Diag. 53 34 53. Bogot, Colombia
1
mmonsalve@sgc.gov.co
Monsalve et al.
Geochemical interpretation was done by graphical methods (Giggenbach & Goguel, 1989; Giggenbach, 1988). Estimation of the
recharge zone was done by the functions found by Rodrguez (2004) to correlate deuterium and oxygen-18 concentration with
elevation for Colombia.
3. GEOLOGICAL FRAME
Geologic history of the northwestern corner of South America relates to accretion of mainly allochthonous terranes (e.g. Jones et al.
1982) on western margin of Guyana Shield. These continuous processes of accretions and subductions have created a series of
regional sutures, recognized as faulted areas trending NNE to NE (Figure 1), and a volcanic belt along old continental margins
uplifted by convergence of tectonic plates movements of Nazca and Caribbean Plates subducting a relatively stable South
American Plate (e.g. Etayo et al., 1983; Toussaint & Restrepo, 1992; Montes et al., 2005; Acosta et al., 2007; Taboada et al., 2000;
Suter et al., 2008). Plates interaction result also in three different tectonic blocks to accommodate the tectonic stress: the Panam
Choc Block, the Maracaibo Block and the North Andes Block, the latter where Colombian cordilleras upraise (Figure 1 left), after
Acosta et al. (2007); Suter et al. (2008); Taboada et al. (2000) & Acosta et al. (2004).
Figure 1: Left: Generalized tectonics and structures of north-western South America. Central Cordillera between Cauca
and Magdalena Rivers; PF: Palestina Fault; SFR: Romeral Fault System; SFC: Cauca fault System; IF:Ibague
Fault. [Sources: Acosta et al. (2007); Suter et al. (2008); Taboada et al. (2000); Acosta et al. (2004)]. Right: Detail of
San Diego Geothermal Area on a geological map generalized from Barrero & Vesga (1976) and this work.
The Central Cordillera of Colombia (CC) is the tallest cordillera on the Colombian Andes (up to 5360 m.a.s.l.), it follows a megasuture known as Romeral Fault System, that divides oceanic type rocks on the west form continental type on east. The San Diego
Area is located east from this fault zone. Rocks are composed of Palaeozoic metamorphites on Greenschist to Amphibolite facies,
these in turn host a series of intermediate to felsic intrusions of ages ranging from Triassic to Early Cenozoic (Gonzles, 1989;
Barrero & Vesga, 1976; Barrero et al, 1969). On the area of interest there is also a graben composed of cretaceous shales with rich
on uranium (Naranjo, 1983) (Figure 1 right).
CC also hosts strato-volcanoes of mainly andesitic composition on a set of active volcanic active segments, south segment on the
Western Cordillera of Colombia followed by the central and northern segment on the CC (Figure 1 left).
4. STRUCTURAL FRAMEWORK OF SAN DIEGO GEOTHERMAL AREA
San Diego Geothermal Area lies on the eastern side of the CC, a zone dominated mainly by a regional NNE trending structure
known as Palestina Fault (Figure 1 left), it ranges for around 500 km and it has an influence zone from 0,5 to 30 km wide
(Feinenger, 1970; Page, 1986). This fault acts a wrench fault with evidences of dextral movement of up to 27,7 km with an inverse
component dipping west. Recent activity of this fault was reported by Collins et al (1970) as sinestral displacements.
Aerial photography used to find morphological lineaments and structures (Figure 2 center), showed the general tendency of N-S to
NNE trending lineaments of 2-5 km long. According to Feninenger (1970) the Palestina fault deflects its trace 5-10 toward east,
on the area of the cretaceous sedimentary graben, lineaments also seem to be controlled by lithology where the relatively soft metasedimentary meets the Sintectonic intrusives (Figure 1 right), thus probably such lithological contacts play a notable role on the
volcanic foci location on San Diego and Escondido.
The axis for the Berlin shales syncline trends NNW with a narrowing tip on the NW where the Palestina fault closes the structure
(Pole diagram, Figure 2 left bottom), this could be the result of a transpressive behavior of the fault.
The Palestina fault trace was followed on field, where evidences of inverse dextral displacements were found on the metasedimentary basement. Fault evidence on volcanic deposits of San Diego volcano may indicate some Neogene- Quaternary activity.
2
Monsalve et al.
Further survey of joints and other structural evidences showed tension on a dominating SEE-NWW axis (Figure 2 right) a possible
result of the dominating transcurrent behavior of Palestina fault.
Syncline bedding
Figure 2: Left: Morphological map with location of volcanic structures. On 30 m DEM at 135 altitude and 180 azimuth
Center: Aerial photography lineaments and structures. Right: map of local jointing (on Schmidth / polar net);
bottom left, clusters for syncline bedding show folding axis (1-3)as a NW to NNW trend, minor (blue) tendency
dipping SE probable on narrowing against Palestina fault
Deposits from San Diego maar, are associated with phreatic to phreato magmatic activity, and consist of pyroclastic breccia,
flows and surges, forming the tuff ring. Powerful and very fine ash deposits with significant variation in facies and abundant
accretionary lapilli associated with intense phreatic activity has been associated with the dacitic ryolithic tuff cone.
Beside the tuff cone deposit, three main eruption phases are recognized, associated to the maar its self. The depositional sequences
are defined by abrupt changes in clast size, dominant facies, block abundance and composition, and apparent water in the system.
Evidence for water in the eruptive system includes, among others, abundant, accretionary lapilli, bedding, and lithics originate from
country rock formations. Changes in dominant facies and facies patterns define the distinct units, which reflect changes in eruptive
dynamics with time. The early eruptions are interpreted to have been dominated by efficient, phreatomagmatic blasts, which caused
3
Monsalve et al.
a progressive deepening of the eruptive center, followed by a transition to dominantly magmatic behavior in the upper sequence,
tuff cone and plug dome growth occurring at the later stage of maar activity.Evidence for water in the eruptive system includes,
among others, abundant, accretionary lapilli, bedding, and lithics originate from country rock formations. Changes in dominant
facies and facies patterns define the distinct units, which reflect changes in eruptive dynamics with time. The early eruptions are
interpreted to have been dominated by efficient, phreatomagmatic blasts, which caused a progressive deepening of the eruptive
center, followed by a transition to dominantly magmatic behavior in the upper sequence, tuff cone and plug dome growth occurring
at the later stage of maar activity.
5.2 Other phreatic structures in Berlin
About 8 km SE of San Diego maar, a circular morphologic feature of about 1 km. of diameter was recognized by satellite imagery.
Field work showed that this structure is formed on metamorphic basement and it has at its interior an isolated accretionary lapilli
rich ash sequence, originated by phreato - magmatic activity (Figure 4). At the moment non pyroclastic ring has been identified
around the circular feature nor the pyroclastic deposits out of it.
Monsalve et al.
Selva de Florencia (a National Park) and El Espritu Santo; the latter characterized for first time more than 40 year ago (Arango et
al., 1970) (Figure 6). Being located in the northernmost Andean Central Cordillera, these springs are located at relatively low
elevations (500-1200 m.a.s.l.) compared with the rest of them located to the South and related to the hydrothermal systems of other
volcanoes as Cerro Bravo, Nevado del Ruiz, Tolima, Huila and Purac (between 1700 and 4520 m.a.s.l). They have low flow rate
(<1L/s), low temperature discharge (<43C) and notable gas emissions with non-detectable concentrations of H2S (by organoleptic
analysis). The highest temperature is registered in some springs from Selva de Florencia and El Espritu Santo (Table 1), in neutral
sodium chloride and neutral mixed sodium chloride bicarbonate sulfate water, respectively. The lowest temperatures (33C) are
found in bicarbonate springs from Laguna de San Diego.
Laguna de
San Diego
El Espritu
Santo
Selva de
Florencia
Figure 6: Location of hot springs in San Diego Geothermal area on geological map (Barrero & Vesga, 1976) . Selva de
Florencia (6 chloride springs), Laguna de San Diego (3 bicarbonate springs), El Espritu Santo (5 Mixed ChlorideBicarbonate-Sulfate springs). Dotted yellow line highlights the preliminary defined geothermal area.
Table 1. Hot springs within San Diego Geothermal Area (data from Ortiz & Alfaro, 2011)
Hot Spring
Town
Department
Locality
Elevation
(msnm)
pH
Flow (L/s)
La Calera
Laguna San Diego
Geyseres I
Geyseres II
Geyseres III
Geyseres IV
Finca Termales
Termal San Antonio
Espiritu Santo I
Espiritu Santo II
Espiritu Santo III
Espiritu Santo IV
Espiritu Santo V
Saman
Saman
Saman
Saman
Saman
Saman
Saman
Saman
Nario
Nario
Nario
Nario
Nario
Caldas
Caldas
Caldas
Caldas
Caldas
Caldas
Caldas
Caldas
Antioquia
Antioquia
Antioquia
Antioquia
Antioquia
540
870
1000
998
983
967
1200
780
770
779
801
790
795
6,5
5,99
6,22
6,26
6,1
6,3
5,99
6,15
6,3
6,0
6,2
6,2
5,92
1
0,5
0,5
0,2
1
1
1,5
0,33
0,25
0,9
0,25
<0,1
<0,1
Temperature
EC (mS/cm)
(C)
32,5
33
33
39
40
42
40
32
43,5
35,5
43,3
36,5
33,4
3530
1345
3200
7770
6500
9190
3420
8740
7410
3920
6720
4960
3870
TDS
(mg/L)
2206
1186
2270
4658
3696
4844
2026
5166
4510
2458
4368
3200
2540
Chemical Classification
Sodium Bicarbonate
Calcium Bicarbonate
Sodium Chloride
Sodium Chloride
Sodium Chloride
Sodium Chloride
Sodium Chloride
Sodium Chloride
Sodium Chloride-Bicarbonate-Sulfate
Sodium Chloride-Bicarbonate-Sulfate
Sodium Chloride-Bicarbonate-Sulfate
The relative chloride-sulfate-bicarbonate diagram (Figure 6, Left) illustrates the differences in composition of these springs. Those
from Selva de Florencia of chloride water type, as it was mentioned before, have the highest probability of reservoir fluid
contribution. This is consistent with their high lithium contents (5-14 ppm, SGC, 2014)). Hot springs from El Espritu Santo show
an important increase in the concentration of sulfate (20% or 320-800 mg/l, SGC, 2014). Although the origin of the sulfate could be
related to mixing with steam heated waters (or sulfate water from sulphur hydrolysis), the pH do not show variation to more acidic
as it would be expected.
On the other hand a linear trend in the variation of calcium and strontium versus sulfate is found, as represented in the Figure 8. It
could suggest a common origin of these dissolved species in El Espiritu Santos hot springs. The Sr is remarkable high (up to 13
mg/l) in front of the springs of this geothermal area and also compared with all the hot springs located in the Central Cordillera
(lower than 2 mg/l in neutral springs and lower than 8 mg/l in acidic springs with high evidence of magmatic fluids contribution
(SGC, 2014).
5
Monsalve et al.
The observed differences in the chemical composition can be attributed to its origin in the case of bicarbonate springs from Laguna
de San Diego, to peripheral waters of the geothermal system rich in CO 2. Differences between the springs from the localities of El
Espritu Santo and Selva de Florencia, could be related to mixing with other water sources of high sulfate and strontium
concentration in an environment of crystalline (metamorphic and igneous) and volcanic rocks (evaporates are not expected in the
zone), an intense interaction with rocks of very different composition along their assent to the surface, or to two different
geothermal reservoirs.
The sodium-potassium-magnesium diagram presented in the Figure 7 (Right) highlights the high dilution of the geothermal fluid
with high relative magnesium concentration. The extrapolation of the dilution trend to the equilibrium line indicate a probably
reservoir equilibrium temperature above 280C.
Na/1000
0
10
10
Cl
0
Equilibrium Line
180
50
28
26
Partial Equilibrium
30
20
75
3
3 40
VO
LC
AN
IC
Immature Waters
0
75
100
10
75
HCO
50
25
12
25
RS
TE
WA
10
LAGUNA DE
SAN DIEGO
25
L
RA
HE
=
SO
4 0
140
80
RIP
PE
EL ESPIRITU
SANTO
160
24
22
10
50
SELVA DE
FLORENCIA
0
20
50
50
25
WA
TE
RS
25
75
75
MATURE
WATERS
K/100
Mg
0
25
50
75
0.5
100
Figure 7. Relative chemical compositions of hot springs from San Diego Geothermal Area. Left: Cl-SO4-HCO3 (after
Giggenbach & Goguel, 1989). Chloride waters are found in Selva de Florencia, assumed as those with a higher fluid
reservoir contribution. Springs from El Espritu Santo show a higher sulfate concentration which origin has not
being determined. Right: Na-K-Mg (After Giggenbach, 1988). Blue: Bicarbonate springs from Laguna de San Diego.
Violet: Chloride-bicarbonate-sulfate springs from Selva de Florencia. Red: Chloride Springs from Selva de
Florencia. Typical high dilution with fresh water is found in hot springs of the three localities. The dilution trend
points out a probable equilibrium temperature higher than 280C for h ot springs from Selva de Florencia.
Figure 8. Calcium (left) and strontium (right) concentration Vs. Sulfate concentration. A directly proportional variation of
calcium and strontium is observed in hot springs from El Espritu Santo, which suggests the possibility of a genetic
relation. The hot springs of Laguna de San Diego and Selva de Florencia do not show such a variation in the case of
Ca Vs SO4. In the case of strontium, although the available results indicate very low strontium level, just few of them
were analysed. Colour conventions as in the Figure 7 .
Other aqueous geothermometers based on silica, sodium /potassium , presented in Table 2, coincide with the possibility of a high
temperature geothermal reservoir in this geothermal area (even above 300C) . The temperatures estimated by the K/Mg
geothermometer is much lower due to the evident dilution process, previously cited.
The water stable isotopes indicate concentrations of deuterium between -57 and -43 and oxygen-18, between -8.2 and -5.2
(Figure 9). Coinciding with the highest discharge temperature of a chloride spring, the highest oxygen-18 enrichment is found in
6
Monsalve et al.
the hottest neutral sodium chloride water from Selva de Florencia, which is coherent with its slightly greater proximity to maturity
character. The recharge zone would be located at the elevation of 1250-1500 (La Calera) and 1920 m.a.s.l. (Laguna San Diego),
between 1800 to 2000 m.a.s.l. for springs from El Espritu Santo and between 1550-1700 m.a.s.l., according to the estimation based
on deuterium contents. From this, there are local effects defining the isotopic composition of the hot springs in each locality.
Table 1. Geochemical temperatures of deep fluids estimated in San Diego geothermal area
Tqtz (C)
Temperature
Tchalcedony T Na/K (C)
no Steam
(C)
(C) (2)
(3)
loss (1)
La Calera
32,5
123
95
237
Laguna San Diego
33
161
137
205
Geyseres I
33
168
146
323
Geyseres II
39
170
147
306
Geyseres III
40
184
164
295
Geyseres IV
42
207
191
289
Finca Termales
40
177
156
279
Termal San Antonio
32
155
130
243
Espiritu Santo I
43,5
130
103
231
Espiritu Santo II
35,5
Espiritu Santo III
43,3
122
94
241
Espiritu Santo IV
36,5
114
85
243
Espiritu Santo V
33,4
(1) TC = (1309/(5,19-LOG(SIO2)))-273. Fournier, 1977. In Nicholson, 1993
(2) TC = (1032/(4,69-LOG(SiO2)))-273. Fournier, 1977. In Nicholson, 1993
(3) TC = (1390/(LOG(AH2/AI2)+1,75))-273. Giggenbach, 1988
(4) TC = (855,6/(LOG(AH2/AI2)+0,8573)) - 273,15. Truesdell (1976)
(5) TC = (4410/(14-LOG(k^2/Mg)))-273. Giggenbach, 1998
Hot Spring
T Na/K
(C)
(4)
T K/Mg (C)
(5)
194
152
322
295
277
267
253
202
186
101
58
108
121
116
123
87
116
106
199
201
104
101
Figure 9. Stable isotopes composition of hot springs from San Diego Geothermal Area. The highest oxygen-18 enrichment
is found in a sodium chloride hot spring from the location of Selva de Florencia. Color conventions as in the Figure 7.
7. CONCLUSIONS
Initial results for geothermal exploration of the San Diego Area narrowed estimates for potentiality of the system, as exploitation
prospect showed favorable evidences for use of this heat source.
Studies in the area of geothermal interest in San Diego, allowed the identification of phreatic phreatomagmatic deposits
associated to San Diego maar. Other possible volcanic foci could exist in the area as several isolated outcrops of pyroclastic rocks
have been mapped in regional geological work, if so the area could be considered as a volcanic province that is offset from the axis
of the Cordillera Central, where the active volcanoes are located.
The predominant type of deposits (phreatic - phreatomagmatic) and its relatively recent age 33,000 years for volcanic deposits
associated to the newly discovered Escondido volcano, indicate from the volcanological point of view, a significant geothermal
potential in area to be verified with the results of studies in other disciplines.
Structural analyses showed evidences for promising geothermal reservoirs, based on the dominant transpressive regime on the area
imposed mainly by the Palestina fault.
San Diego geothermal area probably hosts a high temperature geothermal system as it is inferred by aqueous geothermometers.
Chemical differences between hot springs from Selva de Florencia (neutral sodium chloride waters with low Sr concentration) and
El Espritu Santo (neutral mixed sodium chloride-bicarbonate-sulfate waters with high Sr concentration) could be related to
7
Monsalve et al.
differences in processes of rock water interaction or mixing occurring in the two discharge areas or to their origin from different
reservoirs.
Geochemical analysis of rocks as well as isotope strontium analysis for rock and waters should be done to understand the origin of
these hot springs.
8. REFERENCES
Acosta, J.; Velandia, F.; Osorio, J.; Lonergan, L; Mora; H. 2007. Strike-slip deformation within the Colombian Andes. Captulo de
Libro: Deformation of the Continental Crust: The Legacy of Mike Coward, Vol. 272, pgs:303 - 319, Ed. The Geological
Society of London. UK.
Acosta, J.; Lonergan, L.; Coward, M. P. 2004. Oblique transpression in the western thrust front of the ColombianEastern
Cordillera. Journal of South American Earth Sciences 17 (2004) 181194.
Arango, E., Buitrago, J., Cataldi, R., Ferrara, G. C., Panichi, C. and Villegas, V. 1970. Preliminary Study on the Ruiz Geothermal
Project (Colombia). Geothermica. Special Issue 2. Vol3. 2. Part I:43-59
Barrero & Vesga, 1976. Geologa de la Plancha 188. La Dorada. Map Scale 1:100.000. INGEOMINAS. Bogot.
Barrero, D.; lvarez, J. & Kassem, T. 1969. Actividad gnea Y Tectnica En La Cordillera Central Durante El Meso-Cenozoico.
Documento Tcnico. Publicacin autorizada por el Director del Instituto Nacional de Investigaciones Geolgico - Mineras. 49
Pags. INGEOMINAS. Bogot.
Collins, D. E.; Benarczar, G. & Page, W. D. 1981. Quaternary activity on the Palestina Fault Zone, northwestern Colombia.
Revista CIAF Vol. 6 (1-3), 117. Bogot.
Etayo, F. Barrero, D. Instituto Nacional de Investigaciones Geolgico-Mineras y la participacin de U.S. Geological Survey. 1983.
Mapa de terrenos geolgicos de Colombia. Bogot.
Feininger, T. 1970. The Palestina Fault, Colombia. Geological Society of America Bulletin, v. 81, p. 1201-1216. USA.
Giggenbach, W. F. 1991. Chemical techniques in geothermal exploration. In Application of geochemistry in geothermal reservoir
development Unitar. Edited by Franco DAmore. Unitar/UNDP. Roma. p119-144.
Giggenbach, W. F.Goguel, R. L. 1989. Collection and analysis of geothermal and volcanic water and gas discharges. Department of
Scientific and Industrial Research. Chemistry Division. Report No. CD 2401.Fourth Edition. Petone, NZ. 81 p
Giggenbach, W. F. 1988. Geothermal solute equilibria. Derivation of Na-K-Mg-Ca geoindicators: Geochimica Cosmochimica
Acta, Vol. 52, p. 2749-2765
Gonzles, H. 1989. Anlisis de la nomenclatura estratigrfica de las rocas metamrficas (Litodema A) al este del lmite oriental
de la zona de Falla de Romeral, Cordillera Central, Colombia. Informe Interno. INGEOMINAS. Medelln.
Jones, D. L., Howell, D. G., Coney, P. J. and Monger, J. W. H., 1982. Recognition, character, and analysis of
tectonostratigraphic terranes in western north America: Center for Academic Publication Japan, Oji Seminar Volume
(in press)
Montes, C.; Hatcher, R.; Restrepo-Pace, P. 2005. Tectonic reconstruction of the northern Andean blocks: Oblique convergence and
rotations derived from the kinematics of the PiedrasGirardot area, Colombia. Tectonophysics 399 (2005) 221 250.
Colombia.
Naranjo, J. L. 1983. Research on the uraniferous potential in the black shales of the Berlin-Caldas synclinal. Unpublished
Report, University of Bogota, Colombia (in Spanish).
Nicholson, K. 1993. Geothermal fluids. Chemistry and Exploration Techniques. Springer-Verlag. Berlin. 263 p.
Page, W.D. 1986. Seismic Geology and Seismicity of Northwestern Colombia. Compilacin ISA. Medelln
OLADE, 1982. Estudio de Reconocimiento de los Recursos Geotrmicos de la Repblica de Colombia. Pisa, Italy.
Ortiz, I. and Alfaro, C. 2011. Inventario Nacional de Manantiales Termales. Fase 2011. Caldas, Risaralda, Tolima and Nario.
INGEOMINAS. Informe Tcnico. Bogot. 131 p.
Ortiz, I., 2013. Informe de reconocimiento rea de San Diego Florencia Puente Linda. Unpublished report. SGC. Bogot.
Rodrguez, C., 2004. Lnea meterica isotpica de Colombia. Meteorologa Colombiana. Bogot. pp 43-51
SGC, 2014. Base de datos de Inventario Nacional de Manantiales Termales de Colombia. Grupo Exploracin de Recursos
Geotrmicos.
Suter, F.; Sartori, M.; Neuwerth, R.; Gorin, G. 2008. Structural imprints at the front of the Choc-Panam indenter: Field data from
the North Cauca Valley Basin, Central Colombia. Tectonophysics 460 (2008) 134157. Suiza.
Taboada, A.; Rivera, L. A.; Fuenzalida, A.; Cisternas, A.; Philip, H.; Bijwaard, H.; Olaya, J. & Rivera, C. 2000. Geodynamics of
the northern Andes: Subductions and intracontinental deformation (Colombia). Tectonics, Vol. 19, No. 5. Pags. 787 813.
Toro, G. E. 1989. Caractertzacin del vulcanismo de San Diego y estudios de los depsitos de San Diego (Caldas) y de Nario
(Antioquia), Colombia. Tomo I de las Memorias del Congreso Colombiano de Geologa Pgs. 419-439. Bucaramanga.
Monsalve et al.
Toro 1989. Caracterizacin del vulcanismo de San Diego y estudio de los depsitos de San Diego (Caldas) y de Nario (Antioquia),
Colombia. Memorias del V congreso Colombiano de Geologa, Bucaramanga. Marzo 14 -17, p 419 440.
Toro 1991. Geologa de la zona del volcn de San Diego, departamento de Caldas Colombia. Simposio Magmatismo Andino.
Manizales Colombia.
Monsalve, M. L. y Ortiz, I., 2013. Reconocimiento vulcanolgico rea geotrmica de San Diego. Unpublished report. SGC.
Bogot.
Toussaint, J. F. & Restrepo, J.J. 1992. Cronologa de las acreciones de terrenos alctonos en los Andes Colombianos. Documento
contribucin al Proyecto 279 Terrenos en Amrica Latina del PICG. Pg. 219 221. Medelln.