Abstract

The recently verified superrotation of Earth’s inner core is examined and a new model is presented which is based on the tidal despinning of the mantle and the viscosity of the outer core. The model also takes into account other damping mechanisms arising from the inner core superrotation such as magnetic and gravitational coupling as well as contribution from eddy viscosity in the outer core. The effective viscosity obtained in this model confirms a previously well constrained value of about 103 Pa s. In addition, the model shows that the currently measured superrotation of the inner core must be almost exactly equal to its asymptotic or steady-state value. The effect of extraterrestrial impacts is also investigated, and it is shown that perturbations due to such impacts can only persist over a short geological time.

1. Introduction

In 1996, Song and Richards analyzed two types of seismic waves traveling through Earth’s solid and fluid cores, respectively [1]. They discovered a differential travel time between the two waves, which increased systematically by about 0.3 s from the year 1967 to the year 1995. This temporal change was later attributed to a shift of the lateral velocity gradient in the inner core caused by the inner core rotation [2, 3]. Subsequent studies provided further support for the faster rotation of Earth’s solid inner core compared to the rest of the planet [47].

Because of the anisotropy due to wood-like grain structure of the crystalline iron in the inner core, the speed of seismic waves is different in different directions [8]. If the inner core rotates faster than the rest of Earth, seismic waves that are generated at the same place on Earth but decades apart will have different travel times passing through the inner core and detected diametrically across the globe. Accurate recent investigations, using high-quality waveform doublets generated in the South Sandwich Islands region and detected in and near Alaska, have revealed definite temporal changes over a period of up to 35 years [9]. These experiments confirm that Earth’s inner core rotates faster than the mantle and crust (hereafter simply referred to as mantle), with a superrotation of about 0.27 to 0.53 deg/yr.

Estimates of the inner core differential rotation rates range from low values of 0.2 to 0.3 deg/yr [2] to high values of 3 deg/yr [10]. Deuss [11] presents a compilation of studies on inner core differential rotation rates as a function of publication year that suggests more recent data supports lower inner core differential rotation rates; however, the data is quite scattered. The author cites a recent study [12] that reports no inner core rotation and another study [13] that shows subrotation but the majority of the studies report inner core superrotation with respect to the mantle.

Several Earth-based mechanisms have been suggested for Earth’s inner-core differential rotation. Gubbins [14] attributes the phenomenon to the existence of a large toroidal magnetic field inside Earth’s core which, together with the dipole components, produces a magnetic torque on the inner core that tends to accelerate it relative to the mantle. He then assumes that some equilibrium is reached so that the inner core rotates with constant angular velocity and experiences zero net torque. Glatzmaier and Roberts [15] numerically solve the self-consistent magnetohydrodynamic equations that describe thermal convection and magnetic field generation in a rapidly rotating spherical fluid shell with a solid conducting inner core. Their solution, which serves as an analog for the geodynamo, shows that viscous and magnetic coupling of the outer core with the inner core and the mantle causes time-dependent variations in their respective rotation rates; the inner core usually rotates faster than the mantle. Buffett and Glatzmaier [16] allowed gravitational coupling between the inner core and the mantle by incorporating viscous deformation of the inner core into their numerical simulations of the geodynamo. They discovered that differential rotation between the inner core and the mantle is permitted by allowing the inner core to deform. Numerical calculations by Aurnou et al. [17] showed excess temperature inside the tangent cylinder surrounding the inner core which generates a prograde thermal wind and a strong azimuthal magnetic field inside the tangent cylinder. They conclude that the electromagnetic torque on the inner core resulting from induced azimuthal magnetic fields and the ambient poloidal field equilibrate when certain conditions are met. Dumberry [18] studied the steady and time-dependent rates of inner core rotation based on angular momentum balance between the inner core, fluid core, and mantle. He concluded that the rotational rate of an oscillating inner core is constrained by the changes in mantle rotation induced by gravitational coupling.

Tkalčić and others [19] analyzed earthquake doubles to construct a model for inner core differential rotation rates of 0.25 to 0.48 deg/yr with decadal fluctuations around the mean of 1 deg/yr. The authors suggest that these decadal fluctuations can account for discrepancies between previous core rotation models and agree with recent geodynamo simulations. A three-dimensional model by Livermore and others [20] suggests that axial electromagnetic torque is the dominant influence for inner core differential rotation and that decadal variations of the magnetic field may drive the quasioscillatory nature of the inner core differential rotation.

