Access the full text.

Sign up today, get DeepDyve free for 14 days.

Artificial Satellites
, Volume 44 (3) – Jan 1, 2009

/lp/de-gruyter/solving-the-topographic-potential-bias-as-an-initial-value-problem-mVuETBDaUk

- Publisher
- de Gruyter
- Copyright
- Copyright © 2009 by the
- ISSN
- 0208-841X
- eISSN
- 2083-6104
- DOI
- 10.2478/v10018-009-0021-8
- Publisher site
- See Article on Publisher Site

If the gravitational potential or the disturbing potential of the Earth be downward continued by harmonic continuation inside the Earth's topography, it will be biased, the bias being the difference between the downward continued fictitious, harmonic potential the real potential inside the masses. We use initial value problem techniques to solve for the bias. Fit, the solution is derived for a constant topographic density, in which case the bias can be expressed by a very simple formula related with the topographic height above the computation point. Second, for an arbitrary density distribution the bias becomes an integral along the vertical from the computation point to the Earth's surface. No topographic masses, except those along the vertical through the computation point, affect the bias. (To be exact, only the direct indirect effects of an arbitrarily small but finite volume of mass around the surface point along the radius must be considered.) This implies that the frequently computed terrain effect is not needed (except, possibly, for an arbitrarily small innerzone around the computation point) for computing the geoid by the method of analytical continuation. Keywords: analytical continuation, gravity, potential, topographic bias 1. INTRODUCTION If we disregard or remove the effect of the atmosphere, the gravitational potential of the Earth is harmonic in the exterior of the Earth's surface. As this is not the case below the surface, the analytically continued exterior disturbing potential will be biased. Bjerhammar (1962) introduced a technique for quasigeoid determination with analytical continuation of the surface gravity anomaly to an internal sphere (the Bjerhammar sphere). Once the downward continued gravity anomaly on the sphere has been determined, he showed that any related classical formula that integrates the fictitious gravity anomaly on the internal sphere yields the correct function (e.g., gravity anomaly, gradient of gravity anomaly, height anomaly deflections of the vertical) in the exterior of the Earth. (See also Heiskanen Moritz 1967, Sect. 8-10 Moritz 1980, Sect. 45) On the contrary, if the integral formula is applied with the computational point inside the topography, the result will be a fictitious or biased quantity. For example, if Stokes' formula is applied to determine the geoid height, the result will be biased. However, let us imagine that the disturbing potential (the gravity 76 potential minus the normal potential) can be analytically continued to sea level, that the topographic bias can be determined. Then, by subtracting the bias from the analytically continued potential, the true disturbing potential at sea level is obtained, also the geoid height can be determined by applying Bruns' theorem to the potential. In this way it would be possible to determine the geoid in a very direct way by analytical continuation. It is not proved that the downward continued disturbing potential T will exist in the strong sense within the topographic masses. Nevertheless, if we subtract the effect of a Bouguer shell (V B ) from T [ refer to the approximation theorems of KrarupRunge (Krarup 1969; Moritz 180, Sect. 8) Keldysh-Lavrentieff (Lkof 1972; Bjerhammar 1974)], we may expect that the reduced potential T V B is harmonic can be arbitrarily well approximated all the way down to the geoid. We will not be concerned with such approximations in this article. The interested reader can study, e.g. Moritz (1980, Sect. 45). Usually this problem is studied as that of continuing a finite series of spherical harmonics, representing the external potential of the Earth, analytically to sea level (e.g. Cook 1967, Sjöberg 1977 Ågren 2004), but such a series will suffer not only from the topographic bias, but also from the truncation error of the series , possibly, also from the divergence of the series when applied at sea level. Our main concern in this study is the derivation of the topographic bias. Sjöberg (2007) performed such a study by decomposing the topographic potential into two parts (namely those generated by a Bouguer shell the remaining topography, the terrain, respectively,) compared the analytically continued potential with the true potential. This paper resulted in some critical remarks by Vermeer (2008), which were rather concerned with the analytical continuation error of the spherical harmonic series than with the topographic bias as pointed out by Sjöberg (2008). In two related studies by Sjöberg (2009a) (2009b) the topographic bias was studied by employing digital terrain models as well as a Taylor series. In the present study the approach is different: we study the topographic bias as an initial value problem (IVP). An IVP is the problem of solving an ordinary differential equation related with one or several specific conditions (initial values). In the present case the differential equation is of second order, related with four conditions as specified in Eqs. (2)-(5) below. Related studies for the potential itself were performed by Holota (1994) (1996). Usually a function is defined as being harmoni if it satisfies Laplace's equation in a domain, i.e. in a certaieighbourhood of a computation point. This condition, which may be difficult to prove/meet for the potential close to the boundary of topographic masses, is not important for our study, as our definition of topographic bias simply implies that the function does not satisfy Laplace's equation at the computation point. In this study a spherical template will approximate the Earth's topography with arbitrary number sizes of compartments. The IVP will be studied for a point inside the compartment (i.e. away from its boundary). The limiting case of the topographic bias (when the size of the compartment goes to zero) is studied in Sect. 2.3. The concept of the topographic bias is important in determining the geoid from the external gravity field by the method of analytical continuation. 77 2. FORMULATING SOLVING THE IVP In this section we will fit present some basic formulas in formulating the IVP, then we proceed to solve the problem under different assumptions on the topographic density distribution. 2.1. BASIC FORMULAS The topographic potential bias of the analytically continued disturbing potential T * is defined by the formula b T* T , (1) where T* T are the analytically continued the true disturbing potentials, respectively. A reasonable assumption on the Earth surface will be that it is stahaped. Outside at the Earth's surface the bias its fit radial derivative are both zero, we can therefore introduce the following initial values to warrant continuity at the boundary: b 0 , if r (2) b r 0 , if r , (3) where r are the geocentric radii of the computation point the Earth's surface along the radius vector through the computation point, respectively. Moreover, as the downward continuation (dwc) is based on harmonic functions, T * is harmonic also after dwc inside the topographic masses, it therefore obeys Laplace's equation there, i.e. T* 0 for r R , (4) is the Laplace operator R is the radius of sea level. On the other h, where within the masses the true disturbing potential T obeys Poisson's differential equation, so that b 4 where point. , if R , (5) is the product of the gravitational constant density at the computation Now, let us decompose the topographic potential VPt at the arbitrarily located computation point P (on or above sea level) into the Bouguer shell potential VPB the terrain/ residual topographic potential VPte : VPt VPB VPte . (6) Here we assume (without loss of generality) that the density distribution of the Bouguer shell is radial symmetric determined by the radial density distribution at the computation point P. We also assume a spherical approximation of the shape of sea level with radius R, the upper radius of the Bouguer shell as given by the radius c of topographic height at P, i.e. c P . Note that Eq. (6) decomposes the topographic potential into the contributions from the Bouguer shell VB the terrain Vte, the sphere of radius c is the exterior surface of the Bouguer shell. In the 78 computation of any derivative of VPt we will regard c as a constant (although this leads to an approximation in the solution of the IVP, which we will discuss in Sect. 2.3). This implies that VPt VPB VPB 0 , if rP (7a) VPt 4 ( P) , if R rP (7b) is density. Eqs. (8) where is the Laplace operator, rP is geocentric radius of P (6), (7a) (7b) yield that VPte VPt VPB 0 , if rP c R rP which, by our definition above, implies that the terrain does not contribute to the topographic bias, except possibly for point P located at the surface with radius c. As the second order derivatives of the potential generally do not exist at the boundary, Laplace's Poisson's differential equations of Eqs. (7a) (7b) are not applicable there (Kellog 1954, p. 153). However, this problem is easily circumvented by simply removing the effect dVPte of a small volume of mass of the terrain around the surface point, such that the reduced terrain potential VPte VPte dVPte becomes harmonic along the radius vector, i.e. VPte 0 for r >R . (9) This implies also that the disturbing potential TP (to be downward continued) is d reduced to TP TP dTPte , the indirect effect dVPte must be added to the computed/downward continued potential in order not to introduce an additional bias. It is obvious that the smaller we make the (arbitrary) mass volume generating dVPte , the smaller will be the sum of the direct indirect effects of this volume (but, admittedly, the more critical/difficult becomes the dwc of TP ). This is even more obvious, if we consider the approximation the exterior surface of this volume by a sloping plane centred at the Earth's surface. In case of a constant topographic density, the terrain effect of this innermost mass volume will be exactly zero, because the potential generated by the masses below the surface of the Bouguer shell (to be subtracted for) exactly compensates the potential of the positive mass contribution above the Bouguer shell. Hence, in practice there is no contribution to the topographic bias of TP from the terrain. All these derivations lead to the result that the topographic bias of Eq. (1) is a function of the mass of the Bouguer shell alone. Below we will determine the topographic bias of this potential in the cases of a constant topographic density ( Sect. 2.2) for a general density distribution (Sect. 2.3). In the solutions that follow we will start from the following tentative general model for the topographic bias of the potential T within the topographic masses: bn r n , (10) 79 where the coefficients bn should be consistent with Laplace's Poisson's equations (4) (5) the initial values provided by Eqs. (2) (3). (Note that the model is chosen under the assumption that the coefficients are constants. Cf. Sect. 2.3.) The Laplace operator working on Eq. (10) yields 2 r r bn r n n(n 1)bn r n 2 . ) (11) Hence, it remains to determine the coefficients bn under the initial values of Eqs. (2)(5) by considering Eq. (11). We will consider both a constant radially variable density distribution in the topography (Sects. 2.2-2.3 2.4, respectively). 2.2. SOLUTION FOR A CONSTANT TOPOGRAPHIC DENSITY be constant within the topography, implying that the left h side of Eq. (5) Let must be a constant for any r in the interval R r , i.e. the left h-side must be independent of r. Hence, from Eq. (11) follows that only terms with b 1 , b0 b2 of the series for b in Eq. (10) are different from zero. Thus one obtains b1 b0 b2 r 2 , r (12) from Eq. (5) follows also that b2 2 / 3 . Furthermore, by inserting Eq. (12) into Eqs. (2) (3) we can also solve for b 1 b0 , resulting in the following solution for the bias: 2r 23 r 32 r 2 , R (13) or, (Sjöberg 2007) b 2 H H r (14a) where r , if 0 ,if , if R r r (14b) The topographic bias solved above is that for an Earth model, whose topography is a spherical template, the solution is found for any point inside a compartment of the template. Although this would be the case for any numerical application, the question arises, what happens to the topographic bias in the limit, when the size of the compartment vanishes? The solution to this problem is discussed in the next subsection. 2.3. THE TOPOGRAPHIC BIAS IN THE LIMITING CASE Let us now use the notations b for the topographic bias surface radius of a compartment within the Earth template model. Then it holds that 2r 23 rP 32 rP2 , (15) where k 2 / 3 rP is the radius of the computation point P with R rP c P . Denoting the real Earth radius by , it follows directly from Eq. (15) that lim b 2r 23 rP 32 rP , as . (16) Hence, b approaches smoothly the limiting case. Although, in this case the condition of Eq. (5) is violated, there is no other option for the topographic bias of the real Earth than that of Eq. (16). 2.4. SOLUTION FOR AN ARBITRARY TOPOGRAPHIC DENSITY DISTRIBUTION We now generalize the topographic density distribution along the radial direction at the computation point to the following equation, consistent with Eq. (11): (r ) an r n 2 , R , (17) where we assume that the series is absolutely convergent the coefficients an govern the density distribution. Also in this case we assume that b can be represented by Eq. (10), by inserting Eqs. (11) (17) into Eq. (5) we arrive at n( n(n 1)bn r n ) n 2 an r n 2 , R , (18) where we can identify most of the coefficients bn . The result is bn 4 an , if n 0 (n 1) (19) Hence, the solution caow be written b1 b0 db , r an rn . n(n 1) n( 1,0 (20a) where the known part becomes db 4 (20b) In order to determine the remaining coefficients we reconsider Eqs. (2) (3). From Eqs. (3), (20a) (20b) we obtain b1 4 1,0 a 1 , n 1 (21) Eq. (2) yields a . n 1,0 b0 (22) The solution is thus provided by Eqs. (20a, b)-(22). It requires that the coefficients an be known. As we assumed that the series of Eq. (17) is absolutely convergent, the convergences of Eqs. (20b), (21) (22) follow. A more general form of the solution is provided by the following proposition. Proposition 1: For R the topographic potential bias can be written b 4 ( x) x2 r d x dx . (23) Proof: Let us fit introduce the abbreviation k ( x) 4 ( x) . (24) We want to verify that the possible solution of the bias b( r ) k ( x) x2 r d x dx (25) satisfies Eqs. (2), (3) (5). Eq. (2) follows directly from Eq. (25). Moreover, from b´(r ) db(r ) dr d k ( x) x2 dx r (26) it follows also that Eq. (3) is satisfied. Finally, as b´´(r ) 2 k ( x) x2 dx k (r ) , d r3 (27) it follows that b(r ) b´´(r ) b( 2 b´(r ) r k (r ) , (28) we have thus verified that the cidate solution b satisfies the proposition. However, it remains to show that b is the one only solution. For that reason we assume that there is also another solution g (r ) b(r ) (29) for some function . As function g must satisfy Eqs. (2), (3) (29), it follows '( that f ( ) f '( ) 0 , f must also be harmonic in the interval R r . However, as a harmonic function attains it extreme values at the boundary, it follows that f vanishes in the whole interval. The solution for f can also be obtained as follows: Let us assume that fnr n . (30) Then it holds that 2 f '(r ) ' r f ''(r ) '' n(n 1) f n r n (31) 0 for As f must satisfy Laplace's equation for any r in the interval, it follows that f n all n except for n = 0 n = 2. Hence, f0 f2r 2 (32) f '(r ) 2 f 2 r . Equating f f ' to zero yields that f0 f2 0, (33) (34) which is exactly what we postulated. Note. As with the solution discussed in Subsect. 2.2, the solution derived in this subsection holds strictly for the template model, while for the real Earth topography, we refer to the limiting case of Subsect.2.3. 3. CONCLUDING REMARKS In the case of a constant topographic density we have shown that the topographic potential bias can be expressed by the simple formula of Eq. (13). For an arbitrary density distribution the formula is also simple, being a one-dimensional integral along the vertical through the computation point. Ågren (2004, Sect. 8.4.3) compared Eq. (13) with Rapp's (1997) approach to correct a series representation in external type spherical harmonics of the geoid height, he concluded that the two techniques are practically identical. By adding the topographic potential bias to the downward continued disturbing potential at sea level one can easily use Bruns' formula to determine the geoid height. The downward continuation of the disturbing potential can be performed by gravimetric data alone (e.g., Sjöberg 2003a). Hence, from a theoretical point of view, by applying the technique of analytical continuation the frequently applied terrain effect is not needed to determine the geoid height. (Admittedly, we have removed restored the effect of an arbitrarily small mass of the terrain at the surface point.) On the other h, the terrain effect can be numerically efficient to smooth the observed gravity anomaly or disturbing potential to allow better interpolation Stokes' integration as well in particular to allow a practical solution to the analytical downward continuation problem. However, the removal restoration (r-r) of the terrain effect must therefore not necessarily be carried out with theoretical rigour, but only in such a way that the r-r technique be consistent.) Finally, it should be stated that this article dealt with the topographic bias of the downward continued disturbing potential, we did not dwell upon the common (Molodensky) problem to compute the disturbing potential at the surface from gravity anomalies. A scheme to carry out all these steps is given by the KTH (acronym for Royal Institute of Technology) method for geoid determination, presented, e.g., in Sjöberg (2003b), the origin of the present article was actually to determine the potential bias of this method with clarity. 83 REFERENCES Ågren J. (2004) Regional geoid determination methods for the era of satellite gravimetry, Doctoral dissertation in Geodesy, Royal Institute of Technology, TRITA-INFRA 04-033, Stockholm. www.infra.kth.se/geo/) Bjerhammar A. (1962) Gravity reduction to a spherical surface. Royal Institute of Technology, Division of Geodesy, Stockholm Bjerhammar A (1974) Discrete approaches to the solution of the boundary value problem of physical geodesy. Paper presented at the International School of Geodesy, Erice, Italy Cook A. H. (1967) The determination of the external gravity field of the Earth from observation of artificial satellites. Geophysical Journal of Royal Astronomy Society, Vol. 13. Heiskanen W.A. H. Moritz (1967). Physical geodesy. W H Freeman Co., San Francisco London Holota P. (1994) Two branches of the Newton potential geoid. In H Suenkel I Maon (eds.): Gravity Geoid. IAG Symposia No. 113, Graz, Springer, 205214 Holota P. (1996) Geoid, Cauchy's problem displacement. In J Segawa, H Fujimoto S Okubo (eds.): Gravity, Geoid Marin Geodesy. IAG Symposia No. 117, Tokyo, Springer, 368-375 Kellog O.D.(1953) Foundations of potential theory, Dover Publ., Inc., New York Krarup T. (1969) A contribution to the mathematical foundation of physical geodesy. Geodaetisk Institut, Meddelelse no. 48, Copenhagen Lkof N.S. (1972) Foundations of modern potential theory, Springer, Berlin, Heidelberg, New York Moritz H. (1980) Advanced physical geodesy. Herbert Wichmann Verlag, Karlsruhe. Rapp R.H. (1997) Use of potential coefficient models for geoid undulation determinations using a spherical harmonic representation of the height anomaly/geoid undulation difference, Journal of Geodesy, Vol. 71, 282-289. Sjöberg L. E. (1977) On the erro of spherical harmonic developments of gravity at the surface of the earth, Depth Geod Sci Rep No. 257, The OSU, Columbus, OH Sjöberg L.E. (2003a) A solution to the downward continuation effect on the geoid determined by Stokes's formula, Journal of Geodesy, Vol. 77, 94-100 Sjöberg L.E. (2003b) A computational scheme to model the geoid by the modified Stokes's formula without gravity reductions, Journal of Geodesy, Vol. 77, 423432 Sjöberg L.E. (2007) The topographic bias by analytical continuation in physical geodesy, Journal of Geodesy, Vol. 81, 345-350 Sjöberg L. E.(2008) Answe to the comments by M. Vermeer on L.E. Sjöberg (2007) "The topographic bias by analytical continuation in physical geodesy" J Geod 81:345-350, Journal of Geodesy, Vol. 82, 451-452 84 Sjöberg L.E. (2009a) The terrain correction in gravimetric geoid computation - is it needed? Geophysical Journal International, Vol. 176, 14-18. Sjöberg L.E. (2009b) On the topographic bias in geoid determination by the external gravity field, Journal of Geodesy, Vol. 83, 967-972. Vermeer M. (2008) Comment on Sjöberg (2007 ) "The topographic bias by analytical continuation in physical geodesy "J Geod 81(5),345-350, Journal of Geodesy, Vol. 82, 445-450. Received: 2009-12-03, Reviewed: 2010-01-13, Accepted: 2010-01-25.

Artificial Satellites – de Gruyter

**Published: ** Jan 1, 2009

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.