Bjerhammar Sphere Research Articles

Error bounds for the spectral approximation of the potential of a homogeneous almost spherical body

Several kinds of approximation of the gravitational potential of a homogeneous body by truncated spherical harmonics series are in use in physical geodesy. However, only one of them is capable of a representation converging to the true potential in the whole layer between the Brillouin sphere and the Bjerhammar sphere of the body. We aim at providing various majorizations, namely upper bounds, of the error with the double purpose of proving explicitly the convergence in the sense of different norms and of giving computable bounds, that might be used in numerical studies. The first aim is reached for all the norms. For the second, however, it turns out that among the bounds, when applied to the example of the terrain correction of the Earth, only those referring to the mean absolute error and the mean squared error at the level of Brillouin sphere of minimum radius give significant and useful results. In order to make the computation an easy exercise, a simple approximate formula has been developed requiring only the use of the distribution function of the heights of the surface of the body with respect to the Bjerhammar sphere.

A numerical test of the topographic bias

Abstract In 1962 A. Bjerhammar introduced the method of analytical continuation in physical geodesy, implying that surface gravity anomalies are downward continued into the topographic masses down to an internal sphere (the Bjerhammar sphere). The method also includes analytical upward continuation of the potential to the surface of the Earth to obtain the quasigeoid. One can show that also the common remove-compute-restore technique for geoid determination includes an analytical continuation as long as the complete density distribution of the topography is not known. The analytical continuation implies that the downward continued gravity anomaly and/or potential are/is in error by the so-called topographic bias, which was postulated by a simple formula of L E Sjöberg in 2007. Here we will numerically test the postulated formula by comparing it with the bias obtained by analytical downward continuation of the external potential of a homogeneous ellipsoid to an inner sphere. The result shows that the postulated formula holds: At the equator of the ellipsoid, where the external potential is downward continued 21 km, the computed and postulated topographic biases agree to less than a millimetre (when the potential is scaled to the unit of metre).

Local gravity field modeling using spherical radial basis functions and a genetic algorithm

Spherical Radial Basis Functions (SRBFs) can express the local gravity field model of the Earth if they are parameterized optimally on or below the Bjerhammar sphere. This parameterization is generally defined as the shape of the base functions, their number, center locations, bandwidths, and scale coefficients. The number/location and bandwidths of the base functions are the most important parameters for accurately representing the gravity field; once they are determined, the scale coefficients can then be computed accordingly. In this study, the point-mass kernel, as the simplest shape of SRBFs, is chosen to evaluate the synthesized free-air gravity anomalies over the rough area in Auvergne and GNSS/Leveling points (synthetic height anomalies) are used to validate the results. A two-step automatic approach is proposed to determine the optimum distribution of the base functions. First, the location of the base functions and their bandwidths are found using the genetic algorithm; second, the conjugate gradient least squares method is employed to estimate the scale coefficients. The proposed methodology shows promising results. On the one hand, when using the genetic algorithm, the base functions do not need to be set to a regular grid and they can move according to the roughness of topography. In this way, the models meet the desired accuracy with a low number of base functions. On the other hand, the conjugate gradient method removes the bias between derived quasigeoid heights from the model and from the GNSS/leveling points; this means there is no need for a corrector surface. The numerical test on the area of interest revealed an RMS of 0.48mGal for the differences between predicted and observed gravity anomalies, and a corresponding 9cm for the differences in GNSS/leveling points.

Update of the precision geoid determination in Korea

ABSTRACTFrom 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.

The convergence problem of collocation solutions in the framework of the stochastic interpretation

The convergence of the collocation solution to the true gravity field is a problem defined long ago; some results were derived, in particular by T. Krarup, already in 1981. The problem is taken up again in the context of the stochastic interpretation of collocation theory and some new results are derived, showing that, when the potential T can be really continued down to a Bjerhammar sphere, reasonable convergence results hold true.

The Application of GPS Technique in Determining the Earth's Potential Field

