Geochemical Monitoring of Geothermal Waters (2002–2004) along the North Anatolian Fault Zone, Turkey: Spatial and Temporal Variations and Relationship to Seismic Activity

A total of nine geothermal fields located along an 800-km long E-W transect of the North Anatolian Fault Zone (NAFZ), Turkey were monitored for three years (2002–2004 inclusive; 3-sampling periods per year) to investigate any possible relationship between seismic activity and temporal variations in the chemistry and isotope characteristics of waters in the fields. The geothermal fields monitored in the study were, from west to east, Yalova, Efteni, Bolu, Mudurnu, Seben, Kur?unlu-Çank?r?, Hamamözü, Gözlek and Re?adiye. The chemical (major anion-cation contents) and isotopic (¹⁸O/¹⁶O, D/H, ³H) compositions of hot and cold waters of the geothermal sites were determined in order to both characterize the chemical nature of the individual fields and identify possible temporal variations associated with localized seismic activity. The geothermal waters associated with the NAFZ are dominantly Na-HCO₃, whereas the cold waters are of the Ca-HCO₃ type. The oxygen- and hydrogen-isotope compositions reveal that the hot waters are meteoric in origin as are their cold water counterparts. However, the lower δ¹⁸O, δD and ³H contents of the hot waters point to the fact that they are older than the cold waters, and that their host aquifers are recharged from higher altitudes with virtually no input from recent (post-bomb) precipitation. Although no major earthquakes (e.g., with M ≥ 5) were recorded along the NAFZ during the course of the monitoring period, variations in the chemical and isotopic compositions of some waters were observed. Indeed, the timing of the chemical/isotopic changes seems to correlate with the occurrence of seismic activity of moderate magnitude (3 < M < 5) close to the sampling sites. In this respect, Cl, ³H and Ca seem to be the most sensitive tracers of seismically-induced crustal perturbations, and the Yalova and Efteni fields appear to be the key localities where the effects of seismic activity on the geothermal fluids are most pronounced over the monitoring period. The present study has produced a ‘baseline’ database for future studies directed at characterizing the effects of moderate-major earthquakes on the composition of geothermal waters along the NAFZ. Future work involving longer monitoring periods with more frequent sampling intervals should lead to a better understanding of the underlying mechanism(s) producing the observed chemical and isotopic variations.     

Pan-African magmatism in the Menderes Massif: geochronological data from leucocratic tourmaline orthogneisses in western Turkey

The Menderes Massif, exposed in western Anatolia, is a metamorphic complex cropping out in the Alpine orogenic belt. The metamorphic rock succession of the Massif is made up of a Precambrian basement and overlying Paleozoic-early Tertiary cover series. The Pan-African basement is composed of late Proterozoic metasedimentary rocks consisting of partially migmatized paragneisses and conformably overlying medium- to high-grade mica schists, intruded by orthogneisses and metagabbros. Along the southern flank of the southern submassif, we recognized well-preserved primary contact relationship between biotite and leucocratic tourmaline orthogneisses and country rocks as the orthogneisses represent numerous large plutons, stocks and vein rocks intruded into a basement of garnet mica schists. Based on the radiometric data, the primary deposition age of the precursors of the country rocks, garnet mica schist, can be constrained between 600 and 550?Ma (latest Neoproterozoic). The North Africa–Arabian-Nubian Shield in the Mozambique Belt can be suggested as the possible provenance of these metaclastics. The intrusion ages of the leucocratic tourmaline orthogneisses and biotite orthogneisses were dated at 550–540?Ma (latest Neoproterozoic–earliest Cambrian) by zircon U/Pb and Pb/Pb geochronology. These granitoids represent the products of the widespread Pan-African acidic magmatic activity, which can be attributed to the closure of the Mozambique Ocean during the final collision of East and West Gondwana. Detrital zircon ages at about 550?Ma in the Paleozoic muscovite-quartz schists show that these Pan-African granitoids in the basement form the source rocks of the cover series of the Menderes Massif.     

Pressure–temperature evolution of lawsonite eclogite in Sivrihisar; Tavşanlı Zone–Turkey

