scholarly journals Note on the representation of the earth's surface by means of spherical harmonics of the first three degrees

In my paper on “The Gravitational Stability of the Earth,” dynamical arguments were adduced in favour of the hypothesis that the distribution of density within the earth is such that the surfaces of equal density present, in addition to the inequalities depending upon the diurnal rotation, other inequalities which can be specified by spherical harmonics of the first, second, and third degrees. If this is the case, the surface of the earth, by which I mean the surface of the lithosphere, should present corresponding inequalities, and so also should the equipotential surfaces. Analytically, if the density ρ is given by an equation of the form ρ = f 0 ( r ) + ϵ 1 f 1 ( r )S 1 + ϵ 2 f 2 ( r )S 2 + ϵ 3 f 3 ( r )S 3 , (1) where f 0 ( r ), f 1 ( r ), . . . are functions of the distance r from the centre, S 1 , S 2 , S 3 are spherical surface harmonics of degrees indicated by the suffixes, and ϵ 1 , ϵ 2 , ϵ 3 are small coefficients, then the surface should have an equation of the form r = a + α 1 S 1 + α 2 S 2 + α 3 S 3 , (2) where a and α 1 , α 2 , α 3 are constants, and the α 's are small. The elevations and depressions of the lithosphere should be, at least in their main features, expressible by a formula of this type. The actual elevations and depressions are difficult to determine, because all that can be found by observation is the amount of elevation above, or depression below, a particular equipotential surface, the geoid , or the surface of the ocean, continued beneath the continents. For a first approximation the potential due to such a distribution of density as is expressed by (1) within a surface expressed by (2) would be given by formulæ of the type V = F 0 ( r ) + β 1 F 1 ( r )S 1 + β 2 F 2 ( r )S 2 + β 3 F 3 ( r )S 3 , ( r < a )

1980 ◽  
Vol 78 ◽  
pp. 161-163
Author(s):  
P. Melchior

Note: This short review was given during the Symposium at the request of a number of participants to serve as an elementary introduction to the discussion of the results obtained from different Earth models. It was not prepared in advance and is not an original contribution.The first approximation of the figure of a fluid planet is obtained by assuming hydrostatic equilibrium with respect to its gravitational self-attraction and the centrifugal force. When the speed of rotation is not too fast, the equipotential surfaces may be considered to an excellent approximation as ellipsoids of revolution. It is easy to show that these hydrostatic equipotential surfaces are surfaces of equal density (p).


The equations of motion of a satellite are given in a general form, account being taken of the precession and nutation of the earth. The main part of the paper deals with the motion arising from the gravitational field of the earth, expressed as a general expansion in spherical harmonics. By evaluating the partial derivatives in Lagrange’s planetary equations, • expressions are obtained for the rates of change of the orbital elements. Particular consideration is given to the form of the expressions for the secular terms arising from the first four zonal harmonics.


