Evaluation of the temporal variations of groundwater storage and its interactions with climatic variables using GRACE data and hydrological models: A study from Turkey

Monitoring of the fluctuations of groundwater storage is particularly important in arid and semi-arid regions where water scarcity brings about various challenges. Remote sensing data and techniques play a preponderant role in developing solutions to environmental problems. The launch of Gravity Recovery and Climate Experiment (GRACE) satellites has eased the remote monitoring and evaluation of groundwater resources with an unprecedented precision over large scales. Within the scope of the current study, the latest release (RL06) of GRACE mass concentrations (Mascons) from Jet Propulsion Laboratory (JPL) dataset as well as Global Land Data Assimilation System (GLDAS) models of Noah and Catchment Land Surface Model (CLSM) were used to provide Groundwater Storage Anomalies (GWSA) over Turkey. The temporal interactions of the estimated GWSA with the climatic variables of precipitation and temperature (derived from the reanalysis datasets of CHELSA [Climatologies at High resolution for the Earth's Land Surface Areas] and FLDAS [the Famine Early Warning Systems Network Land Data Assimilation System], respectively) were investigated statistically. The results suggest that there is a descending trend (from 2003 to 2016) for Terrestrial Water Storage Anomalies (TWSA) and GWSA over Turkey with a total loss of 11 and 6 cm of water, respectively. The statistical analysis results also indicate that the monthly variations of GWSA over Turkey are highly correlated with precipitation and temperature at 2-month lag. The analysis of the climatology (long-term) values of monthly GWSA, precipitation and temperature also revealed high agreement between the variables.     

Hydrochemical characteristics and the effects of irrigation on groundwater quality in Harran Plain,GAP Project,Turkey

Improper design, faulty planning, mismanagement and incorrect operation of irrigation schemes are the principle reasons for the deterioration of groundwater quality in a large number of countries, in particular in semi-arid and arid regions. The aim of this study is to determine the dimensions of groundwater quality after surface irrigation was begun in the semi-arid Harran Plain. Physical and chemical parameters of the groundwater including pH, temperature, electrical conductivity (EC), sodium, potassium, calcium, magnesium, chloride, bicarbonate, sulphate, nitrate, nitrite, ammonium, total phosphorus, total organic carbon and turbidity were determined monthly during the 2006 water year. The quality of the groundwater in the study area was assessed hydrochemically in order to determine its suitability for human consumption and agricultural purposes. In the general plain, the EC values measured were considerably above the guide level of 650 μS/cm, while nitrate in particular was found in almost all groundwater samples to be significantly above the maximum admissible concentration of 50 mg/l for the quality of water intended for human consumption as per the international and national standards. Total hardness reveals that a majority of the groundwater samples fall in the very hard water category. Interpretation of analytical data shows that Ca–HCO₃ and Ca–SO₄ are the dominant hydrochemical facies in the study area.     

A geological and geotechnical investigation of the settlement area of Zümrüt Building (Konya, Turkey) which caused 92 fatalities due to its collapse

This study examines the local geological conditions and soil structure as possible causes of the collapse of the Zümrüt Building 2 February 2004. This catastrophe resulted in 92 fatalities and 35 injuries. This study also examines other views which claim weak soil structure, elastic and consolidation settlement of soil and excessive groundwater extraction as well as subsidence resulting from the underground silt erosion as possible factors. Zümrüt Building was constructed on normally consolidated, low plasticity clay. The underground water table was 30 m in depth. The internal friction angle of soil was 8°–30°, its cohesion was between 34 and 127 kN/m² and standard penetration test numbers varied between 11 and 50. The underground water level beneath Zümrüt Building had risen 4.5 m since its construction. Therefore the claim that subsidence resulting from the decrease of underground water level contributed to the collapse is incorrect. Secondly the settlement, resulting from the filling up of the pores created by the silt receding with the underground water, was 4.4 mm in total, and attributing this as the primary cause of the collapse is also incorrect. Soil properties, in situ and laboratory test results showed that the existing and/or expected settlement and the differential ground settlement in the Zümrüt building vicinity had the potential to cause structural damage. The tensile stresses caused by differential settlements recorded here are thought to be an indicator, but not the main cause contributing to the collapse of the building. The Zümrüt Building collapse was due to several compounding mistakes during the construction phase. These were geotechnical and other project faults and the use of low quality construction materials. The resulting catastrophe caused 92 fatalities, 35 injuries and a material loss of approximately US$7 million.     

Analysis of drivers governing temporal salinity and temperature variations in groundwater discharge from Altug Submarine Karst Cave (Kas-Turkey)

Salinity and temperature variations in groundwater discharge from the Altug submarine karst cave have been observed at 28 m below sea level for every 10 min between November 2004 and August 2005 to determine the drivers that govern the salinization. Comparisons between temporal trends of salinity and temperature with those of precipitation, air pressure, sea level and wind velocity revealed an apparent dominance of precipitation regime on the salinity and temperature variations. Spectral analyses applied to observations showed that the air pressure and sea level oscillations are affected by sun and moon tides which do not have an appreciable impact on the salinity and temperature variations. Annual rate of salinization in Altug cave seems inversely related to the inland groundwater head so that the maximum and minimum fresh water contributions occur at mid-spring and late-summer, respectively.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.     

Implementation of reconstructed geomorphologic units in landslide susceptibility mapping: the Melen Gorge (NW Turkey)

In the international literature, although considerable amount of publications on the landslide susceptibility mapping exist, geomorphology as a conditioning factor is still used in limited number of studies. Considering this factor, the purpose of this article paper is to implement the geomorphologic parameters derived by reconstructed topography in landslide susceptibility mapping. According to the method employed in this study, terrain is generalized by the contours passed through the convex slopes of the valleys that were formed by fluvial erosion. Therefore, slope conditions before landsliding can be obtained. The reconstructed morphometric and geomorphologic units are taken into account as a conditioning parameter when assessing landslide susceptibility. Two different data, one of which is obtained from the reconstructed DEM, have been employed to produce two landslide susceptibility maps. The binary logistic regression is used to develop landslide susceptibility maps for the Melen Gorge in the Northwestern part of Turkey. Due to the high correct classification percentages and spatial effectiveness of the maps, the landslide susceptibility map comprised the reconstructed morphometric parameters exhibits a better performance than the other. Five different datasets are selected randomly to apply proper sampling strategy for training. As a consequence of the analyses, the most proper outcomes are obtained from the dataset of the reconstructed topographical parameters and geomorphologic units, and lithological variables that are implemented together. Correct classification percentage and root mean square error (RMSE) values of the validation dataset are calculated as 86.28% and 0.35, respectively. Prediction capacity of the different datasets reveal that the landslide susceptibility map obtained from the reconstructed parameters has a higher prediction capacity than the other. Moreover, the landslide susceptibility map obtained from the reconstructed parameters produces logical results.     

Taxonomical and Biogeochemical Investigation of Elements in Holocene Mollusk Shells in the South and Southwestern Marmara Basin,Turkey

Concentrations of some heavy metals and trace elements such as Cr, Ga, Ni, Zn, Mo, Cu, Pb, Yb, Y, Nb, Ti, Sr, Ba, Mn, Sc, Co, V, Zr, Fe, Al, W, Se, Bi, Sb, As, Cd in recent mollusk shells and factors affecting their distribution and deposits collected from various depths in the southern and southwestern parts of the Marmara Sea are investigated. The distribution of the elements in the shells is categorized into four groups. Of these, concentrations of 12 elements （As, Bi, Cd, Co, Ga, Mo, Nb, Sb, Se, Sc, W and Yb） are below zero [（0.053-0.79）×10^-6]; concentrations of seven elements （Cr, Ni, Pb, V, Y, Zr and Cu） are （1.0-6.0）×10^-6; concentrations of four elements （Ti, Mn, Ba and Zn） are 10- 20×10^-6; and concentrations of five elements （Si, Al, Fe, Mg and Sr） are （47.44-268.11）×10^-6. The taxonomic characteristics of the 29 elements were studied separately in mollusk shells such as Chamalea gallina （Linn6）, Pitar rudis （Poli）, Nassarius reticulatus （Linn6）, Venerupis senescens （Coocconi）, Mytilus galloprovincialis （Lamarck）, Mytilaster lineatus （Gemelin in Linne） and Chlamys glabra. It was found that, in mollusk taxonomy, the elements have unique values. In other words, element concentrations in various mollusk shells depend mainly on the taxonomic characteristics of mollusks. In various bionomic environments different element distributions of the same species are attributed to the different geochemical characters of the each environment. Data obtained in this study indicate that the organisms are the most active and deterministic factors of the environment.     

Determination of potential hazardous elements in Can coals （Canakkale- Turkey）-I

It is known that the increasing use of coal as an energy source led to the growing environmental and health problems. But comprehensive knowledge of coal quality parameters may help to reduce some of these problems. The Canakkale-Can coalfield is located in the western part of Turkey, whose reserves are estimated at 69.3 billion tons, mainly used in the industry in the neighboring areas, specifically for the Can Thermal Power Plant. The aim of this study is to determine the coal quality parameters and examine the origin and distributions of potentially toxic trace elements in lignite which may produce environmental and health hazards in the area. The coal samples were collected from different parts of the coalfield in Can. Proximate and ultimate analyses, sulfur form analyses, X-ray powder diffraction （XRD）, inductively coupled plasma-mass spectrometry （ICP-MS）, and scanning electron microscopy （SEM） were performed on those samples to determine the geochemical profile of hazardous elements.     

Hydrology and hydrogeology of Sakaryabaşı Karstic springs, Çifteler, Turkey

The Sakarya River is one of the largest rivers in Turkey and is fed mainly from Sakaryabaşı springs. The Sakaryabaşı springs are located in the Central Anatolia and issue from confined/semi-confined karst having a thermal component and therefore, having quite different hydrogeological characteristics as compared to the Taurus Karst region, a typical example of the Mediterranean type of karst. The karstic carbonate rocks that form the groundwater reservoir are overlain by a thick semi-pervious overburden of mainly clastics of Neogene age. Tectonics is the major factor controlling the occurrence of the karst springs in the area where topography is rather flat. This study aimed at explaining the occurrence and movement of the karst groundwater within the system by use of hydrogeological, chemical, and isotopic tools. Isotopic composition of the waters revealed that all waters in the region are of meteoric origin and the thermal component is due to deep circulation. The catchment area of the hydrogeological system extends to the south and groundwater movement is towards the outlets, which are in a depression along a major fault. The movement of the groundwater, based on analysis of remotely sensed images, is controlled mainly by structural elements.     

Conceptual hydrodynamic model of the Pamukkale hydrothermal field, southwestern Turkey, based on hydrochemical and isotopic data

The results of study on the hydrochemical and isotope characteristics of shallow and deep waters at Pamukkale hydrothermal field Turkey are described in order to obtain a better understanding of the hydrological circulation. The field can be grouped into two groundwater sub-systems; cold water springs of Ca–HCO₃ type (10–12 °C), and CO₂-rich thermal waters of Ca–HCO₃–SO₄ type (25–58 °C). The occurrence of these water types is closely related to the morphology of the region, where intense tectonism formed horst and graben structures. Hence, two hydrogeological systems were defined: a deep geothermal system which is related to extensive and deep circulation of meteoric water in the regional flow system, and a shallow system which is related to local groundwater flow through sedimentary strata. The meteoric water falling at higher elevations percolates to the local groundwater system at a shallow level and flows to the deep geothermal system. During a deep convection cycle from a recharge to discharge area, the cold water attains heat from the asthenospheric intrusions, causing it to ascend. Variations of chemical and isotopic composition of thermal waters result from their mixing with cool groundwater in a shallow aquifer during their ascent to the surface.

---

Résumé Les résultats d'une étude portant sur les caractéristiques hydrochimiques et isotopiques de puits phréatiques et profonds situés dans le champ hydrothermal de Pamukkale, sont décrits de telle manière à éclairer le fonctionnement des circulations hydrologiques. Le champ peut étre divisé en deux sous-systèmes d'eaux souterraines, l'un avec des eaux de sources froides (10–12 °C) de type Ca–HCO3, et les eaux thermales (25–58 °C) riches en CO₂ et de type Ca–HCO₃–SO4. L'occurrence de ces types d'eaux est fermement liée à la morphologie de la région, oùne tectonique intense a engendré des structures en horsts et en grabens. Dés lors deux systèmes hydrogéologiques ont été définis : un système profond, qui est lié à la circulation extensive et profonde des eaux météoritiques dans le système régional d'écoulement, et un système phréatique lié aux écoulements locaux des eaux souterraines à travers les strates sédimentaires. Les eaux météoritiques aux altitudes élevées, percolent jusqu'aux systèmes locaux phréatiques, puis coulent jusqu'aux systèmes géothermaux plus profonds. Durant le cycle de convection profond des zones de recharge jusqu'aux zones de décharge, l'eau froide atteint les zones chaudes liées aux intrusions athenosphériques, provoquant la remontée. Les variations de la composition chimique et isotopique des eaux thermales résultent dans leurs mélanges avec des eaux souterraines froides dans les aquifères phréatiques durant leur remontée jusqu'à la surface.

---

Resumen Se describen los resultados del estudio de las características isotópicas e hidroquímicas de las aguas someras y profundas para obtener un mejor entendimiento de la circulación hidrológica del campo hidrotermal Pamukkale. El campo puede agruparse en dos sub-sistemas de agua subterránea: manantiales de agua fría del tipo Ca–HCO₃ (10–12°C) y aguas termales ricas en CO₂ del tipo Ca–HCO₃–SO₄(25–58°C). El ambiente de estos tipos de aguas se relaciona estrechamente con la morfología de la región donde el tectonismo intenso ha formado estructuras extensionales tipo graben y horst. Se definieron dos sistemas hidrogeológicos: un sistema geotermal profundo que se relaciona con la circulación profunda y extensa de agua meteórica en el sistema regional de flujo y un sistema somero el cual se relaciona con flujo local de agua subterránea a través de estratos sedimentarios. El agua meteórica que cae en altas elevaciones percola al sistema local de agua subterránea en un nivel somero y fluye hacia el sistema geotermal profundo. Durante un ciclo de convección del área de recarga hacia la zona de descarga, el agua fría se calienta a partir de los intrusivos astenosféricos lo que ocasiona que asciendan. Como resultado de la mezcla de las aguas recalentadas, con agua subterránea fría en un acuífero somero durante el ascenso hacia la superficie, se derivan variaciones en la composición química e isotópica de las aguas termales.