Applications of hybrid genetic algorithms in seismic tomography

Almost all earth sciences inverse problems are nonlinear and involve a large number of unknown parameters, making the application of analytical inversion methods quite restrictive. In practice, most analytical methods are local in nature and rely on a linearized form of the problem equations, adopting an iterative procedure which typically employs partial derivatives in order to optimize the starting (initial) model by minimizing a misfit (penalty) function. Unfortunately, especially for highly non-linear cases, the final model strongly depends on the initial model, hence it is prone to solution-entrapment in local minima of the misfit function, while the derivative calculation is often computationally inefficient and creates instabilities when numerical approximations are used. An alternative is to employ global techniques which do not rely on partial derivatives, are independent of the misfit form and are computationally robust. Such methods employ pseudo-randomly generated models (sampling an appropriately selected section of the model space) which are assessed in terms of their data-fit. A typical example is the class of methods known as genetic algorithms (GA), which achieves the aforementioned approximation through model representation and manipulations, and has attracted the attention of the earth sciences community during the last decade, with several applications already presented for several geophysical problems.In this paper, we examine the efficiency of the combination of the typical regularized least-squares and genetic methods for a typical seismic tomography problem. The proposed approach combines a local (LOM) and a global (GOM) optimization method, in an attempt to overcome the limitations of each individual approach, such as local minima and slow convergence, respectively. The potential of both optimization methods is tested and compared, both independently and jointly, using the several test models and synthetic refraction travel-time date sets that employ the same experimental geometry, wavelength and geometrical characteristics of the model anomalies. Moreover, real data from a crosswell tomographic project for the subsurface mapping of an ancient wall foundation are used for testing the efficiency of the proposed algorithm. The results show that the combined use of both methods can exploit the benefits of each approach, leading to improved final models and producing realistic velocity models, without significantly increasing the required computation time.     

Automated 1D interpretation of resistivity soundings by simultaneous use of the direct and iterative methods*

Noise contamination of measured data greatly affects the final results of inversion. Three types of noise source — random and systematic errors and the uncertainties due to the inadequacy of the mathematical model in representing the actual physical conditions — are discussed in the framework of resistivity sounding data. Two methods are proposed for describing these uncertainties. The first possibility is to smooth the measured data by a combination of simple fitting functions that satisfies the ‘1D smoothness’ criteria and consequently simulates the behaviour of a 1D Schlumberger apparent resistivity curve. The second method is to derive weight coefficients from the differences between the measured and the smoothed data sets. Both methods are carried out under the control of the interpreter. The relative merits and drawbacks of the direct and iterative interpretation methods used for the estimation of the parameters of the layered earth model are summarized. Two variants of the combination of these methods are presented to obtain more powerful and automatic interpretation schemes. In the sequential interpretation, an initial guess supplied by the direct method is improved by the iterative method to obtain a reasonable fit between the measured data and the model response. In the simultaneous interpretation, the successive application of the direct and iterative methods is carried out, starting from the first branch of the apparent resistivity curve. The operation is then shifted to subsequent branches that represent the deeper parts of the geoelectric section. This is similar to the data acquisition applied in direct current sounding in which the depth penetration is increased by expanding the current electrode spacings. The proposed sequential and simultaneous interpretation algorithms require minimum aids and efforts of the interpreter.     

Ambient vibration based seismic evaluation of isolated Gülburnu highway bridge

This paper describes ambient vibration based seismic evaluation procedure of an isolated highway bridge. The procedure includes finite element modeling, ambient vibration testing, finite element model updating and time history analysis. Gülburnu Highway Bridge located on the Giresun–Espiye state highway is selected as a case study. Three dimensional finite element model of the bridge is created by SAP2000 software to determine the dynamic characteristics analytically. Since input force is not measured, Operational Modal Analysis is applied to identify dynamic characteristics. Enhanced Frequency Domain Decomposition and Stochastic Subspace Identification methods are used to obtain experimental dynamic characteristics. Analytical and experimental dynamic characteristic are compared with each other and finite element model of the bridge is updated by changing of material properties to reduce the differences between the results. Analytical model of the bridge after model updating is analyzed using 1992 Erzincan earthquake record to determine the seismic behavior. EW, NS and UP components of the ground motion are applied to the bridge at the longitudinal, transverse and vertical directions, respectively. It is demonstrated that the ambient vibration measurements are enough to identify the most significant modes of highway bridges. Maximum differences between the natural frequencies are reduced averagely from 9% to 2% by model updating. It is seen from the earthquake analyses that friction pendulum isolators are very effective in reducing the displacements and internal forces.     

