An alternative direct method towards mean dynamic topography computations

An alternative “direct method” to “mean dynamic topography” (MDT) computations using satellite altimetry-derived “mean sea surface” (MSS) and “global geopotential model” (GGM), without direct application of the geoid, is devised. The developed approach, which is based on derivation of an equipotential surface of the gravity field of the Earth that fits to global MSS in least squares sense, is formulated via a constrained optimization problem. The validity of our method is numerically tested by computing a global MDT model based on DNSC08 MSS model and EGM2008 GGM as input data.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

Update of the precision geoid determination in Korea

From the late 1990s, many studies on local geoid construction have been made in South Korea. However, the precision of the previous geoid has remained about 15 cm due to distribution and quality problems of gravity and GPS/levelling data. Since 2007, new land gravity data and GPS/levelling data have been obtained through many projects such as the Korean Land Spatilaization, Unified Control Point and Gravity survey on the Benchmark. The newly obtained data are regularly distributed to a certain degree and show much better improvement in their quality. In addition, an airborne gravity survey was conducted in 2008 to cover the Korean peninsula (South Korea only). Therefore, it is expected that the precision of the geoid could be improved. In this study, the new South Korean gravimetric geoid and hybrid geoid are presented based on land, airborne, ship‐borne, altimeter gravity data, geopotential model and topographic data. As for the methodology, the general remove‐restore approach was applied with the best chosen parameters in order to produce a precise local geoid. The global geopotential model EGM08 was used to remove the low‐frequency components using degree and order up to 360 and the short wavelength part of the gravity signal was dealt with by using the Shuttle Radar Topography Mission data. The parameters determined empirically in this study include for Stokes’ integral 0.5° and for Wong‐Gore kernel 110–120°, respectively and 10 km for both the Bjerhammar sphere depth and attenuation factor. The final gravimetric geoid in South Korea ranges from 20–31 m with a precision of 5.45 cm overall compared to 1096 GPS/levelling data. In addition, the South Korean hybrid geoid produces 3.46 cm and 3.92 cm for degrees of fitness and precision, respectively and a better statistic of 2.37 cm for plain and urban areas was achieved. The gravimetric and hybrid geoids are expected to improve further when the refined land gravity data are included in the near future.     

Geostatistical prediction of a local geometric geoid - kriging and cokriging with the use of EGM2008 geopotential model

We investigate prediction abilities of different variants of kriging and different combinations of data in a local geometric (GNSS/leveling based) geoid modeling. In order to generate local geoid models, we have used GNSS/leveling data and EGM2008 geopotential model. EGM2008 has been used twofold. Firstly, it was used as a basic long wave-length trend to be removed from geoid undulation data to generate a residual field of geoid heights modeled later by kriging (remove-restore technique). Secondly, EGM2008-based undulations were used as a secondary variable in a cokriging prediction procedure (as pseudo-observations). Besides the use of EGM2008, the kriging-based local geometric geoid models were generated only on the basis of raw undulations data. Kriging itself was used in two variants, i.e. ordinary kriging and universal kriging for univariate and bivariate cases (cokriging). The quality of kriging-based prediction for all its variants and all data combinations have been investigated on one fixed validation dataset consisting of 86 points and three training data sets characterized by a different density of sampling. Results of this study indicate that incorporation of EGM08 as a long wave-length trend in kriging prediction procedure outperforms cokriging strategy based on incorporation of EGM08 as a secondary spatially correlated variable.     

Investigation of topographic reductions and aliasing effects on gravity and the geoid over Greece based on various digital terrain models

The reduction of gravity-field related quantities (e.g., gravity anomalies, geoid heights) due to the topography plays a crucial role in both geodetic and geophysical applications, since in the former it is an intermediate step towards geoid prediction and in the latter it reveals lateral as well as radial density contrasts and infers the geology of the area under study. The computations are usually carried out by employing a DTM and/or a DBM, which describe the topography and bathymetry, respectively. Errors in these DTMs/DBMs will introduce errors in the computed topographic effects, while poor spatial resolution of the topography and bathymetry models will result in aliasing effects to both gravity anomalies and geoid heights, both influencing the accuracy of the estimated solutions. The scope of this work is twofold. First, a validation and accuracy assessment of the SRTM 3″ (90 m) DTM over Greece is performed through comparisons with existing global models as well as with the Greek 450 m national DTMs. Whenever a misrepresentation of the topography is identified in the SRTM data, it is “corrected” using the local 450 m DTM. This process resulted in an improved SRTM DTM called SRTMGr, which was then used to determine terrain effects to gravity field quantities. From the fine-resolution SRTMGr DTMs, coarser models of 15″, 30″, 1′, 2′ and 5′ have been determined in order to investigate aliasing effects on both gravity anomalies and geoid heights by computing terrain effects at variable spatial resolutions. From the results acquired in two test areas, it was concluded that SRTMGr provides similar results to the local DTM making the use of other older global DTMs obsolete. The study for terrain aliasing effects proved that when high-resolution and accuracy gravity and geoid models are needed, then the highest possible resolution DTM should be employed to compute the respective terrain effects. Based on the results acquired from two the test areas a corrected SRTMGr DTM has been compiled for the entire Greek territory towards the development of a new gravimetric geoid model. Results from that analysis are presented based on the well-known remove-compute-restore method, employing land and marine gravity data, EGM08 as a reference geopotential model and the SRTMGr DTM for the computation of the RTM effects.     