Two approaches to determining the Earth's external potential field by using GPS technique are proposed. The first one is that the relation between the geopotential difference and the light signal's frequency shift, between two separated points, is applied. The second one is that the spherical harmonic expansion series and a new technique dealing with the "downward continuation" problem are applied. Given the boundary value provided by GPS "geopotential frequency shift" on the Earth's surface, the Earth's external field could be determined based on the "fictitious compress recovery" method. Given the boundary value derived by on-board GPS technique on the satellite surface, the Earth's external field could be determined by using a new technique for solving the "downward continuation" problem, which is also based on the "fictitious compress recovery" method. The main idea of the "fictitious compress recovery" is that an iterative procedure of "compress " and "recovery" between the given boundary (the Earth's surface or the satellite surface) and the surface of Bjerhammar sphere is executed and a fictitious field is created, which coincides with the real field in the domain outside the Earth. Simulation tests support the new approach.

A comparison of methods for the inversion of airborne gravity data

Four integral-based methods for the inversion of gravity disturbances, derived from airborne gravity measurements, into the disturbing potential on the Bjerhammar sphere and the Earth’s surface are investigated and compared with least-squares (LS) collocation. The performance of the methods is numerically investigated using noise-free and noisy observations, which have been generated using a synthetic gravity field model. It is found that advanced interpolation of gravity disturbances at the nodes of higher-order numerical integration formulas significantly improves the performance of the integral-based methods. This is preferable to the commonly used one-point composed Newton–Cotes integration formulas, which intrinsically imply a piecewise constant interpolation over a patch centered at the observation point. It is shown that the investigated methods behave similarly for noise-free observations, but differently for noisy observations. The best results in terms of root-mean-square (RMS) height-anomaly errors are obtained when the gravity disturbances are first downward continued (inverse Poisson integral) and then transformed into potential values (Hotine integral). The latter has a strong smoothing effect, which damps high-frequency errors inherent in the downward-continued gravity disturbances. An integral method based on the single-layer representation of the disturbing potential shows a similar performance. This representation has the advantage that it can be used directly on surfaces with non-spherical geometry, whereas classical integral-based methods require an additional step if gravity field functionals have to be computed on non-spherical geometries. It is shown that defining the single-layer density on the Bjerhammar sphere gives results with the same quality as obtained when using the Earth’s topography as support for the single-layer density. A comparison of the four integral-based methods with LS collocation shows that the latter method performs slightly better in terms of RMS height-anomaly errors.

Globally covering a-priori regional gravity covariance models

Abstract. Gravity anomaly data generated using Wenzel’s GPM98A model complete to degree 1800, from which OSU91A has been subtracted, have been used to estimate covariance functions for a set of globally covering equal-area blocks of size 22.5° × 22.5° at Equator, having a 2.5° overlap. For each block an analytic covariance function model was determined. The models are based on 4 parameters: the depth to the Bjerhammar sphere (determines correlation), the free-air gravity anomaly variance, a scale factor of the OSU91A error degree-variances and a maximal summation index, N, of the error degree-variances. The depth of Bjerhammar-sphere varies from -134km to nearly zero, N varies from 360 to 40, the scale factor from 0.03 to 38.0 and the gravity variance from 1081 to 24(10µms-2)2. The parameters are interpreted in terms of the quality of the data used to construct OSU91A and GPM98A and general conditions such as the occurrence of mountain chains. The variation of the parameters show that it is necessary to use regional covariance models in order to obtain a realistic signal to noise ratio in global applications.Key words. GOCE mission, Covariance function, Spacewise approach`

Stability investigations of various representations of the gravity field

A number of problems in physical geodesy are improperly posed in the sense that the approximate solution is not continuously dependent on the given observations. A transformation of the minimum norm least squares estimates for different types of gravity representations to the spectral domain, via a singular value decomposition, allows an easy insight into the stability situation of these representations. For the density layer or point mass model, Bjerhammar’s method, and the least squares collocation model the filter expressions are derived which control their stability. It is shown in which way the shape of the filters is influenced by the a priori variance of the measurement noise, the choice of the radius of the Bjerhammar sphere, and the degree variance model.

On the existence of solutions for the method of Bjerhammar in the continuous case

The method of Bjerhammar is studied in the continuous case for a sphere. By varying the kernel function, different types of unknowns (u*) are obtained at the internal sphere (the Bjerhammar sphere). It is shown that a necessary condition for the existence of u* is that the degree variances (σn2) of the observations are of an order less than n−2. According to Kaula’s rule this condition is not satisfied for the earth’s gravity anomaly field (σn2=n−1) but well for the geopotential (σn2=n−3).