<strong>Harita</strong> <strong>Dergisi</strong> Temmuz 2010 Sayı 144A.İ.KURT vd.Deformasyon Hızlarının İyileştirilmesinde Sabit GPS İstasyonları Zaman SerileriAnalizinden Yararlanılması(Improving Deformation Velocities via Utilizing Time Series Analysis ofContinuous GPS Stations)Ali İhsan KURT 1 , Rasim DENİZ 21 <strong>Harita</strong> <strong>Genel</strong> <strong>Komutanlığı</strong>, Ankara2 İstanbul Teknik Üniversitesi, İstanbulaliihsan.kurt@ hgk.msb.gov.trÖZETDeformasyon ölçmeleri çalışmalarında, kampanyatipi GPS gözlemleri kullanıcıya elindeki alıcı sayısındandaha fazla sayıda noktada gözlem yapma imkânıvererek ölçü bölgesinin daha yüksek konumsalçözünürlükle izlenmesine olanak sağlamaktadır. Diğertaraftan, kampanya tipi veriler zamansal olarakseyrektirler ve küçük sıçramaları tespit edemeyecekkadar yetersizdirler. Sabit GPS (SGPS) noktasındatoplanan sürekli GPS verileri ile yer kabuğuhareketlerine yönelik kinematik modellemeçalışmalarını gerçekleştirmek mümkün olabilmektedir.Bu çalışmada; sabit GPS istasyonları zaman serilerianalizi ile hesaplanan periyodik etkilerin Türkiye UlusalTemel GPS Ağı (TUTGA) hızlarına düzeltme olarakgetirilerek daha gerçekçi hıza ulaşılması için biryaklaşım önerilmiştir. Verileri herkese açık 6 sabit GPSistasyonundan elde edilen veriler GAMIT-GLOBKyazılımı ile değerlendirilerek ITRF2000 referanskoordinat sisteminde zaman serileri analizine hazırhale getirilmiştir. EBOK analizi ile seçilen istasyonlaraait en uygun modelin Beyaz Gürültü ve KırpışmaGürültüsünün birleşimi olduğu görülmüştür. Zamanserileri içindeki anlamlı sinyallerin frekanslarını vegenliklerini istasyon bazında hesaplamak amacıylaGPS zaman serilerine, boşluklu veriler için oldukçauygun bir yöntem olan Lomb-Scargle algoritmasıuygulanmıştır. Tüm istasyonlarda ortak olarak senelikve altı aylık etkilerin baskın olduğu tespit edilmiştir.Böylece, TUTGA hızlarına düzeltme olarak getirilecekperiyodik etkilerin hangilerinin modelleneceğine kararverilmiştir. Senelik ve altı aylık etkilerin, TUTGAhızlarına nasıl ilave edilebileceğine dair bir algoritmageliştirilerek, hız alanının iyileştirilmesi için bir yöntemönerilmiştir. Çalışmadan elde edilen sonuçlara ilişkinbir değerlendirme yapılarak, Türkiye genelinde vebölgesel olarak nasıl uygulanabileceğine ilişkin yatayve düşey hızlar için bir değerlendirme yapılmıştır.Anahtar Kelimeler: Sabit GPS, zaman serileri analizi,spektral analiz, Lomb-Scargle, TUTGA.ABSTRACTCampaign type GPS measurements provide higherresolution spatial monitoring of region by supplying tothe observer more station than of his receiver number.On the other hand, campaign type measurementsreturn baseline measurements with very coarsetemporal sampling. It is possible to make geodetic andgeodynamic calculations, carry out kinematic modelingof crustal movements via using continuous GPS(CGPS) measurements obtained from CGPS stations.In this study, an algorithm is investigated to receivemore accurate velocity field and improve velocityobtained from campaign type GPS measurements byadding corrections to the velocities via utilizing timeseries analysis of continuous GPS stations. Public dataof Ankara, İstanbul, Trabzon, Mersin, Gebze andAntalya stations is used as test measurements. GPSprocesses were done by GAMIT-GLOBK software andtime series are created in ITRF2000. MaximumLikelihood Estimation (MLE) was applied to investigatestochastic noise properties of GPS time series. Timeseries were tested by different noise models and acombination of white noise+flicker noise was found asthe best describing model. Lomb-Scargle algorithmwhich is suitable for unequally spaced time series wasapplied to investigate the amplitudes and frequenciesof dominant signals on a station by station basis.Seasonal signals are a potentially important part of theerror budget for continuous GPS sites. Annual andsemi-annual signals were determined dominantcommon to all sites. Thus, these periodical signalswere chosen to be brought as corrections to TurkishNational Fundamental GPS Network (TNFGN) velocityfield. An algorithm is suggested on how to applyannual and semi-annual effects in order to improveboth deformation measurements and TNFGNvelocities. An assessment is given to evaluate theresults obtained from this study on how to apply bothregional and covering the territory of Turkey forhorizontal and vertical velocity fields.Key Words: Continuous GPS, time series analysis,spectral anaysis, Lomb-Scargle, TNFGN.1. GİRİŞGPS ile yüzey hareketlerinin belirlenmesindehesaplanan ve kullanılan temel parametre,noktaların tektonik vb. nedenlerden kaynaklananyatay ve düşey hızıdır. GPS hızlarını etkileyen enönemli iki parametre ise periyodik etkilerin bağlıolduğu koordinat referans sistemi ve verilerinzamana bağlı hata özellikleridir. GPS noktakoordinatları fiziksel olarak, atmosfer, gel-git,kabuk hareketi, yeraltı su kaynaklarının hareketivb. birçok olaydan etkilenir. Bu etkiler eğer uygunbir modelleme ile düzeltme olarak getirilmez isehata olarak koordinatları etkiler.20
<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ıJeodezik GPS ölçmeleri zaman içindeki yerdeğiştirmeleri belirlemek için, kısa dönemlerlefakat periyodik olarak tekrarlanan kampanya tipigözlemler şeklinde yapılabileceği gibi; uzun birzaman aralığında yer hareketlerini sürekligözlemlemek için tasarlanan bir ağ şeklinde deyapılabilmektedir.Kampanya tipi gözlemler kullanıcıya elindekialıcı sayısından daha fazla sayıda noktadagözlem yapma imkânı vererek ölçü bölgesinindaha yüksek konumsal çözünürlükle izlenmesineolanak sağlamaktadır. Diğer taraftan, kampanyatipi veriler zamansal olarak seyrektirler vedoğrusal olmayan yer hareketlerini veya her yenikampanya gözleminde oluşabilecek küçüksıçramaları tespit edemeyecek kadaryetersizdirler. Bu noktada sabit GPS istasyonlarıuzun süreli olarak kesintisiz nokta konumlarınızamansal çözünürlüğü yüksek düzeydesağlayarak bu tür problemlerin üstesindengelebilmektedirler. Kısacası kampanya tipiölçülerle bölgesel deformasyon hakkındakonumsal olarak daha yüksek çözünürlüklü fakatdüşük frekansta bilgi toplarken, sabit GPSistasyonlarından elde edilen verilerle düşükçözünürlüklü fakat yüksek frekansta verisağlanmaktadır. Sabit GPS istasyonları verianalizleri ile kampanya tipi ölçülerin geliştirilmesiyaklaşımı, yüzey deformasyonları ve doğrusaltektonik hızın daha doğru ve gerçekçihesaplanmasını sağlayacaktır.Hızların bir diğer önemli kullanımı, jeodezikamaçlı ağ sıklaştırmalarıdır. Dinamik jeodezik ağyönteminde noktaların koordinatları ile birliktehızları da belirlenir ve kullanılır. Türkiye’de BüyükÖlçekli <strong>Harita</strong> ve <strong>Harita</strong> Bilgileri Yönetmeliğinegöre, TUTGA’ya dayalı nokta sıklaştırmalarında,C1 ve C2 derece noktaların hızları TUTGAnoktalarının hızlarına bağlı olarak hesaplanmaktave kullanılmaktadır. Farklı zamanlardagerçekleştirilecek nokta sıklaştırmaçalışmalarında, zaman birliği sağlamak içinkoordinatlar hızlarla bir to referans epokunaötelenirler. Bu hızlarla herhangi bir ölçmeepokundaki (t) hızların birbirleriyle uyumuönemlidir. GPS gibi yüksek doğruluklukonumların elde edildiği jeodezik ağlarda, hertürlü dönüşümde kullanılacak gerçek ve güvenilirhızların hesaplanması bu nedenle önemkazanmıştır. TUTGA ve buna bağlı noktahızlarının kampanya çözümlerindenhesaplanması durumunda, yukarıda bahsedilenolgular burada da söz konusu olmaktadır.Bu çalışmada GPS veri değerlendirmeleriGAMIT-GLOBK yazılımı ile yapılmıştır. Önceliklegünlük çözümler oluşturulmuş ve daha sonra yineaynı yazılımla, global IGS istasyonlarındanyararlanarak ITRF2000’de zaman serileriüretilmiştir. Bu aşamada yapılan diğer birçalışma, zaman serileri içerisinden uyuşumsuz vekaba hatalı ölçülerin ayıklanması olmuştur.MATLAB ortamında hazırlanan yazılımla, EnKüçük Kareler Dengelemesi kullanılarakdengeleme sonrası artıklara Pope Kaba Hatatesti uygulanmıştır. Uyuşumsuz ölçü olarak tespitedilen veriler zaman serilerinden çıkarılarakanalize hazır hale getirilmiştir.GPS zaman serilerinin stokastik gürültüözelliklerini araştırmak amacıyla; En BüyükOlabilirlik Analizi (EBOK) uygulanmıştır. EBOKanalizinin başlıca amaçlarından birisi, tektonikhızın hatasının kestiriminde gerçekçi yaklaşımdabulunmaktır. Çünkü yalnızca beyaz gürültününkabul edildiği ve diğer faktörlerin ihmal edildiği birmodelde, tektonik hızın belirsizliğinin daha küçükhesaplanarak yanlış bir kestirimde bulunmaolasılığı bulunmaktadır.Diğer yandan serilerin, zamana bağlı olmayandurağan süreçler (beyaz gürültü) ve zamanabağlı durağan olmayan (renkli) süreçlerden hangiölçüde etkilendikleri ve en uygun gürültü modelibelirlenmiştir.Boşluklu verilerin spektral analizi için buamaca uygunluğu Langbein ve Johnson (1997),Zhang v.d. (1997) ve Mao v.d. (1999)’de de ifadeedilen Lomb-Scargle algoritması kullanılmıştır.Zaman serileri içindeki anlamlı sinyallerinfrekanslarını ve genliklerini istasyon bazındahesaplamak amacıyla GPS zaman serilerineLomb-Scargle algoritması uygulanmıştır.Böylece, kampanya tipi GPS ölçü hızlarınadüzeltme olarak getirilecek periyodik etkilerinhangilerinin modelleneceğine karar verilmiştir.Senelik ve altı aylık etkilerin, TUTGA hızlarınanasıl ilave edilebileceğine dair bir yöntemönerilmiştir. Bunun için öncelikle tüm sabit GPSistasyonlarının her üç bileşeni için hesaplanansenelik ve altı aylık etkiler, enterpolasyonla gridnoktaları için belirlenmiştir. Çalışmadan eldeedilen sonuçlara ilişkin bir değerlendirmeyapılarak, Türkiye genelinde ve bölgesel olaraknasıl uygulanabileceğine ilişkin yatay ve düşeyhızlar için bir değerlendirme örneği geliştirilmiştirve önerilmektedir.21