New Modifications of Stokes’ Integral

2017 ◽  
pp. 45-52
Author(s):  
Lars E. Sjöberg ◽  
Mehdi S. Shafiei Joud
Keyword(s):  
Stokes Integral

Comparison of spectral methods for the evaluation of Stokes integral  

2021 ◽  
Author(s):  
Avadh Bihari Narayan ◽  
Ashutosh Tiwari ◽  
Govind Sharma ◽  
Balaji Devaraju ◽  
Onkar Dikshit
Keyword(s):  
Boundary Value Problems ◽  
Numerical Integration ◽  
Spectral Methods ◽  
Boundary Value ◽  
Fundamental Equation ◽  
Computation Time ◽  
Spherical Approximation ◽  
Integration Technique ◽  
Gravity Measurements ◽  
Stokes Integral

<p>The spherical approximation of the fundamental equation of geodesy defines the boundary value problems. Stokes’s integral provides the solution of boundary value problems that enables the computation of geoid from the properly reduced gravity measurements to the geoid. The stokes integral can be evaluated by brute-force numerical integration, spectral methods, and least-squares collocation. There is a trade-off between computation time and accuracy when we chose numerical integration technique or any spectral method. This research will compare time complexity and the accuracy of different spectral methods (1D-FFT, 2D-FFT, Multi-band FFT) and numerical integration technique for the region in the lower Himalaya, around Nainital, Uttarakhand, India. </p>

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

2011 ◽  
Vol 37 (1) ◽  
pp. 5-14 ◽  
Author(s):  
Ahmed Abdalla ◽  
Robert Tenzer
Keyword(s):  
New Zealand ◽  
Least Squares ◽  
Gravity Data ◽  
Spherical Approximation ◽  
Geopotential Model ◽  
Gravimetric Geoid ◽  
Geoid Model ◽  
Marine Gravity ◽  
Gravity Database ◽  
Stokes Integral

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.

Comparison among Gravimetric, Astrogravimetric and Astrogeodetic Geoids: Case Study for Austria

2021 ◽  
Author(s):  
Hussein Abd-Elmotaal ◽  
Norbert Kühtreiber
Keyword(s):  
Gravity Anomalies ◽  
Geoid Determination ◽  
Data Sets ◽  
Deflections Of The Vertical ◽  
Collocation Technique ◽  
Dense Grid ◽  
Vening Meinesz ◽  
Stokes Integral ◽  
Stokes Kernel

<p>It is used to state that all geoid determination techniques should yield to the same geoid if the indirect effect is properly taken into account (Heiskanen and Moritz, 1967). The current study compares different geoid determination techniques for Austria. The used techniques are the gravimetric, astrogravimetric and astrogeodetic geoid determination techniques. The available data sets (gravity, deflections of the vertical, height, GPS) are described. The window remove-restore technique (Abd-Elmotaal and Kuehtreiber, 2003) has been used. The available gravity anomalies and the deflections of the vertical have been topographically-isostatically reduced using the Airy isostatic hypothesis. The reduced deflections have been used to interpolate deflections on a relatively dense grid covering the data window. These gridded reduced deflections have been used to compute an astrogeodetic geoid for Austria using least-squares collocation technique within the remove-restore scheme. The Vening Meinesz formula has been used to compute an astrogravimetric geoid for Austria. Another gravimetric geoid for Austria has been determined in the framework of the window remove-restore technique using Stokes integral with modified Stokes kernel. All computed geoids have been validated using GNSS/levelling derived geoid. A wide comparison among the derived geoids computed within the current investigation has been carried out.</p>

scholarly journals Accuracy of unmodified Stokes’ integration in the R-C-R procedure for geoid computation

2015 ◽  
Vol 9 (2) ◽  
Author(s):  
Zahra Ismail ◽  
Olivier Jamet
Keyword(s):  
Gravity Anomalies ◽  
Real Data ◽  
Mountain Areas ◽  
General Behaviour ◽  
Spherical Harmonic Degree ◽  
Mountainous Regions ◽  
Optimal Accuracy ◽  
Stokes Integral ◽  
Stokes Kernel ◽  
Geoid Computation

AbstractGeoid determinations by the Remove-Compute-­Restore (R-C-R) technique involves the application of Stokes’ integral on reduced gravity anomalies. Numerical Stokes’ integration produces an error depending on the choice of the integration radius, grid resolution and Stokes’ kernel function.In this work, we aim to evaluate the accuracy of Stokes’ integral through a study on synthetic gravitational signals derived from EGM2008 on three different landscape areas with respect to the size of the integration domain and the resolution of the anomaly grid. The influence of the integration radius was studied earlier by several authors. Using real data, they found that the choice of relatively small radii (less than 1°) enables to reach an optimal accuracy. We observe a general behaviour coherent with these earlier studies. On the other hand, we notice that increasing the integration radius up to 2° or 2.5° might bring significantly better results. We note that, unlike the smallest radius corresponding to a local minimum of the error curve, the optimal radius in the range 0° to 6° depends on the terrain characteristics. We also find that the high frequencies, from degree 600, improve continuously with the integration radius in both semi-­mountainous and mountain areas.Finally, we note that the relative error of the computed geoid heights depends weakly on the anomaly spherical harmonic degree in the range from degree 200 to 2000. It remains greater than 10 % for any integration radii up to 6°. This result tends to prove that a one centimetre accuracy cannot be reached in semi-mountainous and mountainous regions with the unmodified Stokes’ kernel.