Magnetic coupling between the inner and outer cores seems to play a role not only in the superrotation of the inner core but also in the generation of the Earth’s magnetic field. Su et al. [10] reported an anomalous variation in the inner core orientation that temporally coincided with the geomagnetic “jerk” (a sudden change in the strength of Earth’s magnetic field of 1969-1970). This suggests a correlation between Earth’s magnetic field and the inner core superrotation. Glatzmaier and Roberts [21] have suggested that the inner core rotates in response to the magnetic torque and the viscous torque to which it is subjected, with + = 0. The magnetic torque drags the inner core eastward and the viscous torque acts westward. Nevertheless, they state that even though the fluid viscosity in their model is several orders of magnitude greater than is likely for real Earth, the viscous torque on the inner core has little effect.

Although each of these theories provides a reasonable explanation for the differential rotation of the inner core, there is no experimental evidence supporting one over the others, and they are all based on processes that are assumed to be taking place inside Earth. Furthermore, about a decade before the most recent data of 0.27 to 0.53 deg/yr were published, calculations based on magnetic coupling between the inner and outer cores suggested superrotations that were about an order of magnitude higher than these values [10, 21].

Perhaps the most significant and tangible external factor responsible for the inner core superrotation is the tidal forces [22, 23]. Dissipation of tidal energy in oceans and transfer of angular momentum between Earth, Moon, and the Sun results in torques on the mantle causing it to dispin [24], which gradually increases the length of the day continuously. This despinning of the mantle leaves the inner core with a small excess eastward rotational velocity relative to the mantle. Su et al. [10] used the known tidal increase in the length of day of approximately 2 ms per century to extrapolate backward and conclude that the inner core was rotating with the same period as the mantle about 105 years ago. The recent confirmation of superrotation of Earth’s inner core by Zhang and others [9] prompted us to examine this phenomenon in the context of tidal effects from a phenomenological point of view.

In addition to the oceanic tidal effects, yet another external factor could contribute to the differential rotation of the Earth’s inner core which has not been addressed in the literature, Earth impacts. Since the formation of our Solar System some 4.6 billion years ago, collisions and impacts have played a fundamental role in establishing its characteristics, ranging from the accretion of planetesimals and the early formation of planets [25, 26] to the recent series of impacts on Jupiter by the fragments of Shoemaker-Levy 9 comet in July 1994 [27]. These collisions and impacts have affected the dynamics of various components of the Solar System. For example, most planets have obliquity or axial tilt with respect to their orbital planes about the Sun. Earth has an obliquity of about 23.5° while Uranus, the third largest planet in the Solar System, has an obliquity of about 97°. In other words, Uranus is tilted on its side so that its rotation axis is nearly in its orbital plane about the Sun. Yet, giant planets are believed to form with nearly zero obliquity [28]. The axial tilts of planets are believed to have been caused by major impacts [29, 30].

The above discussion is the motivation for examining a second question in this paper. Because Earth’s solid inner core rotates inside the fluid outer core, is it possible for bolide impacts to alter the angular velocity of Earth’s mantle relative to the inner core, resulting in a superrotation, a subrotation or a transition from one to the other? If so, is it possible for these impacts to result in differential rotations that are comparable to the experimentally observed superrotation values, and how long would it take for such perturbations to damp out?

The model presented here could equally be applied for conditions that would generate subrotation or no differential rotation depending on the size and angle of bolide impact. Examination of the probability of direction of differential rotation is beyond the scope of this paper. The focus here is to examine the potential for bolide impact contribution if superrotation of the inner core is present.

2. Tidal Despinning of the Mantle and Viscosity of the Outer Core

As stated earlier, dissipation of tidal energy in the oceans results in gradual despinning of the mantle, which leaves the inner core with a small excess relative eastward differential rotation relative to the mantle called superrotation [24]. This differential rotation is communicated between the mantle and the inner core though a Couette flow in the fluid outer core which tends to damp the superrotation.