The Cretaceous blueschist belt, Tavşanlı Zone, representing the subducted and exhumed northern continental margin of the Anatolide–Tauride platform is exposed in Western Anatolia. The Sivrihisar area east of Tavşanlı is made up of tectonic units consisting of i) metaclastics and conformably overlying massive marbles (coherent blueschist unit), ii) blueschist-eclogite unit, iii) marble–calcschist intercalation and iv) metaperidotite slab. The metaclastics are composed of jadeite–lawsonite–glaucophane and jadeite–glaucophane–chloritoid schists, phengite phyllites, and calcschists with glaucophane–lawsonite metabasite layers. The blueschist-eclogite unit representing strongly sheared, deeply buried and imbricated tectonic slices of accreted uppermost levels of the oceanic crust with minor metamorphosed serpentinite bodies consists of lawsonite-bearing eclogitic metabasites (approximately 90% of the field), lawsonite eclogites, metagabbros, serpentinites, pelagic marbles, omphacite–glaucophane–lawsonite metapelites and metacherts. The mineral assemblage of the lawsonite eclogite (garnet + omphacite > 70%) is omphacite, garnet, lawsonite, glaucophane, phengite and rutile. Lawsonite eclogite lenses are enclosed by garnet–lawsonite blueschist envelopes.Textural evidence from lawsonite eclogites and country rocks reveals that they did not leave the stability field of lawsonite during subduction and exhumation. The widespread preservation of lawsonite in eclogitic metabasites and eclogites can be attributed to rapid subduction and subsequent exhumation in a low geothermal gradient of the oceanic crust material without experiencing a thermal relaxation. Peak P–T conditions of lawsonite eclogites are estimated at 24 ± 1 kbar and 460 ± 25 °C. These P–T conditions indicate a remarkably low geotherm of 6.2 °C/km corresponding to a burial depth of 74 km.     

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high‐pressure metasedimentary rocks

The multiple high‐pressure (HP), low‐temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction‐ and continental accretion‐related evolution of the eastern limb of the long‐lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe‐Mg‐carpholite in three metasedimentary units of the Gondwana‐derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single‐continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica ⁴⁰Ar–³⁹Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi‐equilibrium thermobarometry), on carpholite‐bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite‐bearing assemblages were retrogressed through greenschist‐facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post‐collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tav?anl? Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian‐type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll‐back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.     

Observation of excess heavy metal concentrations in water resources to infer surface water influences on shallow groundwater: a typical example of the Porsuk River (Eskisehir-Turkey)

Water Resources - The Eskisehir province is well-known due to its industrial and agricultural activities, which are a threat for the aquatic environment. Hence, monitoring of water quality in the...     

Susceptibility assessments of shallow earthflows triggered by heavy rainfall at three catchments by logistic regression analyses

Sometimes regional meteorological anomalies trigger different types of mass movements. In May 1998, the western Black Sea region of Turkey experienced such a meteorological anomaly. Numerous residential and agricultural areas and engineering lifelines were buried under the flood waters. Besides the reactivation of many previously delineated landslides, thousands of small-scale landslides (mostly the earthflow type) occurred all over the region. The earthflows were mainly developed in flysch-type units, which have already presented high landslide concentrations. In this study, three different catchments — namely Agustu, Egerci, and Kelemen — were selected because they have the most landslide-prone geological units of the region. The purposes of the present study are to put forward the spatial distributions of the shallow earthflows triggered, to describe the possible factors conditioning the earthflows, and to produce the shallow earthflow susceptibility maps of the three catchments. The unique condition units (UCU) were employed during the production of susceptibility maps and during statistical analyses. The unique condition units numbered 4052 for the Agustu catchment, 13,241 for the Egerci catchment and 12,314 for the Kelemen catchment. The earthflow intensity is the highest in the Agustu catchment (0.038 flow/UCU) and lowest in the Egerci catchment (0.0035 flow/UCU). Logistic regression analyses were also employed. However, during the analyses, some difficulties were encountered. To overcome the difficulties, a series of sensitivity analyses were performed based on some decision rules introduced in the present study. Considering the decision rules, the proper ratios of UCU free from earthflow (0) / UCU including the earthflow (1) for the Agustu, Egerci and Kelemen catchments were obtained as 3, 6, and 5, respectively. Also, a chart for the proper ratio selection was developed. The regression equations from the selected ratios were then applied to the entire catchment and the earthflow susceptibility maps were produced. The landslide susceptibility maps revealed that 15% of the Agustu catchment, 8% of the Egerci catchment, and 7% of the Kelemen catchment have very high earthflow susceptibility; and most of the earthflows triggered by the May 1998 meteorological event were found in the very high susceptibility zones.     

