<strong>Harita</strong> <strong>Dergisi</strong> Temmuz 2010 Sayı 144A.İ.KURT vd.2. GPS VERİLERİNDEN DOĞRUSALTEKTONİK HIZIN HESAPLANMASI<strong>Harita</strong> <strong>Genel</strong> <strong>Komutanlığı</strong> tarafındanoluşturulan TUTGA noktalarında yapılan GPSölçüleri BERNESE V4.0 ve GLOBK V10.03yazılımları ile değerlendirilmiştir. GPSoturumlarındaki ölçülerin değerlendirilmesi veyıllık çözümler BERNESE V4.0 yazılımı, buçözümlerin birleştirilmesi, nokta koordinat vehızlarının hesaplanması ise GLOBK V10.03yazılımı ile yapılmıştır. TUTGA’yı ITRF referanssisteminde tanımlamak amacıyla Türkiye veyakın çevresinde Avrupa ve Asya’da uygundağılımdaki dokuz IGS noktası hesaplamalaradâhil edilmiştir (TUTGA 99A).Doğrusal tektonik hız, levhaların birbirine görehareketi ile oluşan yatay ve düşey koordinatfarklılıklarının zaman (yıl) ile orantılanmasındanoluşmaktadır. Ayrıca bölgesel olarak yapılan birdeformasyon çalışmasında, levhahareketlerinden farklı olarak ölçülen (kayma,çökme vb.) koordinatlar arasındaki farklarınzamanla orantılanmasından, deformasyon hızıbelirlenmektedir. <strong>Genel</strong> olarak GPS verideğerlendirme yazılımlarında koordinatlarıetkileyen bazı fiziksel etkiler otomatik olarak veyaseçime bağlı olarak hesaplamalarda dikkatealınmaktadır. Örneğin; atmosferik modeller, yerdönme parametreleri, gel-git modelleri (yer,kutup, okyanus, atmosfer), yörünge denklemleriiçin radyasyon modeli, okyanus yüklemesi veatmosferik yükleme bu etkileri oluşturmaktadır.Yukarıdaki etkiler GPS veri değerlendirmeaşamasında doğru bir şekilde modellendiğitakdirde, belirlenen koordinatlar arasındaki farklarbüyük ölçüde noktaların zamana bağlı konumsaldeğişikliği (hız) olarak karşımıza çıkmaktadır.Eğer herhangi bir fiziksel etki değerlendirmeaşamasında model ile düzeltilmezse,koordinatların ve dolayısıyla da hızın içine hataolarak dâhil olmaktadır.Henüz modellenmemiş olan bazı periyodiketkiler, zaman serileri analizi ile tespitedilebilmektedir. Bu etkilerin kampanya tipiölçülerle oluşturulan hızlara düzeltme olarakgetirilerek daha gerçek ve güvenilir hızlaraulaşılması, bu çalışmanın amacınıoluşturmaktadır.Tektonik hız vektörü, yatay ve düşey olmaküzere iki düzlemde belirlenmektedir. Her ikibileşen için de iki epok arasındaki koordinatfarkları zamanla ölçeklendirilerek doğrusal hızelde edilmektedir. <strong>Genel</strong> bir kabul olarak GPShızları mm/yıl olarak ifade edilmektedir.VX(t) − X(t)2 1= (1)t2− t1V : Doğrusal hız (mm/yıl),t 1 , t 2 : Ölçü epokları (yıl),X(t 1 ), X(t 2 ) : Koordinat değerleri (mm)Yatay hız vektörü, sırasıyla kuzey-güney vedoğu-batı bileşenleri hesaplandıktan sonrabunların bileşkesi alınarak hesaplanmaktadır.Düşey hız ise, yönü yukarı veya aşağı olacakşekilde tek bir doğrultu üzerinde belirlenmektedir.3. GPS ZAMAN SERİLERİNİN GÜRÜLTÜÖZELLİKLERİYüzey deformasyonları veya gerinimbirikimlerinin hesaplanması için yapılan jeodezikölçülerin kullanıldığı jeofiziksel çalışmalar,sadece ilgili parametrelerin doğru tahminini değilaynı zamanda bu parametrelerin hatalarının dadoğru tahmin edilmesini gerektirmektedir (Maov.d., 1999).Birçok jeofiziksel olayda olduğu gibi GPShataları “Güç Yasası” (Power Law) süreciyletanımlanabilmektedir (Agnew, 1992) veyazamana bağlı olarak aşağıdaki güç spektrumunasahiptirler:P x (f)=P o ( f / f o ) κ (2)f : FrekansP o , f o : Normalizasyon sabitleriκ : Spektral indeksTabiat olayları ile meydana gelen süreçleryüksek frekanslarla karşılaştırıldığında, düşükfrekanslarda daha çok gürültü oranına ve -3
<strong>Harita</strong> <strong>Dergisi</strong> Temmuz 2010 Sayı 144Deformasyon Hızlarının İyileştirilmesinde Sabit GPS İstasyonlarıSabit GPS İstasyonları Zaman Serileri Analizinden Yaralanılmasıdüzensizlikleri, yer altı suları ve atomik saatlerinölçtüğü zamandaki hataları kapsayan genişkapsamlı dinamik süreçler içerisindegözlenmektedir (Mandelbrot, 1983). Zamanserilerindeki gürültülerin spektral indeks ekseniüzerindeki yerleri Şekil 1’de verilmiştir.Rasgele Yürüyüş Gür.(Klasik Brownian Har.)κ (Spk. İndis)-3 -2 -1 01- Durağan Olmayan Sür.- Kesirli Gaussian Har.- Kesirli Rasgele YürüyüşKırpışma Gür.Şekil 1. Zaman serilerindeki gürültülerin spektralindeks ekseni üzerindeki yerleri.Zhang v.d. (1997) Güney Kaliforniya ağından19 aylık sabit GPS istasyonu verisini inceledilerve zaman serilerinin, Beyaz ve KırpışmaGürültüsünün birleşimi (BG+KG) veya κ = -0.4indeks değerindeki kesirli bir gürültüye sahipolduğunu buldular. Mao v.d. (1999) 3 yıllık veriyesahip global olarak dağılmış GPS istasyonlarıverilerini en iyi şekilde Beyaz ve KırpışmaGürültüsünün birleşiminin tanımlayabileceğinigösterdiler. Günlük GPS çözümlerinin analizlerigibi, daha yüksek frekanslı GPS verileri de aynızamanda Beyaz+Kırpışma gürültüsüsergilemektedir (Nikolaidis, 2002).GPS verisinin zamana bağlı gürültüsününiçeriğini anlamak, bu verilerden hesaplananparametrelere gerçekçi yaklaşımlardabulunabilmek için oldukça önemlidir (Nikolaidis,2002). Gürültünün sadece Beyaz Gürültüdenoluştuğunu varsaymak hız belirsizliklerinin eksiktahmin edilmesine yol açacaktır. Örneğin, Zhangv.d. (1997) sadece Beyaz Gürültü modelininyerine Beyaz+Kırpışma Gürültüsü modelikullanıldığı zaman nokta hız belirsizliklerinin 3-6kat daha büyük olduğu sonucuna vardılar. Aynışekilde Mao v.d. (1999), korelasyonlu gürültüyü(Kırpışma Gürültüsü) dahil etmedikleri takdirdehız belirsizliklerinin 5-11 kat daha eksik tahminedileceğini kanıtladılar. Açıkça görüldüğü gibi gözönüne alınan gürültü modeli sonuç hız hatasınıoldukça etkilemektedir.Beyaz Gür.- Durağan Süreçler- Kesirli Brownian Har.- Kesirli BeyazGür.a. En Büyük Olabilirlik Kestirimi (EBOK)YöntemiZaman serilerindeki Beyaz ve Güç Yasası(Kırpışma, Rastgele Yürüyüş) gürültülerininmiktarlarını belirleyebilmek için Langbein veJohnson (1997), Zhang v.d., (1997) ve Mao v.d.(1999) tarafından açıklanan En Büyük OlabilirlikKestirimi (EBOK) yöntemi kullanılmıştır. EBOKyöntemi aynı zamanda veriyi en uyguntanımlayan gürültü modelinin belirlenmesini desağlamaktadır. EBOK yöntemini kullanarakgürültü bileşenlerini kestirmek için olasılıkfonksiyonu, veri kovaryans matrisinindengelenmesiyle maksimize edilmektedir(Williams v.d., 2004).( )( 2π) N/ 2( detC)T −1( −0.5vˆC vˆ )1lik vˆ,C = exp(3)1/ 2lik : Olabilirlik fonksiyonu,N : Epok sayısı,C : Veri kovaryans matrisi,vˆ : Dengeleme sonrası artıklar (orijinal zamanserilerine doğrusal veya doğrusal olmayan birmodel uygulandıktan sonra aynı C kovaryansmatrisi ile ağırlıklandırılmış en küçük kareleryöntemi kullanılarak elde edilen)Zaman serileri modeli genel olarak; serinin x(yatay eksen) eksenini kestiği yer, doğrusal bireğim (hız) ve yıllık ve altı aylık sinyallere aitsinüzoidal terimlerden oluşmaktadır. Model aynızamanda yapay veya ko-sismik atılımların sebepolabileceği kayıklıklar ve post-sismik atılımdankaynaklanabilecek üssel bir fonksiyona aitgenlikleri de ihtiva edebilir. “C” kovaryans matrisi;Beyaz Gürültü, Güç Yasası, Birinci DereceGauss Markov, Özbağlanımlı (autoregressive),Yürüyen Ortalamalı ve Band Geçirgen gibi Gaussolasılıksal (stochastic) modelinin birçok formunutemsil edebilmektedir (Williams v.d., 2004). Buçalışmada; C kovaryans matrisinin aşağıdakiformülde görüldüğü gibi, Beyaz Gürültü ve GüçYasası Gürültüsü olmak üzere iki gürültükaynağının birleşiminden oluştuğu kabuledilmiştir.C = a I + b J(4)2w2κκa w : Beyaz gürültü genliği,b κ : Güç Yasası Gürültüsü genliği,I : Birim matris (NxN boyutlu),Jκ : κ spektral indisine sahip Güç YasasıGürültüsünün kovaryans matrisi.23