Consider two concentric solid spheres separated by a fluid layer between them with a coefficient of viscosity . Let the radius of the inner sphere be and the inside radius of the outer sphere be . Then, if the inner sphere rotates with an angular velocity with respect to the outer sphere, the damping torque on it due to the fluid viscosity is given by [31]Although Dai and others [32] have presented seismic reflection data suggesting significant topography on the inner core boundary, this model can be applied to Earth’s inner core and mantle, with the fluid between them being the outer core. Here the effect of the topography is absorbed in the estimate of effective viscosity.

Using (1), the rotational equation of motion of the inner core,becomeswhere is the rotational inertia of the solid inner sphere. The left hand side in this equation is the net torque on the inner sphere. The first term represents the torque due to the tidal effects, which tends to increase the eastward rotational velocity of the inner core relative to the mantle. is the constant rate at which the rotational speed of the mantle slows down due to the tidal effects. It has a value of 2.3 milliseconds per day per century [33] or 3.45 × 10−5 deg/yr2. The second term on the left hand side of (3) is the viscous damping torque.

Writing in terms of the inner core density and the inner core radius , reduces (3) toFurthermore, defining the characteristic time byreduces (4) toIntegration of this equation using the initial condition givesFrom this equation, we see that is in fact the relaxation time for the process. With the value of obtained from (5), this equation gives the superrotation of the inner core at any time after the dissipation of tidal energy began, about 4 × 109 years ago, since oceans existed as early as the relatively stable Earth. The value of , as mentioned earlier, is 3.45 × 10−5 deg/yr2.

Earth’s inner core has a radius of 1.220 × 106 m and a density of 13000 kg m−3. The liquid outer core has a radius of 3.473 × 106 m [34]. The viscosity of the outer core, however, is highly uncertain depending on its method of determination. In fact, its estimated values from various sources span over 14 orders of magnitude from 10−3 to 1011 Pa s [35]. For example, de Wijs et al. [36] report a value of 1.5 × 10−2 Pa s with an uncertainty of a factor of three through dynamical first principles simulations of liquid iron. Similarly, Rutter et al. [37] report a value of 1.6 × 10−2 Pa s using experimental high pressure study of liquid Fe-S system. On the other hand, using amplitude of forced nutation, Molodenskiy [38] finds a value of about 1 × 106 Pa s.

Although the viscosity of the Earth’s outer core is notoriously uncertain, it is plausibly considered to be bounded by 102 Pa s ≤ ≤ 1011 Pa s [32, 34]. Estimates of the fluid outer cores viscosity show variation from 102 Pa s at the top of the outer core to 1011 Pa s at the bottom near the inner core boundary using Arrhenius extrapolation of pressure dependencies for laboratory measurements on liquid iron [39, 40]. A viscosity value of 1.22 × 1011 Pa s near the solid inner core was found by Smylie [41] using Ekman layer theory to estimate viscous drag forces from Coriolis splitting of the two equatorial translational models of oscillation of the solid inner core. Palmer and Smylie [42] found a viscosity value of 615 Pa s near the top of the outer core from the free decay of free core nutations. Smylie et al. [43] present a viscosity estimate of 2.371 × 103 Pa s for the top of the outer core which relies on the decay of free core nutations method and an estimate for the bottom of the outer core of 1.247 × 1011 Pa s from the method of Smylie [41] and accounting for pressure dependence of the activation volume. The authors found the viscosity profile fits a nearly log-linear trend across the outer core but values still vary by 9 orders of magnitude. Smylie [44] presents this log-linear viscosity profile for the outer core based on the Arrhenius description of temperature and pressure of the viscosity by Brazhkin [39] and also reports a mean viscosity value of 3124 Pa s at the top of the outer core.

Therefore, instead of using a highly variable value of the viscosity to calculate the superrotation of the inner core, we solve the inverse problem. We use the measured value of the superrotation to find the viscosity of the outer core. To do so, we first solve the transcendental equation (7) numerically for the relaxation time , using = 0.4 deg/yr (mean value of the measured superrotation), = 3.45 × 10−5 deg/yr2, and = 4 × 109 yr (age of oceans). This gives a value = 11600 yr (or 3.66 × 1011 s). Then, we use (5) and solve for the viscosity of the outer core . We find = 3.37 × 103 Pa s. This number is in incredible agreement with the value obtained by Bills [24] and by Smylie [44]. Using a well described different model, Bills [24] reports an outer core viscosity of the order of 103 Pa s, assuming a value of superrotation of about 1 deg/yr. It is interesting to note that if we use = 1 deg/yr (instead of 0.4 deg/yr) in our calculations, we find = 1.35 × 103 Pa s. Our calculated viscosity does not agree with the study by Su et al. [10] that implied a viscosity of 10−4 Pa s by ignoring electromagnetic forces and assuming differential rotation of 3 deg/yr.

