Non-Rigid Rotation , Global Geodynamics and Free Nutation

The rotation of the physical Earth not only obeys the law of conservation of angular momentum but follows also the three-finger rule of the right-handed system. The reference frame in the Earth is most essential. The generalized Eulerian equation of motion or the Liouville equation is to represent global geodynamics in general, which is examined directly through the observation of the Earth’s rotation and surface motions, with no additional assumptions after fundamental physics. Rotational acceleration, polar motion, plate tectonics, and seismicity observed on the Earth’s surface indicate that the Earth’s angular momentum is perturbed and adjusting independent of external torques. The internal torque or the rate of change of angular momentum in the absence of external torques represents the Earth’s conservation of perturbed angular momentum for a stable rotation of minimum energy configuration. The internal torques calculated from surface observations exhibit various magnitudes that reflect the overall global geodynamics. The secular components of internal torques are all on the order of 10 ergs, which imply that the Earth’s angular momentum is conserved independent of external torques in a secular sense, and give an order-of-magnitude quantitative constraint to secular geodynamics in the interior. The equation of secular motion in the Earth is generalized for both internal and external perturbations. The equation approximates the global geodynamics that link polar motion with plate motion and mantle flow, while free nutation is induced to conserve the Earth’s angular momentum in the absence of external torques. The gyroscopic effect or gyricity dominates secular motion in the Earth. The kinematics of secular internal torques are consistent with the gross features of plate motion and mantle flow. Secular gyric torques are quadrupolar and continuous to the interior convectively according to the three-finger rule of the right-handed system, while secular rotational torque is toward the north. For a multilayered, deformable, energy-generating and dissipative Earth that allows motion and mass redistribution, stable rotation can be reached only when its major principal, rotation, and instantaneous figure axes are all completely aligned with each other to arrive at the minimum energy configuration of the system. Polar instability consists of incipient polar excitation that gains energy to commence polar instability, as well as polar motion and subsequent plate motion and mantle flow that dissipate energy for polar stability. Angular momentum perturbation via relative angular momentum supplies the excess energy at polar excitation, which can be episodic or continuous, internal and thermal or external and gravitational. Thermal buoyancy perturbs angular momentum, but as a secular motion it is also subject to gyricity, which consequently will guide thermal flow to conserve angular momentum. The Earth’s rotational energy of 2.16 × 10 ergs is the largest energy reservoir that can be directly observed on the Earth’s surface; its slight fluctuation on the order of 10 ergs is sufficient to supply power for any secular geodynamical processes in the Earth. After all, no geodynamics can be independent of the Earth’s perpetual rotation against fundamental physical laws.


Introduction
The Earth is a perpetually rotating heavenly body orbiting in space.The observation of polar motion, plate motion, seismicity, geological deformation, volcanism, as well as the multiple frequency-splits of the Chandler wobble [1,2] indicates motion and mass redistribution is to occur in the solid Earth.The Earth is hence essentially a multilayered non-rigid heavenly body rotating according to the law of conservation of angular momentum, for which an additional physical rule, the three-finger rule of the right-handed system [3], also plays an equally important fundamental role.That is, in a rotating system allowing motion like the Earth, the rate of change of a vector in a reference frame rotating with the system relative to an inertial frame fixed in space differs from that in the inertial frame by a cross product of the angular velocity ω of the rotating frame with the vector.The diffe- rentiation operator d dt ω + ×, therefore, applies to the whole rotating system.This is the gyroscopic effect, which has long been known but ignored in the study of solid Earth geodynamics.
The importance of gyrodynamics comes in light of space dynamics [3].From the observation of the orbiting of space vehicles, it was discovered that the rotation stability of a space vehicle can be significantly affected by the gyroscopic effect if there is energy dissipation in the vehicle [3][4][5][6][7].D'Eleuterio and Hughes [7] hence include the gyroscopic effect, together with that of inertial (mass), dissipative (damping), and stiffness (elastic), as the four mechanical influences on the dynamics of flexible structures, and introduce the concept of gyricity to replace its conventional name under an instrument.Gyricity and gyric are adopted here as convenient shorthand for gyroscopic effect and gyroscopic.As the Earth is a multilayered, deformable, energy-generating, dissipative, and perpetually rotating heavenly body orbiting in space that allows motion and mass redistribution, its rotation dynamics are similar to, if not the same as, those of space vehicles, to obey the fundamental physical laws.The mechanical influences on geodynamics hence also include gyricity.As we shall see below, gyricity indeed dominates the rotation dynamics of a non-rigid Earth.
Any motion in the Earth, fast or slow, instantaneous or secular, is subject to gyricity according to the three-finger rule of the right-handed system, just as that in any other rotating systems orbiting in space.The classic gyricity observed on the Earth is the Coriolis force, which was known to be negligibly small and hence led to the conventional belief that Coriolis effect played a very little role in geodynamics of the solid Earth.This is, analogous to the Earth's rheology or deformability, indeed the case in an instantaneous sense as if the Earth were rigid, but not in a secular sense in a non-rigid Earth.A comparison may help to visualize.The Earth's instanta-neous rotation is at a rate of 1.46 × 10 12 cm/year, while secular motion in the solid Earth is only a few cm/year.As we shall see below, gyricity in the Earth can be of great magnitude, and secular motion in a rotating Earth is not the same as that in an inertial or flat Earth.
Early observations [8][9][10] have already noticed that the torsional deformation of tectonic features on the ocean floors may be attributable to the Coriolis effect; likely are also the rotational inertia of the continents [11], as well as the "vertical axis rotation" of the crustal blocks [12,13].Pan [14,15] has preliminarily studied the (secular) Coriolis and rotational torques in the Earth, and observes that the gross features of the present plate tectonics, such as plate distribution, plate motion, spread of oceanic ridges, and subduction of oceanic slabs, all reflect the quadrupolar nature of gyricity according to the three-finger rule of the right-handed system.This suggests that long-term gyricity may well be the primary driving mechanism for plate motion and mantle flow [15].However, the study was yet preliminary and partial.The present paper concerns an overall examination of the fundamental physics that govern the rotation of the physical Earth, from there global geodynamics starts.It has two foundations: 1) The Liouville equation, as the generalized equation of motion for a rotating Earth orbiting in space that allows motion and mass redistribution, is to represent global geodynamics in general; and 2) direct observation of the Earth's rotation and surface motions, without involving physical properties in the interior.It is practically a study of the Liouville equation directly through the Earth's surface observations according to the three-finger rule of the right-handed system and the law of conservation of angular momentum, while the Earth's non-rigidity is still treated as what the Liouville equation allows.The rotation of a non-rigid Earth is distinctly different from that of a biaxial or slightly triaxial rigid body, which is discussed elsewhere [1,2,14,16,17].

