17 Ağustos 1999 İzmit depremi sonrası Marmara bölgesindeki gravite alanın zaman içindeki değişiminin belirlenmesi

Abstract

Bu çalışmanın amacı; Marmara Bölgesi'nde 17 Ağustos 1999 İzmit Depremi sonucunda oluşan zaman bağımlı gravite alan değişimlerinin ölçülmesi ve analiz edilmesidir.Amaç doğrultusunda kullanılan gravite alan verileri, ?Türkiye'nin Deprem Riski Yüksek Jeo-Stratejik ? ancak tektonik rejimleri farklı ? Bölgelerinde Deprem Davranışının Çok Disiplinli Yaklaşımlarla Araştırılması - TÜRDEP? ve ?Marmara Bölgesi'ndeki Düşey Yerkabuğu Hareketlerinin Mutlak Gravite ve GPS ile Araştırılması? projeleri kapsamında 2006-2011 yılları arasında gerçekleştirilen zamansal bağıl gravite kampanya ölçümleri sonucu elde edilmiştir.Ölçümler, drift hatalarının kontrolü için gidiş-dönüş ölçüm tekniği temel alınarak gerçekleştirilmiştir. Ölçümler de gravite alan değerlerine katılan aletsel ve çevresel etkiler, gravite yöntemi kuramında yer alan ve çalışmanın amacına uygun olarak seçilen düzeltmeler ve indirgemeler ile ölçüm değerlerinden çıkarılmıştır. Düzeltme ve indirgeme işlemleri sonrasında farklı günlerde ölçülen geçkilerin, birbirleriyle birlikteliklerinin sağlanması ve belli bir düzeye indirgenmeleri için dengelenmiştir. Düzeltme, indirgeme ve dengeleme işlemleri; GRAVAP yazılımı kullanılarak gerçekleştirilmiştir. Kampanyalar kapsamında belirlenen gravite alan değerlerinden, her bir ölçüm noktası için gravite zaman serileri oluşturulmuştur. Elde edilen zaman serilerine Ağırlıklı En Küçük Kareler Yöntemi uygulanmış ve yıllık gravite hız değerleri elde edilmiştir. Elde edilen gravite verilerinin, viskoelastik modelleme ile analiz edilmesine karar verilmiştir. Bu amaçla, yeryüzünde oluşan zaman bağımlı deformasyonların kaynağını oluşturan jeoid ve gravite alan değişimlerini hesaplayan bir sayısal model oluşturulmuştur.. Sayısal modelleme yaklaşımı ile istenilen duyarlılıkta ve optimum zamanda gerçekleştirilecek modelleme işlemi için uygun değişkenlerin belirlenmesi ve gerçekleştirilecek modelleme çalışmasının doğruluğu ve güvenilirliğinin tespit edilmesi önemlidir. Bu amaç doğrultusunda; referans analitik çözümler ve sayısal çözüm sonuçları karşılaştırılmıştır. Karşılaştırma sonucunda, sayısal modelin doğruluğu ve güvenilirliği için gerekli çalışma değerleri tespit edilmiştir.Modelleme işleminde öncelikle; Marmara Bölgesi'nin gravite alan değişimleri temel alınarak, sünek yapıya sahip alt kabuk ve üst mantonun viskozite değerinin belirlenmesine yönelik testler yapılmıştır. Testler sonucunda yarısonsuz viskoelastik ortam için viskozite değeri 5x1018 Pa.s civarı olarak belirlenmiştir. Modelleme işlemi sonrasında; Marmara-doğu ve Marmara-batı bölgelerinde model ve ölçülen gravite alan değişimleri arasındaki uyum incelenmiştir. İnceleme sonrasında Marmara-doğu bölgesinde, gravite değişimleri yüksek uyumla modellenebilmiş ve gravite değişime neden olan olaylar tanımlanabilmiştir.. Marmara-batı bölgesinde ise uzun dönem değişimler modellenebilmiş ama modelden elde edilen genlik değerleri ölçülen veriye göre düşük kalmıştır. Marmara-batı bölgesindeki model ve ölçülen veri arasındaki genlik farkının; salt 1999 depremleri sonucu oluşan deformasyonlar ile modellenmesinin yetersiz olmasından, bu bölgedeki kabuk yapısındaki yanal değişimlerden kaynaklandığı öngörülmüştür.

The objective of this study is the measurement and analysis of the time dependent microgravity area changes occurring as a result of the August 17, 1999 İzmit Earthquake in the Marmara Region.The microgravity data for this purpose were obtained by time dependent campaign measurements accomplished between the years 2006-2011 within the scope of the TUBITAK, MRC, Earth and Marine Sciences projects ?TÜRDEP - multi-disciplinary earthquake research in high risk geostrategic regions of Turkey with different tectonic regimes and investigation of vertical displacements of the earth?s crust in the Marmara Region by absolute gravity and GPS?.Round trip measurement technique was employed in order to control the drift errors. Instrumental and environmental effects that may be included in the measurement values are corrected and reduced according to the microgravity area theory and the scope of this study. After the correction and reduction operations, scaling is also performed in order to associate measurements performed at different days. Correction, reduction and scaling operations are accomplished using the GRAVAP software. Microgravity time series were created for each measurement location from the microgravity area values determined during the campaign.Weighted least squares regression was applied to the time series data obtained and annual micro gravity velocity values were obtained. It was decided to analyze the obtained microgravity data by viscoelastic modeling. For this purpose, a numerical model was created to compute the changes in the microgravity area and the geoid, constituting the source of time dependent deformations occurring on the earth. By the numerical modeling approach, it is significant to determine the suitable modeling parameters at the optimum time with the required accuracy. For this purpose, reference analytic solutions were compared against the numerical modeling solutions. The results showed that the numerical model is accurate and reliable.Modeling operation was primarily done in the Marmara Region based on microgravity area changes and tests were performed in order to determine the viscosity values of the low crust of ductile structure and the upper mantle. After the tests, the viscosity value was determined as approximately 5x1018 Pa.s for semi-infinite viscoelastic medium. After the modeling operation, the harmony was investigated between the model and the measured microgravity area changes in the Eastern Marmara and the Western Marmara regions. After this investigation, microgravity changes could be modeled in high harmony in the Eastern Marmara region and the events could be identified inducing the microgravity changes. Long term changes could be modeled in the western Marmara region but the amplitude values obtained from the model were lower than the measured data. It is expected that the amplitude difference between the model and the measured data in the Western Marmara region is either a result of insufficient modeling using only the deformations occurring after the 1999 earthquake or the lateral changes in the crustal structure of this region.