Introduction
Epithermal systems generally occur some distance from a magmatic intrusion; they are mostly produced by intrusions located ~5–6 km below the surface. Epithermal systems can also occur in areas without volcanic activity or volcanic rocks, however (Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000). Known epithermal systems, defined by their characteristic textures and mineral assemblages, are intrusion-related hydrothermal systems below the paleowater table, ranging in depth from 50 to 1000 m, which formed mostly in subaerial volcanic arcs at upper crustal levels which are observed at temperatures of ~50 to 300°C, the maximum temperature limited by the vapor pressure of boiling water (Cooke & Simmons, Reference Cooke and Simmons2000; Hedenquist et al. Reference Hedenquist, Arribas and Gonzalez-Urien2000; Lindgren, Reference Lindgren1933; Sillitoe et al., Reference Sillitoe, Hannington and Thompson1996; Simon et al., Reference Simon, Frank, Pettke, Candela, Piccoli and Heinrich2005; Wade et al., Reference Wade, McAvaney and Gordon2014). Defined epithermal environments are mostly shallow and can host deposits of Au, Ag, and base metals that have also been exploited for a wide range of metals and minerals, including Hg, Sb, Pb, Zn, Cu, Mo, S, halloysite, kaolinite, alunite, and silica (Hedenquist et al. Reference Hedenquist, Arribas and Gonzalez-Urien2000). Volcanic-hydrothermal systems have distinctly different types; these are known as low-sulfidation and high-sulfidation and these reflect the two end-members of the sulfidation states deduced from the sulfide mineral assemblages (Skinner & Barton, Reference Skinner and Barton1973). These two end-member styles of deposit can be distinguished easily based on their hypogene and supergene alteration mineralogy. Kaolin is the most observed mineral group in these alteration processes. If primary kaolin deposits are formed by in situ weathering from granitic, siliceous volcanic, and pegmatitic rocks, however, they may have a supergene origin. If they are formed by hydrothermal activities, they may have hypogene deposits, depending on the origin of the hydrothermal solutions; in some cases, a combination of these two processes may occur (Murray & Keller, Reference Murray, Keller, Murray, Bundy and Harvey1993). According to Dill (Reference Dill2016), supergene kaolin mineralization is often homogeneous and is characterized by oxidized minerals, whereas hypogene kaolin typically has variegated mineralogy with reducing minerals.
The Biga Peninsula, northwest Turkey, is well known for its low- to high-sulfidation ore deposits in which many hydrothermal activities were widely developed (Ağdemir et al., Reference Ağdemir, Kιrιkoğlu, Lehmann and Tietze1994; Çiçek et al., Reference Çiçek, Oyman and Palmer2021; Ece & Schroeder, Reference Ece and Schroeder2007; Ece et al., Reference Ece, Schroeder, Smilley and Wampler2008; Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016; Oyman, Reference Oyman2010; Pirajno, Reference Pirajno1995; Yiğit, Reference Yiğit2012). Much epithermal mineralization including ore minerals, as well as industrial minerals used in the ceramic and cement industries, in different amounts and chemical compositions, was observed in the region. All these mineral associations are the products of acid-sulfate and alkaline fluids of the epithermal system widely developed in the region (Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016; Laçin et al., Reference Laçin, Aysal and Öngen2021; Oyman, Reference Oyman2019; Yiğit, Reference Yiğit2012).
In this region, which has many characteristic features of hydrothermal alteration processes, numerous epithermal mineralizations have occurred along active fault zones. Some of the deposits of the region such as the Sarıbeyli-Sığırlı, Bodurlar, Hamamtepe, Muratlar, and Tepeköy quarries in this site were studied and the mineralogical structures of fossil epithermal sites were revealed, their structural differences were characterized, and origins have been interpreted. All these studies focused firstly on acid-sulfate fluids related to advanced argillic alteration mineral assemblages (Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016) and also alkaline-chloride fluids related to silica sinters (Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022). Thus, attempts have been made to reveal the temporal and spatial changes that the region has undergone from the past to the nearly present in relation to hydrothermal solutions.
The Bahadırlı Kaolin Deposit (BKD), the Duman Kaolin Deposit (DKD), and the Çaltıkara Kaolin Deposit (CKD) are found throughout the Çan-Etili-Bayramiç fault zone in the Biga Peninsula and where kaolinization is intensely developed. These quarries are located approximately halfway between the Sarıbeyli-Sığırlı and Bodurlar kaolin quarries in the north and the Hamamtepe and the Muratlar silica quarries in the south, where siliceous and silicified rocks are found. Identifying the mineralogy, determining the lithological features, and interpreting the environments of formation, based on geochemical and isotopic data of these quarries, has helped to reveal the hydrothermal alteration processes in the region.
The present study aimed to reveal the temporal evolution of the alteration zones in the Etili fault zone, and thus make a detailed investigation of the kaolinite deposits produced from acidic solutions. For this purpose, lithological, mineralogical, and geochemical studies including stable isotopic analyses and age dating were carried out. In addition, biological activity with acid sulfate-rich fluids also was assumed to be effective in the alteration processes of these deposits and, also in support of this hypothesis, micro-morphological images and stable isotopic data were interpreted together.
Geological Characteristics of the Region
Western Turkey is a post-collisional magmatic province with various tectonic sub-units which are: the Rhodope Massif, Sakarya Continent, and Anatolid-Torid Block derived by amalgamation along well defined suture zones. The magmatic activity in the region started in the Biga peninsula with granitic intrusions in the Middle Eocene and ended with the production of similar volcanic and magmatic rocks up to the Late Miocene (Aldanmaz et al., Reference Aldanmaz, Pearce, Thirlwall and Mitchell2000; Altunkaynak & Dilek, Reference Altunkaynak and Dilek2006; Karacık & Yılmaz, Reference Karacιk and Yιlmaz1998). The stratigraphic characteristic of these magmatic rocks changes from calc-alkaline to alkaline and plutonism-dominated to volcanism-dominated type (Altunkaynak & Genç, Reference Altunkaynak and Genç2008; Karacık et al., Reference Karacιk, Yιlmaz, Pearce and Ece2008). The Evciler pluton is an important granodiorite located southeast of the study area. Oligocene–Early Miocene aged felsic-type granodiorite is characterized by high-K, calc-alkaline affinity (Aydın, Reference Aydιn2019; Ercan et al., Reference Ercan, Satir, Steinitz, Dora, Sarifakioglu, Adis, Walter and Yildirim1995; Genç, Reference Genç1998; Özdamar, Reference Özdamar2018; Siyako et al., Reference Siyako, Bürkan and Okay1989; Türkdönmez & Bozcu, Reference Türkdönmez and Bozcu2008; Yılmaz & Karacık, Reference Yιlmaz and Karacιk2001). 40Ar/39Ar radiometric age data from hornblende and biotite from previous studies were determined to be 25.5 ± 0.1–28.3 ± 0.1 Ma and 21.9 ± 0.1–31.1 ± 1.4 Ma, respectively (Altunkaynak et al., Reference Altunkaynak, Sunal, Aldanmaz, Genç, Dilek, Furnes and Yιldιz2012; Birkle & Satir, Reference Birkle, Satιr, Erler, Ercan, Bingöl and Örçen1995; Delaloye & Bingol, Reference Delaloye and Bingol2000; Okay & Satir, Reference Okay and Satir2000; Özdamar et al., Reference Özdamar, Roden, Zou, Billor, Hames, Georgiev and Dunkl2021a). 40Ar/39Ar radiometric dating of alunite obtained from the south of this region indicates ages of 5 ± 0.18 to 32.4 ± 1.2 Ma (Ercan et. al., Reference Ercan, Ece, Schroeder and Gülmez2022) and is in good agreement with the aforementioned sequence, because the inception of hydrothermal alteration in the region dates back to the Oligocene (Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016). Hydrothermal alteration in the region developed along the fault zones and one is the Çan-Etili-Bayramiç fault on which several kaolin deposits occur. The Çan-Etili-Bayramiç Fault Zone existed before the North Anatolian Fault Zone (NAF) and merged with the NAF in Plio-Quaternary times to become the western continuum of the southern branch (Kaymakci et al., Reference Kaymakci, Aldanmaz, Langereis, Spell, Gurer and Zanetti2007; Le Pichon et al., Reference Le Pichon, Imren, Rangin, Şengör and Siyako2014; Zattin et al., Reference Zattin, Okay and Cavazza2005). This tectonic event is the result of the Anatolian plate moving westward and all these tectonic deformations are described as an extensional tectonic regime and are revealed in the form of N–S extension and E–W compression in the Aegean region (Okay & Satir, Reference Okay and Satir2000). Three models were proposed for the extension mechanism in the region. The first, proposed by Dewey and Şengör (Reference Dewey and Şengör1979), is described as a westward escape of the Anatolian plate. The collision between the Arabian platform and Anatolia compressed eastern Anatolia and transferred the stress to the west because Eurasia behaves like a craton. This stress transfer caused extension at the Aegean Sea and western Anatolia. However, the collision between the Arabian platform and Anatolia was at ~15 Ma ago and the detection of Oligocene basins in Western Anatolia invalidated this model (Okay, Reference Okay2001). The second model is gravitational collapse proposed by Seyitoǧlu and Scott (Reference Seyitoğlu and Scott1996). The lithosphere was thickened up to 55 km because of the collision and tended to collapse under its own weight and induce the extension of the Aegean Sea and western Anatolia. The third model, which is the most widely accepted model today, is the extension triggered by the retreating Hellenic trench. It was suggested that the convergence velocity between Africa and Eurasia decreased, and compression in western Anatolia was diminished; this then triggered the slab retreat and started extension in the region (Jolivet & Faccenna, Reference Jolivet and Faccenna2000; Jolivet & Brun, Reference Jolivet and Brun2010; Jolivet et al., Reference Jolivet, Faccenna, Huet, Labrousse, Le Pourhiet, Lacombe and Philippon2013a, Reference Jolivet, Faccenna, Huet, Labrousse, Le Pourhiet, Lacombe and Driussib).
In the study area, extensively observed units are Oligocene–Middle Miocene Çan volcaniclastics in the site consisting of andesite, trachyandesite, basaltic andesite lavas, and rhyolitic tuff. Hydrothermal alteration of Çan volcanics has caused the formation of large kaolin deposits and also silica sinter precipitation in the region (Ercan et al., Reference Ercan, Satir, Steinitz, Dora, Sarifakioglu, Adis, Walter and Yildirim1995; Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022; Erenoğlu, Reference Erenoğlu2017; Türkdönmez & Bozcu, Reference Türkdönmez and Bozcu2008). Volcanic rocks were intercalated with lignite, conglomerate, sand, clay, and carbonate during the sedimentary process in the Miocene; Çan volcanics were covered by all of these volcano-sedimentary units. The kaolin deposits of Bahadırlı, Çaltıkara, and Duman all occur within the Çan volcanics in the Etili area on the NE–SW- and E–W-trending faults that constitute the Çan-Etili-Bayramiç Fault Zone (Fig. 1).
Fig. 1 Geologic map of the study area. The map in the upper left indicates the location of western Turkey and neighboring regions. The regions marked with orange indicate the approximate locations of the Duman, Bahadırlı, and Çaltıkara kaolin deposits (Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022; Türkdönmez & Bozcu, Reference Türkdönmez and Bozcu2012)
Geological and lithological properties of the kaolin deposits
As a result of field work in BKD, DKD, and CKD, argillization was determined to have occurred intensely in the tuffs. The tuffs are observed as block and ash flow type. Various colorations and different secondary mineral developments occurred with the effect of alteration in the tuff. These colorations are characteristic oscillatory zoning observed in the tuffs of the Çan volcanics (Erenoğlu, Reference Erenoğlu2017). The alteration effect is observed in all of the mineral groups in the tuffs, and is more intense in K-feldspars, plagioclase, and the matrix. On the other hand, silicification-related quartz is another typically detected mineral group.
The Bahadırlı kaolin deposit (BKD) is located within the Oligocene Çan volcanics that are east of the study area. This quarry covers ~20,000 m2 (200 m × 100 m) and extends parallel to a fault in the E–W direction which is a branch of the Çan-Etili-Bayramiç fault zone and kaolinization occurred around this fault. The layers in different colors are quite evident at BKD (Fig. 2a, b). The red silicified rocks and less altered maroon iron oxide-rich tuff are observed at the base levels, and beige kaolinite increases at the middle levels of the deposits. Incarnadine, reddish-pink, greenish-gray colored iron oxide- and silica-rich thin bands are interstratified within kaolin-rich layers. The BKD begins with the purple firing zone at the east end and ends with the white-yellow and stockwork-textured opals in the west. The silicified tuff layer covers mostly the upper levels of the deposit. The tuffs affected by silicification are hard and compact, while the tuffs in the deposit unaffected by silicification are porous and brittle (Fig. 2a, b).
Fig. 2 Photographs of the kaolin deposits of the study area: a SW-NE side of the DKD containing coal, silt, and clay layers between the altered tuffs layers, b NE-SW side of the DKD containing clay and altered tuffs, c,d siliceous and iron-rich altered tuff, altered tuffs, and kaolinite layers of the BKD images, e image of CKD containing epithermal vein and iron-rich kaolin and altered tuffs
The Duman kaolin deposit (DKD), with an approximate size of 87,500 m2 (500 m × 175 m), is located within the Çan volcanites. The DKD extends in a NE–SW direction and is located on the same Çan-Etili-Bayramiç fault. In this quarry, kaolinization occurred on two different slopes which are related to faulting, these are the hanging wall block and the footwall block. Altered beige tuffs of the Çan volcanics are found in the upper zones of the hanging wall side, and white-beige colored kaolin is located at the lower levels. The footwall side starts with white altered kaolins at the lower levels and continues with silty, sandy, coaly units of different thicknesses containing lapilli tuff and conglomeratic layers belonging to the Çan Formation. These sedimentary layers start with coal-bearing clay at the base and reach up to 50 cm thick. The overlying coal layer, up to 3 m thick, is intercalated with the upper parts of red, yellow, brown, and beige sediments and volcanics. Slump, fold, and fracture structures were formed from the effect of tectonism in the sedimentary layers (Fig. 2c, d). The upper levels of the footwall are composed of altered tuffs.
The Çaltıkara Kaolin Deposits (CKD), located in the Çan Volcanics to the south of the study area, cover ~35,000 m2 (200 m × 175 m). The deposit extends NE–SW and is a fault-controlled deposit. This fault, which is a branch of the NE–SW trending Çan-Etili fault, is a right-lateral strike-slip normal fault, and kaolinization developed in parallel with this fault. White and beige massive kaolinization is common in the quarry, and an epithermal vein intrusion rich in ore minerals, including pyrite, galena, and bornite, is found in the upper levels of the quarry. In addition, reddish and purplish iron oxide minerals are observed in a wide area developed along the fault zone (Fig. 2e).
Methods of Study
Geological maps at 1/25,000 were made, based on detailed geological observations and representative sampling at open quarries for analysis. The petrological properties of samples collected from kaolins and basement rocks were determined in thin section using a Leica DM4500P polarizing microscope at the Istanbul Technical University (ITU) in İstanbul (Turkey). X-ray diffractometry (XRD) analyses were performed on 15 samples from BKD, 31 samples from DKD, and 12 samples from CKD; this was used to determine the mineralogical compositions. Analysis was carried out using a Bruker D8 Advance model instrument with CuKα radiation and using a LynxEye detector at Istanbul Technical University (ITU) in İstanbul, Turkey. XRD data were collected using CuKα radiation (40 mA, 40 kV), scanning range of 2–72º2θ, and step-size of 0.01940º for bulk random mounts, and 2–42º2θ, and step size of 0.005º2θ for oriented mounts, 0.1 s/step, 0.2 mm divergence slit, Ni filter, and a LynxEye1 solid-state detector. The micro-morphological properties of the kaolins and other minerals were detected using a ZEISS GeminiSEM 500 model Field Emission Scanning Electron Microscopes (FE-SEM) at 15 kV accelerating voltage, 5–15 mA current, and 10–20 s counting time for each element Science and Technology Research and Application Center (BİTAM) at the Necmettin Erbakan University in Konya, Turkey. An SPI-MODULE® sputter coater was used for Au coating under operating conditions of 2 mbar, 15 mA, and 50 s. This provides a coating which is ~150 Å thick.
Major and trace element analyses were performed using Spectro Ciros Vision inductively coupled plasma-optical emission spectrometry (ICP-OES) for Ba, Sc, Mo, Cu, Pb, Zn, Ni, and As at the ACME Laboratory, Vancouver, Canada. The detection limits for major elements, trace elements, and rare-earth elements (REE) were 0.01 to 0.1 wt.%, 0.1 to 5 ppm, and 0.01 to 0.5 ppm, respectively. Some trace and REE analyses were performed using a PerkinElmer Elan 6100 ICP-MS instrument at the ITU-Geochemistry Research Laboratory (JAL) in İstanbul, Turkey. The detection limits for REE analyses were 0.01 ppm. δ34S isotopic analyses were performed on three samples from the BKD, on seven samples from DKD, and on six samples from the CKD, for which purities of alunite were determined by XRD. In addition, one galena and one pyrite sample were taken from the copper mine for δ34S isotopic analysis and the isotopic contents were analyzed; SO2 produced by heating alunite at 1050°C was analyzed using a Finnigan MAT 252 gas-source mass spectrometer at the University of Georgia (UGA) in Georgia, USA. Canyon Diablo Troilite (CDT) was used as the standard and the accuracy is ~ ± 0.2‰. Separation of kaolins (< 2 μm) for isotope analysis: firstly, the kaolin-rich samples were milled for XRD analysis to reduce the particle size to <10–20 mm, and samples were dispersed in distilled water by ultrasonic vibration. Then, the classification of size fractionation of the clays (<2 mm) was performed using a centrifuge, and clay paste homogenized on the smear glass. Finally, oriented clay bases were prepared by air drying the glass slides (Moore & Reynolds, Reference Moore and Reynolds1989). The purity of kaolin samples was determined by XRD; δ18O and δD isotopic analysis were performed on three kaolin samples from BKD, on ten from DKD, and on three from CKD. The δ18O and δD values of kaolinite were determined using the Finnigan MAT 252 or Finnigan Delta E spectrometers at the same university in Georgia, USA. The best available three alunite samples were identified using XRD analysis and 40Ar/39Ar radiometric age analysis was performed by incremental heating at the University of Michigan, Michigan, USA. The incremental heating method utilizes a double vacuum resistance furnace similar to the design of Staudacher et al. (Reference Staudacher, Jessberger, Dorfinger and Kiko1978). These analytical methods and data-treatment procedures are included in detail in the studies by Justet and Spell (Reference Justet and Spell2001) and Spell and McDougall (Reference Spell and McDougall2003).
Results
Mineralogical and Petrographic Studies
In the present study, mineralogical and petrographic investigations were performed to reveal the compositions and textural properties of the Oligo-Miocene lavas of the Çan volcanics, where the intense alteration occurred. These volcanic units are mostly of andesite, basaltic-andesite, and basaltic-trachyandesite types (Türkdönmez & Bozcu, Reference Türkdönmez and Bozcu2008). Plagioclase, K-feldspar, amphibole, biotite, and rarely clinopyroxene crystals with different ratios were detected in these rocks (Fig. 3). The tuffs are observed as block- and ash-flow type and kaolinization is most evident in the tuffs. K-feldspar, plagioclase, and the matrix have been mostly kaolinized while biotite and amphiboles were Fe-oxidized as seen in the photomicrographs (Fig. 3).
Fig. 3 Optical microscopy images of Çan volcanics: a optical photomicrograph (crossed-polarized light) showing twin planes of plagioclase and slightly devitrified (altered) hornblende within the volcanic glass of the andesite; b slightly altered K-feldspars and plagioclases, Fe-oxidized hornblende within the volcanic glass of the andesite; c,d kaolinized plagioclase and K-feldspar phenocryst within the altered volcanic groundmass of the trachyte, the opaque crystals may be pyrite that have been subjected to partial oxidation; e,f kaolinized minerals which have lost their optical properties due to alteration within the tuff (f – plane-polarized light). The light gray patches are K-feldspar and plagioclase (note the local presence of zonation of the Fe-oxide pigment in the photomicrograph)
X-ray diffraction of the BKD samples detected kaolin, quartz, smectite, alunite, and lesser amounts of plagioclase and K-feldspar (Fig. 4a and Table 1). The mineral association representing the lower levels of BKD is kaolinite + quartz ± alunite ± plagioclase ± K-feldspar (Table 1; BH-1 to BH-9-d). However, at middle levels, the quartz content decreases, while smectite increases (Table 1; BH-10-A to BH-10-F). While the kaolin and alunite contents decrease, the amount of quartz increases in the upper levels (Table 1; BH-11 to BH-15-u). Two types of kaolinite were observed during the FE-SEM studies of the kaolinite-rich levels: (1) equidimensional, euhedral with < 6–7 µm kaolinite, and (2) small-grained, anhedral with dimensions of 1–3 µm. The first type is an equidimensional oriented kaolinite form which occurs in partially vermiform booklet stacks (Fig. 5a, b). The second type is randomly oriented crystal platelets (Fig. 3c). Well formed cubic-like alunite crystals which varied from 5 to 10 µm in size were observed. Dissolution pitting was detected in some of the alunite crystals (Fig. 5d). Halloysite has a tubular morphology with outer diameters of 1–2 µm (Fig. 5c, d).
Fig. 4 XRD patterns of the kaolin deposits: a BKD, b CKD, c hanging wall of the DKD, d lower levels of the DKD, e middle levels of the DKD, f upper levels of the DKD. Kln: kaolinite; Alu: alunite; Qz: quartz; Sme: smectite, Pl: plagioclase, Hal: halloysite, K-Fs: K-feldpar (abbreviations from Whitney & Evans, Reference Whitney and Evans2010)
Table 1 Mineralogical contents of samples from the Bahadırlı, Çaltıkara, and Duman kaolin deposits
Fig. 5 FE-SEM images of minerals from the BKD and CKD: a equidimensional and book-shaped kaolinite crystals; b vermiform and randomly oriented kaolinite crystals, BKD; c randomly scattered kaolinite and halloysite minerals, BKD; d well-formed cubic-like alunite crystals and halloysite tubes observed mostly at the surface, BKD; e,f randomly scattered kaolinite on the quartz, BKD; g,h spherical and euhedral quartz with iron-oxide minerals, CKD
In the XRD analysis of CKD samples, kaolin and quartz minerals were detected in all samples and this mineral couple was accompanied by small amounts of alunite. The amount of quartz increases nearer to the fault zone (Fig. 5b and Table 1). A commonly observed (by SEM) clay mineral, equidimensional and randomly oriented, is accompanied mostly by euhedral and anhedral quartz grains. Well crystallized <1–3 µm kaolinite plates look as though they have been sprinkled on the quartz (Fig. 5e, f). Spherical and euhedral quartz coexists, the latter crystals range up to 30 μm in diameter; the spherical quartz ranges between <1 and 30 μm (Fig. 5g, h). Microcrystalline anhedral iron-oxide minerals (<1 µm) were often randomly coated on the surface of the clay minerals and other quartz crystals (Fig. 5h).
Kaolin is the most frequently detected alteration mineral in XRD analysis of DKD samples. In hanging wall block samples, a mineral association of kaolinite + quartz ± alunite ± plagioclase ± K-feldspar ± smectite was detected (Fig. 4c and Table 1). Three different mineral associations were determined by XRD analysis of samples from the footwall block: (1) a kaolinite ± quartz ± plagioclase ± K-feldspar ± smectite association for samples from the lower levels of the footwall block (Fig. 4d); (2) a kaolinite + smectite ± quartz ± plagioclase ± alunite association from the middle levels of the footwall block (Fig. 4e); and (3) a kaolinite + plagioclase + smectite ± quartz ± K-feldspar ± gypsum association for samples from the upper levels of the footwall block (Fig. 4f and Table 1). FE-SEM images of kaolins from the section representing the hanging wall block showed them to be generally equidimensional and book-like but the edges and corners of the kaolin crystals were ragged (Fig. 6a, b). Book-like kaolinites were commonly observed among the altered surface of the coarse and euhedral feldspar minerals. In addition, weathering of halloysite and kaolinite minerals from feldspar minerals was typified quite well in FE-SEM images. Kaolinite to halloysite transformation was observed in the form of bending of kaolinite plates from the edges (Fig. 6c, d). Halloysite was more abundant in the FE-SEM images of samples from the footwall block of the quarry. Halloysite appeared with a hollow open-ended tubular morphology, having internal lumen diameters of 50–80 nm and outer diameters of 100–180 nm which were oriented randomly with varying lengths and thicknesses (Fig. 6e, f). Micro-morphological FE-SEM characterization of DKD samples revealed the ‘halloysitization’ of feldspars (Fig. 6). In addition, from FE-SEM studies, fossil remnants of the microbiota structures which developed intergrowths with clays on the altered mineral surfaces were detected in the altered samples of the hanging wall of the DKD. The structures of these microorganisms are rotiform and each branch was up to <0.5 µm thick (supplementary material Table S1 ESM_1PDF).
Fig. 6 FE-SEM images of altered minerals from the DKD: a vermiform (book-shaped) kaolinite crystals with euhedral altered feldspar minerals; b kaolinized feldspar minerals; c transformed halloysite mineral from feldspar minerals; d alteration of halloysite and kaolinite from feldspar and the bending of the kaolinite; e open-ended tubular halloysite crystals; f advanced altered feldspar with associated halloysite randomly scattered on the feldspar
Geochemistry of the Kaolin Deposits
Chemical analyses of 12 samples from two different rock types was performed: (1) 11 clay-rich samples from the alteration zone; and (2) one tuff sample, representing the host rocks. Major oxide concentrations of the two altered samples from BKD were: Al2O3 10.78–26.23 wt.%, SiO2 54.08–74.90 wt.%, Fe2O3 nearly 0.20 wt.% for the two samples; the Na2O, Ca2O, and K2O contents were small (Table 2). The maximum and minimum major oxide contents of the four samples of the CKD are 12.93–19.38 wt.% for Al2O3, 52.31–66.13 wt.% for SiO2, 0.28–6.79 wt.% for Fe2O3, 0.03–7.28 wt.% for Ca2O, 0.26–4.34 wt.% for Na2O, and 1.96–4.03 wt.% for K2O. In the CKD samples, the Fe2O3 contents of the samples taken from nearer the fault was large (M-26, M-27). Moreover, the CaO and Na2O contents of these samples are also large. The major oxide contents of five DKD samples ranged as follows: Al2O3 from 17.6 to 32.85 wt.%, SiO2 from 40.66 to 68.22 wt.%, Fe2O3 from 0.94 to 8.56 wt.%; the Na2O and Ca2O contents are small but the K2O contents ranged from 0.31 to 5.78 wt.% (Table 2).
Table 2 Chemical analyses of the BKD, CKD, DKD, and parent rocks (major elements in wt.%, trace and rare earth elements in ppm). B represents BKD samples, M represents CKD samples, and D represents DKD samples. T represents the tuff taken from the region and least affected by alteration
The Ba contents of the BKD, CKD, and DKD samples are close and range as follows: 833–1430, 608–1071, and 742–1503 ppm, respectively. The Ba content of the host rock was 921 ppm. The Sr concentrations of the BKD, CKD, and DKD samples range 737.5–834.5, 567.7–1459.6, and 50.7–93.4 ppm, respectively. The Sr content of the DKD samples is small compared to other deposits and host rocks (205.6 ppm) (Table 2). The Rb concentration in some samples of DKD was much greater (13.2–184.3 ppm) than those of the BKD and CKD deposits and host rock (Table 2). The Y, La, and Ce contents of the BKD samples ranged: 14.4–14.7, 17.6–33.3, and 33.2–59.6 ppm, respectively. The Y, La, and Ce concentrations of the CKD samples ranged: 1.4–35.7, 7.8–33.6, and 9.6–64.5 ppm, respectively. The Y, La, and Ce contents of the DKD were greater than the other deposits, ranging from 20.9 to 32.85, 40.1 to 61.2, and 79.0 to 190.9 ppm, respectively. The Y, La, and Ce contents of the host rock were small at 19.3, 0.6 and 0.5 ppm, respectively (Table 2). As can be seen from Table 2, the REE contents of the DKD are generally greater than the other two deposits and host rock.
Stable Isotope Geochemistry of the Kaolin Deposits
δ34S isotope analysis was performed on alunite minerals, the purity of which was determined by XRD analysis. The δ34S values of the BKD samples were enriched at 19.8 and 20.2‰ and the samples were slightly depleted at –6.3‰. The δ34S isotopic DKD alunite values were slightly depleted at –0.9 and –1.7‰ (DD-1 and DD-2) whereas the DC-1, DC-2, DC-3, and DC-4 samples δ34S isotopic contents were heavily depleted at –15.3 to –17.2‰. The DU-1 sample δ34S isotopic content was highly enriched at 10.5‰. Three δ34S isotopic values of the CKD samples ranged from 5.7 to 16.8‰; these samples were enriched in δ34S. The δ34S isotopic contents of galena and pyrite were slightly depleted at –0.3 and –0.1‰, respectively (Table 3).
Table 3 δ18O and δD isotopic content of kaolinite and calculated model formation temperature, assuming the δ18O value of meteoric water to be –8‰ (Baba & Gunduz, Reference Baba and Gunduz2010). δ34S isotopic content of alunite and 40Ar/39Ar dating analyses from the alunite minerals in the study area
δ18O and δD analyses were performed on the purified and enriched kaolin samples. The δ18O values of the BKD kaolin samples ranged from 11.8 to 17.0‰ and δD values ranged from –99 to –93‰. The δ18O values of the DKD kaolin samples ranged from 18.3 to 10.3‰ and δD values ranged from –70 to –102‰. The δ18O ranged from 10.3 to 12.2‰ and δD from –49 to –82‰ for the CKD (Table 3).
Geochronology
Three pure alunite samples taken from the Bahadırlı and Duman kaolinite deposits were analysed for 40Ar/39Ar age dating (Table 3). Two age dates for CKD were taken from Ercan et al. (Reference Ercan, Ece, Schroeder and Gülmez2022). The timing of hydrothermal alteration of the kaolin deposits was Oligocene-based on the data obtained from 40Ar/39Ar dating: (1) 32.2 ± 1.0 to 32.4 ± 1.2 Ma in the BKD; (2) 29.7 ± 0.1 Ma in the DKD; and (3) 22.6 ± 0.2 to 28.4 ± 0.4 Ma in the CKD.
Discussion
Mineralogical Characterization of Quarries
Detailed assessment of kaolin and alunite group minerals in the alteration zones of the Bahadırlı, Duman, and Çaltıkara kaolin deposits provided more detailed information on acid sulfate-rich geothermal fluids. Photomicrograph investigations indicated that secondary kaolinization occurred from primary K-feldspar, plagioclase phenocrysts, and the matrix of the porphyritic textured lavas and tuffs (Fig. 3).
The mineralogical association of BKD was kaolinite + quartz ± alunite ± smectite ± plagioclase ± K-feldspar (Fig. 4a). The increasing alunite content and the small smectite content in the lower zones where kaolinization is significant indicated formation in a low-pH environment; the increasing amount of smectite in the samples taken from the upper zones was relatable to the increase in pH, however. A kaolinite + quartz + smectite mineral assemblage in the upper zone represented the areas in which the hydrothermal solution pH is neutralized vertically. This may be interpreted as meteoric water inflow. The mineralogical composition of CKD is quite homogeneous and consists of kaolinite + quartz ± alunite (Fig. 4b), which is the most typical mineral association of hydrothermal alteration zones indicating an acidic environment (Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000). The mineralogical composition of the DKD hanging wall block samples was similar to that of the BKD, but the footwall block exhibited a different and heterogeneous mineralogical distribution. The hanging wall block was raised by tectonic activity; the Miocene lacustrine sediments above this wall were probably eroded completely and the consequent hanging wall block represents the Oligocene. Miocene lacustrine deposits in the footwall block were not affected by erosion; the mineralogical association of this block is more complex than that of the hanging wall block, therefore. Footwall blocks were considered as lower levels, middle levels, and upper levels in terms of mineralogical variations. The mineralogical composition of the lower levels of the footwall block is similar to that of the hanging wall block; during the deposition of Miocene lacustrine deposits, the middle and upper levels were affected by freshwater flow into the lake environment and the solution acidity decreased, therefore. In conjunction with increasing pH, smectite became the dominant mineral. Primary plagioclase minerals were also observed with the decrease in solution acidity at these levels. The gypsum minerals detected in some samples of the upper levels of the footwall block were associated with lacustrine deposits rather than hydrothermal alteration because gypsum was not detected in any of the other hydrothermal kaolin deposits in the region.
Kaolinization is well developed in the kaolinite-rich zone of the BKD, and semi-booklet and randomly scattered kaolinite crystals were detected in the FE-SEM images (Fig. 5). Kaolinite crystal dimensions vary generally between 2 and 5 µm. Randomly scattered kaolinite has similar dimensions to the book-shaped kaolinite. Book-shaped and vermiform kaolinites were observed with the altered feldspar minerals in the hanging wall samples of the DKD (Fig. 6). When the increase in crystallization index in kaolinization occurred, it proceeded by gradual structural changes concomitant with crystal coarsening and changes from booklet to blocky (vermiform) morphology (Beaufort et al., Reference Beaufort, Cassagnabere, Petit, Lanson, Berger, Lacharpagne and Johansen1998). This morphological progress indicates that the advanced stage of kaolinization is usually possible with increasing temperature and decreasing pH. Randomly oriented kaolinite crystals were observed commonly in the FE-SEM images of the CKD and this crystal structure is usually formed in alteration zones where low temperatures were effective (Bailey, Reference Bailey, Murray, Bunday and Harvey1993).
In XRD analysis of the BKD, CKD, and DKD, alunite was detected in all deposits, but only in FE-SEM analyses of BDK samples were alunite crystals of <10 µm observed (Fig. 5d). Alunite found typically in low-temperature and steam-heated environments have grain sizes of 20–50 µm, whereas grain sizes of hypogene alunite are coarse and in the range of 50 to 100 µm (Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000). This suggests that alunite in the BKD formed at relatively low temperatures in a steam-heated environment. Dissolution pits were observed on alunite crystal surfaces and they were interpreted to be the result of either secondary surface weathering or of dissolution under highly acidic conditions (Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016). In Fig. 5d, alunite intergrowth with halloysite was observed, and one or two halloysite tubes were noted on the alunite crystals. The mineral association of halloysite and alunite was interpreted as representing an acidic (< 4 pH) environment of formation (Ece et al., Reference Ece, Schroeder, Smilley and Wampler2008). Halloysite tubes were observed in some samples from BKD and in footwall block samples from DKD (Fig. 5c, d and 6e, f). The halloysite tubes are quite long (1000–2000 nm), but the tube diameters (< 100 nm) are small. The morphology and dimensions of the halloysite can be related to the chemical compositions of the tubes. A suggestion by Bailey (Reference Bailey, Murray, Bunday and Harvey1993) was that long tubes are generally rich in A12O3 but poor in Fe2O3. The weathering of feldspar crystals and the formation of short (<500 nm) and thick (>100 nm) halloysite crystals were observed in the upper-level samples of the hanging wall block in DKD (Fig. 6c, d). These images revealed the halloysitization of the feldspars quite clearly.
The coexistence of euhedral (<30 μm) and spherical quartz (<30 μm) crystals was observed in the FE-SEM analysis of CKD (Fig. 5g, h). This morphological difference in the quartz crystals indicates that quartz polymorphs occurred under different physicochemical conditions and may be non-contemporary (Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022). Spherical opal morphologies are cited as immature quartz. Below the vadose zone, opal transforms to chalcedony over time, and again with the effect of time, chalcedony transforms into coarse crystalline quartz (Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000; Sillitoe, Reference Sillitoe, Kirkham, Sinclair, Thorpe and Duke1993, Reference Sillitoe2015).
Evaluation of the regional environment from the general mineralogical composition of the quarries is very important in order to understand the effective hydrothermal alteration processes of the study area and surrounding Etili. The XRD and FE-SEM analyses indicated that BKD, CKD, and DKD are composed mainly of kaolinite, halloysite, alunite, and quartz. All these minerals are commonly associated with advanced argillic alteration assemblages at lower temperatures and these are considered the basic mineral assemblage related to acid-sulfate geothermal fluids. The mineral groups in all of the kaolin deposits indicate a high-sulfidation environment under acidic (pH 2–4) conditions (Dill et al., Reference Dill, Botz and Luppold2005; Ece et al., Reference Ece, Schroeder, Smilley and Wampler2008, Reference Ece, Ekinci, Schroeder, Crowe and Esenli2013; Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000; Pirajno, Reference Pirajno2012). Kaolinite and alunite-group minerals are acid-stable in low-temperature environments above the groundwater table (Sillitoe, Reference Sillitoe, Kirkham, Sinclair, Thorpe and Duke1993). Alunite and kaolin in association occur typically under highly acidic and oxidizing conditions in a high-sulfidation environment (Berger & Henley, Reference Berger and Henley1989; Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000). In this environment of formation, the pH of the steam-heated acid-sulfate fluids is ~2–3 and the presence of hydrochloric acid appears to be a source of this acidity (Schoen et al., Reference Schoen, White and Hemley1974). Fluids in the pH 2–3 range can dissolve minerals such as volcanic glass, feldspar, and mica but not secondary kaolinite and alunite group minerals because Al is insoluble over these pH ranges (Schoen et al., Reference Schoen, White and Hemley1974; Sillitoe, Reference Sillitoe, Kirkham, Sinclair, Thorpe and Duke1993). All the above-mentioned research indicates that the kaolin deposits with associations rich in kaolin + alunite minerals are formed by the effect of steam-heated, acid sulfate fluids with low pH. Smectite + halloysite + feldspar associations are determined at the upper levels of the BKD and upper and middle levels of the DKD. These minerals represent neutral alkali-chloride conditions overlain by a steam-heated, acid-sulfate environment (Simon et al., Reference Simon, Frank, Pettke, Candela, Piccoli and Heinrich2005). Halloysite, smectite, and illite assemblages indicate pH 5 to 6 in the modern geothermal systems; such mineral associations suggest, therefore, that hydrothermal fluids had been neutralized vertically through interaction with the host rock (Reyes, Reference Reyes1990).
Geochemical Characterization of Quarries
The major oxide contents of the samples taken from the BKD, CKD, and DKD are different from each other, and this difference was evident in the SiO2 contents (Table 2). Silica-rich fluids, which are the end members of alkali-chlorine hydrothermal fluids in epithermal systems, precipitate as SiO2 in the pores and fractures of the rocks when they reach high saturation while moving through the rocks or along fracture zones. Previous research in this region revealed that an intense silicification generally occurs in the region (Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022). Therefore, the high SiO2 content in the samples can be associated with silicification. The SiO2/Al2O3 ratio is between 2.69 and 5.11 for the BKD, 2.7 to 3.2 for the CKD, and 1.23 to 3.35 for the DKD; these ratios are far from an ideal kaolinite value (0.992–1.082). Because the kaolinite crystal forms matched the ideal kaolinite form in FE-SEM and XRD analysis, these high rates were interpreted as a silicification effect. The enrichment in Fe2O3 in the samples was interpreted as the precipitation of Fe-rich minerals from the epithermal solution enriched with FeO as a result of alteration of iron-rich minerals when suitable physicochemical conditions were encountered. The amounts of the other major and trace elements in the total mass decreased with the increase in the amount of SiO2 related to silicification.
The main fractionation in the major oxide composition in and around the study area was influenced by the altered mineral associations of the host rock. The alteration of K-feldspar, plagioclase, biotite, hornblende, and matrix resulted in enrichment of Al, Si, Fe, Mg, and K in the solution and kaolinization occurred in the acidic open hydrological system, while in an immobile and slightly alkaline environment ‘smectitization’ took place (Erkoyun et al., Reference Erkoyun, Kadir and Huggett2019). Kaolinite and alunite minerals are compatible with the characteristic alteration products of the BKD, CKD, and DKD hydrothermal system. These indicate that the chemistry and pH of hydrothermal fluids changed spatially and temporally from acid-sulfate (kaolin group minerals) and slightly alkaline (smectite) to alkali-chloride solutions (silica sinter and silicification processes). The increment of K/(Ca + Na) ratio associated with S and Al indicates the formation of alunite under acidic environmental conditions, compatible with those reported in western and central Anatolia (Çelik Karakaya et al., Reference Çelik Karakaya, Karakaya, Temel and Yavuz2021; Ece & Schroeder, Reference Ece and Schroeder2007; Ece et al., Reference Ece, Schroeder, Smilley and Wampler2008; Kadir et al., Reference Kadir, Külah, Eren, Onalgil and Gürel2014, Reference Kadir, Ateş, Erkoyun, Külah and Esenli2022; Mutlu et al., Reference Mutlu, Sarιiz and Kadir2005; Sayın, Reference Sayιn2007). The alunite minerals are typical hypogene minerals that formed during the initial phase of kaolinite formations and were observed frequently in all kaolin deposits in the study area.
The Ba, Sr, Rb, and K contents are enriched in all three kaolin deposits. Ba and Sr are affected quite rapidly by acidic hydrothermal fluids where alteration processes were effective. The role of feldspars and mica in the alteration processes causes significant enrichment in Ba and Sr elements in the alteration product. Most altered volcanics display an apparent enrichment in large ion lithophile elements (LILE) and that they retain the magmatic-arc volcanic origin (Ece & Schroeder, Reference Ece and Schroeder2007). For LILE enrichment of the altered rocks, K-feldspar is recommended as a source mineral, if the geothermal solution is insufficiently saturated in deep reservoirs, K-feldspar hydrolyzes and as a result, the dissolution of K-feldspar changes the chemical composition of geothermal fluids, especially for the enrichment of LILE (Dill et al., Reference Dill, Bosse, Henning and Fricke1997).The depletion of Ce, Y, and La, together with the negative Eu anomaly and enrichment in Sr and Ba, shows that the fractionation by feldspar and mica originated from volcano-sedimentary rocks during the magmatic and hydrothermal alteration processes and was the primary source for the formation of kaolinite and alunite (Kadir et al., Reference Kadir, Ateş, Erkoyun, Külah and Esenli2022). Trace element values for the altered minerals at BKD, CKD, and DKD were determined using the studies of Dill et al. (Reference Dill, Bosse, Henning and Fricke1997). Most of the kaolin samples at BKD and CKD fell in the hypogene area, but the DKD samples fell in the supergene area (Fig. 7a). According to MacLean and Barrett (Reference MacLean and Barrett1993), in a Zr–TiO2 binary diagram used to clarify the hypogene-supergene effects of kaolin deposits, all samples indicate a mixed-origin type (Fig. 7b). Projections, based on elemental data from these two binary graphs, indicate that all the altered samples from three kaolin deposits occurred under the effects of supergene and hypogene solutions. In addition, the DKD kaolins were more affected by the supergene solutions which can be related to waters of lacustrine origin (Fig. 7).
Fig. 7 a Ba + Sr/Ce + Y + La binary diagram and b Zr/TiO2 binary diagram used for the determination of the origin of the acidic waters causing the formation of kaolinization in the Bahadırlı, Çaltıkara, and Duman kaolin deposits (modified from Dill et al., Reference Dill, Bosse, Henning and Fricke1997)
The REE contents of DKD are enriched compared to CKD and BKD (Fig. 8). The REEs are relatively insoluble; many studies have shown that REEs can be mobilized and can fractionate during the formation, deposition, and diagenesis of sediments and rocks, thereby providing a valuable geochemical tool for diagenetic studies of highly weathered materials (Cheshire et al., Reference Cheshire, Bish, Cahill, Kertesz and Stack2018). Ion adsorption-style deposits occur in the weathered crusts of granitic rocks where REE behaving as positively charged trivalent ions during weathering are considered to be adsorbed on negatively charged surfaces of clay minerals such as kaolinite and halloysite (Murakami & Ishihara, Reference Murakami and Ishihara2008). Under low ionic strength, as the REEs are adsorbed through ion exchange and surface complexation; the REE sorption capacity of kaolin and halloysite increases at high pH levels and they also adsorb REE in a similar ratio (Yang et al., Reference Yang, Liang, Ma, Huang, He and Zhu2019). For this reason, it is suggested that the probable formation environment of DKD had a higher pH than other kaolin deposits and, therefore, REE adsorption occurred more here than in the other two deposits. The high pH of the solution effective in this deposit was also suggested and supported by its mineralogical composition.
Fig. 8 Graphic showing REE contents of kaolin samples and tuff taken from Bahadırlı, Çaltıkara, and Duman quarries. Blue lines: Bahadırlı kaolinite, gray lines: Çaltıkara kaolinite, green lines: Duman kaolinite, red line: Tuff. Vertical axis: concentration (ppm)
Stable Isotope Characterization of Quarries
The alunite δ34S isotope values from all of the deposits are quite different from each other and the sulfur isotope data indicate at least three different environments of formation. The δ34S isotope values of the galena and pyrite samples ranged from –0.1 to –0.3‰ and two alunite samples from the hanging wall block of the DKD ranged from –1.7 to 0.9‰ (Table 3). These pointed to a magmatic-hydrothermal sulfur source (–3 ± 3%; Rye et al., Reference Rye, Bethke and Wasserman1992). However, the δ34S values in alunite were depleted at the footwall blocks samples of the DKD. The δ34S isotope values ranged from –15.3 to –17.2‰ for this block, indicating the presence of reduced sulfur produced by a microbial reducing environment or a sedimentary source for sulfur (Cunningham et al., Reference Cunningham, Rye, Steven and Mehnert1984, Reference Cunningham, Rye, Rockwell, Kunk and Councell2005) (Table 3). In the case of oxidation of the starting material in sulfur redox processes, the product formed becomes enriched in δ34S, while reduction causes the formation of products depleted in δ34S (Seal et al., Reference Seal, Alpers and Rye2000). Sulfate-reducing bacteria are found mostly in anoxic media such as below the sediment–water interface and in anoxic water bodies. These bacteria can stay alive in the temperature range 0 to 110°C and pH conditions of 5 to 9.5 (Postgate, Reference Postgate1984). The data indicating this high-pH environment are consistent with the interpretations about the environment obtained from the mineralogical and chemical data of the DKD. The BKD, CKD, and DKD (DU-1, representing the upper levels of the footwall block) samples were enriched in δ34S, with very large positive values being attained (Table 3). During pyrite precipitation by bacterial reduction of sulfate, δ32S was enriched in the sulfide, which caused enrichment of δ34S in the residual sulfate-rich solution; therefore, secondary sulfate minerals are enriched in δ34S (Seal et al., Reference Seal, Alpers and Rye2000). Ercan et al. (Reference Ercan, Ece, Schroeder and Gülmez2022) suggested that the organic acid input may result from the interaction of lacustrine coal deposits with groundwater/geothermal fluids, observed commonly in the study area. In addition, the detection of retiform structures of microfossils in micromorphological studies ((supplementary material Table S1, ESM-1pdf)) also supported the possibility of the existence of bio-organism activities in this quarry.
δ18O versus δD isotopic values of the BKD fall on the kaolinite line (two samples) and halloysite line (one sample) and also kaolinite δ18O versus δD isotopic values of the DKD fall between the kaolinite and halloysite lines (Fig. 7). The δ18O and δD isotopic values of CKD are very close to the supergene/hypogene line proposed by Sheppard and Gilg (Reference Sheppard and Gilg1996) (Fig. 9) and only one kaolinite sample falls near the magmatic water area. Although δ18O and δD isotope values of all CKD kaolins indicated a hypogene origin, kaolin samples from BKD and CKD indicate a supergene origin. All of these isotope data are very consistent with and confirm the formation environment results obtained from trace-element data. In addition, the isotopic data revealed clearly that DKD was more affected by hypogene origin solutions, but DKD was more exposed to supergene effects. The calculated model temperatures of the CKD range from 69 to 84ºC, these values reflecting temperatures of the fossil geothermal fields. The CKD probably formed in a hypogene environment and isotopic data plot very close to the line representing kaolinites in equilibrium with meteoric waters at 35ºC. These lower temperatures were then possibly related to a steam-heated process, which condensed waters in the vadose zone, mixed with meteoric water, and then isotopically re-equilibrated (Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016).
Fig. 9 The δ18O and δD binary diagram. The δ18O and δD isotopic values of kaolins from the Bahadırlı, Çaltıkara, and Duman kaolin deposits, hydrothermal kaolins from active geothermal fields, and thermal waters near kaolin locations in Japan (Hayba et al., Reference Hayba, Bethke, Heald, Foley, Berger and Bethke1985). The meteoric water line is δD = 6.8*δ18O + 10.5 (Harris et al., Reference Harris, Compton and Bevington1999). Supergene vs. hypogene (S/H) lines are taken from Sheppard and Gilg (Reference Sheppard and Gilg1996). Equilibrium lines for kaolinite and halloysite at various temperatures were calculated assuming meteoric water of δ18O = –8‰ for the Kirazlı region of Çanakkale; this value is taken from Baba and Gunduz et al. (Reference Baba and Gunduz2010)
The area to the left of the S/H line indicates a hypogene origin for the CKD, whereas the areas to the right-hand side indicate a supergene environment for the BKD and DKD. Component-calculated model formation temperatures of the BKD and DKD deposit are between 40–72ºC and 32–61ºC, respectively, thus reflecting isotopic equilibrium with meteoric water in the supergene environment (Table 3).
The consensus is that calc-alkaline plutonism was active from Early Eocene to Early Miocene and hydrothermal alterations were active from Late Eocene through Miocene. The findings of 40Ar/39Ar alunite dates (22.6 ± 0.2 to 32.4 ± 1.2 Ma) (Table 3) are consistent with the alteration dates of primary volcanic rocks (29.7 to 32.4 Ma) (Altunkaynak & Yılmaz, Reference Altunkaynak and Yιlmaz1998; Altunkaynak et al., Reference Altunkaynak, Sunal, Aldanmaz, Genç, Dilek, Furnes and Yιldιz2012; Altunkaynak & Dilek, Reference Altunkaynak and Dilek2013; Buyukkahraman & Çoban, Reference Buyukkahraman and Coban2013; Ercan et al., Reference Ercan, Ece, Schroeder and Karacik2016; Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022; Özdamar, Reference Özdamar, Zou, Billor and Hames2021b; Yılmaz & Karacık, Reference Yιlmaz and Karacιk2001). Consequently, radiometric dating in the region revealed that the active period of calc-alkaline plutonism and the alteration processes due to hydrothermal activity were highly correlated and occurred at approximately the same time. These dates indicate that the alunite at BKD and CKD probably formed during igneous hydrothermal activity, but the isotopic studies of DKD revealed that, besides the hydrothermal activity, groundwater was mixed with organic acid derived from the coal deposits of lacustrine origin. Because the amount of alunite collected from the CKD was insufficient, 40Ar/39Ar dating was not carried out, but the dating obtained from the south of this deposit (Ercan et al., Reference Ercan, Ece, Schroeder and Gülmez2022) was in the range 22.6–28.5 Ma. These data are compatible with the dating of DKD and BKD quarries which are located along the same fault.
Conclusions
Data from the mineral assemblages and chemical and isotopic compositions of the Bahadırlı, Çalıkara, and Duman kaolin deposits can be used to help distinguish between supergene and hypogene origins in hydrothermal kaolin deposits. Based on two different mineral assemblages of alteration groups, two different hydrothermal stages were proposed for the region. In the first stage, representing BKD, CKD, and the hanging-wall side of DKD, highly acidic and oxidizing conditions prevailed, as shown by the presence of alunite and kaolinite, an indicator of high-temperature mineralization, formed in relatively deep epithermal zones (Berger & Henley, Reference Berger and Henley1989; Hedenquist et al., Reference Hedenquist, Arribas and Gonzalez-Urien2000; Sillitoe, Reference Sillitoe, Kirkham, Sinclair, Thorpe and Duke1993). During the second stage, which is represented by the middle and upper levels of the DKD, kaolinite, smectite, and halloysite were formed and some of the original textures of the tuffs were destroyed due to lower intensity chemical reactions of geothermal acids, organic humic acids, and meteoric waters. The mixture of groundwater and/or geothermal waters and humic acid fluids during the alteration in DKD was also supported by isotopic data. In this study, humic acids are suggested teo originate from lacustrine coal deposits, which are commonly observed in the Etili region (Table 4).
Table 4 Determination of hypogene and supergene alteration criteria of kaolin deposits in the Etili region
During the Oligocene–Early Miocene to Middle and Upper Miocene, the evolution of alteration processes and hydrothermal fluids in the region is shown clearly in the footwall block of the DKD. The footwall blocks of the DKD reflect the mineralogical, chemical, and isotopic properties dominated by the effect of acid-sulfate hydrothermal fluids at the base. However, the mineralogical composition, geochemical properties, and 34S isotope characterization of the samples of the upper zones of the footwall block revealed the hypogene-supergene transition with the effect of the alkali-chloride fluids.
This study revealed that three kaolin deposits located in the same host rock and very close to each other in the Çanakkale Region experienced different stages and degrees of hydrothermal alteration processes. Each deposit was presumably affected by the groundwater and meteoric waters differentially heated by hydrothermal processes. The BKD, CKD, and DKD quarries are very close to each other, offer many different mineralogical, geochemical, and isotopic data, and can be used as a guide in solving and understanding the origin of similar formations throughout the world.
Supplementary Information
The online version contains supplementary material available at https://doi.org/10.1007/s42860-022-00214-4.
Acknowledgements
The authors are grateful to ITU-JAL labs for providing analytical support during geochemical analysis in ITU, Turkey. They also thank the Kalemaden Mining Corp. for logistical support and accommodation during the fieldwork.
Funding
This research was supported by TUBITAK Project No: 112Y369 (Project coordinator: Prof. Dr. Ömer Işık Ece).
Data Availability
All data related to the manuscript are presented in tables.
Code availability
Not applicable.
Credit authorship contribution statement: Hatice Ünal Ercan: Conceptualization, investigation, writing—original draft, writing—review and editing, visualization. Ömer Işık Ece: Conceptualization, investigation, writing—review and editing, project administration, funding acquisition. Emin Çiftçi: Conceptualization, investigation, writing—review and editing. Ayça Aydın: Investigation, formal analysis, review and editing.
Declarations
Conflicts of Interest
The authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.