Mapping of total nitrogen,available phosphorous and potassium in Amik Plain,Turkey

Soil nitrogen, phosphorous, and potassium concentrations accurately revealed spatial distribution maps and site-specific management-prone areas through inverse distance weighting (IDW) method in the Amik Plain, Turkey. Spatial mapping of soil nitrogen, phosphorous, and potassium is a very severe need to develop an economically and environmentally sound soil management plans. The objectives of this study were (a) to map spatial variability of total N, available P, and exchangeable-K content of Amik Plain’s soils and (b) to locate problematic areas requiring site specific management strategies for the nutrient elements. Spatial analyses of Kjeldhal-N, Olsen-P, and exchangeable-K concentrations of the soils were performed by the IDW method. Mean N content for surface soils (0–20 cm) was 1.38 g kg⁻¹, available P was 28.19 kg ha⁻¹ and exchangeable-K was 690 kg ha⁻¹ with the differences between maximum and minimum being 7.63 g N kg⁻¹, 242 kg P ha⁻¹, and 2,082 kg K ha⁻¹. For the surface soil, site-specific management-prone areas of Kjeldahl-N, Olsen-P, and exchangeable-K for “low and high + very high” classes were found to be 20.1–17.8%, 24.7–10.0%, and 4.1–39.6%, respectively. Consequently, lands with excessive nutrient elements require preventive-leaching practices, whereas nutrient-poor areas need fertilizer applications in favor of increasing plant production.     

Fossil findings from the Sıcak Çermik fissure ridge-type travertines and possible hominid tracks,Sivas, Central Turkey

S?cak Çermik (Sivas) is an important geothermal and recent travertine formation area in Central Anatolia. The majority of travertines found in the region comprise fissure-ridge type travertines according to morphological classification. At the location called Tepe Çermik within the travertine area, fill containing fossil bone fragments of Equus sp., Bovidae and other abundant animals formed within the fracture axis of a N–S striking fissure-ridge travertine developed under control of tectonic forces. The finds of these fossils in fissure-ridge travertines linked to tectonic forces indicates formation of a unique fossil environment created under the control of these forces. The Accelerator Mass Spectrometry Radiocarbon Dating analyses of fossils from the study area determined the fills were older than 43,000 years. The U/Th age of a sample from the most recently-formed banded travertine in the axis of the fracture was identified as 278,540 ± 18,436 years. As a result, the ages of fossils found within this fill are thought to be between 43,000 and 278,540 ± 18,436 years old. The high amount of perissodactyla and artiodactyla fossils found within fill in the axis of the fissure-ridge travertine probably indicates the presence of hominids who chose the region for hunting or settlement. The Equus sp. and Bovidae fossil samples found in the axis of the fracture indicate that in the dry and cold glacial period the paleogeography in a large portion of Anatolia comprised desert-like steppe.     

Use of interactive multisensor snow and ice mapping system snow cover maps (IMS) and artificial neural networks for simulating river discharges in Eastern Turkey

Basins located in Eastern Turkey are largely fed by snowmelt runoff during spring and early summer seasons. This study investigates the efficiency of artificial neural networks (ANNs) in snowmelt runoff generation. Although ANNs have been used for streamflow simulating/forecasting in the last two decades, using satellite-based snow-covered area (SCA) maps and meteorological observations as inputs to ANN provides a novel basis for estimating streamflow. The proposed methodology is implemented over Upper Euphrates River Basin in Eastern Turkey. SCA data was acquired from Interactive Multisensor Snow and Ice Mapping System (IMS) for an 8-year period from February 2004 to September 2011. Meteorological observations including daily cumulative precipitation and daily average air temperatures were obtained from Turkish State Meteorological Services. The simulation results are promising with coefficient of correlation varying from 0.67 to 0.98 among proposed models. Past days discharge was found to substantially improve the forecast accuracy. The paper presents the expected basin discharge for 2011 water year based on meteorological observations and SCA input.