توموگرافی دوبعدی سرعت امواج حجمی فشارشی در پوستة فوقانی با استفاده از زمین‌لرزه‌های محلی شمال باختر ایران

نوع مقاله : مقاله پژوهشی‌

نویسندگان

1 دانشکده علوم زمین، دانشگاه تحصیلات تکمیلی علوم پایه، زنجان، ایران

2 پژوهشکده زلزله‌شناسی، پژوهشگاه بین‌المللی زلزله‌شناسی و مهندسی زلزله، تهران، ایران

3 گروه فیزیک زمین، موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران

چکیده

در این پژوهش با استفاده از زمان‌سیر امواج حجمی P ثبت شده از زمین‌لرزه‌های محلی، توزیع دوبعدی ساختار سرعتی در پوسته بالایی شمال باختر ایران تعیین شده است. وجود خطای زیاد در تعیین عمق رویدادها سبب می‌شود تا تنها تغییرات جانبی سرعت بررسی شود و از تغییرات عمقی آنها صرف‌نظر گردد. برای این مطالعه، از چندین مجموعه داده شامل داده‌های ایستگاه‌های دائمی شبکه تبریز، وابسته به مؤسسه ژئوفیزیک دانشگاه تهران (1996 تا 2013 میلادی)، داده‌های ایستگاه‌های موقت نصب شده توسط پژوهشگاه بین‌المللی زلزله‌شناسی و مهندسی زلزله (در بازه زمانی آوریل تا جولای 2004 میلادی) و داده‌های ایستگاه‌های موقت نصب شده توسط دانشگاه تحصیلات تکمیلی علوم پایه زنجان (2009 تا 2011 میلادی) استفاده شده است. از کل زمین‌لرزه‌های ثبت شده، تنها 940 زمین­لرزه برای توموگرافی انتخاب شدند. در این پژوهش، ابتدا یک مدل اولیه سرعت به‌همراه یک برآورد اولیه برای ضریبهای زمین لرزه‌ها (شامل مکان و زمان وقوع زمین‌لرزه) تعیین و سپس مقدار تأخر و یا تقدم زمانی فازهای دریافتی توسط ایستگاه‌ها نسبت به‌مقدار پیش‌بینی شده توسط مدل اولیه محاسبه شده است. این مقادیر توسط برنامه SIMULPS14 به مدل سرعتی دوبعدی برگردانده شده و یک تفاوت آشکار در ساختار سرعتی در دو طرف گسل شمال تبریز را آشکار نموده است. علاوه بر آن، یک بی‌هنجاری کم‌سرعت در قسمت مرکزی گسل شمال تبریز مشاهده می‌شود که با ضخامت زیاد رسوبات در آن ناحیه متناسب است.

کلیدواژه‌ها


عنوان مقاله [English]

2D local earthquake tomography of P-wave velocity in the upper crust of NW Iran

نویسندگان [English]

  • Farzam Feizaghaee 1
  • Khalil Motaghi 1
  • Mohammad Tatar 2
  • Abdolreza Ghods 1
  • Ali Moradi 3
1 Earth Sciences Department, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan, Iran
2 Seismological Research Center, International Institute of Earthquakes Engineering and Seismology (IIEES), Tehran, Iran
3 Institute of Geophysics, University of Tehran, Tehran, Iran
چکیده [English]

NW Iran is part of the complex tectonic system caused by the interaction between Arabian plate, Anatolia and Eurasia. The North Tabriz Fault (NTF) is one of the main structural features of the region and is considered  to be the eastern termination of the Gailatu-Siah-Chesmeh-Khoy fault (Karakhanian et al., 2004), which merges with the Maku and the Nakhichevan dextral strike-slip faults and continues to move farther east. Part of the northward motion of Arabia is transferred to Anatolia by this complex system of faults and, the oblique orientation of the motion relative to the Zagros mountain range,results in the partitioning of the motion between shortening and thickening in the Caucasus and right-lateral strike-slip motion along the NTF.
In this research, we investigated the laterally two-dimensional velocity structure in the upper crust of NW Iran (mainly around the NTF) using local earthquake P-waves tomography. Several data sets were utilized, including Pg phase pickings of the Tabriz Network permanent stations governed by Institute of Geophysics, University of Tehran (1996 to 2013), temporary seismic stations installed around the North Tabriz Fault by International Institute of Earthquake Engineering and Seismology (IIEES) (April to July 2004) and temporary seismic stations installed by Institute for Advanced Studies in Basic Sciences (IASBS) (2009 to 2011) that merged  with our data set so as to improve ray coverage in the eastern parts of the study area. The merged data set, recorded by 72 stations, consisted of more than 20,000 local earthquakes out of which, only 940 earthquakes were good enough to be selected for the local earthquake tomography. The velocity structures were resolved via a simultaneous solution of the coupled hypocenter and velocity model programmed in SIMULPS14. The tomographic images obtained from the linearized inversion are dependent on the initial velocity models and hypocenter locations. We primarily calculated the initial velocity model through the use of the 940 earthquake selected datasets. The time difference between the observed phase arrival time and predicted arrival time was then calculated and called travel time residuals. The residuals were further used as inputs for SIMULPS14 simulator to be converted into velocity model, which would in turn be used to adjust earthquake location parameters. Following four iterations for our inversion process, we obtained a 2D velocity tomogram that clearly showed different velocity structures on the two sides of the NTF. The velocity contrast across the NTF might have been caused by existence of different kinds of rocks on the two side of the fault trace. The North Tabriz Fault is an active and steep strike-slip fault generating strong structural differences around its surficial trace. It is a WNW–ESE trending fault in which the motion is concentrated on the fault at a rate of 7 mm/year. Such a strong rate of sliding explains the clear structure difference on the two sides of the fault. An anomalous low velocity feature can be seen in the central part of NTF. Comparing the velocity tomogram with the geological map of the region, one can observe that there exists a thick sediment basin in the same area. The low velocity anomaly is probably related to the thick, low velocity sediments deposited in that area.