scholarly journals Generalized geoidal estimators for deterministic modifications of spherical Stokes' function

2010 ◽  
Vol 40 (1) ◽  
pp. 45-64 ◽  
Author(s):  
Michal Šprlák
Keyword(s):  
Gravity Model ◽  
Spherical Harmonic ◽  
Gravity Data ◽  
Surface Integral ◽  
Terrestrial Gravity ◽  
Global Gravity ◽  
Higher Derivatives ◽  
Stokes Integral ◽  
Global Gravity Model ◽  
Local Geoid

Generalized geoidal estimators for deterministic modifications of spherical Stokes' function Stokes' integral, representing a surface integral from the product of terrestrial gravity data and spherical Stokes' function, is the theoretical basis for the modelling of the local geoid. For the practical determination of the local geoid, due to restricted knowledge and availability of terrestrial gravity data, this has to be combined with the global gravity model. In addition, the maximum degree and order of spherical harmonic coefficients in the global gravity model is finite. Therefore, modifications of spherical Stokes' function are used to obtain faster convergence of the spherical harmonic expansion. Decomposition of Stokes' integral and modifications of Stokes' function have been studied by many geodesists. In this paper, the proposed deterministic modifications of spherical Stokes' function are generalized. Moreover, generalized geoidal estimators, when the Stokes' integral is decomposed in to spectral and frequency domains, are introduced. Higher derivatives of spherical Stokes' function and their numerical stability are discussed. Filtering and convergence properties for deterministic modifications of the spherical Stokes' function in the form of a remainder of the Taylor polynomial are studied as well.

scholarly journals Applications of Voronoi and Delaunay Diagrams in the solution of the geodetic boundary value problem

2012 ◽  
Vol 18 (3) ◽  
pp. 378-390
Author(s):  
C. A. B. Quintero ◽  
I. P. Escobar ◽  
C. F. Ponte-Neto
Keyword(s):  
Boundary Value Problem ◽  
Gravity Anomaly ◽  
Indirect Effect ◽  
Boundary Value ◽  
Original Data ◽  
Target Area ◽  
Surface Integration ◽  
A Cell ◽  
Stokes Integral ◽  
Original Information

Voronoi and Delaunay structures are presented as discretization tools to be used in numerical surface integration aiming the computation of geodetic problems solutions, when under the integral there is a non-analytical function (e. g., gravity anomaly and height). In the Voronoi approach, the target area is partitioned into polygons which contain the observed point and no interpolation is necessary, only the original data is used. In the Delaunay approach, the observed points are vertices of triangular cells and the value for a cell is interpolated for its barycenter. If the amount and distribution of the observed points are adequate, gridding operation is not required and the numerical surface integration is carried out by point-wise. Even when the amount and distribution of the observed points are not enough, the structures of Voronoi and Delaunay can combine grid with observed points in order to preserve the integrity of the original information. Both schemes are applied to the computation of the Stokes' integral, the terrain correction, the indirect effect and the gradient of the gravity anomaly, in the State of Rio de Janeiro, Brazil area.

scholarly journals The Newtonian Potential and the Demagnetizing Factors of the General Ellipsoid

2017 ◽  
Author(s):  
Giovanni Di Fratta
Keyword(s):  
Integral Representation ◽  
Explicit Expression ◽  
Representation Formula ◽  
Newtonian Potential ◽  
Bounded Domains ◽  
Dimensional Problem ◽  
Smooth Functions ◽  
Integral Representation Formula ◽  
Transport Theorem ◽  
Stokes Integral

The objective of this paper is to present a modern and concise new derivation for the explicit expression of the interior and exterior Newtonian potential generated by homogeneous ellipsoidal domains in $\mathbb{R}^N$ (with $N \geqslant 3$). The very short argument is essentially based on the application of Reynolds transport theorem in connection with Green-Stokes integral representation formula for smooth functions on bounded domains of$\mathbb{R}^N$, which permits to reduce the N-dimensional problem to a 1-dimensional one. Due to its physical relevance, a separate section is devoted to the derivation of the demagnetizing factors of the general ellipsoid which are one of the most fundamental quantities in ferromagnetism.

Sign in / Sign up

Export Citation Format

Share Document