Admittedly, our values of the viscosity are calculated at the inner core boundary as our theoretical model considers rotation of the inner core. However, the extremely high viscosity of 1011 Pa s near the inner core strongly couples the lower zone of the outer core to the inner core. Thus, the differential rotation effectively takes place further out in the outer core where the viscosity is considerably lower. Since we are solving the inverse problem, our calculated viscosity corresponds to this region where the differential rotation effectively takes place. Ironically, due to the functional form of (5), the values of the viscosity calculated for a superrotation of 0.4 deg/yr with = 12000–13000 kg/m3 and with ranging from the inner core radius all the way to 98% of the outer core radius vary in a narrow range of 103-104 Pa s.

Let us now return to the notorious uncertainty in the reported values of outer core viscosity. Generally these numbers fall into two distinct categories. The high values are based on seismologic, geodetic, and geomagnetic observations of Earth [32, 34, 41], whereas the low values are based on theory and laboratory investigations of liquid metals [24, 45, 46]. It has been suggested that this dichotomy is possibly due to contribution of eddy viscosity caused by fluid motion, whereas the liquid metal investigations only account for intrinsic or molecular viscosity [24]. In addition, it is possible that electromagnetic and gravitational couplings also play roles. Since all these effects arise from the differential rotation of the inner core; therefore, they should each be a function of . Furthermore, according to Lenz’ Law and LeChatelier Principle, they should all act in such a way to reduce or eliminate the cause, which is the inner core superrotation.

Let the combined torque resulting from all these damping effects, which for the sake of simplicity from now on we refer to as core coupling, be = . Expanding this function in a Taylor series, we haveBecause is small, we can neglect the second- and higher-order terms in the expansion. In addition, since we must have = 0, we obtainAdding this retarding torque to the left hand side of (3), we obtainRepeating the calculations as before, instead of (6) we getwhere is a constant with the dimension of time. Integration of this equation with the initial condition = 0 giveswhere is defined bywhich is the new relaxation time, taking into account both the viscosity of the outer core as well as all other damping effects which we called core coupling.

Equation (12) has exactly the same functional form as (7) except that is replaced by . Therefore, if we solve this equation for using the known value of , we find the same value for that we found for before; that is, = 11600 yr. Thus, the relaxation time of 11600 yr corresponds not only to the viscosity of the outer core, but also to all other damping effects as well. Using this value in (5) and solving for ; therefore, gives an effective value for the outer core viscosity that could be greater than the actual value. It is interesting to note that since the relaxation time is 11600 yr, the value of today is so close to its saturation value that it would be impossible to measure the difference.

3. Extraterrestrial Impacts

Because geomagnetism and its reversals are associated with the differential rotation of Earth’s inner core, therefore, any changes in the Earth’s magnetic field can be attributed to a change in the inner core’s superrotation. Muller and Morris [47] have suggested that the impact of a large extraterrestrial object on Earth can produce a geomagnetic reversal through a mechanism involving a sequence of events. If so, it would be possible for such an impact to alter the superrotation of the inner core. In what follows, we study the dynamics of a bolide impact and calculate the change in the steady-state value of of the differential rotation of Earth’s inner core resulting from such impacts. We then investigate the time scale over which such a perturbation in the inner core superrotation would damp out.

Consider an asteroid or comet of mass and velocity v relative to the center of Earth just before impact, as shown in Figure 1. For simplicity, we assume that the velocity vector of the asteroid or comet is in the equatorial plane of Earth. The center of mass of the bolide-Earth system is located at a distance offrom Earth’s center, where is the mass of Earth and is its radius. For a typical bolide, and therefore is very small. For example, for a bolide of radius  km and a density comparable to that of Earth, we haveThen, with = 6.378 × 106 m, we find < 0.2 m, which is less than 20 cm from Earth’s center. Therefore, we can safely take the center of mass of the bolide-Earth system at the moment of impact to be simply Earth’s center.