Assessment of landslide susceptibility for a landslide-prone area (north of Yenice, NW Turkey) by fuzzy approach

. Regional landslide susceptibility assessments pose complex problems. To solve these problems, numerous approaches, such as statistical analysis, geotechnical engineering approach, geomorphologic approach and fuzzy logic, have been employed. However, all the available methods for regional landslide susceptibility assessments have some uncertainties due to a lack of knowledge and variability. Minimizing these uncertainties provides realistic approaches. Use of the fuzzy logic approach to produce a landslide susceptibility map of a landslide-prone area in NW Turkey is the main purpose of the present study. For this purpose, the study includes five main stages, these being the preparation of a landslide inventory of the study area, the application of factor analysis, the extraction of fuzzy if-then rules, the use of a geographical information system, and the control of the reliability of the resulting landslide susceptibility map. Slope angle, slope aspect, land use, weathering depth, water conditions and topographical elevation were considered as landslide conditioning factors for the study area. A total of 23 if-then rules was extracted from the field data. Employing these rules, fuzzified index maps representing each parameter were obtained. Finally, combining these maps, the landslide susceptibility map of the area was prepared. When compared with the landslide susceptibility map, the landslides identified in the area were found to be located in the very high- and high-susceptibility zones. As far as the performance of the fuzzy approach for processing is concerned, the images appear to be quite satisfactory, the zones determined on the map being zones of relative susceptibility.     

Prediction of earthquake hazard by hidden Markov model (around Bilecik,NW Turkey)

Earthquakes are one of the most important natural hazards to be evaluated carefully in engineering projects, due to the severely damaging effects on human-life and human-made structures. The hazard of an earthquake is defined by several approaches and consequently earthquake parameters such as peak ground acceleration occurring on the focused area can be determined. In an earthquake prone area, the identification of the seismicity patterns is an important task to assess the seismic activities and evaluate the risk of damage and loss along with an earthquake occurrence. As a powerful and flexible framework to characterize the temporal seismicity changes and reveal unexpected patterns, Poisson hidden Markov model provides a better understanding of the nature of earthquakes. In this paper, Poisson hidden Markov model is used to predict the earthquake hazard in Bilecik (NW Turkey) as a result of its important geographic location. Bilecik is in close proximity to the North Anatolian Fault Zone and situated between Ankara and Istanbul, the two biggest cites of Turkey. Consequently, there are major highways, railroads and many engineering structures are being constructed in this area. The annual frequencies of earthquakes occurred within a radius of 100 km area centered on Bilecik, from January 1900 to December 2012, with magnitudes (M) at least 4.0 are modeled by using Poisson-HMM. The hazards for the next 35 years from 2013 to 2047 around the area are obtained from the model by forecasting the annual frequencies of M ≥ 4 earthquakes.     

Causes, mechanism and environmental impacts of instabilities at Himmetoğlu coal mine and possible remedial measures

 Since the commencement of mining at the Himmetoğlu coal mine, northwest Turkey, serious stability problems have led to interruptions in mining and some environmental effects. A geotechnical investigation was initiated in 1997 and the significant factors that influence the stability have been defined. This paper outlines the results of the field and laboratory studies associated with the causes and mechanism of the slope instabilities and their environmental impacts. The possible remedial measures to improve the stability and to minimize the environmental problems are also described. Back-analyses and data from long-term monitoring indicate that the failures occur along two or three planar surfaces by combination of faults and localized strata steepening adjacent to the faults. The stability is sensitive to changes in length of the lower part of the basal sliding surface and shear strength of the bedding surfaces in the overburden. Suitable remedial measures include slope flattening (i.e. staged bench stripping), proper drainage and spreading of a rock blanket on the pit floor to increase spoil pile stability. Received: 18 April 2000 · Accepted: 15 August 2000