Geoid Model Determination Research Articles

Assessments of Gravity Data Gridding Using Various Interpolation Approaches for High-Resolution Geoid Computations

This article investigates the role of different approaches and interpolation methods in gridding terrestrial gravity anomalies. In this regard, first of all, simple and complete Bouguer anomalies are considered in gravity data gridding. In the comparison results of gridding these two Bouguer anomaly datasets, the effect of the high-frequency contribution of topographic gravitation (by means of the terrain correction) is clarified. After that, the role of the used interpolation algorithm on the resulting grid of mean gravity anomalies and hence on the geoid modeling accuracy is inspected. For this purpose, four different interpolation methods including geostatistical Kriging, nearest neighbor, inverse distance to a power (IDP), and artificial neural networks (ANNs) are applied. Here, the IDP and nearest neighbor methods represent simple-structured algorithms among the interpolation methods tested in this study. The ANN method, on the other hand, is preferred as a complex, optimization-based soft computing method that has been applied in recent years. In addition, the geostatistical Kriging method is one of the conventional methods that is mostly applied for gridding gravity data in geodesy and geophysics. The calculated gravity anomalies in grids are employed in high-resolution geoid model computations using the least squares modifications of Stokes formula with additive corrections (LSMSA) technique. The investigations are carried out using the test datasets of Auvergne, France that are provided by the International Service for the Geoid for scientific research. It is concluded that the interpolation algorithms affect the gravity gridding results and hence the geoid model determination. The ANN method does not provide superior results compared to the conventional algorithms in gravity gridding. The geoid model with 4.1 cm accuracy is computed in the test area.

Geoid model determination for the Hellenic area “Hellas Geoid 2023”

The latest geoid model "HELLAS GEOID 2023" (HG2023) derived by the Hellenic Military Geographical Service is the most comprehensive model for the entire Hellenic area. Long-term gravity data, orthometric and geometric heights, seabed topography and a high-resolution digital terrain model were implemented in the calculations. Data evaluation and their accuracy estimation were of major importance to ascertain compatibility among data sources. In this study, data from neighboring countries were used as they were essential for the completeness of the model developed. The technique remove–compute–restore was adopted for the separation of the gravity signal in the heterogeneous data utilizing EIGEN 6C4 (full degree and order 2190), since the global geodynamic model was found to fit best in the Greek region. Mean gravity residual contribution to the local geoid model was calculated using Stokes' theorem in the frequency domain and implementing the Fourier transform using the Wang and Core modification. The final surface of the resulting gravimetric geoid was adapted to the existing height system of the State utilizing the Least Squares Collocation method by fitting points of known orthometric and geometric heights, distributed throughout Greece. The external accuracy of the estimated geoid model was estimated at the level of 6 cm.

Gravimetric geoid model determination of Central Ukraine area using combination of LSC and truncated SHA methods

Gravimetric geoid model determination of Central Ukraine area using combination of LSC and truncated SHA methods

Utilizing Airborne LiDAR and UAV Photogrammetry Techniques in Local Geoid Model Determination and Validation

This investigation evaluates the performance of digital terrain models (DTMs) generated in different vertical datums by aerial LiDAR and unmanned aerial vehicle (UAV) photogrammetry techniques, for the determination and validation of local geoid models. Many engineering projects require the point heights referring to a physical surface, i.e., geoid, rather than an ellipsoid. When a high-accuracy local geoid model is available in the study area, the physical heights are practically obtained with the transformation of global navigation satellite system (GNSS) ellipsoidal heights of the points. Besides the commonly used geodetic methods, this study introduces a novel approach for the determination and validation of the local geoid surface models using photogrammetry. The numeric tests were carried out in the Bergama region, in the west of Turkey. Using direct georeferenced airborne LiDAR and indirect georeferenced UAV photogrammetry-derived point clouds, DTMs were generated in ellipsoidal and geoidal vertical datums, respectively. After this, the local geoid models were calculated as differences between the generated DTMs. Generated local geoid models in the grid and pointwise formats were tested and compared with the regional gravimetric geoid model (TG03) and a high-resolution global geoid model (EIGEN6C4), respectively. In conclusion, the applied approach provided sufficient performance for modeling and validating the geoid heights with centimeter-level accuracy.

