Outline
Comptes Rendus

Hydrology, hydrogeology
Water–rock interaction in the geothermal systems related to post-collision zone volcanism: A case study based on multivariate statistical analysis from the Kavak geothermal field (Konya, Turkey)
Comptes Rendus. Géoscience, Volume 355 (2023), pp. 311-329.

Abstracts

Water–rock interaction is the focus of geothermal energy studies and can be documented by traditional geochemical methods such as ion ratio method and hydrogeochemical modelling etc. Statistical approaches are also vital for the quantitative models, and mainly combined with the traditional methods. In this study, we re-evaluate the published data, including water chemistry and volcanic and metamorphic whole-rock geochemistry from the Kavak geothermal field (Konya, Turkey) by using multivariate statistical analysis (e.g. factor analysis and clustering analysis) to research possible interaction between the thermal waters and surrounding rocks.

The Kavak geothermal field (KGF) overlies a metamorphic basement composed of the Paleozoic metacarbonates and metaclastic rocks, yet is located near the Erenlerdağ–Alacadağ volcanic complex (ErAVC). An example of unimodal orogenic volcanism in an extensional geodynamic setting, the Neogene ErAVC is composed of widespread high-K calcalkaline andesite to rhyodacite lavas with enclaves and their pyroclastic counterparts. The Kavak geothermal field covers a small area (7.5 km 2 ) and lies along the Seydişehir fault zone in the southeast of the ErAVC. The Kavak thermal waters are meteoric in origin and peripheral waters (Ca–Na–HCO 3 ) in the geothermal system related to the orogenic volcanism. The Kavak thermal waters are characterised by high K + and Na + cations, and low pH (between 6.4–6.9 pH) values relative to the cold waters around the KGF. Two types of thermal waters were identified in the KGF based on the slight difference in their hydrochemistry and discharge temperature. The first type thermal water (22 °C) has higher TDS and Cl/Br ratio and lower dissolved silica and Br content relative to the second type of water (up to 45 °C). The chemical relationship between the KGF and high-K ErAVC is clearly seen in the cation-based diagrams. Multivariate statistical analysis confirms that first type and second type thermal waters identified based on their hydrochemistry formed two separate statistical groups, and suggests that the chemistry of the KGF waters was mainly controlled by the composition of the ErAVC rather than those of the basement metamorphic rocks as a result of water–rock interaction.

L’interaction eau–roche est au centre des études sur l’énergie géothermique et peut être documentée par des méthodes géochimiques traditionnelles telles que la méthode du rapport d’ions et la modélisation hydrogéochimique, etc. Les approches statistiques sont également essentielles pour les modèles quantitatifs, et principalement combinées aux méthodes traditionnelles. Dans cette étude, nous réévaluons les données publiées, y compris la chimie de l’eau et la géochimie des roches volcaniques et métamorphiques du Champ géothermique de Kavak (Konya, Turquie) à l’aide d’une analyse statistique multivariée (telle que l’analyse factorielle, et l’analyse de regroupement) afin de rechercher des indices d’interaction possible entre les eaux thermales et les roches environnantes.

Le champ géothermique de Kavak (KGF) repose sur un socle métamorphique composé de métacarbonates paléozoïques et de roches métaclastiques, situé près du complexe volcanique d’Erenlerdağ–Alacadağ (ErAVC). En tant qu’exemple de volcanisme orogénique unimodal dans un cadre géodynamique extensionnel, l’ErAVC du Néogène est majoritairement composé d’andésite calcoalcaline à haute teneur en potassium, ainsi que de rhyodacite à enclaves et leurs produits pyroclastiques associés. Le champ géothermique de Kavak (KGF) couvre une petite superficie (7,5 km 2 ) et se situe le long de la zone de faille de Seydişehir au sud-est de l’ErAVC. Dans ce système géothermique lié au volcanisme orogénique, les eaux thermales de Kavak (T22–45 °C) sont d’origine météorique et périphérique (Ca–Na–HCO 3 ). Ces eaux sont caractérisées par des teneurs en cations K + et Na + élevées et un pH bas (entre 6.4–6.9) relatif aux eaux froides autour du KGF. Deux types d’eaux thermales ont été identifiés dans le KGF, grâce à de légères différences hydrogéochimiques et de températures de sortie. Le premier type d’eau thermale (T22 °C) possède une TDS et un rapport Cl/Br relativement élevé ainsi que des teneurs plus faibles en silice dissoute et en Br (T>45 °C), ceci vis-à-vis du second type d’eau. La relation géochimique liant le KGF et les formations riches en potassium de l’ErAVC sont clairement mises en évidence par des diagrammes basés sur les cations. L’analyse statistique multivariée confirme que les deux types d’eaux thermales forment deux groupes séparés statistiquement. Ceci suggère que la chimie des eaux du KGF est principalement contrôlée par le biais de l’interaction eau-roche avec l’ErAVC plutôt qu’avec les roches métamorphiques du socle.

Metadata
Received:
Revised:
Accepted:
Published online:
DOI: 10.5802/crgeos.249
Keywords: Kavak geothermal field, Water–rock interaction, Hydrogeochemistry, Data analysis of the thermal waters, Silicate weathering
Mot clés : Champ géothermique de Kavak, Interaction eau–roche, Hydrogéochimie, Analyse des données des eaux thermales, Altération des silicates

Mesut Gündüz 1; Ayla Bozdağ 2; Ali Ferat Bayram 2; Ali Bozdağ 2; Kürşad Asan 2; Paul Sardini 3