Reference Frame, Liouville Equation and Global Geodynamics
In the study of the rotation of a multilayered non-rigid Earth orbiting in space, the most essential and difficult problem is the choice of reference frame; otherwise it may end up not truly representing the physical Earth.So we repeat a brief review here to emphasize.Munk and MacDonald [18] have already pointed out that it is unlikely to find a truly body-fixed frame in the Earth.Chao [19] further points out that there is an inconsistency between the reference frame used for observation and those for theoretical calculation.Pan [17] observes that all conventional theoretical frames are idealized systems that are not truly located in the physical Earth.The Earth's rotation must refer to an inertial frame fixed in

OPEN ACCESS IJG
space, but which can be assumed to be instantaneously coinciding with a frame located in the Earth without loss of generality [17,18,20].Pan [14][15][16][17] demonstrates that a physically consistent theoretical frame has to be unique, physically located in the Earth, consistent with observation, and always associated with polar motion.A frame close to such conditions is adopted here.As depicted in Figure 1, the ( ) , , a b c frame is the Earth's principal axes prior to polar excitation, referring to which the inertia tensor was diagonal; whereas, the ( ) , , x y z frame is the axes of the Earth's inertia tensor that will appear simultaneously with polar excitation, referring to which the inertia tensor is not diagonal.The ( ) , , x y z frame is the reference frame, which differs from the ( ) , , a b c frame by the axial near-symmetry angle pair ( ) , θ φ , where θ is the deviation angle between the c-and z-axes, and φ is the azimuth angle between corres- ponding equatorial axes of the two systems.The frame is geocentric.Its z-axis or instantaneous figure axis is aligned with the axis of reference [18] or the geographic axis [21] around which the rotation axis physically revolves, but which is, as Pan [2,[14][15][16][17] points out, not a principal axis.Its y-axis is along the direction of secular polar drift, while x-axis is perpendicular to the y-and z-axes in the right-handed system.In this frame, the Earth is axially near-symmetrical and slightly triaxial even it was originally biaxial [2,[14][15][16].The Liouville equation is fully described in the frame.
The Liouville equation or the generalized Eulerian equation of motion [2,[14][15][16][17][18]20,22] governs the rotation of a system that allows motion in the system.The Earth is such a rotating system; this is fundamental physics.The equation states that according to the law of conservation of angular momentum, the rate of change of total angular momentum or the internal torque [14] of a rotating system is equal to the external torque exerted on the system, while motion in the system follows the threefinger rule of the right-handed system.The Liouville equation can either be linearized to examine the polar excitation mechanism or rotation irregularities excited by motion and mass redistribution [1,2,14,[16][17][18]20,22], or otherwise directly adopted as the Earth's equation of motion for what after polar excitation or global geodynamics such as the driving mechanism of polar motion, plate motion and mantle flow [14,15], as well as free nutation in the absence of external torques.An incipient solution of the linearized Liouville equation [14,16,17] predicts the multiple splits of the Chandler frequency, as is confirmed by observation [1,2].Pan [14,15] has also preliminarily demonstrated that two secular internal torques, the Coriolis and rotational torques, are likely the primary driving mechanisms for plate motion and mantle flow.The present paper is an introductory examination of the internal torques in the Liouville equation that represent the overall global geodynamics, as well as an integration of global geodynamics with rotation dynamics.
Let I be the Earth's inertia tensor, ω rotation ve- locity, h relative angular momentum arising from motion and mass redistribution [2,14,15,17], and L external torque.The Liouville equation in the ( ) , , x y z frame rotating relative to an instantaneously coinciding inertial frame fixed in space is, then, where the overdot designates d dt in the ( ) , , x y z frame.Equation (1) represents global geodynamics in general.The left-side terms in the equation are the rate of change of angular momentum, a reaction to external torque that reflects the central-force dynamics in the Earth [23].The terms can thus be called the internal torques of the system [14,15,24], though their dimension is not exactly in dyne-cm as that of a torque.Throughout this paper, we shall use erg as their unit for the convenience of direct conversion to energy, and for a more refined orderof-magnitude resolution.Observation of non-tidal rotational acceleration, polar motion, plate motion, and seismicity on the Earth's surface indicates the Earth's angular momentum is changing or adjusting independent of external torques.Internal torques therefore must exist in the Earth independent of external torques that can be calculated directly from surface observations, without involving physical properties in the interior [14,15].The calculation is reliable at least on an order-of-magnitude sense, which allows us to quantitatively examine the in-

OPEN ACCESS IJG
C. PAN 66 ternal torques from the Earth's surface according to the three-finger rule of the right-handed system and the law of conservation of angular momentum.Order-of-magnitude is not likely to change with new observations or the inclusion of the Earth's rheology.
The first term I ω ⋅  in Equation ( 1) arises from inertia changes due to motion and mass redistribution, notably secular polar drift, plate motion, and mantle flow, in a rotating Earth.According to the three-finger rule of the right-handed system, it belongs to gyricity but is not a typical gyric torque as that can be directly represented by the cross product of ω with motion, and it will disappear if the Earth were rigid.Pan [14,15] names it the Coriolis torque and has quantitatively analyzed it in rather detail.Its magnitude is on the order of 10 23 ergs, which, if converted to force, is about 12 to 80 times greater than the extensional force exerted at the bottom of the lithosphere by upper mantle thermal convection as calculated by Richter [25].However, thermal flow, as a secular motion, can also induce such a gyric torque according to the three-finger rule of the right-handed system.The torque will then, for the sake of polar stability, guide the flow convectively to offset the perturbation to angular momentum from thermal buoyancy, which will be further explored below.
The second term I ω ⋅  in Equation (1) arises from rotational acceleration and is non-gyric [15].The term exists both in a rigid and a non-rigid Earth.Using the Earth's moments and products of inertia calculated by Pan [16] and observed rotational acceleration [26][27][28][29][30][31][32], the different components of the torque are calculated and listed in Table 1, where α is the angular distance between the rotation and instantaneous figure axes, and β is the azimuth angle of the rotation axis.The calculation procedure is in Pan [15].This torque has diverse components.Its periodic and irregular components are on the order of 10 26 ergs, while its secular components, tidal and non-tidal, are on the order of 10 23 ergs, the same as that of I ω ⋅  .Its secular non-tidal component is called the secular rotational torque, which has also been quantitatively analyzed in rather detail [15].
The third term h  in Equation ( 1) is the rate of change of relative angular momentum; the term will disappear if the Earth were rigid.Relative angular momentum has two components; one is due to motion and the other due to the residual products of inertia induced by the gyricity from rotation to motion [2,[14][15][16][17].Consequently, its rate of change consists also of two components, that due to acceleration and that due to the rate of change of the residual products of inertia.Pan [14] ignores the second component and calls the first the relative disturbing torque for its transient role at the onset of polar excitation.Calculation gives it a magnitude on the order of 10 23 ergs, but can also reach a much higher magnitude transiently.At α = 6.7 × 10 −7 rad.

OPEN ACCESS IJG
For instance, a h  from the impact of a giant asteroid can be as high as 10 28 ergs [17], if the asteroid is about 10 km in diameter as suggested by Alvarez and Asaro [33].On the other hand, the second component of h  behaves exactly like I ω ⋅  ; i.e., secular and gyric.Pan [14] discusses this secular gyric component under the rate of change of excitation function, and contends that even a secular h  arising from the motion of a single plate will affect polar stability.
The fourth term 1) is a typical gyric torque.It arises from the slight separation of the rotation axis from the invariant angular momentum axis, and is the only gyric torque that exists both in a rigid and a non-rigid Earth.It is the largest internal torque in the Earth, at least six orders greater than any other internal torques.Its magnitude depends on the Earth's axial nearsymmetry and triaxiality [15].If the Earth is biaxial, then its z-component becomes zero, but its x-and y-components are still greater than all the others.Gyricity in the Earth can thus be of large magnitude and not necessary all negligibly small as is conventionally believed.This torque is periodic and is primarily a reaction to the precession due to external attractions from the Moon and the Sun.However, nutation such as that due to mantle-core coupling [34][35][36][37][38] is observationally indistinguishable from precession [18,39].As we shall see below, free nutation arising from the Earth's unstable rotation will play a critical role in the Earth's equation of secular motion for the conservation of angular momentum of a perturbed Earth in the absence of external torques.
The fifth term h ω × in Equation ( 1) arises from mo- tion and mass redistribution in a rotating Earth; it hence does not exist in a rigid Earth.Any motion and mass redistribution in the Earth will induce such a gyric torque.Pan [15] calculates that a h ω × arising from the motion of a single plate is on the order of 10 22 ergs, but can reach 10 23 ergs if the accompanying flow underneath is also included.On the other hand, if h is a relative angular momentum induced by motion in a part of the Earth such as the flows in the outer core, h ω × then can represent the nutation arising from mantle-core coupling, and so on.
The above examination manifests that internal torques are of various magnitudes and are not all balanced by external torques, while the known reaction to external torques is order-of-magnitude greater than all the other internal torques.This indicates perturbation to the Earth's angular momentum can be independent of external torques.Among the five internal torques in Equation (1) ω ω × ⋅ is primarily a reaction to precession due to external torques, but is observationally indistinguishable from the free nutation in the absence of external torques, but free nutation, as we shall see below, has much smaller magnitude.Whereas, the periodic and irregular components of I ω ⋅  , listed in Table 1, are diverse and have their individual causes, internal or external, known or yet unknown.Nevertheless, the secular components of the internal torques are all on the order of 10 23 ergs, to which no external counterparts are found, as secular tidal torque is already counted for.This indicates that the Earth's angular momentum is conserved independent of external torques in a secular sense, which links secular global geodynamics to the Earth's non-rigid rotation via the equation of secular motion [15].

Generalization of the Earth's Equation of Secular Motion
Pan [15] demonstrates that polar motion, plate motion, and mantle flow are interrelated and their interaction can be traced back to the evolution of polar instability.A special equation of secular motion in the Earth for internal perturbations [15] was derived according to the law of conservation of angular momentum and the three-finger rule of the right-handed system.The equation can be further generalized to account for both internal and external perturbations.Let 0 I and 0 ω still be the constant inertia tensor and rotation velocity of an Earth in complete stable rotation of minimum energy configuration, and I and ω their variable counter- parts at perturbation [15].The relative angular momentum h now becomes due either to motion and mass redistribution within, mass depletion from, or external disturbance to the Earth; i.e., it represents either internal or external perturbation to the Earth's angular momentum.Then, at the very onset moment of perturbation, the instantaneous figure axis of I appears with h and shifts away from the major principal axis of 0 I to a new position at the excitation axis [2,17], while the rotation axis ω revolves away from its original position at 0 ω , which was aligned with the major principal axis of 0 I , around the instantaneous figure axis of I according to the three-finger rule of the right-handed system.However, while the instantaneous figure axis of I shifts its direction in space with the appearance of h , conserva- tion of angular momentum will keep the direction of the rotation axis ω still nearly fixed in space besides pre- cession and nutation around the invariant angular momentum axis.This is polar excitation [2,15,17].At this very moment, the law of conservation of angular momentum requires that a transfer of angular momentum must occur independent of external torques; i.e.,

OPEN ACCESS IJG
C. PAN 68 0 0 Equation ( 2) differs from Equation (3) in Pan [15] with the approximately equal sign  , which needs explanation.If the relative angular momentum h is due to mo- tion and mass redistribution within the Earth, the rightside of Equation (2) will equal exactly to the left-side.However, if h is from an external disturbance such as a planetsimal impact to the proto-Earth [40,41], the rightside can become greater than the left-side because of the transient impact and its adding mass to the Earth.On the other hand, if h is due to mass depletion such as the escape of the Moon from the Earth [42], the right-side will be smaller than the left-side.So Equation ( 2) implies that at the onset moment of perturbation, the Earth's total angular momentum may change via moving mass in or out the Earth; the internal torques thus may not be totally in balance at that time.Nevertheless, h must supply energy to the Earth for the transfer from a stable rotation of minimum energy configuration to an excited configuration of higher energy budget.Applying the differentiation operator d dt ω + × to Equation (2), and with the same physical reasoning, algebraic manipulations, and notations as those in Pan [15], the generalized equation of secular motion in the Earth is, ( ) Equation ( 3) approximates the linkage between the Earth's non-rigid rotation, secular global geodynamics as well as free nutation in the absence of external torques; it differs from the special equation of secular motion only for internal perturbations [15] by the approximately equal sign  as well as the secular gyric torque s h  .The equ- ation differs from Equation (1) by that it consists of only the secular internal torques and responding free nutation in the absence of external torques, while Equation ( 1) is the generalized equation of motion to represent the overall global geodynamics that consist of all internal torques, secular, periodic, irregular, and transient in the presence of external torques.One distinction in Equation ( 3) is the right-side gyric torque, which is not a reaction to precession but appears independent of external torques in response to the free nutation induced by secular motion.It will disappear only after the Earth reaches its complete stable rotation of minimum energy configuration at 0 0 s s

I I
ω ω ⋅ = ⋅ .The presence of this torque implies that secular motion will also induce free nutation in the absence of external torques for the sake of conservation of angular momentum.However, since Equation ( 3) is derived from Equation (2) in a general form before specified to represent only the secular motion, it cannot rule out that other apparently unbalanced internal torques of different magnitudes, such as those shown in Table 1, may also be able to induce free nutation in the absence of external torques for the sake of conservation of angular momentum.Free nutation shall be the reaction of a rotating Earth to internal perturbation for the sake of conservation of angular momentum, just like precession is a reaction to external attractions.Free nutation of different magnitudes is then buried in the much greater precession to become observationally indistinguishable, which needs further investigation.As exhibited in Equation ( 3), the secular internal torques are balanced by responding free nutation, the Earth's angular momentum is hence conserved in the absence of external torques in a secular sense.Now 0 0 I ω ⋅ in Equation ( 3) represents the ultimate angular momentum for the Earth to reach its complete stable rotation of minimum energy configuration, so it may not necessarily be identical to its original "fossil" at incipient polar excitation as that in Equation ( 2).This is because not only redistributing mass within, but also adding mass to, or depleting mass from the Earth may excite polar instability.All terms in Equation ( 3) belong to gyric torques except the second term at the leftside, s s I ω ⋅  .Gyricity therefore dominates secular mo- tion in the Earth, and is the major effort for a perturbed Earth to search for its stable rotation of minimum energy configuration.
The left-side secular internal torques in Equation (3) exhibit the driving mechanism of secular motion in the Earth.Among  h  in a rotating Earth, which enhances secular mo- tion accordingly.Concurrently, the rheological equatorial bulge will yield in response to secular motion and rotational deformation in a quadrupolar migration according to the three-finger rule of the right-handed system, which may appear as a wave-like secular propagation similar to the propagation of diurnal tides [18].The secular gyric torques and rotational deformation will disappear if the Earth were rigid.On the other hand, the right-side term in Equation ( 3) is a correction for polar stability [15]; stable rotation can be reached only when 0 0 s s

I I
ω ω ⋅ = ⋅ , and then the accompanying free nutation ceases.The magnitude of this right-side term in Equation ( 3) is equivalent to the excess energy supplied by angular momentum perturbation and gained at polar excitation.It is on the order of 10 23 ergs at the present time, negligibly small comparing to the Earth's total rotational energy of 2.16 × 10 36 ergs.
Pan [14,15] observes that since ω is always toward the north, according to the three-finger rule of the right-handed system, kinematically the secular gyric torque possesses two components, lateral and radial, and is quadrupolar; i.e., they are opposite in opposite hemispheres, depending on the direction of secular polar drift or lateral motion in the hemisphere.Surface lateral component will induce a subduction toward the center of the Earth, while subduction induces a lateral returning flow in the interior.A returning lateral flow will then induce an upward flow convectively, while an uprising plume induces a lateral motion on the surface.The migration of hot spots observed along the Hawaiian island chain can belong to such an effect; i.e., the hot plume is guided by gyricity.As exhibited in Figure 2 [43] and also in Gordon [44], the secular gyric torques can be responsible for the breakup of the lithosphere into quadrupolar plates more or less around 90˚-by-90˚ in dimension.Observation [14,15,45,46] shows secular polar drift is away from the Pacific Ocean at the north pole and towards the Pacific Ocean at the south pole, which thus induce westward subduction of the Pacific plate in the northern hemisphere and eastward subduction of the same plate in southern hemisphere.The quadrupolar effect of gyricity becomes a little chaotic near the equator and not confronted with a stronger continental plate, so the confronted oceanic plates may go down in both directions, such as that near the Hebridges, Fiji, and Tonga Islands [15].On the other hand, the only non-gyric secular torque s s I ω ⋅  is toward the north, which can be responsi- ble for northward plate motion, as those also displayed in Figure 2, the collision of the India plate with the Eurasia plate, the subduction of the Pacific plate underneath the Aleutian arcs and Alaska, as well as the subduction of the Cocos subplate underneath Middle America.No southward subduction is observed in the Earth.

Polar Stability
The polar or rotation stability implied in Equation ( 3) is different from that of a biaxial rigid body [2,3,17].A biaxial rigid body can keep a stable rotation even its rotation axis is misaligned with its major principal axis, since it lacks either energy dissipation or motion and mass redistribution, while the equatorial bulge acts as a stabilizer.However, for a deformable, energy-generating and dissipative Earth that allows motion and mass redistribution, a stable rotation can be reached only when its major principal, rotation, and instantaneous figure axes are all completely aligned with each other [2,15,17]; i.e., only when At that moment all secular internal torques will disappear or completely balance with each other, and then the accompanying free nutation ceases.This is the minimum energy configuration.Any separation of the axes is polar excitation that needs additional energy to attain, and the Earth has to dissipate the excess energy gained at polar excitation through subsequent polar motion, accordingly plate motion and mantle flow, for the "ass-and-carrot" pursuit of the axes for realignment [47].The direction of the rotation axis is nearly fixed in space besides precession and nutation, while the principal axes shift responding to mass redistribution [2,17].In a deformable, energy-generating and dissipative Earth that allows motion and mass redistribution, a realignment of the axes will, therefore, only be able to revert back to their original direction in space but not original position in the Earth.This is an unstable rotation.In such an Earth, the rheological equatorial bulge is no longer as much a stabilizer as that in a biaxial rigid body, because it will migrate with the unstable rotation via quadrupolar deformation [18].The incipient solution of the Liouville equation [14] confirms the rotation of a perturbed non-rigid Earth is unstable, for it gives an exponentially increasing secular polar drift as well as an exponentially damping wobble, as an unstable rotation should.
Polar instability consists of incipient polar excitation that gains energy, and subsequent polar motion that dissipates energy.Polar excitation, as we shall see below, can either be episodic or continuous, but polar motion is always continuous until a stable rotation is reached.The polar wandering mechanism discussed by Goldreich and Toomre [48] is a simplified account of polar instability, because it does not differentiate polar excitation from polar motion.The incipient solution of the Liouville equation [14,16] gives a static term for incipient polar excitation, as well as a secular polar drift term and a wobble term for subsequent polar motion.However, there is no clear distinction between the perturbation that excites polar instability and the subsequent adjustment that damps polar instability; both belong to mass redistribution in the Earth.One criterion for their physical division is that the mass redistribution that induces motion and products of inertia excites polar instability, and that reduces motion and products of inertia damps polar instability.Equation (3) expresses the latter.Nevertheless, if a motion, such as seasonal fluctuations of atmosphere or the atmospheric and oceanic excitation of the Earth's wobbles [49,50], that induces little mass redistribution in the solid Earth, then there will be little change in the position of the principal axes in the solid Earth.However, according to the three-finger rule of the right-handed system, the relative angular momentum induced by the motion will force the rotation axis to revolve away from the major principal axis around an instantaneous figure axis that has been shifted to a new position at the excitation axis [2,17,18], about which the inertia tensor is no longer diagonal due to its axial near-symmetry [2,14,17].Polar excitation by motion alone thus can also induce products of inertia in the solid Earth.In such a case, polar excitation is due to continuous motion as well as the products of inertia arising from the near-symmetry of the instantaneous figure axis at its new position, the excitation axis.The instantaneous figure axis thus will shift responding to either motion or mass redistribution, while the principal axes shift responding only to mass redistribution.

Perturbation and Energy Supply
We have so far demonstrated that secular internal torques are not the angular momentum perturbation that will excite polar instability, but represent the effort of an already perturbed Earth in search for a stable rotation through self-adjustment according to the three-finger rule of right-handed system for the sake of conservation of angular momentum.They thus dissipate energy and will appear in the Earth independent of external torques only after unstable rotation of higher energy configuration has already been excited.Pan [14,16] predicts the Chandler wobble possesses multiple frequency splits and is slow damping or quasi-permanent, as is confirmed by observation [1,2].Polar motion hence does not need a simultaneous polar excitation to maintain.However, a quasipermanent polar motion will still be damped out eventually for a stable rotation, if there is no repeated angular momentum perturbation or episodic polar excitation to suspend it.Ahrens [40] and Wood [41] contend that the proto-Earth was formed through rapid and violent impact accretion of large planetsimals, which could repeatedly excite large polar motion or even a change of the Earth's obliquity [51].Pines and Shaham [42] suggest that some violent events in geological history, such as the escape of the Moon from the Earth or the impact of a giant meteorite, might have excited polar instability.Pan [17] demonstrates that the magnitude of polar excitation owing to meteorite or asteroid impact depends on the products of inertia that can be induced by the impact.For instance, if the impact asteroid was around 10 km in diameter [33], its polar excitation was about 1.15 × 10 −5 , slightly greater than the present Chandler wobble.In the above external perturbations, the energy supply is gravitational and mechanical.On the other hand, geological observation has long established that geological history was intersected by mountain buildings or diastrophisms between geological epochs.Pan [14,[52][53][54][55][56] calculates that the diastrophisms are able to excite polar instability episodically on the order of 10 7 years, initiated by explosive release of excess thermal energy accumulated in the interior through gigantic failures along some susceptible zones in

OPEN ACCESS IJG
C. PAN 71 the lithosphere.The energy supply for such internal perturbations is thermal.The excess thermal energy that can be accumulated in the Earth in 10 7 years is on the order of 10 35 ergs [52,53].Pan [14] has quantitatively analyzed ten cases of such polar excitation, which can reach a magnitude as large as 10 −2 , four orders greater than the present Chandler wobble.All above perturbations, internal or external, are episodic.Spada et al. [57] consider the excitation of polar wandering by subduction, while Steinberger and O'Connell [58] discuss changes of the Earth's rotation axis owing to advection of mantle density heterogeneities.Gross [49] and Gross et al. [50] investigate the atmospheric and oceanic excitation of the Earth's wobbles.For such cases, perturbation and consequent polar excitation are continuous, while energy supply is internal and thermal or gravitational.Polar excitation will occur simultaneously with polar motion until thermal convection or atmospheric and oceanic processes stop, but damping of polar motion will continue until a complete stable rotation is reached.
The Earth's dissipation of excess energy for a stable rotation is independent of the loss of kinetic rotational energy through tidal friction.This is because, as we have noted above, secular rotational torque comes from nontidal secular rotational acceleration.The energy dissipation through the Earth's self-adjustment will stop as soon as the Earth has reached its complete stable rotation, but the loss of rotational energy through tidal friction will continue as long as external attractions are there.It is interesting that the magnitude of lunar tidal torque, according to Stacey [59], is 4.4 × 10 23 ergs, on the same order of magnitude as that of secular internal torques.If the viscosity of the low-velocity layer is sufficiently low or the mantle is non-Newtonian [15], it cannot rule out that tidal torque might cause a secular lag of the lithosphere to induce apparent quadrupolar gyric motion of the tectonic plates over the low-velocity layer according to the three-finger rule of the right-handed system.If this were the case, then the driving mechanism of plate motion would involve also the exchange of angular momentum with the Moon, of which the energy supply is external and gravitational.

Convection
As has been pointed out above, secular gyric torques are continuous to the interior convectively according to the three-finger rule of the right-handed system.Pan [14,15] has preliminarily studied the (secular) Coriolis and rotational torques in the mantle and named the internal-torque-induced flow transvection, which is a form of convection induced by gyricity and is thus different from thermal convection, though both belong to secular motion in the mantle.Equation (3) and related discussions define transvection.It is a reaction of a deformable and dissipative Earth to angular momentum perturbation; i.e., it is the Earth's effort to damp its products of inertia for a stable rotation to conserve angular momentum.Transvection is, therefore, a mass redistribution that occurs only after polar instability has been excited.On the other hand, thermal convection is an action of thermal buoyancy in the interior; i.e., it is a mass redistribution that induces products of inertia.In other words, thermal convection belongs to the angular momentum perturbation that supplies energy to excite polar instability, while transvection is the subsequent adjustment through gyricity that dissipates energy for polar stability.Transvection hence will guide thermal flow according to the quadrupolar rule of gyricity to offset the perturbation from thermal buoyancy for the sake of angular momentum conservation.
Mantle convection has long been believed to be the driving mechanism for plate motion.Convection models are generally divided into four categories: Thermal plumes [60,61], upper mantle convection [25,[62][63][64][65][66], layered convection [67][68][69][70][71][72], and whole mantle convection [68,70,[73][74][75][76].A common problem for the various models is the physical properties in the interior; it is not yet certain whether the mantle is Newtonian or non-Newtonian.Newtonian mantle prevails because we know more of Newtonian surface flows, and it is also easier to handle mathematically.However, observations and experiments also suggest that the rheological behavior of the upper mantle may be non-Newtonian [62,63,[77][78][79][80].Under the order-of-magnitude quantitative constraint from surface secular motion observations, a preliminary study [15] shows a non-Newtonian mantle is more compatible with secular internal torques.If the mantle is Newtonian, then the viscosity of the upper mantle used in that study was either too high or non-secular; otherwise, the viscous torque adopted for the upper mantle was inadequate [15].Since the physical properties in the interior are deduced indirectly on the surface and are non-secular, how will they behave in the true physical conditions in the mantle in a geological time-frame is a good question.Calculation based on surface observations [15] shows that the driving effect of upper mantle thermal convection on plate motion is one order smaller than that of secular internal torques; only whole mantle convection is able to induce a torque that is on the order of 10 23 ergs.This implies if thermal convection is to drive plate motion, it has to be whole mantle convection.
There is another concern for thermal convection, energy source or power supply.Heat from the core and radioactivity in the interior are believed to be the energy reservoirs.Nevertheless, the energy sources in the interior are indirectly estimated with assumptions [59], not like the rotational energy and secular internal torques that

OPEN ACCESS IJG
are directly observed on the Earth's surface.The various thermal convection models are therefore more kinematical than dynamical.As a perturbation to angular momentum, thermal convection is a mass redistribution that will excite polar instability.However, as a secular motion, it has also to follow the three-finger rule of the right-handed system just like any other motions in a rotating Earth.This is to say, if thermal convection is to drive mantle flow, it will cause an inertia change to induce the secular Coriolis torque s s I ω ⋅  , which consequently will follow the quadrupolar rule of gyricity to guide the flow to offset the perturbing effect from thermal buoyancy for polar stability.So no matter how thermal convection behaves, whether it is to drive plate motion, gyricity will always appear to dominate secular motion in the interior, in order to search for a stable rotation for the conservation of angular momentum.Thermal convection can be either integrated with the Earth's rotation as a perturbation through relative angular momentum s h [17] or incorporated into s s I ω ⋅  as a mantle flow, and then secular motion in the mantle will become much more complex.Thermal convection in a rotating Earth is not the same as that in an inertial Earth or flat Earth.
Secular internal torques are depth-dependent; i.e., they decrease with depth.In a non-Newtonian mantle, they are likely balanced by the viscoplastic torque about the bottom of the upper asthenosphere [15], which gives transvection an appearance analogous to the upper-mantle or layered convection [72]; the 660-km discontinuity is thus not a barrier but a likely lower boundary for transvection.It is, therefore, not unreasonable to perceive that hot plumes may originate in the asthenosphere.The interaction between mantle thermal convection and the Earth's rotation is worth an investigation.The fundamental physics that governs the flows in the Earth's interior should not be much different from that for the more mobile layers, the atmosphere and oceans.Mantle dynamics cannot be independent of the Earth's perpetual rotation.

More Discussions
Movements in the more rigid layer of the Earth, the lithosphere, are not in flows but in quadrupolar plates more or less 90˚-by-90˚ in dimension moving or subducting against each other according to the three-finger rule of the right-handed system, reflecting the quadrupolar nature of gyricity [15].Secular internal torques therefore furnish dynamics to plate tectonics.The stress distribution in the lithosphere is expected to concentrate more along the confronting boundaries.This helps to explain the differences between intraplate and interplate seismicities; why the former is less frequent than the lat-ter.The temporal variation of global seismicity can be correlated with secular internal torques, and such a correlation may help to learn about the global trend of seismic activity.On the other hand, neither plate motion nor flows in the different layers in the Earth will have much effect on the Earth's principal gravity field, which is governed by Newton's gravitational law [81].Nevertheless, secular motion can affect the part of gravity that is due to rotation; i.e., it is responsible for the adjustment of the slight imbalances of the Earth's gravity caused by its unstable rotation.
The energy required for secular motion in the Earth, as has been calculated above, is on the order of 10 23 ergs, negligible in comparison to the Earth's total rotational energy of 2.16 × 10 36 ergs, which is the largest energy reservoir that can be directly observed on the Earth's surface.This implies a slight fluctuation of the rotational energy, or an angular momentum perturbation that can provide to the Earth an unnoticeable amount of energy on the order of 10 23 ergs, is sufficient to supply power for any secular dynamical processes in the Earth and maintain the presently observed polar motion, plate motion, and mantle flow.Similar to shaping a pottery on a wobbling turning table, the deformation and adjustment for the Earth to keep its near-spherical shape after gravitation reflect more of the effect of rotation rather than that of a dynamical process in the interior.We may have long overlooked this fact because we are too much used to our instantaneous life of twenty-four hours for a day and take it for granted.Dynamical processes in the interior can dominate mass redistribution only when the Earth is in a stable rotation of minimum energy configuration.Yet, whenever a dynamical process in the interior perturbs the Earth's angular momentum, polar instability will be excited immediately.Secular internal torques will then appear accordingly, and rotation, subsequently gyricity, will again dominate secular motion in the Earth.

Conclusion
The physical Earth is a multilayered, deformable, energygenerating, dissipative, and perpetually rotating heavenly body orbiting in space according to the law of conservation of angular momentum, while motion in a rotating Earth follows the three-finger rule of the right-handed system.The Earth's rotation dynamics is thus integrated with global geodynamics, and the generalized Eulerian equation of motion or the Liouville equation represents global geodynamics in general.An examination of the Liouville equation via observation of the Earth's rotation and surface motions manifests that the Earth's angular momentum is perturbed and adjusting independent of external torques.Secular internal torques represent the Earth's effort to conserve its angular momentum for a

OPEN ACCESS IJG
stable rotation of minimum energy configuration, which appear in the Earth independent of external torques only after polar instability has been excited.Secular internal torques calculated from surface observations are all on the order of 10 23 ergs, which gives a quantitative orderof-magnitude constraint to secular geodynamics in the interior.The Earth can be excited from a stable rotation of minimum energy configuration to an unstable rotation of higher energy budget by either internal or external perturbations.For the sake of conservation of angular momentum, internal torques in a perturbed Earth are balanced by free nutation in the absence of external torques.The Earth's generalized equation of secular motion is an approximation of the global geodynamics that links polar motion, plate motion and mantle flow to free nutation.Gyricity dominates secular motion in the Earth, which, together with secular rotational torque, may well be the primary driving mechanism for plate motion and mantle flow.A deformable, energy-generating and dissipative Earth can reach a stable rotation of minimum energy configuration only when its major principal, rotation, and instantaneous figure axes are all completely aligned with each other; separation of the axes excites polar instability.Polar instability consists of incipient polar excitation that gains energy to commence polar instability, and subsequent polar motion, plate motion and mantle flow that dissipate energy for polar stability.The energy gained at polar excitation is supplied by angular momentum perturbation, which can be either episodic or continuous, internal and thermal or external and gravitational.According to the three-finger rule of the right-handed system, secular internal torques will induce transvection, a passive form of convection arising from gyricity that is different from active thermal convection.Thermal buoyancy supplies energy to excite polar instability, while transvection is the gyric effort of a rotating Earth to dissipate energy for polar stability.However, as a secular motion, thermal buoyancy is also subject to gyricity, which then guides the flow to follow the three-finger rule of the right-handed system for angular momentum conservation.The rotational energy of 2.16 × 10 36 ergs is the largest energy reservoir that can be directly observed on the Earth's surface; in comparison the energy supply for secular motion in the Earth is insignificant.An unnoticeable fluctuation of the Earth's rotational energy on the order of 10 23 ergs is sufficient to supply for any secular geodynamical processes in the Earth.No geodynamics in a deformable, energy-generating and dissipative Earth orbiting in space that allows motion and mass redistribution can be independent of its perpetual rotation, as long as the Liouville equation is its generalized equation of motion that follows the law of conservation of angular momentum and the three-finger rule of the right-handed system.This paper presents an introduction to the unifi-cation of global geodynamics with the Earth's rotation dynamics through the Liouville equation under fundamental physical laws.

Figure 1 .
Figure 1.The coordinate systems of the Earth prior to and at polar excitation.

Table 1 . Different components of I
ω ⋅  calculated from direct surface observations.