Outline
Comptes Rendus

Internal Geophysics (Space Physics)
Earth's oblateness and its temporal variations
Comptes Rendus. Géoscience, Volume 338 (2006) no. 14-15, pp. 1123-1129.

Abstracts

The oblateness is a fundamental property of the Earth under stable rotation. Different oblateness parameters can be defined, each having its physical significance. We retrace the main developments in concept as well as in measurement toward a qualitative understanding of the oblateness and particularly of J2, the so-called oblateness coefficient that ‘glues’ together distinct concepts of oblateness. We examine the numerical values and discuss the dynamical and geometrical consequences of the oblateness. Finally, we report the recent findings of the variations in J2, including an outstanding multi-year anomaly that occurred since 1998, leading to a discussion of possible cause(s) in terms of large-scale mass transports in the Earth's atmosphere–hydrosphere–cryosphere–solid Earth–core system.

L'aplatissement de la Terre en rotation est un paramètre important, car il renseigne sur les propriétés globales de la planète. Dans cet article, nous décrivons la signification de ce paramètre et présentons des résultats de mesure de l'aplatissement terrestre et de son évolution temporelle par les satellites géodésiques depuis le milieu des années 1970. Nous discutons l'interprétation des fluctuations observées, en particulier l'anomalie de grande ampleur observée vers 1998, dont la cause est très probablement liée à une redistribution de masse à grande échelle dans le système océans–atmosphère–cryosphère–terre solide.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crte.2006.09.014
Keywords: Earth's oblateness, $ {J}_{2}$, Mass redistribution, Post-glacial rebound
Mot clés : Aplatissement de la terre, $ {J}_{2}$, Redistribution de masse, Rebond postglaciaire

Benjamin Fong Chao 1