High Resolution Constituents of the Earth’s Gravitational Field

During the General Assembly of the European Geosciences Union in April 2008, the new Earth Gravitational Model 2008 (EGM08) was released with fully-normalized coefficients in the spherical harmonic expansion of the Earth’s gravitational potential complete to degree and order 2159 (for selected degrees up to 2190). EGM08 was derived through combination of a satellite-based geopotential model and 5 arcmin mean ground gravity data. Spherical harmonic coefficients of the global height function, that describes the surface of the solid Earth with the same angular resolution as EGM08, became available at the same time. This global topographical model can be used for estimation of selected constituents of EGM08, namely the gravitational potentials of the Earth’s atmosphere, ocean water (fluid masses below the geoid) and topographical masses (solid masses above the geoid), which can be evaluated numerically through spherical harmonic expansions. The spectral properties of the respective potential coefficients are studied in terms of power spectra and their relation to the EGM08 potential coefficients is analyzed by using correlation coefficients. The power spectra of the topographical and sea water potentials exceed the power of the EGM08 potential over substantial parts of the considered spectrum indicating large effects of global isostasy. The correlation analysis reveals significant correlations of all three potentials with the EGM08 potential. The potential constituents (namely their functionals such as directional derivatives) can be used for a step known in geodesy and geophysics as the gravity field reduction or stripping. Removing from EGM08 known constituents will help to analyze the internal structure of the Earth (geophysics) as well as to derive the Earth’s gravitational field harmonic outside the geoid (geodesy).     

Ability of the EGM2008 high degree geopotential model to calculate a local geoid model in Valencia,Eastern Spain

A new generation of global geopotential models (GGM) is being developed. These solutions offer a file with fully-normalized spherical harmonic coefficients of the Earth’s gravitational potential up to a degree greater than 2000 with very low commission errors. This paper analyses the recent Earth Gravitational Model EGM2008, developed up to degree and order 2159 with additional coefficients to degree 2190 and order 2159, which means recovering the gravitational field up to approximately 20 km wavelengths. 223 GPS/levelling points of the new Spanish High Precision Levelling Network in the Valencia region (Eastern Spain) are used as external tool for evaluation in that particular region. The same evaluation has been performed to other different global (EGM96 and EIGENCG03C), continental (EGG97), regional (IGG2005 and IBERGEO2006) and local (GCV07) geoid models for comparison purposes only. These comparisons show that EGM2008 is the geoid model that best fits to the GPS/levelling data in that region.     

Detection of Local Geoid Deformations by Gravity Disturbances from GPS/Gravimetric Observations

This paper deals with a method for detection of local geoid deformations; as a consequence, the methods main application concerns geoid adjustment to GPS/levelling points. This is based on the fact that these points should present no local geoid deformation to avoid errors in the adjustments. These type of miscalculations would lead to an incorrect adjustment and result in further errors in subsequent studies with GPS in the proximity at the point with local deformation.The method proposed is based on predictions of gravity disturbance from geoid undulations using Poisson integral with modified kernel, and its comparison with the gravity disturbance from GPS and gravimetric observations.The use of gravity disturbance instead of gravity anomalies has been chosen since gravity disturbance is a quantity derived from GPS and not from levelling. The loss of accuracy arising with a local height reference system is therefore theoretically avoided as far as the differences in geodetic reference systems regarding positions of gravity measurements and coefficients of the global models are accounted for.Extended numerical tests using computed geoidal undulations and the corresponding gravity disturbances obtained from the geopotential model GPM98cr computed up to degree 720 illustrate the validity of the proposed method and its usefulness as local geoid deformations detection tool.Finally, the method is tested using real GPS/Gravimetric data and geoid models IBERGEO95 and EGG97 with good results.     

Validation of recent geopotential models in Tierra Del Fuego