1. Any estimate of the rigidity of the Earth must be based partly on some observations from which a deformation of the Earth’s surface can be inferred, and partly on some hypothesis as to the internal constitution of the Earth. The observations may be concerned with tides of long period, variations of the vertical, variations of latitude, and so on. The hypothesis must relate to the arrangement of the matter as regards density in different parts, and to the state of the parts in respect of solidity, compressibility, and so on. In the simplest hypothesis, the one on which Lord Kelvin’s well-known, estimate was based, the Earth is treated as absolutely incompressible and of uniform density and rigidity. This hypothesis was adopted to simplify the problem, not because it is a true one. No matter is absolutely incompressible, and, the Earth is not a body of uniform density. It cannot be held to be probable that it is a body of uniform rigidity. But when any part of the hypothesis, e. g ., the assumption of uniform density, is discarded, the estimate of rigidity is affected. Different estimates are obtained when different laws of density are assumed. Again, whatever hypothesis we adopt as regards the arrangement of the matter, so long as we consider the Earth to be absolutely incompressible and of uniform rigidity, different estimates of this rigidity are obtained by using observations of different phenomena. Variations of the vertical may give one value, variations of latitude a notably different value. It follows that “the rigidity of the Earth” is not a definite physical constant. But there are two determinate constant numbers related to the methods that have been used for obtaining estimates of the rigidity of the Earth. One of these numbers specifies the amount by which the surface of the Earth yields to forces of the type of the tide-generating attractions of the Sun and Moon. The other number specifies the amount by which the potential of the Earth is altered through the rearrangement of the matter within it when this matter is displaced by the deforming influence of the Sun and Moon. If we adopt the ordinarily-accepted theory of the Figure of the Earth, the so-called theory of “fluid equilibrium,” and if we make the very probable assumption that the physical constants of the matter within the Earth, such as the density or the incompressibility, are nearly uniform over any spherical surface having its centre at the Earth’s centre, we can determine both these numbers without introducing any additional hypothesis as to the law of density or the state of the matter. We shall find, in fact, that observations of variations of latitude lead to a determination of the number related to the inequality of potential, and that, when this number is known, observations of variations of the vertical lead to a determination of the number related to the inequality of figure. [ Note added , December 15, 1908.—This statement needs, perhaps, some additional qualification. It is assumed that, in calculating the two numbers from the two kinds of observations, we may adopt an equilibrium theory of the deformations produced in the Earth by the corresponding forces. If the constitution of the Earth is really such that an equilibrium theory of the effects produced in it by these forces is inadequate, we should expect a marked discordance of phase between the inequality of figure produced and the force producing it. Now Hecker’s observations, cited in § 6 below, show that, in the case of the semidiurnal term in the variation of the vertical due to the lunar deflexion of gravity, the agreement of phase is close. If, however, an equilibrium theory is adequate, as it appears to be, for the semidiurnal corporeal tide, a similar theory must be adequate for the corporeal tides of long period and for the variations of latitude.]


The magnetic field generated in the core of the Earth is often represented by spherical harmonics of the magnetic potential. It has been found from looking at the equations of spherical harmonics, and from studying the values of the spherical harmonic coefficients derived from data from Magsat, that this is an unsatisfactory way of representing the core field. Harmonics of high degree are characterized by generally shorter wavelength expressions on the surface of the Earth, but also contain very long wavelength features as well. Thus if it is thought that the higher degree harmonics are produced by magnetizations within the crust of the Earth, these magnetizations have to be capable of producing very long wavelength signals. Since it is impossible to produce very long wavelength signals of sufficient amplitude by using crustal magnetizations of reasonable intensity, the separation of core and crustal sources by using spherical harmonics is not ideal. We suggest that a better way is to use radial off-centre dipoles located within the core of the Earth. These have several advantages. Firstly, they can be thought of as modelling real physical current systems within the core of the Earth. Secondly, it can be shown that off-centred dipoles, if located deep within the core, are more effective at removing long wavelength signals of potential or field than can be achieved by using spherical harmonics. The disadvantage is that it is much more difficult to compute the positions and strengths of the off-centred dipole fields, and much less easy to manipulate their effects (such as upward and downward continuation). But we believe, along with Cox and Alldredge & Hurwitz, that the understanding that we might obtain of the Earth’s magnetic field by using physically reasonable models rather than mathematically convenient models is very important. We discuss some of the radial dipole models that have been proposed for the nondipole portion of the Earth’s field to arrive at a model that agrees with observations of secular variation and excursions.


The propagation of radio waves of very low frequency to great distances is conveniently treated by regarding the space between the earth and the ionosphere as a wave-guide. Several authors have found that the least attenuated modes are profoundly affected by the earth’s curvature. This effect is investigated for several models of the ionosphere. It is found, in particular, that for frequencies greater than about 30 kc/s some modes are possible for which the energy is concentrated in a region near the base of the ionosphere, and the field strength near the ground is small. It is useful to think of such modes as being composed of waves repeatedly reflected at the inside spherical surface of the ionosphere, the rays being chords of this sphere. By analogy with sound waves these modes are called ‘whispering gallery modes’. The theory uses wave admittance and reflexion coefficient variables because these satisfy differential equations which are convenient for integration using a digital computer. The curvature of the earth is allowed for by using the method of the modified refractive index, but the earth’s magnetic field is neglected. Formulae for the m ode condition and the excitation of the various modes by a transmitter are given and discussed. A new way of dealing with an ionosphere having a continuous electron density profile is presented. The results of some numerical calculations are given both for a sharply bounded homogeneous ionosphere and for an exponential profile of electron density.