The time rate of change of the total angular momentum of a system of particles is equal to the net external torque on the system, all measured with respect to some inertial coordinate system:However, if (16) is written relative to the center of mass of the system, then it is valid even if the center of mass is accelerating [48].

For the bolide-Earth system, the net external torque about Earth’s center is zero and the angular momentum of the system about that point is conserved during the impact. Therefore, we havewhere is the speed (magnitude of the velocity vector) of the bolide and is the angle between the bolide’s velocity vector and the tangent to Earth in the direction of Earth’s rotation, as shown in Figure 1. The first term on the left hand side is the angular momentum of Earth before the impact, in which is the rotational inertia of Earth and Ω is its angular velocity. The second term is the angular momentum of the bolide before the impact. The term on the right is the angular momentum of the combination after the impact, in which and are very small compared to and Ω, respectively.

Expanding the right hand side of (17), neglecting the second-order infinitesimal term   Ω, and rearranging the remaining terms, we findwhich may also be written asEquation (19) simply states that the change of the angular momentum of Earth is that imparted to it by the bolide.

Dividing both sides of (18) Ω, we obtainNow, because the bolide lodges itself near Earth’s surface. Furthermore, assuming Earth to be a solid sphere of uniform density for the moment (we address this point later on), we have . Substituting in (20), we obtainFinally, usingwhere and are the mean densities of the bolide and Earth, respectively, we obtainEquation (23) measures the fractional change in Earth’s angular velocity as a result of a bolide impact. This equation assumes that the velocity vector of the bolide is in the equatorial plane of Earth, and that Earth is a solid sphere of uniform density. However, since the outer core is fluid, it does not immediately transfer the angular momentum resulting from the bolide impact to the solid inner core. Therefore, we need to only consider the rotation of the mantle in our calculations. The moment of inertia of the mantle, , is the difference between the moment of inertia of the entire Earth and that of the inner and outer cores:where is the radius of the outer core, is the combined mean density of the inner and outer cores, and is the mean density of the entire Earth. If (24) is written aswhere(23) becomeswhere now and refer, respectively, to the angular velocity and the change in angular velocity of the mantle.

Because Earth’s mass and mean radius are, respectively, = 5.976 × 1024 kg and = 6.371 × 106 m, we calculate a mean density of = 5517 kg m−3 for Earth. Furthermore, the radius of the outer core and the combined mean density of the inner and outer cores are, respectively, = 3.473 × 106 m and = 11770 kg m−3 [34]. These data, together with (26), give = 0.8973. Thus, the correction to the fractional change in the angular velocity of the mantle as compared to entire Earth is only about 10.3%.

The change in the angular velocity of the mantle as a result of asteroid or comet impact can be positive, negative, or zero depending on the parameters involved in (27). There are tens of thousands of asteroids in the Solar System and about 220 have diameters greater than 100 km. Ceres, the largest known asteroid, has a diameter of about 950 km [49]. The orbit of some of these asteroids is such that they could potentially be perturbed into a collision course with Earth.

An object attracted by Earth’s gravitational field from large distances will impact Earth with a minimum speed of about 11 km s−1. A more likely collision, however, is when Earth intercepts an object that is being attracted by the Sun which has a much more powerful gravitational attraction and the impact speed will be much higher. In fact, it has been estimated that if the bolide responsible for the Chicxulub crater was a comet, it must have struck Earth with a speed as high as 70 or 80 km s−1 [50].

To obtain a numerical estimate for the change of the angular velocity of the mantle as a result of an asteroid or comet impact, we consider a bolide of radius 15 km and a density comparable to that of Earth, impacting Earth with a speed of 40 km s−1. The angle of impact can be anywhere from 0 to 180° as shown in Figure 1. For between 0 and 90°, the impact speeds up the rotation of the mantle relative to the inner core, and for between 90° and 180° the impact slows it down. Let us consider the special cases of near glancing angles = 0 and = 180°. Because Earth’s angular velocity isfrom (27), with = 0, we find = 2.257 × 10−10 rad/s = 0.41 deg/yr. On the other hand, if the impact is in the direction = 180°, we find a change in the mantle’s angular velocity of  deg/yr. Because the change of the differential rotation of the inner core relative to the mantle as a result of an impact is , the impact in the former case decreases the superrotation of the inner core by 0.41 deg/yr and that in the latter case increases it by 0.42 deg/yr. A strong enough impact in the direction can impart a change in that exceeds the value of , resulting in a subrotation of the inner core.