This work presents a validation study of global geopotential models (GGM) in the region of Fagnano Lake, located in the southern Andes. This is an excellent area for this type of validation because it is surrounded by the Andes Mountains, and there is no terrestrial gravity or GNSS/levelling data. However, there are mean lake level (MLL) observations, and its surface is assumed to be almost equipotential. Furthermore, in this article, we propose improved geoid solutions through the Residual Terrain Modelling (RTM) approach. Using a global geopotential model, the results achieved allow us to conclude that it is possible to use this technique to extend an existing geoid model to those regions that lack any information (neither gravimetric nor GNSS/levelling observations). As GGMs have evolved, our results have improved progressively. While the validation of EGM2008 with MLL data shows a standard deviation of 35 cm, GOCO05C shows a deviation of 13 cm, similar to the results obtained on land.     

Comparison of the qualities of recent global and local gravimetric geoid models in Iran

A number of regional gravimetric geoid models have recently been determined for the Iran area, and a common problem is to select the best model, e.g. for engineering applications. A related problem is that in order to improve the local geoid models, the selection of the best Global Geopotential Model (GGM) model for the region is essential, to be used in a combined solution from GGM and local gravimetric data. We discuss these problems by taking advantage of 260 GPS/levelling points as an external tool for validation of different global and local geoid models in the absolute and relative senses. By using relative comparisons of the height differences between precise levelling and GPS/geoid models we avoid possible unknown systematic effects between the different types of observables.The study shows that the combination of the newly released GRACE model (GGM02C) with EGM96 geoid model fits the GPS/levelling data in Iran with the best absolute and relative accuracy among the GGMs. Among the local geoid models, the newly gravimetric geoid model IRG04 agrees considerably better with GPS/levelling than any of the other recent local geoid models. Its rms fit with GPS/levelling is 55 cm. Hence, we strongly recommend the use of this new model in any surveying engineering or GPS/levelling projects in the area.     

Use of EGM08 model and shuttle radar topography mission data for geoid computation in the State of Rio de Janeiro,Brazil: A case study with Voronoi/Delaunay discretisations

The EGM08 geopotential model complete to degree and order 2159 was used in a remove-compute-restore (RCR) method for the geoid computation in the State of Rio de Janeiro, Brazil. Terrain and indirect effect corrections were computed using a 6-arcsec resolution DTE, derived from the TOPODATA Project (Shuttle Radar Topography Mission data) raised by the National Institute for Space Research. INPE, Brazil. We applied Voronoi/Delaunay discretisations for discrete Stokes integration. In these schemes, target area is partitioned into polygons/triangles, respectively, and the computation is carried out by point-wise numerical integration and no gridding is mandatory. For both procedures, the cells were produced using either observed gravity data combined with gridded Bouguer derived information. Particularly in Delaunay scheme, as the gravity anomalies are interpolated into the triangular cells, and geoid undulations are computed for their vertices, Stokes function singularity was gone. Externally estimated errors resulting from a comparison with GPS/leveling data were presented for both the schemes and classical ones, as well as for the EGM08 undulations. They yielded RMS differences equal to 0.105 m, 0.110 m, 0.110 m, 0.115 m and 0.228 m, respectively, for Voronoi, Delaunay, Voronoi/Delaunay with gridded-data alone and EGM08, computed between 32 GPS/leveling points.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

Transformation between gravimetric and GPS/levelling-derived geoids using additional gravity information

The transformation from the gravimetric to the GPS/levelling-derived geoid using additional gravity information for the covariance function of geoid height differences has been investigated in a test area in south-western Canada. A “corrector surface” model, which accounts for datum inconsistencies, long-wavelength geoid errors, vertical network distortions and GPS errors, has been constructed using least-squares collocation. The local covariance function of geoid height differences is usually obtained from residual values between the GPS/levelling and gravimetric geoid heights after the elimination of all known systematic distortions. If additional gravity data (in the form of gravity anomalies) are available, the covariance function of geoid height differences can be determined by the following steps: (1) transforming the GPS/levelling-derived geoid heights into gravity anomalies; (2) forming differences between the computed in step 1 and given gravity anomalies; (3) determining the parameters of the local covariance function of the gravity anomaly differences; (4) constructing an analytical covariance model for the geoid height differences from the covariance function of the gravity anomaly differences using the parameters derived in step 3. The advantage of the proposed method stems from the great number of gravity data used to derive the empirical covariance function. A comparison with the least-squares adjustment shows that the standard deviation of the residuals of the predicted geoid height differences with respect to the control point values decreases by 2.4 cm.     

海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.     

