This study gives a processing strategy with variance component estimation aiming to performa proper stochastic model in a gravity monitoring network.First, a program which is able to make some necessary reductions and to adjust gravityobservations was written in Matlab. The basic properties of the program called GrAd asfollows: 1) Gravimetic tide reductions are made according to some tidal potential models, 2)robust estimation with Denmark weight function and median methods are used for outlierdetection, 3) the observations belonging to the lines of a network are combined properly in theadjustment, 4) gravity values of the points are estimated with a trace-minimum adjustment fordeformation analysis, 5) deformation analysis is performed for different epochs with theglobal test and localization procedures and the temporal gravity changes are obtainedaccording to the statistically-verified reference points in the network. In order to estimate thevariance components, the method of BIQUE (Best Invariant Quadratic Unbiased Estimation)was used in the study. An algorithm in which the components of different meters areestimated iteratively (iterative-BIQUE) was adapted to GrAd.For numerical applications, Marmara Region Gravity Network (MAGRANET) with 25gravity points was used. The network was measured with 2 LaCoste&Romberg relativegravity meters between 2003-2005 semi-annually as 5 periods. For free adjustments ofdifferent periods, the weigths of the meters? readings, a) were selected equally (=P) as usual,b) were obtained according to the estimated variance components. As a result of theestimation of the components, it is shown that the meters? weigths are different(approximately P and P/2, respectively). Comparing the estimated gravity values for the casesof a and b, the differences between ?7 and 11 ?Gal were observed. Moreover, using theestimated gravity values and their cofactor matrices from the two cases separately, the periodsof October 2004-October 2003, October 2005-October 2004 and June 2005-June 2004 wereanalysed with localization procedure. According to the deformation analyses? results, thetemporal gravity changes computed for the two cases vary each other significantly.Especially, for the June 2005-June 2004 period, the variations reach up to 17 ?Gal. As a resultof the statistically significant differences, it is concluded that the improper stochastic part of aGauss-Markoff model for gravity networks affects adversely the interpretation of the temporalgravity changes.With the temporal gravity changes obtained for the case b, it is seen clearly that the regionbetween Istanbul and Izmit as well as the southern part of the Izmit Bay have differentanomaly patterns related to the changes on the other points of the MAGRANET. Themagnitudes of temporal gravity changes in the corresponding parts of the network reach up to50 ?Gal per year.Key Words: Temporal gravity changes?Variance component estimation?BIQUE method?Gravity networks?Deformation analysis?Marmara Region.

Bu çalışmada, zamansal gravite değişimlerinin izlendiği gravite ağları için, stokastik modelinvaryans bileşen kestirimiyle belirlendiği bir değerlendirme stratejisi ortaya konmaktadır.Öncelikle Matlab'de, ölçülerin değerlendirme ve dengelemesini yapabilen bir programoluşturulmuştur. GrAd olarak adlandırılan bu programın temel özellikleri şunlardır: 1) Gelgitindirgemesi çeşitli gelgit potansiyel modellerine göre yapılmaktadır. 2) Uyuşumsuz ölçüaraştırması için medyan ölçütü ve Danimarka ağırlık fonksiyonunun kullanıldığı robustyöntem uygulanmaktadır. 3) Gravite ağını oluşturan geçkiler uygun biçimdebirleştirilmektedir. 4) Deformasyon analizi için gerekli gravite bilgisi tüm iz minimumyöntemine göre dengeleme ile üretilmektedir. 5) Global test ve yerelleştirme işlemlerine göredeformasyon analizi yapılmakta, istatistiksel olarak tanımlanan referans noktalarına göregravite değişimleri hesaplanmaktadır. Varyans bileşen kestirimi için BIQUE (Best InvariantQuadratic Unbiased Estimation) yöntemi seçilmiştir. Bileşenlerin iteratif olarak kestirildiği biralgoritma (iteratif-BIQUE) GrAd programına eklenerek farklı gravite ölçerler için varyansbileşen kestirimi yapılmaktadır. Elde edilen varyans bileşenleri ile yeni ağırlık matrisioluşturulmakta ve gravite ağı dengeleme işlemi yinelenmektedir.Sayısal uygulama için 25 noktadan oluşan Marmara Bölgesi Gravite Ağı (MAGRANET)kullanılmıştır. 2 adet LaCoste&Romberg bağıl gravite ölçerle 2003-2005 yılları arasında 5periyot ölçülen bu ağ, GrAd programı ile serbest dengelenmiştir. Periyot dengelemelerindegravite ölçer okumalarının ağırlıkları, a) alışıldık biçimde eşit (P) öngörülmüş ve b) varyansbileşen kestirimine göre oluşturulmuştur. Kestirim sonrası gravite ölçer ağırlıklarının farklıolması gerektiği (yaklaşık P ve P/2) ortaya çıkmıştır. Her iki durumda serbest dengelemelerleelde edilen gravite değerleri karşılaştırıldığında, ?7 ve 11 ?Gal arasında değişen farklargözlenmiştir. Bununla birlikte, her iki durum için elde edilen gravite değerleri ve bunlarınağırlık katsayıları matrisleri kullanılarak deformasyon analizleri yapılmıştır. Ekim 2004-Ekim2003, Ekim 2005-Ekim 2003 ve Haziran 2005-Haziran 2004 dönemlerine ilişkin analizlerinfarklı sonuçlar verdiği gözlenmiştir. Haziran 2005-Haziran 2004 dönemi için iki durumakarşılık yerelleştirmeler sonrası elde edilen gravite değişimleri arasında 17 ?Gal'e varanfarklar çıkmıştır. Böylece, gerçeği yansıtmayan bir stokastik modelin deformasyon analizisonuçlarını önemli ölçüde etkilediği gösterilmiştir.(b) durumu için elde edilen gravite değişimlerine göre, Marmara Bölgesinde özellikle Gebze-Tuzla arasında ve İzmit Körfezinin güneyinde yer alan Gemlik bölgesinde farklı anomalilerinbulunduğu saptanmıştır. Bu bölgedeki noktalarda zamansal gravite değişimleri 50 ?Gal'eulaşmaktadır.Anahtar Kelimeler: Zamansal gravite değişimi?Varyans bileşen kestirimi?BIQUE yöntemi?Gravite ağları?Deformasyon analizi?Marmara Bölgesi.