Figure 2 shows plots of the change in superrotation of Earth’s inner core, , as a result of bolide impacts of various diameters and speeds. Bolide densities are assumed to be comparable to that of Earth and the impacts are in the direction = 180°. The limits of the measured values of Earth’s inner core superrotation are also shown in the same figure for comparison.

Immediately after impact, the value of the differential rotation of the inner core relative to the mantle changes to . To find out how this new value decays, we solve (11) using the initial condition . This givesBecause = 11600 yr, (29) shows that in a time scale of a few tens of thousand years the effect of bolide impact damps out. Therefore, impacts only contribute to geologically short-term perturbations in the steady-state value of the inner-core superrotation.

4. Discussion and Conclusion

There are a number of existing theories explaining the superrotation of Earth’s inner core. Most of these theories are based on internal Earth processes. Although each theory provides a reasonable account for the phenomenon, there is not sufficient experimental evidence to support one theory over the others. The external tidal effect, on the other hand, appears to be one of the most reasonable explanations, which is also supported by experimental evidence and measurements.

The frictional torques generated by the tidal effects tend to slow down the rotational motion of the mantle, thereby increasing the superrotation of the inner core. This effect is counterbalanced by the viscous forces of the Earth’s fluid outer core as well as other effects such as the electromagnetic and gravitational couplings and eddy viscosity caused by fluid motion in the outer core. The combined effect gives rise to an effective viscosity of the outer core on the order of 103 Pa s, in agreement with a previously reported value based on a different model [24]. The viscosity of the molten material of the outer core without these additional contributions, however, could be considerably less. Because the relaxation time for the inner core superrotation is of the order of 104 yr, the superrotation today is almost exactly equal to its saturated or steady-state value.

Bolide impacts, yet another external factor, are possibly responsible for sudden changes in the differential rotation of Earth’s inner core from time to time and can render changes comparable in magnitude to the presently observed values. Although the existing geological evidence indicates few major bolide impacts with Earth in recent geological times, the large number of craters on the Moon strongly suggests that Earth, with a considerably larger cross-sectional area, has also been subject to a large number of impacts since its formation. Most of the resulting craters, however, have been hidden or destroyed by Earth’s weathering and tectonic processes. Therefore, the number of bolide impacts on Earth is not limited to the number of existing craters but rather is much higher. Furthermore, as our science and technology advances, we discover more impact craters on Earth of which we were previously unaware. For example, it was originally believed that the Vredefort dome in South Africa was formed by a volcanic explosion. But in the mid-1990s, evidence revealed that it was in fact the site of a huge bolide impact about two billion years ago. Evidence of four impacts older than Vedefort that occurred between 3.2 to 3.5 billion years ago has also been found in South Africa. These sites, however, are not easily recognizable as impact structures on Earth’s surface.

Asteroids approaching Earth and the Moon have an average speed of about 17 km s−1 and, depending on their composition, can have densities ranging in the extreme from about 3 kg m−3 (rock) to about 8 kg m−3 (iron). The diameters can be as high as tens of kilometers. As shown in Figure 2, there is a wide range of diameters and velocities for an impacter which can cause superrotations of Earth’s inner core comparable to the observed values. In fact, several asteroids among the near Earth objects fall in this category. For example, 1036 Ganymed, the largest asteroid in this class is about 32 km across. According to Figure 2, with an impact speed of about 10 km s−1, an asteroid of this size can cause a superrotation in Earth’s inner core that is in the range of the observed values. Furthermore, cumulative effects from smaller impacts may also result in such superrotations.

Nonetheless, as stated earlier, these impacts that can change the steady-state differential rotation of the Earth’s inner core can be viewed as short lived perturbations in geological terms, with decaying times of the order of a few tens of thousands of years. The discrepancy between our damping time scale of a few tens of thousands of years and those obtained by Tkalčić et al. [19] and by Livermore et al. [20] of the order of decades is the result of invoking different mechanisms and models for the inner-core superrotation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.