کلیدواژه‌ها [English]

  • Local Earthquake Tomography (LET)
  • velocity structure
  • SIMULPS14
  • North Tabriz Fault (NTF)
  • Vp/Vs Ratio
Al-Lazki, A. I., Sandvol, E., Seber, D., Barazangi, M., Turkelli, N., and Mohamad, R., 2004, Pn tomographic imaging of mantle lid velocity and anisotropy at the junction of the Arabian, Eurasian and African plates: Geophysical Journal International, 158, 1024–1040.
Berberian, M., 1976, Contribution to the seismotectonics of Iran (partІІ): p. 518, Geological Survey of Iran, Tehran.
Berberian, M. and Yeats, R. S., 1999, Pattern of historical earthquake rupture in the Iranian Plateau: Bulletin of the Seismological Society of America, 89, 120-139.
Djamour, Y., Vernant, P., Nankali, H.R., and Tavakoli, F., 2011, NW Iran-eastern Turkey present-day kinematics: results from the Iranian permanent GPS network: Earth and Planetary Science Letters, 307 (1), 27–34.
Eberhart-Phillips, D., 1986, 3-D velocity structure in Northern California Coast Ranges from inversion of local earthquake arrival times: Bulletin of the Seismological Society of America, 76, 1025-52.
Eberhart-Phillips, D., 1990, Three dimensional P and S velocity structure in Coalinga region, California: Journal of Geophysical Research, 95:15, 15343-15363.
Geiger, L., 1910, Herdbestimmung bei Erdbeben aus den Ankunftszeiten. Nachr. Koninglichen Gesell: Wiss. Gottingen Mathematical Physics, 4, 331.
Hasslinger, F., 1999, Velocity structure, seismicity and seismotctonics of Northwestern Greece, between the Golf of Arta and Zakinthos: Ph. D. Thesis, diss. ETH. 12966, Zurich, P.159.
Hessami, K., Tabassi, H., Abbassi, M. R., Azuma, T., Okumura, K., Echigo, T., and Kondo., H., 2003, Surface expression of the Bam Fault Zone in southeastern Iran: Causative Fault of the 26 December 2003 Bam Earthquake: Journal of Seismology and Earthquake Engineering, 4, 5-14.
Maggi, A. and Priestley, K., 2005, Surface waveform tomography of the Turkish-Iranian plateau: Geophysical Journal International, 160, 1068-1080.
Menke, W., 1989, Geophysical data analysis: Discrete inverse theory: Academic Press, New York.
Moradi, A. S., Hatzfeld, D. and Tatar, M., 2011, Microseismicity and seismotectonics of the North Tabriz fault (Iran): Tectonophysics, 504, 22–30.
Sandvol, E., Al-Damegh, K., Calvert, A., Seber, D., Barazangi, M., Mohamad, R., Gok, R., Turkelli, N., and Gurbuz, C., 2001, Tomographic imaging of Lg and Sn propagation in the Middle East: Pure and Applied Geophysics 158, 1121–1163.
Stein, S., and Wysession, M., 2003, Introduction to Seismology, Earthquakes, and Earth Structure: Blackwell, Boston.
Thurber, C. H., 1983, Earthquake locations and three-dimensional crustal structure in the Coyote Lake area, central California: Ph. D. Thesis, Massachusetts Institute of Technology, Cambridge, MA.
Thurber, C. H., 2003, Seismic Tomography of the Lithosphere with Body Waves: Pure and Applied Geophysics, 160, 717-737.
Um, J., and Thurber, C. H., 1987, A fast algorithm for two-point seismic ray tracing: Bulletin of the Seismological Society of America, 77, 792-86.