Numerical aspects of EGM08-based geoid computations in Fennoscandia regarding the applied reference surface and error propagation

So far the recent Earth's gravity model, EGM08, has been successfully applied for different geophysical and geodetic purposes. In this paper, we show that the computation of geoid and gravity anomaly on the reference ellipsoid is of essential importance but error propagation of EGM08 on this surface is not successful due to downward continuation of the errors. Also we illustrate that some artefacts appear in the computed geoid and gravity anomaly to lower degree and order than 2190. This means that the role of higher degree harmonics than 2160 is to remove these artefacts from the results. Consequently, EGM08 must be always used to degree and order 2190 to avoid the numerical problems.     

An estimation of the height system bias parameter <Emphasis Type="Italic">N</Emphasis><Subscript>0</Subscript> using least squares collocation from observed gravity and GPS-levelling data

This paper deals with the analysis of gravity anomaly and precise levelling in conjunction with GPS-Levelling data for the computation of a gravimetric geoid and an estimate of the height system bias parameter N_o for the vertical datum in Pakistan by means of least squares collocation technique. The long term objective is to obtain a regional geoid (or quasi-geoid) modeling using a combination of local data with a high degree and order Earth gravity model (EGM) and to determine a bias (if there is one) with respect to a global mean sea surface. An application of collocation with the optimal covariance parameters has facilitated to achieve gravimetric height anomalies in a global geocentric datum. Residual terrain modeling (RTM) technique has been used in combination with the EGM96 for the reduction and smoothing of the gravity data. A value for the bias parameter N_o has been estimated with reference to the local GPS-Levelling datum that appears to be 0.705 m with 0.07 m mean square error. The gravimetric height anomalies were compared with height anomalies obtained from GPS-Levelling stations using least square collocation with and without bias adjustment. The bias adjustment minimizes the difference between the gravimetric height anomalies with respect to residual GPS-Levelling data and the standard deviation of the differences drops from 35 cm to 2.6 cm. The results of this study suggest that N_o adjustment may be a good alternative for the fitting of the final gravimetric geoid as is generally done when using FFT methods.     

Improving airborne strapdown vector gravimetry using stabilized horizontal components

Integrating the deflections of the vertical along the flight line can yield geoid profiles which are valuable in the study of geodesy and geophysics, fortunately, the deflections can be measured directly by vector gravimetry. Airborne vector gravimetry using a Strapdown Inertial Navigation System and the Global Navigation Satellite System (SINS/GNSS) has shown promising results in previous studies. However, the quality of the SINS and GNSS is a major limitation; in particular, the attitude errors induced by the gyros will result in large measurement errors to the horizontal components of the gravity disturbance, and these measurement errors represent the behavior of low-frequency trend. An airborne vector gravimetry method used to remove the bias and low-frequency trends in the gravity disturbance estimated for each survey line has been developed. This method uses the horizontal components of the gravity disturbance computed from EGM2008 (Earth Gravitational Model 2008) as a reference. Firstly, the horizontal measurement results obtained from the gravimeter are divided into high- and low-frequency components according to the resolution of the EGM2008, and then, the bias and low-frequency trends of the low-frequency components are corrected using a linear fit to the EGM2008 reference data. Finally, the ultimate results can be acquired after combining the high-frequency components and the corrected low-frequency components. The data used was obtained from the SGA-WZ, which is the first strapdown airborne gravimeter developed in China. The results of this method are promising. The internal accuracy of the gravity disturbance's horizontal components for repeated survey lines exceeds 3.5 mGal, and the corresponding resolution is approximately 4.8 km based on 160-s data smoothing and an airplane averaging speed of approximately 216 km/h. After applying the WCF (Wavenumber Correlation Filter), the internal accuracy of the horizontal components exceeds 2 mGal. This can satisfy the requirement of the application in geodesy and solid earth geophysics.     

珠峰地区似大地水准面精化与珠峰顶正高的确定

本文使用珠峰及其周边地区的重力数据与SRTM3、1∶50000 DEM、GTOPO30地形数据,以该地区的44个GPS水准点为控制,选择国内外的EGM96、WDM94、IGG05B、DQM2000D和CG03C作为参考重力场模型,采用移去-恢复技术,首次完成了珠峰地区分辨率为2.5′×2.5′高精度似大地水准面,其精度达到±9 cm,并据此推算了珠峰顶高程异常值. 再利用登山线路上的新测重力点与珠峰地形数据,依据严格的重力归算理论及移去-恢复技术,完成了珠峰顶似大地水准面与大地水准面差值计算,即珠峰正常高与正高的换算.     

Regional gravity field modeling based on rectangular harmonic analysis

Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.