Investigations on the accuracy concerning the derivation of absolute gravimetric height anomalies and deflections of the vertical based on the theory of Molodensky

The investigations refer to the compartment method by using mean terrestrial free air anomalies only. Three main error influences of remote areas (distance from the fixed point >9°) on height anomalies and deflections of the vertical are being regarded:

1. The prediction errors of mean terrestrial free air anomalies have the greatest influence and amount to about ±0″.2 in each component for deflections of the vertical and to ±3 m for height anomalies;
2. The error of the compartment method, which originates from converting the integral formulas of Stokes and Vening-Meinesz into summation formulas, can be neglected if the anomalies for points and gravity profiles are compiled to 5°×5° mean values.
3. The influences of the mean gravimetric correction terms of Arnold—estimated for important mountains of the Earth by means of an approximate formula—on height anomalies may amount to 1–2 m and on deflections of the vertical to 0″0.5–0″.1, and, therefore, they have to be taken into account for exact calculations.

The computations of errors are carried out using a global covariance function of point free air anomalies.     

Comparison of mean gravity anomalies at elevation with corresponding ground anomalies

In support of requirements for the U.S. Air Force Cambridge Research Laboratories, gravity anomalies have been upward continued to several elevations in different areas of the United States. One area was 34⁰ to 40⁰ N in latitude and 96⁰ to 103⁰ W in longitude, generally called the Oklahoma area. The computations proceeded from 26, 032 point anomalies to the prediction of mean anomalies in 14, 704, 2.5′×2.5′ blocks and 9,284, 5′×5′ blocks. These anomalies were upward continued along 28 profiles at 5′ intervals for every 30′ in latitude and longitude. These anomalies at elevations were meaned in various patterns to form mean 30′×30″, 1⁰×1⁰, 5⁰×5⁰ blocks. Comparisons were then made to the corresponding ground values. The results of these comparisons lead to practical recommendations on the arrangement of flight profiles in airborne gravimetry.     

Assessment of ZTD derived from ECMWF/NCEP data with GPS ZTD over China

The accuracy and feasibility of computing the zenith tropospheric delays (ZTDs) from data of the European Center for Medium-Range Weather Forecasts (ECMWF) and the United States National Centers for Environmental Prediction (NCEP) are studied. The ZTDs are calculated from ECMWF/NCEP pressure-level data by integration and from the surface data with the Saastamoinen model method and then compared with the solutions measured from 28 global positioning system (GPS) stations of the Crustal Movement Observation Network of China (CMONOC) for 1 year. The results are as follows: (1) the error of the integration method is 1–3 cm less than that of the Saastamoinen model method. The agreement between the ECMWF ZTD and GPS ZTD is better than that between NCEP ZTD and GPS ZTD; (2) the bias and root mean square difference (RMSD), especially the latter, have a seasonal variation, and the RMSD decreases with increasing altitude while the variation with latitude is not obvious; and (3) when using the full horizontal resolution of 0.5° × 0.5° of the ECMWF meteorological data in place of a reduced 2.5° × 2.5° grid, the mean RMSD between GPS and ECMWF ZTD decreases by 4.5 mm. These results illuminated the accuracy and feasibility of computing the tropospheric delays and establishing the ZTD prediction model over China for navigation and positioning with ECMWF and NCEP data.     

Determination of Indian absolute geodetic datum by least-square solution technique

The Everest spheroid, 1830, in general use in the Survey of India, was finally oriented in an arbitrary manner at the Indian geodetic datum in 1840; while the international spheroid, 1924, in use for scientific purposes; was locally fitted to the Indian geoid in 1927. An attempt is here made to obtain the initial values for the Indian geodetic datum in absolute terms on GRS 67 by least-square solution technique, making use of the available astro-geodetic data in India, and the corresponding generalised gravimetric values at the considered astro-geodetic points, as derived from the mean gravity anomalies over1°×1° squares of latitude and longitude in and around the Indian sub-continent, and over5° equal area blocks covering the rest of the earth’s surface. The values obtained independently by gravimetric method, were also considered before actual finalization of the results of the present determination.     

A comparison of altimeter and gravimetric geoids in the Tonga Trench and Indian Ocean areas

A gravimetric geoid computed using different techniques has been compared to a geoid derived from Geos-3 altimeter data in two 30°×30° areas: one in the Tonga Trench area and one in the Indian Ocean. The specific techniques used were the usual Stokes integration (using 1°×1° mean anomalies) with the Molodenskii truncation procedure; a modified Stokes integration with a modified truncation method; and computations using three sets of potential coefficients including one complete to degree 180. In the Tonga Trench area the standard deviation of the difference between the modified Stokes’ procedure and the altimeter geoid was ±1.1 m while in the Indian Ocean area the difference was ±0.6 m. Similar results were found from the 180×180 potential coefficient field. However, the differences in using the usual Stokes integration procedure were about a factor of two greater as was predicted from an error analysis. We conclude that there is good agreement at the ±1 m level between the two types of geoids. In addition, systematic differences are at the half-meter level. The modified Stokes procedure clearly is superior to the usual Stokes method although the 180×180 solution is of comparable accuracy with the computational effort six times less than the integration procedures.     

Detailed gravimetric geoid for the North Atlantic

