Environmental hazard analysis of a gypsum karst depression area with geophysical methods: a case study in Sivas (Turkey)

Because the mixture of seawater and freshwater in the Gyeongin-Ara Waterway in South Korea can lead to the intrusion of saline water into surrounding aquifers, systematic management through the establishment of a groundwater protection area is required. The analytic hierarchy process (AHP) model is used to delineate this protection area based on two primary factors and five secondary factors related to saline water intrusion. The study area is divided into 987 gridded cells with a unit size of 100 × 100 m, and the final evaluation score for each cell is calculated using the AHP model. Consequently, several artificial neural network models based on a multilayer perceptron are developed using the AHP’s secondary criteria and the evaluation score. Comparing the evaluation scores of ANN and AHP, more than 180 samples are required in the ANN model to insure high R2 between the original and estimated values. The ANN model is more consistent than the AHP model when determining groundwater protection area, because it can be re-constructed due to the changes in some secondary criteria and also changed due to a standardization process. The final evaluation score by the ANN model based on 300 samples, with the highest R2, is calculated and the regions with a score higher than 2.0 are selected as the groundwater protection area, accounting for 15% of the total cells. This area is similar to the range within approximately 200 m of the GA Waterway and also includes some changing sites in hydrogeochemistry and electric conductivity, which is produced by saline water intrusion. If the land-use type, groundwater levels, and some other criteria change at any cell, the ANN model can be re-executed to verify whether the cell belongs to a groundwater protection area. Considering that salinity of groundwater near the waterway can be affected by various factors including well depth, pumping conditions, and groundwater levels, the ANN model, which is a non-linear model, can be more effective for prediction than the AHP model.     

Hydrogeophysical modelling of Hisarcik (Kütahya) geothermal field,western Turkey

Western Anatolia hosts many low-to-moderate and high-temperature geothermal sources in which active faults play a dominant role to control the recharge and the discharge of geothermal fluid. In this study, we used the two-dimensional geoelectric structure of Kütahya Hisarcık geothermal field, and created a conceptual hydrogeophysical model that includes faults, real topographical variations and geological units. The temperature distribution and fluid flow pattern are also investigated. The depth extension of Hisarcık Fault, electrical basement and low resistivity anomalies related to the presence of geothermal fluid are determined by using resistivity studies in the area. Numerical simulations suggest that Hisarcık fault functioning as a fluid conduit primarily enables hot fluid to be transported from depth to the surface. It is shown that the locations of predicted outflow vents coincide with those of hot springs in the area.     

An improved input energy spectrum verified by the shake table tests

Input energy is the principal component of the energy balance equation. It is beneficial to determine, through its components, how the recoverable and irrecoverable energies are distributed within the structural elements. Several equations and attenuation relations to define mass-normalized input energy spectra exist in the literature. They are mainly proposed for elastic systems subjected to far-fault EQs. There is a lack of experimental verification of these proposed spectra. In this paper, experimental assessment was performed to the existing spectra, and further improvements were accomplished. For this purpose, steel cantilever columns were tested on the shake table for two specific historical EQs coincidently having similar spectral acceleration values. Based on the experimental results, a three-part mass-normalized relative input energy spectrum was formulated including soil type, EQ (corner period, intensity, duration, spectral acceleration, and velocity), and structural behavioral characteristics (period and structural damping). The proposed input energy spectrum was experimentally calibrated and numerically validated for various EQs featuring near- and far-field types. Analytical and experimental comparisons were made between the previously developed spectrum and the newly proposed one. The validation studies and the statistical evaluations exposed that the proposed spectrum yielded better agreement with the experimental and numerical results.     

Evaluation of the physico-mechanical parameters affecting the deterioration rate of Ahlat ignimbrites (Bitlis,Turkey)

The paper principally focuses on the durability assessment of various stratigraphic levels of Ahlat ignimbrites collected from the eastern region of Turkey. A total of four different ignimbrite types with dissimilar color, texture and particularly welding degree were tested in laboratory. The laboratory tests performed on the ignimbrite specimens indicate that the welding degree as well as the lithic material content mainly controls the strength and capillarity properties of the ignimbrites. In addition, the durability of highly porous ignimbrites strongly depends upon the degree of welding. The effect of several weathering agents on the ignimbrites was evaluated on the basis of decay constant parameter. Accordingly, salt and ice crystallization pressures are a couple of major destructive agents acting within the micropores of the ignimbrites. Conversely, the investigated specimens are relatively durable against cyclic wetting–drying. Statistical evaluations reveal that the pore diameter is the major controlling factor on the deterioration rate of the ignimbrites after specifically recurrent freeze–thaw cycles. Moreover, the dry unit weight of the ignimbrites is more significant than the uniaxial compressive strength considering the deterioration rates during wetting–drying and salt crystallization. A less significant relationship was obtained between pore diameter and salt crystallization decay constant.     