Creation of the Gravimetric Network and Determination of the GEOÏD in the Mediterranean Basin (Region-Tunisia)

The determination of geoid models with great precision (centimeter) was always at the center of interest in geodesy research [1] [2]. One of the most used methods to calculate the geoid model is the method called Remove-Compute-Restore (R-C-R). This method applies the stokes’ integral formula by the use of short, medium and long wavelength information via the three main stages R-C-R. The GravSoft software implemented by [3] was used for this study. Geodetic heights, a digital terrain model (SRTM) and leveled GPS points were used as inputs. The geoid modeling was carried out on the North region of Tunisia (Grand Tunis and Bizerte) on an area of 83 × 83 km. The accuracy of the quasi-geoid provisional reached 3.1 cm.

Modelling Local Gravity Anomalies from Processed Observed Gravity Measurements for Geodetic Applications

As the application of gravity data in applied sciences such as geodesy, geodynamics, astronomy, physics and geophysics for earth shape determination, geoid model determination, computation of terrestrial mass displacement, orbit computation of natural and artificial celestial bodies, realization of force standards and derived quantities and density distribution in the different layers in the upper crust and having considered the cost of direct gravity survey, the study presents modelling local gravity anomalies from processed observed gravity measurements for geodetic application in Benin City. A total of 22 points were used. The points were respectively observed with CHC900 dual frequency GNSS receivers and SCINTREX CG-5 Autograv to obtain their coordinates and absolute gravity values. The theoretical gravity values of the points were computed on the Clarke 1880 ellipsoid to obtain their local gravity anomalies. The free air and the Bouguer corrections were applied to the computed gravity anomalies to obtain the free air and the Bouguer gravity anomalies of the points. Least squares adjustment technique was applied to obtain the model variables coefficient/parameters, as well as to fit the fifth-degree polynomial interpolation surface to the computed free air and the Bouguer gravity anomalies. Kriging method was applied using Surfer 12 software to plot the computed and the models' free air and Bouguer gravity anomalies. Microsoft Excel programs were developed for the application of the models in the study area. The Root Mean Square Errors (RMSEs) and the standard errors of the two models were computed to obtain the dependability, as well as reliability of the models. It is recommended that whenever either free air or Bouguer gravity anomalies of points within Benin City are to be obtained for application in applied sciences, the determined models should be applied.

Relationship between Polynomial Geometric Surfaces Terms and Observation Points Numbers and Effect in the Accuracy of Geometric Geoid Models

The application of the geometric method of local geoid model determination which requires the fitting of geometric surfaces to known geoid heights to enable geoid heights of new points to be interpolated involves the use of least squares technique for computation of the models' parameters. The selection of polynomial geometric surfaces depends on the size of the study area, the variation of the geoid heights and the number of measurement points. The accuracy of the geometric geoid model increases as the number of observation points approximates the number of geometric surface terms. But in most cases, the number of observation points is not considered. To this effect, this paper presents the relationship between geometric surfaces terms and observation points numbers and effect in the accuracy of geometric geoid models. A total of 23 points of known local gravimetric geoid heights were used. Two polynomial geometric (third and fifth degrees) surfaces were fitted to the geoid heights at various observation point numbers and compared to determine the relationship between the number of model terms and that of observation points and effect in the accuracy of the models. Least squares adjustment technique was applied to obtain the model parameters. The differences between the models and the known geoid heights of the points were computed and used to obtain the RMSEs as well as the accuracy of the models. The obtained results showed that the accuracy of the polynomial geometric geoid models tends to the highest as the number of measurement points approximates the number of the model terms and in a unique solution where the number of observation points is equal to the number of the polynomial geometric model terms, the model accuracy is highest. The paper recommends that the geometric method of local geoid model determination should be strictly applied in small areas. Where the method will be applied in considerable large areas, higher degrees polynomial geometric surfaces with a larger number of terms approximating the number of observation points should be applied. This will enable a proper fit of the polynomial surface to the known geoid heights, as well as high accuracy to be obtained.