1735 ◽  
Vol 39 (438) ◽  
pp. 98-105 ◽  

The centrifugal force arising from the diurnal rotation of the Earth, depresseth it at the poles, and renders it protuberant at the Equator


1894 ◽  
Vol 55 (331-335) ◽  
pp. 210-217 ◽  

In a paper which was read before the Royal Society in June, 1890, I showed that the principal phenomena of terrestrial magnetism and the secular changes in its horizontal and vertical components could be explained on the assumption of an electro-dynamic substance (presumably liquid or gaseous) rotating within the crust of the earth in the plane of the ecliptic, and a little slower than the diurnal rotation. By means of some electro-mechanism, new to experimental science, which I termed a magnetarium, the period of backward rotation of the internal electro-dynamic sphere required for the secular variations of the magnetic elements on different parts of the earth’s surface was found to be 960 years, or 22.5 minutes of a degree annually. It was also demonstrated that the inclination of the axes of the electro-dynamic and terrestrial globes to each other of 20° 30', was the cause of the inequality of the declination periods about the same meridian in the northern and southern hemispheres; as instanced in the short period of outward westerly declination at London, and the long period of outward westerly declination at the Cape of Good Hope and St. Helena.


2020 ◽  
Author(s):  
Hadi Amin ◽  
Lars E. Sjöberg ◽  
Mohammad Bagherbandi

&lt;p&gt;According to the classical Gauss&amp;#8211;Listing definition, the geoid is the equipotential surface of the Earth&amp;#8217;s gravity field that in a least-squares sense best fits the undisturbed mean sea level. This equipotential surface, except for its zero-degree harmonic, can be characterized using the Earth&amp;#8217;s Global Gravity Models (GGM). Although nowadays, the satellite altimetry technique provides the absolute geoid height over oceans that can be used to calibrate the unknown zero-degree harmonic of the gravimetric geoid models, this technique cannot be utilized to estimate the geometric parameters of the Mean Earth Ellipsoid (MEE). In this study, we perform joint estimation of W&lt;sub&gt;0&lt;/sub&gt;, which defines the zero datum of vertical coordinates, and the MEE parameters relying on a new approach and on the newest gravity field, mean sea surface, and mean dynamic topography models. As our approach utilizes both satellite altimetry observations and a GGM model, we consider different aspects of the input data to evaluate the sensitivity of our estimations to the input data. Unlike previous studies, our results show that it is not sufficient to use only the satellite-component of a quasi-stationary GGM to estimate W&lt;sub&gt;0&lt;/sub&gt;. In addition, our results confirm a high sensitivity of the applied approach to the altimetry-based geoid heights, i.e. mean sea surface and mean dynamic topography models. Moreover, as W&lt;sub&gt;0&lt;/sub&gt; should be considered a quasi-stationary parameter, we quantify the effect of time-dependent Earth&amp;#8217;s gravity field changes as well as the time-dependent sea-level changes on the estimation of W&lt;sub&gt;0&lt;/sub&gt;. Our computations resulted in the geoid potential W&lt;sub&gt;0 &lt;/sub&gt;= 62636848.102 &amp;#177; 0.004 m&lt;sup&gt;2&lt;/sup&gt;s&lt;sup&gt;-2&lt;/sup&gt; and the semi-major and &amp;#8211;minor axes of the MEE, a = 6378137.678 &amp;#177; 0.0003 m and b = 6356752.964 &amp;#177; 0.0005 m, which are 0.678 and 0.650 m larger than those axes of the GRS80 reference ellipsoid, respectively. Moreover, a new estimation for the geocentric gravitational constant was obtained as GM = (398600460.55 &amp;#177; 0.03) &amp;#215; 10&lt;sup&gt;6&lt;/sup&gt; m&lt;sup&gt;3&lt;/sup&gt;s&lt;sup&gt;-2&lt;/sup&gt;.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document