Hydrogeochemical evaluation of Umut geothermal field (SW Turkey)

Tekkehamam geothermal field is located in the South of Menderes Graben (Aegean region) and is one of the most important geothermal sites of Western Anatolia. Umut geothermal field is a part of the Tekkehamam field. This study was conducted in order to determine the origin and hydrogeochemical properties of the geothermal waters. For this purpose, sampling was done in order to check the chemistry of the water, and ¹⁸O, ²H isotope analyses done at four wells, nine natural springs and three cold water sources. According to the results of the chemical analysis, the geothermal waters were determined to be of Na + K-SO₄ type. Additionally, ¹⁴C and ³H analyses were done in selected well and spring waters for the purpose of age determination of groundwater; most of the waters were determined to be submodern. Geothermometer calculations show that the reservoir temperature for the Umut geothermal field ranges between 148 and 180 °C. Stable isotope results indicate that Umut geothermal waters are meteoric in origin. Mixing between shallow and deep waters is the dominant subsurface process that determines the physical and chemical character of the waters.     

Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

Evaluation of a recently proposed record selection and scaling procedure for low‐rise to mid‐rise reinforced concrete buildings and its use for probabilistic risk assessment studies

This paper evaluates a recent record selection and scaling procedure of the authors that can determine the probabilistic structural response of buildings behaving either in the elastic or post‐elastic range. This feature marks a significant strength on the procedure as the probabilistic structural response distribution conveys important information on probability‐based damage assessment. The paper presents case studies that show the utilization of the proposed record selection and scaling procedure as a tool for the estimation of damage states and derivation of site‐specific and region‐specific fragility functions. The method can be used to describe exceedance probabilities of damage limits under a certain target hazard level with known annual exceedance rate (via probabilistic seismic hazard assessment). Thus, the resulting fragility models can relate the seismicity of the region (or a site) with the resulting building performance in a more accurate manner. Under this context, this simple and computationally efficient record selection and scaling procedure can be benefitted significantly by probability‐based risk assessment methods that have started to be considered as indispensable for developing robust earthquake loss models. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of river channel morphology and bank characteristics on water surface boundary delineation using high‐resolution passive remote sensing and template matching

Accurate mapping of water surface boundaries in rivers is an important step for monitoring water stages, estimating discharge, flood extent, and geomorphic response to changing hydrologic conditions, and assessing riverine habitat. Nonetheless, it is a challenging task in spatially and spectrally heterogeneous river environments, commonly characterized by high spatiotemporal variations in morphology, bed material, and bank cover. In this study, we investigate the influence of channel morphology and bank characteristics on the delineation of water surface boundaries in rivers using high spatial resolution passive remote sensing and a template‐matching (object‐based) algorithm, and compare its efficacy with that of Support Vector Machine (SVM) (pixel‐based) algorithm. We perform a detailed quantitative evaluation of boundary‐delineation accuracy using spatially explicit error maps in tandem with the spatial maps of geomorphic and bank classes. Results show that template matching is more successful than SVM in delineating water surface boundaries in river sections with spatially challenging geomorphic landforms (e.g. sediment bar structures, partially submerged sediment deposits) and shallow water conditions. However, overall delineation accuracy by SVM is higher than that of template matching (without iterative hierarchical learning). Vegetation and water indices, especially when combined with texture information, improve the accuracy of template matching, for example, in river sections with overhanging trees and shadows – the two most problematic conditions in water surface boundary delineation. By identifying the influence of channel morphology and bank characteristics on water surface boundary mapping, this study helps determine river sections with higher uncertainty in delineation. In turn, the most suitable methods and data sets can be selectively utilized to improve geomorphic/hydraulic characterization. The methodology developed here can also be applied to similar studies on other geomorphic landforms including floodplains, wetlands, lakes, and coastlines. Copyright © 2014 John Wiley & Sons, Ltd.     

An investigation of seismicity for the Central Anatolia region,Turkey

The aim of this study is to investigate the seismicity of Central Anatolia, within the area restricted to coordinates 30–35° longitude and 38–41° latitude, by determining the “a” and “b” parameters in a Gutenberg–Richter magnitude–frequency relationship using data from earthquakes of moment magnitude (Mw)?≥?4.0 that occurred between 1900 and 2010. Based on these parameters and a Poisson model, we aim to predict the probability of other earthquakes of different magnitudes and return periods (recurrence intervals). To achieve this, the study area is divided into six seismogenic zones, using spatial distributions of earthquakes greater than Mw?≥?4.0 with active faults. For each seismogenic zone, the a and b parameters in the Gutenberg–Richter magnitude–frequency relationship were calculated by the least squares method. The probability of occurrence and return periods of various magnitude earthquakes were calculated from these statistics using the Poisson method.