Derivation of Theoretical Gravity Model on the Clarke 1880 Ellipsoid for Practical Local Geoid Model Determination

Derivation of Theoretical Gravity Model on the Clarke 1880 Ellipsoid for Practical Local Geoid Model Determination

Practical Local Geoid Model Determination for Mean Sea Level Heights of Surveys and Stable Building Projects

A local geoid model for Evboriaria, Benin City using the geometric (GPS/Levelling) method was determined for calculation of mean sea level heights. Fifty points were established for the model and ten points were used for interpolation. The geoid heights were determined by finding the difference between the observed orthometric heights and the ellipsoidal heights. The polynomial regression model D was used for the interpolation of the orthometric heights. The computed mean standard deviation between the observed orthometric heights and the interpolated orthometric heights was ± 21cm. A mean geoidal undulation of 28.410m was computed using the gravimetric method. The computed orthometric heights using the gravimetry mean geoidal undulation were compared with the observed orthometric heights and seen to be identical. It is recommended that orthometric heights of project areas should be determined from GPS observations with the local geoid model of the area also determined.

Software development and its description for Geoid determination based on Spherical-Cap-Harmonics Modelling using digital-zenith camera and gravimetric measurements hybrid data

Over several years the Institute of Geodesy and Geoinformatics (GGI) was engaged in the design and development of a digital zenith camera. At the moment the camera developments are finished and tests by field measurements are done. In order to check these data and to use them for geoid model determination DFHRS (Digital Finite element Height reference surface (HRS)) v4.3. software is used. It is based on parametric modelling of the HRS as a continous polynomial surface. The HRS, providing the local Geoid height N, is a necessary geodetic infrastructure for a GNSS-based determination of physcial heights H from ellipsoidal GNSS heights h, by H=h-N. The research and this publication is dealing with the inclusion of the data of observed vertical deflections from digital zenith camera into the mathematical model of the DFHRS approach and software v4.3. A first target was to test out and validate the mathematical model and software, using additionally real data of the above mentioned zenith camera observations of deflections of the vertical. A second concern of the research was to analyze the results and the improvement of the Latvian quasi-geoid computation compared to the previous version HRS computed without zenith camera based deflections of the vertical. The further development of the mathematical model and software concerns the use of spherical-cap-harmonics as the designed carrier function for the DFHRS v.5. It enables - in the sense of the strict integrated geodesy approach, holding also for geodetic network adjustment - both a full gravity field and a geoid and quasi-geoid determination. In addition, it allows the inclusion of gravimetric measurements, together with deflections of the vertical from digital-zenith cameras, and all other types of observations. The theoretical description of the updated version of DFHRS software and methods are discussed in this publication.

Testing Stokes-Helmert geoid model computation on a synthetic gravity field: experiences and shortcomings

We report on testing the UNB (University of New Brunswick) software suite for accurate regional geoid model determination by use of Stokes-Helmert’s method against an Australian Synthetic Field (ASF) as “ground truth”. This testing has taken several years and has led to discoveries of several significant errors (larger than 5mm in the resulting geoid models) both in the UNB software as well as the ASF. It was our hope that, after correcting the errors in UNB software, we would be able to come up with some definite numbers as far as the achievable accuracy for a geoid model computed by the UNB software. Unfortunately, it turned out that the ASF contained errors, some of as yet unknown origin, that will have to be removed before that ultimate goal can be reached. Regardless, the testing has taught us some valuable lessons, which we describe in this paper. As matters stand now, it seems that given errorless gravity data on 1′ by 1′ grid, a digital elevation model of a reasonable accuracy and no topographical density variations, the Stokes-Helmert approach as realised in the UNB software suite is capable of delivering an accuracy of the geoid model of no constant bias, standard deviation of about 25 mm and a maximum range of about 200 mm. We note that the UNB software suite does not use any corrective measures, such as biases and tilts or surface fitting, so the resulting errors reflect only the errors in modelling the geoid.

