THE EVALUATION OF THE NEW ZEALAND'S GEOID MODEL USING THE KTH METHOD / NAUJOSIOS ZELANDIJOS GEOIDO MODELIO VERTINIMAS KTH METODU / ОЦЕНКА МОДЕЛИ ГЕОИДА НОВОЙ ЗЕЛАНДИИ С ПРИМЕНЕНИЕМ МЕТОДА КТН

Abstract

We compile a new geoid model at the computation area of New Zealand and its continental shelf using the method developed at the Royal Institute of Technology (KTH) in Stockholm. This method utilizes the least-squares modification of the Stokes integral for the biased, unbiased, and optimum stochastic solutions. The modified Bruns-Stokes integral combines the regional terrestrial gravity data with a global geopotential model (GGM). Four additive corrections are calculated and applied to the approximate geoid heights in order to obtain the gravimetric geoid. These four additive corrections account for the combined direct and indirect effects of topography and atmosphere, the contribution of the downward continuation reduction, and the formulation of the Stokes problem in the spherical approximation. The gravimetric geoid model is computed using two heterogonous gravity data sets: the altimetry-derived gravity anomalies from the DNSC08 marine gravity database (offshore) and the ground gravity measurements from the GNS Science gravity database (onshore). The GGM coefficients are taken from EIGEN-GRACE02S complete to degree 65 of spherical harmonics. The topographic heights are generated from the 1×1 arc-sec detailed digital terrain model (DTM) of New Zealand and from the 30×30 arc-sec global elevation data of SRTM30_PLUS V5.0. The least-squares analysis is applied to combine the gravity and GPS-levelling data using a 7-parameter model. The fit of the KTH geoid model with GPS-levelling data in New Zealand is 7 cm in terms of the standard deviation (STD) of differences. This STD fit is the same as the STD fit of the NZGeoid2009, which is the currently adopted official quasigeoid model for New Zealand. Santrauka Stokholmo Karališkajame technologijos institute (KTH) sukurtu metodu apskaičiuotas naujas Naujosios Zelandijos ir kontinentinio šelfo geoido modelis. Taikoma Stokso integralo mažiausiųjų kvadratų modifikacija, įvertinant paklaidas ir jų nevertinant bei ieškant optimalių stochastinių sprendinių. Modifikuotas Bruno ir Stokso integralas sieja regioninius žemyninius gravimetrinius duomenis su globaliuoju geopotencialo modeliu (GGM). Gravimetriniam geoidui gauti skaičiuojamos keturios papildomos pataisos: topografinės situacijos ir atmosferos tiesioginės ir netiesioginės įtakos, redukcijos įtakos ir Stokso integralo taikymo sferiniam paviršiui. Gravimetrinis geoido modelis apskaičiuotas pagal du duomenų rinkinius: DNSC08 jūrinių gravimetrinių duomenų bazėje (šelfas) esančias altimetriniu metodu nustatytas sunkio pagreičio anomalijas ir žemyninės dalies gravimetrinių matavimų duomenis iš GNS gravimetrinės duomenų bazės (pakrantė). GGM koeficientai imti iš EIGEN-GRACE02S modelio sferinių iki 65 laipsnio harmonikų. Topografiniai aukščiai sugeneruoti iš Naujosios Zelandijos 1×1 sekundės detaliojo skaitmeninio reljefo modelio ir iš 30×30 sekundžių globaliojo aukščių modelio SRTM30_PLUS V5.0. Gravimetriniams ir GPS niveliacijos duomenims sujungti taikytas mažiausiųjų kvadratų 7 parametrų metodas. KTH metodu sudaryto geoido modelio vidutinė kvadratinė paklaida 7 cm. Tai sutampa su NZGeoid 2009 geoido modelio, taikomo Naujoje Zelandijoje, tikslumu. Резюме Модель геоида континентального шельфа Новой Зеландии построена с применением метода, созданного в Королевском технологическом институте Стокгольма. Данный метод основан на модификации решения интеграла Стокса методом наименьших квадратов с оценкой или без оценки погрешностей и поиском оптимальных статистических решений. Модифицированный интеграл БрунаСтокса объединяет региональные надземные гравиметрические данные с глобальной геопотенциальной моделью (GGM). Для определения гравиметрического геоида вычисляются дополнительные поправки прямого и косвенного влияния топографии и атмосферы, редукции и применения проблемы Стокса для сферической поверхности. Гравиметрическая модель геоида вычисляется на основе двух баз данных: альтиметрическим методом определенных аномалий силы тяжести в базе морских гравиметрических данных DNSC08 (шельф) и надземной части гравиметрических измерений из базы данных GNS. Коэффициенты GGM взяты из сферических гармоник до 65 степени модели EIGENGRACEO2S. Топографические высоты сгенерированы из детальной цифровой модели рельефа Новой Зеландии с сеткой 1×1 секунду и из глобальной модели высот SRTM30_PLUSv5.0 с сеткой 30×30 секунд. Для объединения гравиметрических и GPSнивелирных данных применялся метод наименьших квадратов с 7 параметрами. Среднеквадратическая погрешность модели геоида, созданной по методу КТН, равна 7 см. Точность аналогична точности применяемой в Новой Зеландии модели геоида NZGeoid2009.