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.     

Long-chain alkenones and U37′ variability along a south–north transect in the Western Pacific Ocean

In this paper, we present 50 surface water samples collected during the IMAGES III cruise (June–July 1997) along a transect from New Zealand to the China Sea (42°S–178°E, 21°N–120°E) covering a temperature range from 13.3 to 30.4 °C. A very worthwhile aspect of this study is a coupling of both biomarker (alkenone) and coccolithophorid counting. We show that the U₃₇^k′–temperature relationship is very similar to the Prahl et al. [Geochim. Cosmochim. Acta 52 (1988) 2203] culture calibration and to the global core top calibration of Müller et al. [Geochim. Cosmochim. Acta 62 (1998) 1757]. However, in the warmest surface waters of the Western Pacific ocean (>26.4 °C) where Gephyrocapsa oceanica is likely the most widespread species, the associated U₃₇^k′ has a constant value of 1.0. The consequence is that above this temperature threshold, U₃₇^k′ cannot be used as an accurate paleothermometer.     

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.     

Estimating potential landslide sites of an upland sub-watershed in Western Ghat’s of Kerala (India) through frequency ratio and GIS

The purpose of this study is to assess the susceptibility of landslides in parts of Western Ghats, Kerala, India, using a geographical information system (GIS). Landslide inventory of the area was made by detailed field surveys and the analysis of the topographical maps. The landslide triggering factors are considered to be slope angle, slope aspect, slope curvature, slope length, distance from drainage, distance from lineaments, lithology, land use and geomorphology. ArcGIS version 8.3 was used to manipulate and analyse all the collected data. Probabilistic-likelihood ratio was used to create a landslide susceptibility map for the study area. The result was validated using the Area under Curve (AUC) method and temporal data of landslide occurrences. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations. As the result, the success rate of the model was (84.46%) and the prediction rate of the model was (82.38%) shows high prediction accuracy. In the reclassified final landslide susceptibility zone map, 5.68% of the total area is classified as critical in nature. The landslide susceptibility map thus produced can be used to reduce hazards associated with landslides and to land cover planning.     

Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Geomorphological hazards related to deep dissolution phenomena in the Western Italian Alps: Distribution, assessment and interaction with human activities

Deep dissolution affects great part of soluble rocks (e.g. gypsum and anhydrite) of the Western Italian Alps. The related superficial phenomena (sinkholes, gravity-induced processes and a local worsening of geomechanical rock properties) are not limited to typical karsts landscape and cause slope instability also affecting populated sites and infrastructures. The paper aims to describe general characteristic of dissolution phenomena, to interpret their conditioning factors and evolutionary stages and to assess possible hazards due to their superficial effects.The search for evidences of deep dissolution leads to the selection of representative sites in the central part of the Western Italian Alps (Piemonte and Valle d'Aosta Region). Detailed geological and geomorphological studies have been used to classify the selected sites by type, size and variable state of activity. Very different evolutionary stages of dissolution phenomena have been interpreted by comparison of case-studies: some are early “embryonic”; others are more evolved, up to typical sinkholes, or even remodelled by other phenomena. Some cases show an extreme complexity in the interactions between corrosion phenomena and other geomorphic processes: slope deformations, from one side, and karst, fluvial and glacial phenomena, to the other. A wide range of movement rates on slope instabilities induced by deep dissolution have been estimated by topographic and geomorphic data. Geochemical data on removed rocks by dissolution indicate 0.4 mm/year values for local subsidence. Historical and technical data indicate low frequency of major dissolution-induced collapses, but highlight widespread damages to tunnels, roads and buildings, especially along slopes.     

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.