The use of minimum curvature surface technique in geoid computation processing of Egypt

According to the wide spread use of satellite-based positioning techniques, especially Global Navigation Satellite Systems (GNSS), a greater attention has been paid to the precise determination of geoid models. As it is known, leveling measurements require high cost and long time in observation process that make it not convenient for the practical geodetic purposes. Thus obtaining the orthometric heights by GNSS is the most conventional way of determining these heights. Verifying this goal was the main objective behind the current research. The current research introduces a numerical solution of geoid modeling by applying a surface fitting for a few sparse data points of geoid undulation using minimum curvature surface (MCS). The MCS is presented for deriving a system of linear equations from boundary integral equations. To emphasize the precise applicability of the MCS as a tool for modeling the geoid in an area using GPS/leveling data, a comparison study between EGM2008 and MCS geoid models, is performed. The obtained results showed that MCS technique is a precise tool for determining the geoid in Egypt either on regional and/or local scale with law distortion at check points.

Derivation of Orthometric Heights from GPS Measured Heights Using Geometric Technique and EGM96 Model

As a result of wide spread use of satellite based positioning techniques, especially Global Positioning System (GPS), a greater attention has been focused on precise determination of geoid models with an aim to replace the classical leveling with Global Navigation Satellite System (GNSS) measurements. In this study, geometric technique of deriving orthometric height from GPS survey along a profile and the use of EGM 96 geoid model for deriving orthometric height from GPS data (using GNSS solution software) are discussed. The main focus of the research is to critically examine the potentials of these methods with a view to establishing the optimum technique as an alternative to classical differential levelling. From the results obtained, the standard errors are 1.453m and 1.450m for EGM 96 model and the geometrical approach respectively. From the graphical representation of the residuals from the two methods, it was observed that the two curves suddenly became sinusoidal from station 9 (corresponding to SB08 in the tables). This similarity pattern of the residuals makes it difficult to draw a conclusive judgment between the two methods examined; however, from the standard errors, it could be inferred that the geometrical technique gave a better result over EGM 96 model. Key Words: Geometrical Interpolation, EGM 96 Model, Orthometric Heights, Ellipsoidal Heights

Improved geoid model determination for Japan from GRACE and a regional gravity field model

Abstract A highly improved gravimetric geoid model for Japan, JGEOID2008, is developed on a 1×1.5 arc-minute grid by combining a GRACE (the Gravity Recovery And Climate Experiment)-based global geopotential model, GGM02C, surface (land and ship-borne) gravity measurements, and an altimetry-derived marine gravity model, KMS2002. In the combination, a semidiscrete two-dimensional wavelet analysis/reconstruction method is employed, selecting the spatial wavelength signals of the highest quality out of the respective data sets. Intercomparison with GPS/leveling geoid undulations shows substantial improvement of JGEOID2008 over the previous model, JGEOID2004, and reveals that the systematic errors at long wavelengths contained in JGEOID2004 have been effectively removed. Deviations of JGEOID2008 from the mean sea surface height at tidal stations on isolated islands were comparable to the differences in the sea surface dynamic heights (SSDH) in the Japan Sea, the Nansei Islands and the Izu Island chain from that in Tokyo Bay. The deviations show good agreement with SSDH features estimated from oceanographic observation, indicating that JGEOID2008 has an accuracy within 10 cm. The geoid model is strongly expected to serve as a reference in ocean dynamics studies.