1 Space Geodesy Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
@article{CRGEOS_2006__338_14-15_1123_0,
     author = {Benjamin Fong Chao},
     title = {Earth's oblateness and its temporal variations},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {1123--1129},
     publisher = {Elsevier},
     volume = {338},
     number = {14-15},
     year = {2006},
     doi = {10.1016/j.crte.2006.09.014},
     language = {en},
}
TY  - JOUR
AU  - Benjamin Fong Chao
TI  - Earth's oblateness and its temporal variations
JO  - Comptes Rendus. Géoscience
PY  - 2006
SP  - 1123
EP  - 1129
VL  - 338
IS  - 14-15
PB  - Elsevier
DO  - 10.1016/j.crte.2006.09.014
LA  - en
ID  - CRGEOS_2006__338_14-15_1123_0
ER  - 
%0 Journal Article
%A Benjamin Fong Chao
%T Earth's oblateness and its temporal variations
%J Comptes Rendus. Géoscience
%D 2006
%P 1123-1129
%V 338
%N 14-15
%I Elsevier
%R 10.1016/j.crte.2006.09.014
%G en
%F CRGEOS_2006__338_14-15_1123_0
Benjamin Fong Chao. Earth's oblateness and its temporal variations. Comptes Rendus. Géoscience, Volume 338 (2006) no. 14-15, pp. 1123-1129. doi : 10.1016/j.crte.2006.09.014. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2006.09.014/

Version originale du texte intégral

1 Introduction

The rotating Earth is oblate, that is, it is slightly ‘flat’ in the North Pole–South Pole direction, compared to the slightly ‘bulging’ Equator. This is the result of the hydrostatic balance between the dominant gravitational force, which wants to pull the Earth into a spherically symmetric configuration, and the centrifugal force due to Earth's rotation, which wants to expel mass away from the rotating axis but in the end only manages to modify the Earth into an slightly oblate body.

Quantitatively, this oblateness is about 1 part in 300, which is very close (but see below) the ratio of the centrifugal force on the equator to the gravitational force. This is by far the Earth's largest deviation from a spherically symmetric body. There are certain thermodynamic, but secondary, processes that cause departures from the rotational-hydrostatic equilibrium. Sustained by the internal heat engine and manifested as external gravity anomalies reflecting lateral heterogeneity of internal mass distribution, these deviations are in general no more than parts per million in relative terms.

Yet all these deviations are not static or constant; they change with time. The rotation of the Earth is itself changing over geological time, and the aforementioned mass heterogeneities also vary on timescales upwards from millions of years. On more human timescales, there operates a myriad of dynamic processes that involve mass redistribution in or on the Earth, from tides to atmosphere–ocean circulations, to internal phenomena like earthquakes, post-glacial rebound and core flows. For the Earth, these changes are typically on the order of parts per billion at the largest [9,24]. The present article is a story about the oblateness in particular and how and why it changes with time, where we examine the geophysical implications.

2 The Earth's oblateness parameters and their inter-relationships

As long as the Earth is a 3-D body, we shall use the word oblateness to describe its off-spherical shape. Traditionally, the term ‘flatness’, or ‘ellipticity’, has been used; these names are imprecise because the Earth, of course, is not ‘flat’, and it is a 2-D geometric object only when we try to draw it on paper.

There are several parameters in use to describe the oblateness; each one has its significance depending on the application in question. For simplicity let's for the moment assume the Earth is axially symmetric, or a body of revolution and so essentially a 2-D body, which is a good approximation.

By equating the general expression of the spherical harmonic expansion of the external gravitational potential field V with that of the mass distribution of a body, one concludes that the spherical harmonic coefficients, or Stokes coefficients, of V are simply normalized multi-poles of the density function of the body [4,5]. When specialized, the degree-2 Stokes coefficients are related to the body's inertia tensor elements through a set of equations known as generalized MacCullagh's formulas. In particular, the degree-2 zonal (order 0) Stokes coefficient is given by:

J2=(CA)/Ma2(1)
Symbol namesake in honor of Sir Harold Jeffreys [19] and called the oblateness coefficient for the Earth, J2 has the physical meaning of the difference of the axial or polar (greatest) moment of inertia C with the equatorial (least) moment of inertia A, normalized by Ma2 where M and a are the Earth's mass and mean Equatorial axis, respectively. Its corresponding term in the harmonic expansion of V is the dominant term next to the ‘monopole’ term representing the total mass [4].

We can express J2 in the following form:

J2=[C/Ma2][(CA)/C]ηH(2)
where ηC/Ma2 is a fundamental functional of the Earth's internal structure, and H(CA)/C is called the dynamic oblateness, which can be determined from the observation of the astronomical precession of the Earth [27]. The Earth's η can be readily determined by knowing the values of J2 and H. However, this has not been feasible for other planets because their H's are generally unresolved.

The ‘geopotential’ field is V, modified by the centrifugal potential, i.e. V12a2ω2, where ω is the angular speed of Earth's rotation. If one approximates the equipotential surface, known as the geoid, to an oblate spheroid of revolution, then one obtains the geoid oblateness for the Earth which can be given by the Clairaut's first relationship (for a review, see [27]):

f=(ac)/a=32J2+12m(3)
where c is the mean axial or polar axis of the (ellipsoidal) Earth, and ma3ω2/GM is the ratio of the centrifugal force aω2 on the equator to gravity GM/a2 (where G is the gravitational constant). Eq. (3) is based on the first-order theory, which is sufficient for the present purpose for the Earth (for high-order formulas see, e.g., [23]). For planets with much larger m (for example, Jupiter, see below), the second-order effects become rather significant [29].

Two more concepts of oblateness can be defined at this point: suppose the Earth is under rotational-thermo-hydrostatic equilibrium. A hypothetical hydrostatic geoid oblateness fH representing an idealized Earth can be defined; to first order fH can be found by [18]:

fH=(5/2)m/[1+(25/4)(11.5η)2](4)
Under such equilibrium, the Earth's geometric surface would conform to and coincide with the geoid, so the geometric oblateness, similarly defines as Eq. (3), is simply equal to f. That is why f is often referred to as defining the ‘figure’ of the Earth. In reality, due to its heterogeneities, the Earth's geometric oblateness (as properly defined or approximated) would depart slightly from f. Furthermore, the Earth's true f also departs slightly from fH, where the departure signifies interesting geophysical dynamics sustaining non-equilibrium. The lateral heterogeneity and non-equilibrium configuration of the Earth also manifest themselves in the (relatively small) difference between the two equatorial principal moments of inertia A and B. In the above we have assumed that the Earth is axially symmetric where A=B, which would be the case if the Earth were an otherwise spherically symmetric body subject only to an axial rotation. The real Earth, of course, is not so (we will further discuss this below.) In fact, the strict definition of J2 is [C(A+B)/2]/Ma2, to which Eq. (2) is only an approximation, or valid only for axial-symmetric Earth.

Now let us examine the numerical values. m is known to be 3.46775×10−3, or 1/288.371, close to 1 part in 300 or 1/300. According to Eq. (3), half of it contributes directly to the geoid oblateness f. For f, the remaining contribution comes from 32J2, which, of course, shares the same dynamic origin as m, i.e. Earth's rotation. J2 is measured from satellite geodesy (see below) to a high accuracy, 1.082626×10−3, about one third of 1/300. So, the two terms in (3) contribute almost the same amount to f, i.e., 50% each, and f itself becomes close to 1/300, at f=3.35281×10−3 or 1/298.257. Finally, the dynamic oblateness H in Eq. (2) is observed to be 3.27379×10−3 or 1/305.456, again close to 1/300.

Are these matchings in values just fortuitous? From dynamical considerations, one can rightly ‘guess’ that all parameters should be on the order of m, which they indeed are. However, upon closer examination as follows, they do not necessarily have to have such similar values, so in a sense the latter is fortuitous.

For a reasonable Earth configuration, we should have Hf because of its moderate sensitivity to the internal density profile, although it is well recognized that H would be somewhat less than f because of the smaller oblateness of the interior layers (due to smaller centrifugal force) and the higher density toward the center of the Earth (and hence proportionally lower importance in contributing to the moment of inertia) [27]. Putting this condition into Eqs. (2) and (3), we see the following: the two terms would contribute near-equal shares to f in Eq. (3) and hence all the values would be close to m or 1/300, only if η=1/3.

The interesting, but certainly not out of the ordinary, fact is that, knowing J2 and H in Eq. (2), the Earth in reality has η=0.33069, indeed almost exactly 1/3! This of course does not have to be the case, but one does expect a η value somewhat less than 0.4, that of a uniform-density sphere, for a ‘reasonable’ centrally-heavy, terrestrial planet body such as the Earth.

Based on the PREM Earth model (Preliminary Reference Earth Model [15]) derived from seismological data, the Earth should have an estimated hydrostatic fH of 1/299.66 [29], about 0.5% smaller than the observed. This corresponds to a hydrostatic J2 of 1.0722×10−3, about 1% smaller than the observed. On the other hand, Liu and Chao [21] formulated the relation between A, B and the two Stokes coefficients of degree 2 and order 2. Using the gravity-observed values for the latter, they get BA=7.260×10−6Ma2, amounting to 69.4-m difference in the equivalent geoidal semi-major axis and semi-minor axis on the Equator. Although only 1/150 that of CA, this amount is comparable to the non-hydrostatic portion in CA, as pointed out by [17]. They concluded that the non-hydrostatic portion of the three principal moments of inertia A, B, C only describes a triaxial body and appears to have no preference in orientation. As far as the aforementioned excess oblateness over the hydrostatic value is concerned, this does not favor the notion that this excess oblateness is a remnant, lagging ‘memory’ of the past, as the Earth slows down due to the tidal braking.

3 Comparative planetology

For a contrast, let us compare the Earth with the giant planet Jupiter. Jupiter has a faster rotation and a much larger mass, and hence larger radius and gravity. Its m=0.0892=1/11.2. We can expect that the geoid oblateness f and the dynamic oblateness H to be similar to m, but not necessarily very close in value. The observed J2=0.01469. Adopting second-order formulas [29], which are more accurate than Eq. (2), f=0.0649=1/15.4. Assuming rotational-hydrostatic equilibrium, η=0.254 (cf. Eq. (4)), indicating that, not surprisingly, Jupiter is somewhat more centrally-heavy than the Earth. The derived H=(1/η)J2=0.0578=1/17.3. We further expect that the geometric oblateness is the same as f, except possibly for some small departures from rotational-hydrostatic equilibrium.

In another extreme example, let us consider a non-rotating, uniform-density body not under hydrostatic equilibrium (hence the shape sustained by its internal material strength), such as an asteroid. Then there exists an analytical, but complex, relationship between the spherical harmonic coefficients of gravity and geometrical shape [7]. For the present discussion, let us further assume a special case where the body is a slightly oblate spheroid. Then, letting m=0 in Eq. (2), we have the geoid oblateness f=32J2. Since the body is not under hydrostatic equilibrium, the geometric oblateness is not equal to f; rather, according to Eq. (9) of [7], it equals 53f. Finally, the dynamic oblateness H=(1/η)J2, when η=25 (for a uniform spherical body), equals 53f, the same as the geometric oblateness, as expected.

4 Consequences of oblateness

We live in the Earth's gravity field, controlled by the dominant monopole term GM/a2. We hardly notice any consequence of the oblateness of the Earth (or for that matter the rotational centrifugal force). However, dynamically, the Earth's oblateness is an essential element in our livelihood – it stabilizes our Earth's rotation. The Earth is ‘bombarded’ all the time by countless geophysical agents exerting external torques as well as internal torques or mass transports that exchange angular momenta. Yet its rotation axis hardly changes relative to the Earth-fixed geography. This is not true if it were spherically symmetric: then the crawl of a bug or a firing of a canon, for instance, would completely ‘tumble’ the Earth relative to the (spatially stationary) rotation axis [16,22]. On the other hand, it is well known from classical mechanics that the rotation of a body about its principal axis of the greatest moment of inertia (C) is a stable one. What prevents large shifts of the rotation axis from happening is the extra oblateness in the form of CA with the associated extra angular momentum, which is to be overcome by any geophysical agent that tries to shift the Earth's rotational pole positions. Since the oblateness itself is a consequence of the rotation in the first place, it can be stated that the rotating planet is self-stabilized.

A corollary of the above, but on a less dramatic scale, the Earth's dynamic oblateness under the tidal torques exerted by the Moon and Sun gives rise to the astronomical precession of the Earth's rotation axis in space, and hence is the deciding factor for the precessional period. That in fact is how the dynamic oblateness H is determined. On the same token, H acts as the restoring factor that prescribes the free wobble, known as the Chandler wobble, of the Earth's polar motion. The period of the Chandler wobble would be 1/H days if the Earth were a rigid body, but was found to be significantly lengthened by the Earth's non-rigidity, or finite elasticity [22].

As stated, the geometrical shape of the Earth largely conforms to the oblate geoid. Therefore, the mean equatorial radius and the mean polar radius of the geoid differ by as much as ac=fc21 km. In particular, the global sea level follows closely this oblate geoid, only undulating on top of the geoid geographically no more than 200 m peak-to-peak and temporally less than 10 m or so. The land topography undulates up to 10 km, but largely supported isostatically.

As such, the oblateness also affects various geophysical quantities. For example, in the space geodesy enterprise using near-Earth satellites, the oblateness term resides in all Earth surface geometry that locates the geodetic observatories and altimetric targets. Similarly, the oblateness prevails in the external gravity field that significantly affects the satellite orbits from which geodetic measurements are made. On the Earth surface, together with the centrifugal force field, the oblateness gives the surface gravity a slight latitudinal dependence which is actually the largest term in the surface gravity anomaly on the global scale. In another example, the Earth's elastic free oscillation modes (often excited by large earthquakes) see splitting in their otherwise degenerate characteristic periods due to Earth's oblateness and rotation, completely analogous in the atomic world to the Stark splitting and Zeeman splitting, respectively, as such splitting is determined by the symmetry properties common to different dynamic systems [2].

5 Historical Notes

Sir Isaac Newton, based on his law of gravitation and force laws, was the first to realize that the Earth under rotational equilibrium should possess a non-vanishing oblateness. The value of 1/f favored by him and given in the Principia was 230. Cassini subsequently came up with a negative value, −95, presumably owing to certain systematic errors. The value had evolved [19], since the Peru/Lapland expedition in the 1740s from a value between 179 and 266, to 301, 295, 297.0, and finally in the early 1950s to Sir Harold Jeffreys' 297.1±0.4 which is within 0.4% of the modern value (298.257).

Then came the space age, ushered in by the launch of USSR's Sputnik I spacecraft in October 1957. A month later Sputnik II was launched, and within a few weeks, by monitoring the nodal precession of its orbit in space, our knowledge of J2 grew almost an order of magnitude, to about 0.1% of the modern value. This measurement was arguably one of the very first scientific triumphs of space enterprise.

Today, after nearly half a century of precise orbit determination of dozens of geodesy-quality satellite orbits around the Earth, the Earth's global gravity field has been solved to harmonic degrees as high as 120, among which the average J2 coefficient has been determined to the accuracy of seven significant figures (1.082627×10−3) [20].

Since the 1980s, thanks to the advent of the technique of satellite laser ranging [1], tiny temporal variations around the average value of J2 began to be noted. The variation occurs in the last digit of the above-quoted number and beyond, typically no more than one part in a billion! This will be discussed next.

6 How and why does Earth's J2 change?

Mass transports in the atmosphere–hydrosphere–cryosphere–solid Earth–core system (the ‘Earth system’) occur on all temporal and spatial scales for a host of geophysical and climatic reasons [9,24]. According to Newton's gravitational law, such mass transport will cause the gravity field to change with time, producing time-variable gravity signals.

Increasingly refined models for the Earth's static gravity field in terms of spherical harmonic components have been determined by means of decades of precise orbit tracking data of many geodetic satellites. On top of that, low-degree components of Earth's time-variable gravity have been clearly observed by the space geodetic technique of satellite laser ranging (SLR) [1]. Although tiny in relative terms (no more than 1 part per billion), these variations signify global-scale mass redistribution in the Earth system.

In particular, the lowest-degree zonal harmonic is Earth's oblateness coefficient J2, whose temporal variation was the first to be detected among all gravity components. A ‘secular’ decrease in J2 (over the observed quarter century) was first identified from the SLR satellite nodal precession acceleration. Its main excitation source has since been attributed to the post-glacial rebound (PGR) of the solid Earth [26,30], and for additional secondary causes [3]. Subsequently, many studies reported strong seasonal as well as weaker non-seasonal signals, primarily in J2, but of late also in the next-lowest harmonics [13] and geocenter. The prominent seasonal J2 signals (with primarily annual amplitude 3×10−10) have been correlated with mass transports in the atmosphere, oceans, and land hydrology [8,11,25].

Such was the case until around the turn of the century beginning in 1998, when the SLR data began to reveal that Earth's J2 had suddenly deviated significantly from the PGR secular decreasing trend (at about 2.8×10−11 yr−1). This ‘1998 anomaly’ embarked on a reverse, increasing trend over the following years, before quieting back down to the ‘normal’ decreasing trend. This was reported by [10,12]. Fig. 1a shows an updated time series of the SLR-observed J2, using SLR data from up to nine satellites, with more satellites becoming available with time [12]. Note that the relevant 18.6-yr lunar-driven ocean tide amplitude was set to the value recovered in a 21-yr comprehensive solution for the secular zonal rates, low-degree static and annual terms, and the 18.6-yr and the much smaller 9.3-yr lunar tides. Fig. 1(b) is the same time series but after the removal of (i) the atmospheric contribution calculated according to the global NCEP reanalysis data assuming an inverted-barometer effect, and (ii) the least-squares fit of the remaining seasonal signals, which are attributable to (the poorly known) seasonal mass redistribution in the oceans and land hydrology. The PGR slope (the solid line) and the 1998 anomaly are clearly evident.

Fig. 1

(a) Time series of the SLR-observed, non-tidal J2. (b) Same as (a), but after the removal of NCEP-calculated atmospheric contribution (with inverted-barometer effect) and any remaining seasonal signals.

(a) Variations temporelles de J2 mesurées par les satellites laser (après correction des effets de marées). (b) Identique à (a), après retrait des termes saisonniers et de la contribution atmosphérique (calculée à partir des données de pression atmosphérique NCEP, dans le cas d'une réponse de type « baromètre inverse » de l'océan).

A number of possible causes for the 1998 anomaly was speculated by Cox and Chao [12], including oceanic water mass redistribution, melting of polar ice sheets and high-latitude glaciers, global sea level rise, and material flow in the fluid core. Dickey et al. [14] emphasized and demonstrated the importance of the melting of high-latitude glaciers. Chao et al. [10] report an oceanographic event that took place in the extratropic North + South Pacific basins that was found to match remarkably well with the time evolution of the J2 anomaly; the phenomenon appears to be part of the Pacific Decadal Oscillation immediately following the episode of the 1997–1998 El Niño.

The difficulty in identifying the definite cause(s) in the above w.r.t. J2 stems from the extremely low geographical resolution of the zonal harmonic function in question, namely the degree-2 Legendre function. Thus, a positive J2 anomaly only tells us that a net transport of mass from higher latitudes to lower latitudes (across the nodal latitude of the degree-2 Legendre function, namely ±35.3°) has occurred, in either or both Northern and Southern Hemispheres. For example, an equivalent of as much as 3000 km3 of water melted from Greenland and spread into the oceans would be needed to produce the first half of the J2 anomaly where the relative change is +7×10−11 per year; but we have no way of telling without other ancillary evidences or observations. On the other hand, the space gravity mission of GRACE (launched in March 2002, with an expected lifetime of over 5 yr), using the satellite-to-satellite tracking (SST) technique, is yielding gravity information at much higher geographical resolution than the SLR-based information. For example, GRACE is able to detect centimeter level water-height-equivalent mass changes over an area of about 1000 km across from month to month [28]. However, considering the relatively weak sensitivity of the SST to the longest-wavelength (lowest-degree) gravity components, GRACE's utility in measuring the variation of J2 in particular awaits to be seen. The same applies to future gravity missions of GOCE (using gravity gradiometer) and other SST measurements in ‘follow-on’ gravity missions under planning.

7 Relationship between Earth's rotation and J2 change

As stated, the Earth's oblateness arises from its rotation; the rotational-hydrostatic relationship, to first order, is given in Eq. (4), where the oblateness is proportional to m, which is in turn proportional to ω2. Therefore,

J˙2J2=2ω˙ω(5)

For example, the Earth's secular spin-down due to the tidal braking would lead to a secular ‘rounding’ of the Earth (barring possible temporal retardation under viscosity), thus decreasing J2. Numerically, at the tidal-braking rate of ω˙=6.5×10−22 rads−2, that decreasing rate of J2 is about 6.1×10−13 yr−1, contributing only 2% of the observed decreasing rate of J2 (see above).

On the other hand, any change in J2 will cause ω to change, as dictated by the conservation of angular momentum for the Earth. For instance, a decreasing J2 means a faster spin (analogous to a spinning skater pulling arms closer to the body), and vice versa. That effect can be shown to be [6]:

ω˙ω=2.01J˙2(6)
where the coefficient 2.01 is evaluated from Earth parameters. For example, the decrease in J2, at the rate of 6.1×10−13 yr−1 due to the tidal braking of ω given above, will in turn feedback to cause ω to increase, but only by as little as 2.8×10−24 rads−2, or 0.053 μs in the equivalent length-of-day per year. That is completely negligible in today's measurement. On the other hand, the observed J2 rate of change 2.8×10−11 yr−1 (see above) presumably speeds up the Earth rotation by 2.4 μs in the equivalent length-of-day per year, which is still negligible.

8 Epilogue

Although numerically small, the oblateness is a fundamental property of the Earth under stable rotation. Its existence and cause, its dynamical and geometrical consequences, its values and departures from idealized models, and its temporal evolution due to mass transports in the Earth system are all fascinating topics in geophysics, which reveal insights towards the understanding of the structure and dynamical behavior of the Earth. The measurement and monitoring of the Earth's oblateness have been a triumph as well as a scientific target of the modern space geodesy. As one sees deeper and finer into the Earth's oblateness, there is little doubt that the Earth will surprise and further fascinate us with a continuing story unfolding with time.

Acknowledgement

This paper is completed under the support of the NASA Solid Earth program. I am grateful to Christopher Cox for providing Fig. 1.


References

[1] Contributions of Space Geodesy to Geodynamics: Technology (AGU; D.E. Smith; D.L. Turcotte, eds.), Geodynamics Series, vol. 25, Am. Geophys. Union, 1993

[2] B.F. Chao Symmetry and terrestrial spectroscopy, Geophys. J. R. Astron. Soc., Volume 66 (1981), pp. 285-312

[3] B.F. Chao Geodesy is not just for static measurements any more, EOS Trans. Amer. Geophys. Union, Volume 84 (2003), pp. 145-156

[4] B.F. Chao On inversion for mass distribution from global (time-variable) gravity field, J. Geodyn., Volume 39 (2005), pp. 223-230 | DOI

[5] B.F. Chao; R.S. Gross Changes in the Earth's rotation and low-degree gravitational field induced by earthquakes, Geophys. J. R. Astron. Soc., Volume 91 (1987), pp. 569-596

[6] B.F. Chao; W.P. O'Connor Effect of a uniform sea level change on the Earth's rotation and gravitational field, Geophys. J., Volume 93 (1988), pp. 191-193

[7] B.F. Chao; D.P. Rubincam The gravitational field of Phobos, Geophys. Res. Lett., Volume 16 (1989), pp. 859-862

[8] B.F. Chao; A.Y. Au Temporal variation of Earth's zonal gravitational field caused by atmospheric mass redistribution: 1980–1988, J. Geophys. Res., Volume 96 (1991), pp. 6569-6575

[9] B.F. Chao; V. Dehant; R.S. Gross; R.D. Ray; D.A. Salstein; M.M. Watkins; C.R. Wilson Space geodesy monitors mass transports in global geophysical fluids, EOS Trans. Am. Geophys. Union, Volume 81 (2002), pp. 47-250

[10] B.F. Chao; A.Y. Au; J.-P. Boy; C. Cox Time-variable gravity signal of an anomalous redistribution of water mass in extratropic Pacific during 1998–2002, Geochem. Geophys. Geosyst., Volume 4 (2003), p. 1096 | DOI

[11] M.K. Cheng; B.D. Tapley Seasonal variations in low degree zonal harmonics of the Earth's gravity field from satellite laser ranging observations, J. Geophys. Res., Volume 104 (1999), pp. 2667-2681

[12] C. Cox; B.F. Chao Detection of a large-scale mass redistribution in the terrestrial system since 1998, Science, Volume 297 (2002), pp. 831-833

[13] C.M. Cox, B.F. Chao, A. Au, Interannual and Annual Variations in the Geopotential Observed Using SLR, in: 14th International Workshop on Laser Ranging, San Fernando, CA, USA, June 2004

[14] J.O. Dickey; S. Marcus; O. de Viron; I. Fukumori Recent Earth oblateness variations: Unraveling climate and postglacial rebound effects, Science, Volume 298 (2002), pp. 1975-1977

[15] A. Dziewonski; D. Anderson Preliminary reference earth model, Phys. Earth. Planet. Inter., Volume 25 (1981), pp. 297-356

[16] T. Gold Instability of the Earth's axis of rotation, Nature, Volume 175 (1955), pp. 526-529

[17] P. Goldreich; A. Toomre Some remarks on polar wandering, J. Geophys. Res., Volume 74 (1969), pp. 2555-2567

[18] H. Jeffreys The Earth, its Origin, History and Physical Constitution, Cambridge University Press, New York, 1959

[19] D. King-Hele A Tapestry of Orbits, Cambridge University Press, New York, 1992

[20] F. Lemoine et al., The development of the joint NASA and NIMA geopotential model EGM96, NASA/TP 1998-206861, 1998

[21] H.S. Liu; B.F. Chao The Earth's equatorial principal axes and moments of inertia, Geophys. J. Int., Volume 106 (1991), pp. 699-702

[22] W.H. Munk; G.J.F. MacDonald The Rotation of the Earth, Cambridge University Press, New York, 1960

[23] S.M. Nakiboglu Hydrostatic figure and related properties of the Earth, Geophys. J. R. Astron. Soc., Volume 57 (1979), pp. 639-648

[24] NAS, J.O. Dickey (Eds.), Satellite Gravity and the Geosphere, National Research Council, Washington, D.C., 1997

[25] R.S. Nerem; R.J. Eanes; P. Thompson; J.L. Chen Observations of annual variations of the Earth's gravitational field using satellite laser ranging and geophysical models, Geophys. Res. Lett., Volume 27 (2000), pp. 1783-1786

[26] D.P. Rubincam Postglacial rebound observed by Lageos and the effective viscosity of the lower mantle, J. Geophys. Res., Volume 89 (1984), pp. 1077-1087

[27] F.D. Stacey Physics of the Earth, Brookfield Press, Brisbane, 1992

[28] B.D. Tapley; S. Battadpur; M. Watkins; C. Reigber The gravity and climate recovery experiment: Mission overview and early results, Geophys. Res. Lett., Volume 31 (2004), p. L09607 | DOI

[29] C. Yoder Astrometric and geodetic properties of Earth and the solar system (T.J. Ahrens, ed.), Global Earth Physics: A Handbook of Physical Constants, AGU Ref. Shelf 1, Am. Geophys. Union, Washington, DC, 1995

[30] C.K. Yoder; J.G. Williams; J.O. Dickey; B.E. Schutz; R.J. Eanes; B.D. Tapley Secular variations of Earth's gravitational harmonic J2 coefficient from Lageos and nontidal acceleration of Earth rotation, Nature, Volume 303 (1983), p. 757


Comments - Policy