A detailed gravimetric geoid in the North Atlantic Ocean, named DGGNA-77, has been computed, based on a satellite and gravimetry derived earth potential model (consisting in spherical harmonic coefficients up to degree and order 30) and mean free air surface gravity anomalies (35180 1°×1° mean values and 245000 4′×4′ mean values). The long wavelength undulations were computed from the spherical harmonics of the reference potential model and the details were obtained by integrating the residual gravity anomalies through the Stokes formula: from 0 to 5° with the 4′×4′ data, and from 5° to 20° with the 1°×1° data. For computer time reasons the final grid was computed with half a degree spacing only. This grid extends from the Gulf of Mexico to the European and African coasts. Comparisons have been made with Geos 3 altimetry derived geoid heights and with the 5′×5′ gravimetric geoid derived byMarsh andChang [8] in the northwestern part of the Atlantic Ocean, which show a good agreement in most places apart from some tilts which porbably come from the satellite orbit recovery.     

Determination of geoidal heights and deflections of the vertical for the hellenic area using heterogeneous data

Mean gravity anomalies, deflections of the vertical, and a geopotential model complete to degree and order180 are combined in order to determine geoidal heights in the area bounded by [34°≦ϕ≤42°, 18°≦λ≦28°]. Moreover, employing point gravity anomalies simultaneously with the above data, an attempt is made to predict deflections of the vertical in the same area. The method used in the computations is least squares collocation. Using empirical covariance functions for the data, the suitable errors for the different sources of observations, and the optimum cap radius around each point of evaluation, an accuracy better than±0.60m for geoidal heights and±1″.5 for deflections of the vertical is obtained taking into account existing systematic effects. This accuracy refers to the comparison between observed and predicted values.     

The gains of small circular,square and rectangular filters for surface waves on a sphere

Meissl has derived weighting functions for converting point gravity anomaly degree variances into mean anomaly variances over a circular cap on a sphere. If the cap is sufficiently small so that the cap on a sphere degenerates into a circle on a plane, the problem may be considered that of the gain of a circular filter for a surface wave whose wave number depends on the spherical harmonic degree. The Meissl weights then become replaced by diffraction integrals of optical physics. The expected gain for a square filter for waves coming from random directions is derived and shown to be close to the gain of a circular filter with the same area. The expected gains and cross-gains for rectangular filters are also derived. When weighted by an anomaly degree variance model, these gains and cross-gains can be used to determine rectangular anomaly variances and covariances for arbitrary bandwidths. Using the Tscherning-Rapp model, analytic gravity anomaly variances and covariances are calculated for 1°×1° blocks.     

A comparison of recent Earth gravitational models with emphasis on their contribution in refining the gravity and geoid at continental or regional scale

Since the publication of the Earth gravitational model (EGM)96 considerable improvements in the observation techniques resulted in the development of new improved models. The improvements are due to the availability of data from dedicated gravity mapping missions (CHAMP, GRACE) and to the use of 5′ × 5′ terrestrial and altimetry derived gravity anomalies. It is expected that the use of new EGMs will further contribute to the improvement of the resolution and accuracy of the gravity and geoid modeling in continental and regional scale. To prove this numerically, three representative Earth gravitational models are used for the reduction of several kinds of data related to the gravity field in different places of the Earth. The results of the reduction are discussed regarding the corresponding covariance functions which might be used for modeling using the least squares collocation method. The contribution of the EIGEN-GL04C model in most cases is comparable to that of EGM96. However, the big difference is shown in the case of EGM2008, due not only to its quality but obviously to its high degree of expansion. Almost in all cases the variance and the correlation length of the covariance functions of data reduced to this model up to its maximum degree are only a few percentages of corresponding quantities of the same data reduced up to degree 360. Furthermore, the mean value and the standard deviation of the reduced gravity anomalies in extended areas of the Earth such as Australia, Arctic region, Scandinavia or the Canadian plains, vary between −1 and +1 and between 5 and 10 × 10⁻⁵ ms⁻², respectively, reflecting the homogenization of the gravity field on a regional scale. This is very important in using least squares collocation for regional applications. However, the distance to the first zero-value was in several cases much longer than warranted by the high degree of the expansion. This is attributed to errors of medium wavelengths stemming from the lack of, e.g., high-quality data in some area.     

Improving Orthometric Heights Using Precise GPS Measurements

The vertical component obtained from the Global Positioning System (GPS) observations is from the ellipsoid (a mathematical surface), and therefore needs to be converted to the orthometric height, which is from the geoid (represented by the mean sea level). The common practice is to use existing bench marks (around the four corners of a project area and interpolate for the rest of the area), but in many areas bench marks may not be available, in which case an existing geoid undulation is used. Present available global geoid undulation values are not generally as detailed as needed, and in many areas they are not known better than ±1 to ±5 m, because of many limitations. This article explains the difficulties encountered in obtaining precise geoid undulation with some example computations, and proposes a technique of applying corrections to the best available global geoid undulations using detailed free-air gravity anomalies (within a 2° × 2° area) to get relative centimeter accuracy. Several test computations have been performed to decide the optimal block sizes and the effective spherical distances to compute the regional and the local effects of gravity anomalies on geoid undulations by using the Stokes integral. In one test computation a 2° × 2° area was subdivided into smaller surface elements. A difference of 37.34 ± 1.6 cm in geoid undulation was obtained over the same 2° × 2° area when 1° × 1° block sizes were replaced by a combination of 5' × 5' and 1' × 1' subdivision integration elements (block sizes).