1 Graduate School of Natural and Applied Sciences, Muğla Sıtkı Koçman University, TR-48000, Muğla, Turkey
2 Geological Engineering Department, Konya Technical University, TR-42250, Konya, Turkey
3 Institut de Chimie des Milieux et Matériaux de Poitiers (IC2MP), CNRS, (E2-HydrASA), Université de Poitiers, F-86073, Poitiers, France
License: CC-BY 4.0
Copyrights: The authors retain unrestricted copyrights and publishing rights
@article{CRGEOS_2023__355_G2_311_0,
     author = {Mesut G\"und\"uz and Ayla Bozda\u{g} and Ali Ferat Bayram and Ali Bozda\u{g} and K\"ur\c{s}ad Asan and Paul Sardini},
     title = {Water{\textendash}rock interaction in the geothermal systems related to post-collision zone volcanism: {A} case study based on multivariate statistical analysis from the {Kavak} geothermal field {(Konya,} {Turkey)}},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {311--329},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {355},
     year = {2023},
     doi = {10.5802/crgeos.249},
     language = {en},
}
TY  - JOUR
AU  - Mesut Gündüz
AU  - Ayla Bozdağ
AU  - Ali Ferat Bayram
AU  - Ali Bozdağ
AU  - Kürşad Asan
AU  - Paul Sardini
TI  - Water–rock interaction in the geothermal systems related to post-collision zone volcanism: A case study based on multivariate statistical analysis from the Kavak geothermal field (Konya, Turkey)
JO  - Comptes Rendus. Géoscience
PY  - 2023
SP  - 311
EP  - 329
VL  - 355
PB  - Académie des sciences, Paris
DO  - 10.5802/crgeos.249
LA  - en
ID  - CRGEOS_2023__355_G2_311_0
ER  - 
%0 Journal Article
%A Mesut Gündüz
%A Ayla Bozdağ
%A Ali Ferat Bayram
%A Ali Bozdağ
%A Kürşad Asan
%A Paul Sardini
%T Water–rock interaction in the geothermal systems related to post-collision zone volcanism: A case study based on multivariate statistical analysis from the Kavak geothermal field (Konya, Turkey)
%J Comptes Rendus. Géoscience
%D 2023
%P 311-329
%V 355
%I Académie des sciences, Paris
%R 10.5802/crgeos.249
%G en
%F CRGEOS_2023__355_G2_311_0
Mesut Gündüz; Ayla Bozdağ; Ali Ferat Bayram; Ali Bozdağ; Kürşad Asan; Paul Sardini. Water–rock interaction in the geothermal systems related to post-collision zone volcanism: A case study based on multivariate statistical analysis from the Kavak geothermal field (Konya, Turkey). Comptes Rendus. Géoscience, Volume 355 (2023), pp. 311-329. doi : 10.5802/crgeos.249. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.5802/crgeos.249/

Version originale du texte intégral (Propose a translation )

1. Introduction

Water–rock interaction (WRI), one of the most important factors that deeply control the geochemistry of thermal waters is a prime interest in geothermal energy studies. Although this phenomenon has been known for many years by earth scientists, its nature has always been debated [Wohletz and Heiken 1992]. WRI encompasses primarily three types of reactions, namely adsorption, dissolution, and precipitation that take place at the mineral-fluid interfaces [Teng 2005]. WRI can be documented by theoretical, experimental, field, and laboratory studies, based mainly on the geochemical and isotopic data. After the acquisition data, WRI can be studied by traditional methods such as ion ratio method and hydrogeochemical modelling software (e.g. PHREEQC). Multivariate statistical analysis techniques have recently been combined with these methods. Factor analysis and principal component analysis (PCA) are commonly used multivariate statistical analysis techniques [Poznanović-Spahić et al. 2023; Zhang et al. 2022]. These statistical methods are useful to analyse all parameters in a large dataset as a whole, which allows practically to screen, reduce and classify the data. Compared to the traditional methods, they serve as a more powerful tool on the deciphering origin, interaction with the rocks and circulation pathway of thermal waters [Hao et al. 2020; Yidana et al. 2012].

The Kavak geothermal field (KGF) is a fault-controlled thermal system of meteoric origin located in the southwest of the Erenlerdağ–Alacadağ volcanic complex [Bayram 1992; Bozdağ 2016; Burçak et al. 2002; Davraz et al. 2022a, b; Göçmez et al. 2005; Göçmez and Şen 1998] (Figures 1 and 2). Using traditional methods, previous studies emphasised that water–rock interactions and decomposition processes, and argued the interaction between the thermal waters and metamorphic basement rocks played an important role in the KGF [Bozdağ 2016; Davraz et al. 2022a, b; Karaisaoğlu and Orhan 2018]. However, the possible interaction between the thermal waters and volcanic rocks of the Erenlerdağ–Alacadağ volcanic complex (e.g. ErAVC) as a geothermal heat source in the KGF was not evaluated in the previous studies. Hence, we here reprocess the published geochemical data from the thermal waters, volcanic and metamorphic rocks in the KGF using multivariate statistical analysis (e.g. factor analysis and clustering analysis) and evaluate whether interactions between the Kavak thermal waters and volcanic rocks in the ErAVC play a decisive role in water chemistry.

Figure 1.

Simplified tectonic map of Turkey showing the major suture zones [Okay and Tüysüz 1999], arc systems (red lines), continental blocks (black lines), and distribution of the main Post-Collision (Late Cenozoic) volcanic fields [MTA 2013]. İAESZ: İzmir-Ankara-Erzincan Suture Zone, BSZ: Bitlis Suture Zone, ASZ: Antalya Suture Zone, NAF: North Anatolian Fault, EAF: East Anatolian Fault, DSF: Dead Sea Fault, AF: Afyon Fault, EAVP: East Anatolian volcanic provence, KV: Karacadağ volcanics, CAVP: central Anatolian volcanic provence, GVP: Galatian volcanic provence, ErAVC: Erenlerdağ–Alacadağ volcanic provence, KAIVP: Kırka-Afyon-Isparta volcanic provence, WAVP: West Anatolian volcanic provence.

Figure 2.

Simplified geological map of the west of Konya [Keller et al. 1977; General Directorate of Mineral Research and Explorations 1/100,000 geology map] from Gündüz [2023].

2. Geological and hydrogeological setting

Located in the Central Taurus, the KGF is in the Beyşehir Basin characterised by a graben-like structure, including several units with different stratigraphic, lithologic, tectonic, volcanic, and metamorphic features (Figure 2). The Beyşehir Basin, which controls the morpho-tectonic structure of the Kavak geothermal field, is mostly related to the northwest–southeast trending (sub) parallel normal faults. The Mio-Pliocene fluvio-lacustrine deposits around the KGF are proof of the fault-controlled intermontane basin in the region [Koç et al. 2017, 2018]. In addition, this graben system (i.e. Beyşehir Basin) is surrounded by the Taurus belt from the SW and the Afyon zone from the NE. It is also separated from the Konya graben in the north. The morpho-tectonic structures of the study area and seismic activities in recent years indicate that the faults in the region (i.e. Abazdağı, Beyşehir, and Seydişehir fault zones, etc.) are a part of the active-neotectonic structure [Aksoy 2019]. The Kavak geothermal field is an example of the fault-controlled spring system [Keegan-Treloar et al. 2022] and thermal water springs of the KGF discharge along a fault in the north (e.g. TS4 and TS5 thermal waters) and fault in the south (e.g. TS1, TS2 and TS3 thermal waters) parallel to the Seydişehir Fault Zone (Figures 23). The southern fault acts actually as a structural (semi) barrier, leading to the existence of two types of thermal waters in terms of their chemical compositions, surface temperatures and location.

Figure 3.

Geology map of the study area and location of thermal springs and drilling wells [Bozdağ 2016].

The Central Tauride is an assemblage of various autochthonous and allochthonous crystalline-metamorphosed basement rocks [Dean and Monod 1970; Göncüoğlu et al. 2007; Göncüoğlu and Kozlu 2000; Gürsu et al. 2003; Moix et al. 2008; Özgül 1976; Robertson et al. 2013; Şengör et al. 2019; Turan 2010]. The basement rocks around the KGF consist of Paleozoic metacarbonates and metaclastic rocks intercalated with metamarls [Karadağ 2014]. These basement rocks are unconformably overlain by the Neogene-Quaternary sedimentary and volcanic rocks [Eren 1996; Hakyemez et al. 1992; Koç et al. 2012; Özkan and Söğüt 1998; Turan 2020].

The Neogene volcanic rocks (ErAVC) cover approximately 1500 km2 area (SW-Konya) and are separated from the Sulutas volcanic complex (SVC) by the Anatolide (Afyon zone) [Asan et al. 2021] (Figure 2). The ErAVC rocks are characterised by a large compositional range (53 wt% to 72 wt% SiO2) and a unimodal distribution [Asan 2017; Asan et al. 2021] in the TAS diagram (Figure 4), and they were plotted into basaltic andesite represented by enclaves to rhyolite fields [Le Bas et al. 1986]. The ErAVC rocks are calc-alkaline in the AFM diagram of Irvine and Baragar [1971], and have high-K content (Figures 5 and 6). These volcanic rocks were first described as calc-alkaline andesite and dacite, rarely basaltic andesite and rhyolitic composition, and they were suggested to be related to subduction-related magmatism by Keller et al. [1977]. The authors presented the first geological, geochemical (major oxides, very limited trace element, and Sr isotopes), and geochronological (K/Ar radiometric age) data of these rocks by Besang et al. [1977]. In these studies, it was shown that volcanic products are composed of lava domes/lava flow, pyroclastic fall and flow (e.g. block-and ash flows and ignimbrites) deposits [Keller et al. 1977]. It was stated by the researchers that the volcanism prolonged a wide period from Miocene (11.95 ± 0.02 Ma) to Pliocene (3.35 ± 0.08 Ma). Sr–Nd isotopic data of the ErAVC presented by Temel et al. [1998] show that the intermediate and felsic volcanic rocks have high Sr and low Nd isotopes. The authors conclude that the rocks were related to the subduction of the African plate under the Anatolian plate. On the other hand, Koçak and Zedef [2016] divided Konya volcanics into eight different lithological units and stated that these lithological units mainly contain mafic minerals (Mg#0.69–0.71) (e.g. pyroxene “augite and enstatite”, calcic amphibole, biotite, opaques) and felsic minerals (e.g. plagioclase “An30–86” and fewer quartz phenocryst) (Figure 7).

Figure 4.

Geochemical classification diagrams of the ErAVC; pointed [Keller et al. 1977], gray dashed-line [Temel et al. 1998], pink-colored [Asan 2017; Asan et al. 2023]: Total-Alkali Silica (TAS) diagram of [Le Bas et al. 1986]. Alkaline-sub alkaline division (dashed line) from Irvine and Baragar [1971]. All analyses were recalculated to 100% on a water-free basis for TAS classification.

Figure 5.

Sub-alkaline ErAVC rocks; pointed [Keller et al. 1977], gray dashed-line [Temel et al. 1998], pink-colored [Asan 2017; Asan et al. 2023] on the AFM diagram [Irvine and Baragar 1971].

Figure 6.

Plot of SiO2 vs. K2O (wt%) illustrating andesite types of Gill [1981]; pointed [Keller et al. 1977], pink-colored arrow [Asan 2017; Asan et al. 2023].

Figure 7.

Microphotos of main mineralogical and textural features of the ErAVC rocks; (a) orthopyroxene, (b) clinopyroxene, (c) amphibole (hornblende), (d) biotite, (e) plagioclase, (f) K-feldspar (sanidine), (g) quartz and (h) opaque minerals (ilmenite or magnetite).

The Paleozoic metacarbonates that seem to be the main reservoir in the KGF have high permeability because of their highly fractured and karstified characteristics. Fractured quartzites of the Paleozoic metaclastic rocks are permeable whereas phyllite, clay, and metasiltstone within the Paleozoic metaclastics are relatively impermeable. Basin-filling limestones, conglomerates, and sandstones of the Neogene-Quaternary sedimentary deposits are permeable, but silty and clayey layers are impermeable [Bozdağ 2016]. Andesitic to rhyodacitic lavas from the ErAVC rocks acting as possible caprocks is impermeable whereas their pyroclastics (e.g. tuff, lapilli tuff and breccia, etc.) are permeable. The ErAVC pyroclastics serving as possible geothermal reservoir or aquifer are widespread from the central to distal facies where they are interbedded with basin-filling fluvio-lacustrine deposits (Figures 2 and 3).

3. Material and methods

The published geochemical data from water and whole-rock samples in the KGF were used in this study. The water samples (e.g. thermal wells, springs, and cold waters) were analysed by Bozdağ [2016] for major anion–cation, trace elements, and isotope compositions in different geochemical laboratories (see the given reference for the detailed analytical methods of water samples). Metamorphic rock samples (e.g. metaclastics and metacarbonate) from the basement [Karadağ 2014] and volcanic rock samples (e.g. lava flows, domes, enclaves and their pyroclastic equivalents) from the ErAVC [Asan 2017; Asan et al. 2023; Karakaya 2009; Keller et al. 1977; Temel et al. 1998] were analysed for major oxide and trace elements using XRF, ICP-OES or ICP-MS (see the given references for the detailed analytical methods of whole-rock samples).

In this study, cluster analysis and factor analysis were chosen as data analysis methods using GEOstats [Gündüz and Asan 2022], a statistical data analysis program to perform these complex multivariate analyses. Using a class of cluster analysis designed to find groups of similar items within a data set, hierarchical cluster method produces a hierarchy of clusters, ranging from small clusters of very similar items to larger clusters of increasingly dissimilar items. The data matrix of cluster analysis is to be in standard form, n × p matrix with n rows of samples (e.g. sampling locations) and p columns of variables (e.g. geochemical data, discharge temperature, EC etc.). Computation may be in “Q-mode analysis (e.g. matrix size n × n)” or “R-mode analysis (e.g. matrix size p × p)”. R-mode analysis calculates distances (similarities) among all pairs of variables whereas Q-mode analysis calculates a matrix of distances between all pairs of samples. The hierarchical cluster analysis produces a graph known as a dendrogram. In dendrogram, the most similar two samples (e.g. A and B) are joined to form the first cluster, followed by another similar sample(s) (e.g. C or D) forming the second cluster, and so on. In this study, we first used “R-mode” to classify and characterise the KGF thermal waters based on the geochemical data, leading to the identification of homogeneous groups or subgroups in the data set. This grouping has the advantage to compare between water chemistry and rock or mineral chemistry, simplifying to interpret water–rock interactions. For example, if K, Na and Ba analysis results from thermal waters form a cluster, this can be interpreted to the interaction of thermal waters with a rock including alkali feldspar. On the other hand, Q-mode analysis was used to classify sampling locations into similar groups, and to identify thermal water types contributing to these groups. In this study, water samples taken along the KGF from 12 locations of which 4 locations are from cold waters and 8 locations from thermal waters (e.g. Cold Waters: CW1, CW2, CW3, CW4; Thermal Springs: TS1, TS2, TS3, TS4, TS5; Thermal Wells: TW1, TW2, TW3). These waters were evaluated by Q-mode analysis to identify water origins and tectonic (e.g. fault) control (e.g. spatial distribution) on their chemistry. For example, if springs are aligned along a lineament, sampling locations are to form a cluster in dendogram as expected from fault-controlled springs. To test the grouping in the cluster analysis, we finally used factor analysis which is designed to reduce a large number of variables into fewer numbers of factors. Factor analysis collects variables into groups termed “factors”, which seem to behave similarly. Factor analysis is computed by “R-mode analysis” or “Q-mode analysis”. Therefore, it is anticipated that “factor analysis” and “cluster analysis” produce similar groups as they do in this study.

4. Results

4.1. Hydrochemical characteristics of the Kavak thermal waters

In this study, the Piper diagram was used to assess the main chemical composition of the Kavak thermal waters. In addition, thermal water samples from the Pamukkale geothermal field [e.g. an example of interactions mainly with metamorphic rocks, Alçiçek et al. 2019] and Varto geothermal field [e.g. an example of interactions mainly with volcanic rocks, Karaoğlu et al. 2019] were selected to compare with the Kavak thermal waters (Figure 8).

Figure 8.

Piper plot showing the chemical composition of thermal waters in the Kavak geothermal field [Bozdağ 2016]. Pamukkale thermals from Alçiçek et al. [2019] and Varto themals from Karaoğlu et al. [2019].

The KGF can be characterised by two types of thermal waters based on the slight difference observed in their hydrochemistry. The first type thermal waters including TS1, TS2 and TS3 have higher TDS (up to 3200 mg/l) and Cl/Br ratio (950–1215), and lower dissolved silica (22–23 mg/l) and Br content (234–268 mg/l) relative to the second type water including TS4, TS5, TW1, TW2 and TW3 (TDS: 1850–2490 mg/l; Cl/Br ratio: 300–465; dissolved silica: 39–44 mg/l and Br content: 408–614 mg/l). Li, Cs, Rb, Mn, Fe, Cu, and As are also slightly higher, but Cl, Na, K, Ca, Mg, B, Ba, Sr and Zn concentrations are lower in second type thermal water than those of the first type. In addition, outlet temperatures of these thermal waters discharging along two different faults in north and south are also slightly different (e.g. ∼22 °C for the first type of thermal waters and 25 to 46 °C for the second ones), and discharge along two different faults in north and south (see Figure 3).

Both type thermal waters in the KGF are chemically of Ca–Na–HCO3 type although they show slight difference in their hydrochemistry. Also, they have higher total dissolved solids (TDS) values (average of 2566 mg/l) than those of cold waters (average of 359 mg/l) which probably offer long flow pathways, long residence times and intensive water–rock interactions [Bozdağ 2016; Erbaş and Bozdağ 2022]. Dissolution of the Paleozoic carbonates forming the main reservoir of thermal waters in the KGF is the reason for the high Ca2+ (283.1–549.4 mg/l) and HCO3 (1445–2687 mg/l) concentrations. Although the source of Mg2+ in Kavak thermal waters is mainly the dissolution of dolomite and dolomitic limestone in the main aquifer in KGF, it is not the dominant cation. The reason for the low Mg values (average 72 mg/l) in thermal waters can be explained by the Mg depletion in the formation of phyllosilicates (e.g. clay minerals, chlorite and micas) in the Paleozoic metamorphics forming the secondary aquifer and/or ion exchange process. Na+ (232–387 mg/l) is the second common cation in the thermal waters and K+ (average 86 mg/l) is reached considerable values. Accordingly, the high Na+ and K+ in the Kavak thermal waters may be assumed to be derived from the decomposition mostly of intrusions of the ErAVC and regional metacrystalline rocks, progressive reactions with feldspar in the Paleozoic metamorphics and ion exchange process. As also seen in the diagram, the similarity of the main ion concentrations of the Kavak thermal waters with that of the Varto thermal waters located in the south of the Varto caldera “alkaline and calc-alkaline volcanic rocks” [Buket and Temel 1998] may be important for a better understanding of the relationship between regional geology and water–rock interaction.

The differences of major cations between Pamukkale and Varto thermal waters and the Kavak thermal waters can also be seen in the Schoeller diagram (Figure 9). The Kavak thermal waters have very high cation values compared to cold waters as a result of relatively deep and long-term circulation and contact with mainly volcanic rocks and partially crystalline basement or reservoir rocks. Additionally, Giggenbach’s triangular diagram is often used for water classification in geochemical investigations performed in convergent-plate settings. In the diagram (Figure 10), the Kavak thermal water samples are located in the area of peripheral waters and this type water generally occurs at relatively shallow depths at some distance from the geothermal field.

Figure 9.

Semi-logarithmic Schoeller diagram, comparison of the concentrations of the Kavak geothermal (red line), the Pamukkale geothermal (darkblue line), and Varto (green line) samples. Pamukkale thermals from Alçiçek et al. [2019] and Varto themals from Karaoğlu et al. [2019].

Figure 10.

A triangular diagram of major anions has been suggested by Giggenbach [1988] for thermal waters, symbols as Figure 8. Pamukkale thermals from Alçiçek et al. [2019] and Varto themals from Karaoğlu et al. [2019].

It is known that silicate weathering is effective in terms of SiO2, Ca2+, Mg2+, Na+, and K+ concentrations in groundwater. On the other hand, carbonate dissolution is effective in terms of Ca2+, Mg2+, and HCO3 concentrations. Accordingly, it is seen in the Ca/Na-Mg/Na diagram that the Kavak thermal waters fall close to the silicate weathering area and to the Varto thermal waters (Figure 11). The (Na++K+)∕(Na++K++Ca2+) ratios are higher than 0.45 indicating that the chemistry of the Kavak thermal waters is mainly controlled by the water–rock interaction. (Na++K+ − Cl)∕(Na++K+ − Cl + Ca2+) ratio < 0.20, and Na+∕(Na++Cl) ratio > 0.50 indicate that the source of Na+ ion is plagioclase decomposition and ion exchange process [Hounslow 1995]. Additionally, SiO2 content, which is generally between 1–30 mg/l in natural waters, is higher than 30 mg/l in the Kavak thermal waters, indicating silicate weathering. The ion exchange process can be indicated by chloro-alkaline indices (CAI-1 and CAI-2) [Schoeller 1965], and negative CAI indices reflect the presence of ion exchange in the system. The CAIs are calculated by the following equations (all ionic concentrations are expressed by meq/L);

CAI-1 = [ Cl ( Na + + K + ) ] Cl (1)
CAI-2 = [ Cl ( Na + + K + ) ] ( HCO 3 + SO 4 2 + NO 3 ) .   (2)
The negative of CAI-1 and CAI-2 indicates that Ca2+ or Mg2+ ions in the groundwater have been replaced by Na+ ion in a hydrodynamic medium (Figure 12). However, Şahinci [1991] specified that generally, the chloro-alkaline indices of water emerging from igneous and metamorphic rocks are also negative because the contribution of the alkali ions releasing from the decomposition of silicates is much more than the chloride ion. Accordingly, the hydrolysis of alkali feldspars in the Paleozoic metamorphics may be main process for especially high Na+ ion in the Kavak thermal waters. However, the low Na2O contents (<1 wt%) and the low modal feldspar proportions of the Paleozoic metamorphites [Karadağ 2014] indicate that neither ion exchange nor hydrolysis of alkali feldspars could not provide high sodium ions to thermal waters.

Figure 11.

Ca2+/Na+ vs. Mg2+/Na+ diagram showing end-members; Evaporite dissolution, Silicate weathering, and Carbonate dissolution, symbols as Figure 8. Pamukkale thermals from Alçiçek et al. [2019] and Varto themals from Karaoğlu et al. [2019].

Figure 12.

Chloro-alkaline indices (CAI) of the Kavak geothermal waters, symbols as Figure 8.

4.2. Data analysis and geochemical origin of the Kavak thermal waters

To better understand the details of the water–rock interaction and to classify sampling locations into similar groups, and to identify two types of thermal water contributing to these groups in the KGF, data analysis methods were performed by using the published hydrochemistry data from Bozdağ [2016]. Firstly, cluster analysis was performed by the method of shortest distance on the Kavak thermal waters to reveal the relationship between thermal water and rocks/minerals based on the R-mode analysis and to research spatial distribution of two types of thermal waters (e.g. their relationship with the faults) based on the Q-mode analysis. The clustering analysis and factor weights were calculated for variables from thermal water taken from 8 different locations (TS1, TS2, TS3, TS4, TS5; TW1, TW2, TW3) in the Kavak geothermal field by using GEOstats [Gündüz and Asan 2022].

The hierarchical cluster analysis with R-mode is based on correlation coefficients (Table 1), producing dendrogram (Figure 13). The dendrogram consists of two main groups (Group I: K, Ba, Sr, Ca, Na, Mg, Al, P; Group II: Cs, Rb, Li, Fe, SiO2) and 4 sub-groups. These groups may be attributed to rock-forming minerals especially in volcanic rocks represented by the ErAVC and metamorphic rocks from the basement. For example, Subgroup (1) Sr, Ba and K are important components for biotite and K-feldspar: sanidine is modally rare in the ErAVC rocks, but they have high normative orthoclase (Or; KAlSi3O8 end-member) content (Table 2). K may also be associated with sericite, and other phyllosilicate minerals in the crystalline basement rocks (e.g. phyllite, schist, etc.). On the other hand, Sr is an important element for plagioclases as a result of Ca-Sr substitution. Subgroup (2) Na, Ca and Mg may be interpreted to plagioclase and mafic silicates (e.g. amphibole, pyroxene) in volcanic rocks or calcite-dolomite in meta-carbonate rocks. Also, Al and P can be accepted in the same group as plagioclase and apatite minerals, respectively. Subroups (3) and (4) Li, Rb, Cs, Fe, and Si are important in the crystal structure of amphiboles and biotites.

Figure 13.

R-mode cluster analysis (Dendrogram) of the Kavak thermal waters (single-linkage method).

Table 1.

R-mode correlation matrix of the Kavak thermal waters (for April 2015)

SiO2 Ca Mg Na K Fe Li P Rb Sr Ba Al Cs
SiO2 1                        
Ca −0.86 1                      
Mg −0.86 0.96 1                    
Na −0.76 0.96 0.98 1                  
K −0.84 0.82 0.94 0.86 1                
Fe 0.72 −0.28 −0.31 −0.11 −0.44 1              
Li 0.79 −0.67 −0.59 −0.47 −0.56 0.72 1            
P −0.62 0.74 0.78 0.81 0.79 −0.07 −0.26 1          
Rb 0.70 −0.48 −0.43 −0.29 −0.47 0.78 0.95 −0.18 1        
Sr −0.94 0.92 0.96 0.90 0.94 −0.48 −0.63 0.84 −0.52 1      
Ba −0.96 0.93 0.95 0.89 0.90 −0.50 −0.67 0.79 −0.55 0.99 1    
Al 0.12 −0.18 −0.09 −0.16 0.12 −0.12 −0.25 −0.01 −0.36 −0.08 −0.14 1  
Cs 0.40 −0.19 −0.08 0.05 −0.06 0.69 0.81 0.34 0.80 −0.12 −0.17 −0.23 1

Table 2.

Results of CIPW norm calculation (%) of the ErAVC [Asan 2017; Asan et al. 2023]

  Min Max Mean
Quartz 0.99 27.75 13.27
Anorthite 13.02 28.52 21.21
Albite 20.90 34.19 27.80
Orthoclase 9.87 24.29 18.26
Pyroxene 6.88 30.43 16.24
Ilmenite 0.63 1.98 1.29
Magnetite 0.64 2.45 1.31

The Cluster analysis based on the Q-mode forms two main groups (Group I: TS1, TS2, TS3; Group II: TW1, TW2, TW3, TS4, TS5) and three subgroups (Figure 14, Table 3). The Group I, including subgroup 1, is represented by the first type of thermal waters (e.g. higher TDS and Cl/Br ratio, and lower dissolved silica and Br content with discharge temperature of 22 °C) discharged along the southern fault whereas the Group II includes the second type thermal waters (e.g. lower TDS and Cl/Br ratio, and higher dissolved silica and Br content with discharge temperature of 25 to 46 °C) discharged along the northern fault, which is close to the ErAVC (see Figures 23). In the Group II, well samples (e.g. Subgroup 2; TW1, TW2, TW3) and spring samples (e.g. Subgroup 3; TS4, TS5) also form two separate sub-groups.

Figure 14.

Q-mode cluster analysis (Dendrogram) of the Kavak thermal waters (single-linkage method).

Table 3.

Q-mode correlation matrix of the Kavak thermal waters (for April 2015)

  TS1 TS2 TS3 TS4 TS5 TW1 TW2 TW3
TS1 1
TS2 0.997 1
TS3 0.999 0.999 1
TS4 0.763 0.808 0.785 1
TS5 0.801 0.844 0.822 0.995 1
TW1 0.813 0.855 0.830 0.973 0.986 1
TW2 0.849 0.884 0.865 0.975 0.987 0.990 1
TW3 0.845 0.881 0.861 0.977 0.989 0.993 0.999 1

The grouping of elements resulted from R-mode and Q-mode factor analysis is similar to hierarchical clustering analysis. In R-mode factor analysis (Figure 15, Table 4), K, Sr, Ba, Ca, Na, Mg and P form a group similar to those of the Group I in the R-mode cluster analysis, whereas SiO2, Li, Rb, Fe and Cs form another group similar to the Group II in the R-mode cluster analysis. On the other hand, in Q-mode factor analysis (Figure 16, Table 5), TS1 TS2, TS3 samples from the first type thermal waters and TW1, TW2, TW3, TS4, TS5 samples from the second type thermal waters form two separate groups as in Q-mode cluster analysis.

Figure 15.

R-mode factor weights of the Kavak thermal waters on the diagram.

Table 4.

R-mode factor analysis and eigenvalue results of the Kavak thermal waters for 13 variables

Component  Factor 1 Factor 2 Factor 3 Factor 4 Factor 5
SiO2 −0.965 0.116 0.181 0.146 0.006
Ca 0.934 0.202 −0.107 0.270 −0.037
Mg 0.946 0.268 0.025 0.109 0.142
Na 0.872 0.420 0.017 0.234 0.090
K 0.917 0.177 0.223 −0.149 0.219
Fe −0.561 0.661 0.168 0.464 −0.064
Li −0.780 0.583 −0.006 −0.205 0.089
P 0.735 0.529 0.262 −0.147 −0.287
Rb −0.668 0.697 −0.093 −0.041 0.165
Sr 0.977 0.183 0.009 −0.105 −0.004
Ba 0.980 0.144 −0.083 −0.073 −0.032
Al −0.004 −0.427 0.900 0.025 0.056
Cs −0.311 0.896 0.180 −0.237 −0.058
Eigenvalue 8.22 2.92 1.05 0.54 0.21
Cumulative eigenvalue 8.22 11.13 12.18 12.72 12.93
Variability (%) 63.22 22.43 8.08 4.13 1.59
Cumulative variability (%) 63.22 85.65 93.73 97.85 99.45

Figure 16.

Q-mode factor weights of the Kavak thermal waters on the diagram.

Table 5.

Q-mode factor analysis and eigenvalue results of the Kavak thermal waters for 8 variables

Component Factor 1 Factor 2 Factor 3 Factor 4 Factor 5
TS1 0.920 0.392 0.004 −0.004 −0.001
TS2 0.947 0.322 0.007 0.020 0.004
TS3 0.933 0.360 0.019 −0.004 −0.001
TS4 0.951 −0.288 0.110 −0.004 0.021
TS5 0.970 −0.233 0.061 0.026 −0.027
TW1 0.972 −0.205 −0.096 0.059 0.009
TW2 0.986 −0.148 −0.049 −0.056 −0.002
TW3 0.985 −0.158 −0.051 −0.035 −0.002
Eigenvalue 7.35 0.61 0.03 0.01 0.01
Cumulative eigenvalue 7.35 7.96 7.99 8.00 8.00
Variability (%) 91.82 7.67 0.38 0.11 0.02
Cumulative variability (%) 91.82 99.49 99.87 99.98 100.00

5. Water–rock interaction and conceptual model of the Kavak geothermal field

Whole rock- and mineral chemistry of volcanic rocks are key features of better understanding the water–rock interaction in geothermal fields such as KGF where volcanic rocks act as a heat source. The relative cation chemistry of waters is compared with rocks in Ca–Mg–Na+K triangular plot to research water–rock interaction in the KGF. In the Ca–Mg–Na+K triangular diagram, both the ErAVC volcanic rocks and the Kavak thermal waters plotting close to the Na+K field suggest that high Na–K concentrations of the Kavak thermal waters can be related to the ErAVC volcanic rocks and their intrusions. However, the cold waters show a high tendency towards Ca and Ca+Mg around the KGF (Figures 17 and 18).

Figure 17.

The relative cation chemistry of the Kavak thermal and cold water [Bozdağ 2016] and comparison with the ErAVC rocks; pointed [Keller et al. 1977], gray dashed-line [Temel et al. 1998], pink-colored [Asan 2017; Asan et al. 2023] in the (Na+K)–Ca–Mg triangular plot. Altered volcanics are from Karakaya [2009]. Metamorphic rocks (metapelites, metapsammites, and metamarls) are from Karadağ [2014].

Figure 18.

The relative cation chemistry of the Kavak thermal water and cold water [Bozdağ 2016] and comparison with the ErAVC rocks; pointed [Keller et al. 1977], gray dashed-line [Temel et al. 1998], pink-colored [Asan 2017; Asan et al. 2023] in the Na–K–(Ca+Mg) triangular plot. Altered volcanics are from Karakaya [2009]. Metamorphic rocks (metapelites, metapsammites, and metamarls) are from Karadağ [2014].

The Kavak geothermal field is recharged mainly by infiltrating of meteoric waters, which go down through the faults and fracture/cracks systems acting as hydrothermal conduits [Bozdağ 2016]. The infiltrated waters are heated by the high thermal gradient of the Neogene (Miocene–Pliocene) Erenlerdağ–Alacadağ volcanic rocks in the deep, and then the density of these waters heated at the deep decreases as a result of increasing temperature and pressure, and so the thermal waters move up to the surface and/or shallower depths of the hydrothermal pathways (Figure 19). The low tritium values [<1 TU; Bozdağ 2016] of the Kavak thermal waters also suggest deep circulation and longer residence times. The following geological model of the region (Figure 19) is also supported by the δ13C values in the Kavak travertines and thermal waters because high δ13C is revealed as a result of reaching very deep the faults fractures in the region [Karaisaoğlu and Orhan 2018].

Figure 19.

Unscaled schematic illustration of the Erenlerdağ–Alacadağ volcanic complex (ErAVC) and Kavak Travertines modified from Bozdağ [2016]. The geochemical framework of the geothermal systems that are related to subduction zone volcanism [Henley and Ellis 1983].

Fractured and karstified Palaeozoic carbonates constitute the main reservoir rocks for the Kavak geothermal waters and, pyroclastic rocks of the ErAVC can be accepted as the secondary aquifer. Fine-grained metasandstone, metasiltstone and phyllite in the Paleozoic metamorphites and the clayey and silty levels of the Neogene units act as cap rocks for the Kavak geothermal system. The chemical composition of the Kavak thermal waters is mainly controlled water–rock interactions including dissolution of the Paleozoic carbonates, and the decomposition mostly of intrusions of the ErAVC.

Based on the chemical and statistical data, it has been identified that there are two distinct types of the thermal waters in the KGF. The reason for the occurrence of two different types of thermal waters is thought to be that the fault in the south, where the first type of thermal waters (TS1, TS2 and TS3) are located, acts as a barrier or semi-barrier. This is because faults act as a barrier and/or channel for flow, depending on their composition and the stage of fault evolution. In other words, a fault may act as a conduit immediately after deformation and it can later act as a barrier due to the precipitation of minerals [Caine et al. 1996]. Therefore, the partial interruption of the relationship between the first and second types of the thermal waters due to the southern fault resulted in slight differences, particularly in the minor ion contents and TDS values. In addition, the mixing of the first type of thermal water with colder water caused a greater decrease in temperature compared to the second type of thermal water.

6. Conclusions

In this study, the water–rock interaction mechanism was investigated based on the multivariate statistical analysis (e.g. factor analysis and clustering analysis) using the published water and whole rock geochemical data from the Kavak Geothermal Field (Konya, Turkey). The Kavak Geothermal Field (KGF) is controlled by the Seydişehir fault that is a part of a graben-like structure in an extentional basin. The KGF overlies a metamorphic basement, yet is located near the Erenlerdağ–Alacadağ volcanic complex (ErAVC). The metamorphic basement is represented by low-degree metaclastic and metacarbonate rocks whereas the volcanic rocks are composed of orogenic-type high-K calc-alkaline andesite and rhyodacite. The Kavak thermal waters are meteoric in origin and peripheral waters (Ca–Na–HCO3), and are characterised by high K+ and Na+ cations, and low pH (between 6.4–6.9 pH) values relative to the cold waters. Two types of thermal waters were identified in the KGF based on the slight difference in their hydrochemistry. The first type thermal water has higher TDS and Cl/Br ratio, and lower dissolved silica and Br content relative to the second type water. These thermal waters lying along two different faults allso have a slightly different discharge temperature. Q-mode cluster and factor analysis confirmed that first type and second type thermal waters formed two separate statistical groups, which is compatible with the observed slight geochemical differences between two types of thermal waters. The reason for the occurrence of two different types of thermal waters is thought to be that the fault in the south, where the first type of thermal waters (TS1, TS2 and TS3) are located, acts as a barrier or semi-barrier. Mixing of the first type of thermal water with colder water also caused a greater decrease in temperature compared to the second type of thermal water.

When compared, hydrochemical properties of the Kavak thermal waters are similar to the ErAVC relative to the metamorphic basement rocks. R-mode cluster and factor analysis performed in this study suggests that the chemistry of the KGF waters was mainly controlled by the composition of the ErAVC rather than those of the basement metamorphic rocks as a result of water–rock interaction based on the quantitative of cation exchanges. In this regard, the Kavak geothermal system was revised to a new conceptual model considering the influence of the heat source on the chemistry of the regional thermal waters as a result of water–rock interaction.

Declaration of interests

The authors do not work for, advise, own shares in, or receive funds from any organization that could benefit from this article, and have declared no affiliations other than their research organizations.

Acknowledgments

This article is based on the first author’s study at University of Poitiers. The authors thank to the editor Dr. François Chabaux and an anonymous reviewer for their critical and constructive comments that improved the quality of this paper. The first author also thanks to Campus France for scholarships to foreign students in France.


References

[Aksoy, 2019] R. Aksoy Extensional neotectonic regime in West–Southwest Konya. Central Anatolia, Turkey, Int. Geol. Rev., Volume 61 (2019) no. 14, pp. 1803-1821 | DOI

[Alçiçek et al., 2019] H. Alçiçek; A. Bülbül; I. Yavuzer; M. C. Alçiçek Origin and evolution of the thermal waters from the Pamukkale Geothermal Field (Denizli Basin, SW Anatolia, Turkey): insights from hydrogeochemistry and geothermometry, J. Volcanol. Geotherm. Res., Volume 372 (2019), pp. 48-70 | DOI

[Asan et al., 2021] K. Asan; H. Kurt; M. Gündüz; G. Gençoğlu Korkmaz; M. Ganerød Geology, geochronology and geochemistry of the miocene sulutas volcanic complex, Konya-Central Anatolia: genesis of orogenic and anorogenic rock associations in an extensional geodynamic setting, Int. Geol. Rev., Volume 63 (2021) no. 2, pp. 161-192 | DOI

[Asan et al., 2023] K. Asan; H. Kurt; M. Gündüz; G. Gençoğlu Korkmaz Geochemistry and Petrology of the Unimodal Versus Bimodal Neogene Volcanism in the Konya Volcanic Field, 2023 Central Anatolia,Turkey. (under review)

[Asan, 2017] K. Asan Çarpışma Sonrası Ortamda Ultra-potasikten Kalk-alkalen Volkanizmaya Geçişin Jeokimyasal ve Sr–Nd–Pb Izotopik Özellikleri, Konya-Türkiye, Tübitak, Konya, 2017 (Proje No: 113Y415)

[Bayram, 1992] A. F. Bayram Seydişehir ve Kaşaklı sıcaksu kaynaklarının hidrojeoloji incelenmesi, Selçuk University, Konya, 1992, 83 pages

[Besang et al., 1977] C. Besang; F. J. Eckhardt; W. Harre; H. Kreuzer; P. Müller Radiometrische Altersbestimmungen an neogenen Eruptivgesteinen der Turkei, Geol. Jahrb., Volume B25 (1977), pp. 3-36

[Bozdağ, 2016] A. Bozdağ Hydrogeochemical and isotopic characteristics of Kavak (Seydişehir-Konya) geothermal field, Turkey, J. Afr. Earth Sci., Volume 121 (2016), pp. 72-83 | DOI

[Buket and Temel, 1998] E. Buket; A. Temel Major-element, trace-element, and Sr-Nd isotopic geochemistry and genesis of Varto (Muş) volcanic rocks, Eastern Turkey, J. Volcanol. Geotherm. Res., Volume 85 (1998) no. 1–4, pp. 405-422 | DOI

[Burçak et al., 2002] M. Burçak; O. Duman; K. Bekar Seydişehir (Konya)-Kavakköy Jeotermal SahasıJeotermal-Jeofizik Etüt Raporu, Raport No: 10666, General Directorate of Mineral Research and Exploration (MTA), Ankara, 2002

[Caine et al., 1996] J. S. Caine; J. P. Evans; C. B. Forster Fault zone architecture and permeability structure, Geol. Soc. Am., Volume 24 (1996), pp. 1025-1028

[Davraz et al., 2022a] A. Davraz; M. T. Nalbantçılar; İ. Önden Hydrogeochemical characteristics and trace element of geothermal systems in Central Anatolia, Turkey, J. Afr. Earth Sci., Volume 195 (2022), 104666 | DOI

[Davraz et al., 2022b] A. Davraz; M. T. Nalbantçılar; S. Varol; İ. Önden Hydrogeochemistry and reservoir characterization of the Konya geothermal fields, Central Anatolia/Turkey, Geochem. Geophys. Geosyst., Volume 82 (2022) no. 2, 125867

[Dean and Monod, 1970] W. T. Dean; O. Monod The Lower Paleozoic stratigraphy and faunas of the Taurus Mountains near Beyşehir, Turkey. I. Stratigraphy, Bull. Br. Mus. (Nat. Hist.) Geol., Volume 19 (1970), pp. 411-426

[Erbaş and Bozdağ, 2022] H. A. Erbaş; A. Bozdağ Hydrogeochemical characteristics and evaluation of the geothermal fluids in the Gazlıgöl geothermal field (Afyonkarahisar), Western Anatolia, Turkey, Geothermics, Volume 105 (2022), 102543 | DOI

[Eren, 1996] Y. Eren Sille-Tatköy (Bozdağlar masifi-Konya) kuzeyinde Alpin öncesi bindirmeler, Türkiye Jeoloji Kurultayı Bülteni, Volume 11 (1996), pp. 163-169

[Giggenbach, 1988] W. F. Giggenbach Geothermal solute equilibria. Derivation of Na–K–Ca–Mg geoindicators, Geochim. Cosmochim. Acta, Volume 52 (1988) no. 12, pp. 2749-2765 | DOI

[Gill, 1981] J. B. Gill Orogenic Andesites and Plate Tectonics, Springer Verlag, Berlin, 1981, p. 401

[Göncüoğlu and Kozlu, 2000] M. C. Göncüoğlu; H. Kozlu Early Paleozoic evolution of the NW Gondwanaland: data from southern Turkey and surrounding regions, Gondwana Res., Volume 3 (2000) no. 3, pp. 315-323 | DOI

[Göncüoğlu et al., 2007] M. C. Göncüoğlu; Ş. Çapkınoğlu; Ş. Gürsu; P. Noble; N. Turhan; U. K. Tekin; C. Okuyucu; Y. Göncüoğlu The Mississippian in the Central and Eastern Taurides (Turkey): constraints on the tectonic setting of the Tauride-Anatolide Platform, Geol. Carpathica, Volume 58 (2007) no. 5, pp. 427-442

[Göçmez and Şen, 1998] G. Göçmez; O. Şen Kavakköy Seydişehir çevresinin jeoloji ve hidrojeoloji incelemesi, F.Ü. Jeoloji eğitiminin 20, yıl sempozyum bildirileri, Elazığ, 1998, pp. 607-613

[Göçmez et al., 2005] G. Göçmez; İ. Kara; M. T. Nalbantçılar; A. Güzel Konya’daki Sıcak ve Mineralli suların Hidrokimyasal Özellikleri, II, Ulusal Hidrolojide Izotop Teknikleri Sempozyumu, Izmir, 2005, pp. 121-132

[Gündüz and Asan, 2022] M. Gündüz; K. Asan GEOstats: an excel-based data analysis program applying basic principles of statistics for geological studies, Earth Sci. Inf., Volume 15 (2022), pp. 705-712 | DOI

[Gündüz, 2023] M. Gündüz The origin and eruption mechanism of the Kilistra ignimbrites deduced from chemostratigraphic, geochronologic, remote sensing, and GIS methods, Erenlerdağ-Alacadağ volcanic complex SW Konya-Central Anatolia, PhD Thesis, Muğla Sıtkı Koçman University, Muğla (2023) (194)

[Gürsu et al., 2003] S. Gürsu; H. Kozlu; M. C. Göncüoğlu; N. Turhan Orta Torosların Batı Kesimindeki Temel Kayaları ve Alt Paleozoyik Örtülerinin Korelasyonu, TPJD Bülteni, Volume C.15 (2003) no. s.2, pp. 129-153

[Hakyemez et al., 1992] Y. Hakyemez; E. Elibol; M. Umut; B. Bakırhan; İ. Kara; H. Dağıstan; T. Metin; N. Erdoğan Konya-Çumra-Akören dolayının jeolojisi, 1992 (Rapor No: 9449, MTA Ankara)

[Hao et al., 2020] C. Hao; Y. Huang; D. Ma; X. Fan Hydro-geochemistry evolution in Ordovician limestone water induced by mountainous coal mining: A case study from North China, J. Mt. Sci., Volume 17 (2020), pp. 614-623 | DOI

[Henley and Ellis, 1983] R. W. Henley; A. J. Ellis Geothermal systems ancient and modern: a geochemical review, Earth-Sci. Rev., Volume 19 (1983) no. 1, pp. 1-50 | DOI

[Hounslow, 1995] A. W. Hounslow Water Quality Data: Analysis and Interpretation, CRC Press, Boca Raton, 1995, p. 416

[Irvine and Baragar, 1971] T. N. Irvine; W. R. A. Baragar A guide to the chemical classification of the common volcanic rocks, Can. J. Earth Sci., Volume 8 (1971) no. 5, pp. 523-548 | DOI

[Karadağ, 2014] M. M. Karadağ Geochemistry, provenance and tectonic setting of the Late Cambrian-Early Ordovician Seydişehir Formation in the Çaltepe and Fele areas, SE Turkey, Geochemistry, Volume 74 (2014) no. 2, pp. 205-224 | DOI

[Karaisaoğlu and Orhan, 2018] S. Karaisaoğlu; H. Orhan Sedimentology and geochemistry of the Kavakköy Travertine (Konya, central Turkey), Carbonates Evaporites, Volume 33 (2018), pp. 783-800 | DOI

[Karakaya, 2009] N. Karakaya REE and HFS element behaviour in the alteration facies of the Erenler Dağı Volcanics (Konya, Turkey) and kaolinite occurrence, J. Geochem. Explor., Volume 101 (2009) no. 2, pp. 185-208 | DOI

[Karaoğlu et al., 2019] Ö. Karaoğlu; M. Bazargan; A. Baba; J. Browning Thermal fluid circulation around the Karliova triple junction: Geochemical features and volcano-tectonic implications (Eastern Turkey), Geothermics, Volume 81 (2019), pp. 168-184 | DOI

[Keegan-Treloar et al., 2022] R. Keegan-Treloar; D. J. Irvine; S. C. Solorzano-Rivas; A. D. Werner; E. W. Banks; M. J. Currell Fault-controlled springs: A review, Earth-Sci. Rev., Volume 230 (2022), 104058 | DOI

[Keller et al., 1977] J. Keller; D. Jung; K. Burgath; F. Wolff Geologie und Petrologie des neogenen Kalkalkali Vulkanismus von Konya (Erenler Dağ-Alaca Dağ-Massiv, Zentral-Anatolien), Geol. Jahrb., Volume B25 (1977), pp. 37-117

[Koç et al., 2012] A. Koç; N. Kaymakci; D. J. J. van Hinsbergen; K. F. Kuiper; R. L. M. Vissers Tectono-Sedimentary evolution and geochronology of the Middle Miocene Altınapa Basin, and implications for the Late Cenozoic uplift history of the Taurides, southern Turkey, Tectonophysics, Volume 532–535 (2012), pp. 134-155 | DOI

[Koç et al., 2017] A. Koç; N. Kaymakci; D. J. J. van Hinsbergen; K. F. Kuiper Miocene tectonic history of the Central Tauride intramontane basins, and the paleogeographic evolution of the Central Anatolian Plateau, Global Planet. Change, Volume 158 (2017), pp. 83-102 | DOI

[Koç et al., 2018] A. Koç; D. J. J. van Hinsbergen; C. G. Langereis Rotations of normal fault blocks quantify extension in the Central Tauride intramontane basins, SW Turkey, Tectonics, Volume 37 (2018) no. 8, pp. 2307-2327 | DOI

[Koçak and Zedef, 2016] K. Koçak; V. Zedef Geochemical Characteristics of the Lava Domes in Yatağan Village and Sağlık Town, From Erenlerdağı (Konya, Central Turkey) Volcanites, Acta Geobalcanica, Volume 2 (2016) no. 1, pp. 7-19 | DOI

[Le Bas et al., 1986] M. J. Le Bas; R. W. Le Maitre; A. Streckeisen; B. Zanettin A chemical classification of volcanic rocks based on total alkali-silica diagram, J. Petrol., Volume 27 (1986) no. 3, pp. 745-750

[Moix et al., 2008] P. Moix; L. Beccaletto; H. W. Kozur; C. Hochard; F. Rosselet; G. M. Stampfli A new classifcation of the Turkish terranes and sutures and its implication for the paleotectonic history of the region, Tectonophysics, Volume 451 (2008) no. 1–4, pp. 7-39 | DOI

[MTA, 2013] MTA Magmatic Rocks Map of Turkey, General Directorate of Mineral Research and Explorations, Ankara, 2013

[Okay and Tüysüz, 1999] A. I. Okay; O. Tüysüz Tethyan sutures of Northern Turkey, The Mediterranean Basins: Tertiary Extension within the Alpine Orogen (B. Durand; L. Jolivet; F. Hovarth; M. Séranne, eds.), Geological Society, London, 1999, pp. 475-515

[Poznanović-Spahić et al., 2023] M. Poznanović-Spahić; G. Marinković; D. Spahić; S. Sakan; I. Jovanić; M. Magazinović; N. Obradović Water-Rock interactions across volcanic aquifers of the lece andesite complex (Southern Serbia): geochemistry and environmental impact, Water, Volume 15 (2023), 3653 | DOI

[Robertson et al., 2013] A. H. F. Robertson; O. Parlak; T. Ustaömer Late Palaeozoic-early Cenozoic tectonic development of southern Turkey and easternmost Mediterranean region: evidence from the interrelations of continental and oceanic units, Geological Development of Anatolia and the Easternmost Mediterranean Region (O. Parlak; U. C. Ünlügenç, eds.), Geological Society, Special Publications, London, 2013, pp. 9-48

[Schoeller, 1965] H. Schoeller Qualitative evaluation of groundwater resources, Methods and Techniques of Groundwater Investigations and Development (H. Schoeller, ed.), The United Nations Educational, Scientific and Cultural Organization, Paris, 1965, pp. 54-83

[Temel et al., 1998] A. Temel; M. N. Gündoğdu; A. Gourgaud Petrological and geochemical characteristics of Cenozoic high-K calc-alkaline volcanism in Konya, Central Anatolia, Turkey, J. Volcanol. Geotherm. Res., Volume 85 (1998) no. 1–4, pp. 327-354 | DOI

[Teng, 2005] H. H. Teng Water–rock interactions, Water Encyclopedia (J. H. Lehr; J. Keeley, eds.), John Wiley & Sons, New Jersey, 2005, p. 846

[Turan, 2010] A. Turan Akören (Konya, Orta Toroslar) Çevresinin Jeolojik Özellikleri, S.Ü. Müh.-Mim. Fak. Derg., Volume 25 (2010) no. 4, pp. 17-36

[Turan, 2020] A. Turan Akkise-Yalıhüyük (Konya) Arasının Stratigrafisi-Stratigraphy of Between Akkise and Yalıhüyük (Konya), DEÜ Mühendislik Fakültesi Fen ve Mühendislik Dergisi, Volume 22 (2020) no. 64, pp. 369-382 | DOI

[Wohletz and Heiken, 1992] K. Wohletz; G. Heiken Volcanology and Geothermal Energy, University of California Press, Berkeley, 1992, 432 pages

[Yidana et al., 2012] S. M. Yidana; B. Banoeng-Yakubo; A. Aliou; T. M. Akabzaa Groundwater quality in some Voltaian and Birimian aquifers in Northern Ghana-application of mulitvariate statistical methods and geographic information systems, Hydrol. Sci. J., Volume 57 (2012) no. 6, pp. 1168-1183 | DOI

[Zhang et al., 2022] K. Zhang; X. Deng; J. Gao; S. Liu; F. Wang; J. Han Insight into the process and mechanism of water-rock interaction in underground coal mine reservoirs based on indoor static simulation experiments, ACS Omega, Volume 7 (2022) no. 41, pp. 36387-36402 | DOI

[Özgül, 1976] N. Özgül Torosların Bazı Temel Jeoloji Özellikleri, Türkiye Jeol. Kur. Bült, Volume 19 (1976), pp. 65-78

[Özkan and Söğüt, 1998] A. M. Özkan; A. R. Söğüt Dilekçi (Konya Batısı) Çevresindeki Neojen Çökellerinin stratigrafisi, Mühendislik Bilimleri Dergisi, Volume 15 (1998) no. 2–3, pp. 1131-1138

[Şahinci, 1991] A. Şahinci Doğal suların jeokimyası, Reform Matbaası, Izmir, 1991, 548 pages

[Şengör et al., 2019] A. M. C. Şengör; L. Nalan; S. Gürsel; Z. Cengiz; S. Taylan The phanerozoic palaeotectonics of Turkey. Part I: an inventory, Mediterr. Geosci. Rev., Volume 1 (2019), pp. 91-161 | DOI


Comments - Policy