3D Tomography of Isfahan and Sharekord regions using local earthquakes


1 M.Sc. Student, Earth Physics Department, Institute of Geophysics, University of Tehran, Iran

2 Assistant Professor, Earth Physics Department, Institute of Geophysics, University of Tehran, Iran

3 Professor, Earth Physics Department, Institute of Geophysics, University of Tehran, Iran


*نگارنده رابط:          تلفن: 61118230-021         دورنگار: 88630479-021                               E-mail: ebayram@ut.ac.ir


One of the most important purposes in seismology is determination of crustal velocity using earthquakes data that have been recorded by regional and local seismic stations. The more precise the process is carried out, the better shall be the results reached in various studies in the area including earthquake locating, seismicity of area, seismic zone mapping or the determination of the plane of faults causing earthquakes Imaging velocity structure of Earth's interior using travel times inversion commonly called seismic tomography which is usually done two or three-dimensionally. . This method has extensively been used in recent decades by researchers. Since the seismic waves are associated with valuable information about direction of propagation and environmental properties, using earthquakes as natural seismic sources are very useful in seismic tomography as well as the artificial sources such as limited and controlled explosions, air guns and bore-hole sources.. The seismic tomography characterizes the size, geometry and extent of velocity anomalies. In this method subsurface structure is modeled initially by several parameters and improved by the inversion of seismic travel times data.
In this study the crustal velocity structure was determined using three-dimensional inversion of local earthquakes travel times recorded by seismic networks of Institute of Geophysics University of Tehran (IGUT) occurred within the period from 2000 and 2012 in the study area. The study area is bounded within 31oE to 34oE and 50oN to 53oN. We used the VELEST software in one-dimensional modeling section. The procedure of study is minimizing the differences between observed and calculated travel time by applying the initial obtained model. This software simultaneously optimizes the earthquake locations, crustal velocity model and station corrections using the Joint-Hypocenter-Determination (JHD) method. 
The initial estimates for P waves velocity and crust thickness of the region are achieved using the travel-time curve of primary phases of all earthquakes occurred in the area. Then the larger relative earthquakes are selected and the best crustal one-dimensional model was derived by simultaneously inverse modeling method using this data set and VELEST algorithm. This method can be considered as one of the useful methods in study of 1D crust structure. The results proposed a 4 layers model of crust in which the P wave velocity is equal to 5.4 km/s for depths less than 5 km,   6.0 km/s for depths from 5 km to 20 km,  6.2 km/s for depths of 20 km to 32 km and finally  6.9 km/s for depths of 32 km to 47 km. The thickness of crust and Pn velocity are respectively obtained 47 km and 7.9 km/s.  The aim of this work is obtaining an optimal crust model that can aid to improve seismic data and can be used to determine the next earthquake locating. Then the obtained crustal model is used as an initial model to study of three-dimensional inverse modeling of crust in the region by using FAST algorithm. All the earthquakes relocated using new obtained model. In this study a data set recorded by the 8 seismic stations of Isfahan and Shahrekord networks were used. The resolution of final solution of 3D model was investigated using synthetic dataset (checkerboard model) that shows fair resolution for various depths. The lateral variations of the main resolved structures in the model obtained are highly correlated with the faulting systems in the region.