On the Mechanism of the Magnetic Dynamo of the Planets

Sh. Sh. Dolginov Academy of Sciences, Moscow, U.S.S.R.

Results of testing the effectiveness of the theory of precessional dynamos in the genera- tion of the magnetic fields of the planets are presented. It is shown that the magnetic state of and of the planets , Jupiter, and can be satisfactorily described by the formula

17 — IT v I J 3 "i fli — fl3 -^j- -fff — sin«3 where H, v, T, 0, and a are the dipole fields, volumes of liquid cores, periods of rotation, rates of precession, and angles between precession vector and angular rotation, respec- tively, for the planets and Earth. The v, corresponds to known models of the internal structure. It is shown that the magnetic state of Mercury satisfies this formula if the dynamic flattening of the planet t = 5.7 x 10~3 - 8.3 x 10s.

The idea that the magnetic field of the With a model of a purely iron-nickel core, Earth is related to the operation of some the radioactive elements U and Th are forced sort of dynamo mechanism in its highly con- out because of their chemical and physical ductive liquid core has now been largely con- incompatibility with iron at the pressures firmed. present in the core. Thermal sources in the Modern kinematic models of the terrestrial core have been related to processes of con- magnetic dynamo have been found capable tinuing differentiation of matter in the of describing the basic peculiarities of the Earth: the settling of heated iron, melted terrestrial magnetic field (refs. 1 and 2). At from the mantle (ref. 3) and the latent heat the same time, no physical theory of the ter- of melting, liberated upon crystallization of restrial field as yet presented considers the the outer core at the boundary with the inner actual parameters and processes in the core solid core (ref. 4). of the Earth. This is particularly true con- In order for thermal convection to occur in cerning the uncertainty of the mechanism the liquid core of the Earth, the temperature generating the terrestrial dynamo. Three gradient must exceed the adiabatic gradient. mechanisms are known: In a number of publications during 1971- 1973 (ref. 8), concepts of the isothermal 1. Convection in the core under the in- state of the core were discussed. A stable fluence of thermal sources (refs. 3 and thermal state is incompatible with thermal 4). convection in the core, to which the genera- 2. Convection under the influences of tion of the magnetic field is related. Ideas of forces of gravity (ref. 5). the thermal conditions of the Earth's core 3. Convection in the core caused by pre- have changed, on the basis of assumptions of cession of the axis of rotation of the the presence of radioactive potassium K40 in Earth (refs. 6 and 7). the core (ref. 9).

493 494 COSMOCHEMISTRY OF THE AND PLANETS

The effectiveness of thermal sources in any liquid in the core, leading to the motion nec- mechanism of generation of the magnetic essary for the operation of the dynamo. In field has been assumed to be low, because of the final analysis, the energy of the magnetic the low efficiency of thermal machines. Ac- dynamo is taken from the kinetic energy in cording to one hypothesis (ref. 5), the pro- the Earth's rotation. cesses of differentiation of material indicated The details of the precessional mechanism by Urey and Verhoogen are actually in- generating the magnetic field have not been volved in the generation of the magnetic field, worked out. There are contradictory opinions but only because of the gravitational energy concerning its effectiveness. Braginskiy (ref. released during these processes. Convection 12) turned his attention to the fact that the in the core is caused by direct upwelling of speed of the liquid under the influence of the impurities of light silicon, accompanying the Poincare force fluctuates with the frequency process of crystallization and sinking of of the main rotation <,> and therefore cannot heavy iron. directly influence slow processes in the dy- The hypothesis of the gravitational nature namo. On the other hand Dolginov (ref. 13) of the "motor" driving the Earth's dynamo notes that the precession of a liquid in a encounters great difficulties in light of cur- nonspherical envelope results in the appear- rent hypotheses of the nonhomogeneous for- ance of flows of a spiral nature that are sim- mation of the planets (refs. 10 and 11). As ilar to tidal flows. Flows of a spiral nature concerns thermal sources, their role is at facilitate the generation of magnetic fields least as great in the creation of the necessaiy (ref. 14). conductivity and viscosity in the core. At the present time, the magnetic fields of Malkus suggested that the precession of the Earth, Venus, Mars, Jupiter, and Mer- the Earth was the motive force of the mag- cury have been studied. For the first four netic dynamo. The precession of the planets planets, certain models of the internal struc- is caused by the action of the gravitational ture are known (refs. 15-20), as well as the fields of the Sun and other celestial bodies on parameters of rotation and dynamic flatten- the of the rotating planets. ing (refs. 21 and 22), and magnetic fields The presence of precession, at an angular (refs. 23-27). velocity n, results in an additional angular This paper presents a model and results of acceleration n = [0 X «] and an additional a test of the effectiveness of the mechanism inertial force Fr, — — />[n X «] X f. Malkus of the precessional dynamo in the generation called this the Poincare force, after the man of the magnetic fields of the Earth, Mars, and who first solved analytically the problem of Jupiter. Furthermore, conditions are studied the motion of a fluid in a precessing spheroi- under which the magnetic states of Venus dal container. The experiments performed by and Mercury can be explained, on the basis Malkus showed that with certain relation- of the hypothesis of precession as the motive ships of w and n in the precessing liquid, tur- force of the planetary dynamo. bulent flow arises, which is not predicted by The approach is based on the following as- classical theory. sumptions : The rate of precession of the planet is di- rectly proportional to its dynamic flattening, 1. The processes of generation of the field and independent of its radius and mass. Since in the planets are similar in terms of model- the core and mantle of the Earth have some- ing and are determined by a number of di- what differing dynamic flattening, the grav- mensionless parameters. ity fields of the Sun and the Moon create 2. The magnetic fields are proportional to different torques on the core and mantle, the volumes of the liquid cores of the planets causing stress both in the core and in the (ref. 28). mantle, which tend to equalize the rates of 3. The rates of convection of matter are the two precessions. These stresses act on the proportional to the Poincare force Frj MAGNETIC DYNAMO OF THE PLANETS 495

—' [ is the angular veloc- the metallic state. According to Zharkov et ity of rotation, and n is the angular velocity al. (ref. 19), this occurs at depths of of precession. r/Rj = 0.8. The magnetically active-area ap- 4. Under these conditions, a simple depen- parently begins at a point ranging from dence between the dipole fields of two planets r/Rj = 0.7 to r/Rj = 0.15. The volume of the 1 and 2 is assumed to exist: magnetically active area of Jupiter, with

V, T2 a sin these dimensions, exceeds the volume of the Yl-2 sin , (1) liquid core of the Earth by some 1600 times. Now, formula (1) and the data of table 1 where y is the ratio of the dipole fields -~- lead to the following values of calculated and on the surface at the magnetic equator; « is measured field strength ratios between the the angle between vectors , the volume of the liquid core of the second Venus planet. The dimensions of the core of Mars The planet Venus apparently does not have have been very approximately determined. a magnetic field of its own which exceeds 10 The measured values of radius, mass, and gammas (refs. 24 and 27). It is natural to moment of inertia are satisfied in various expect that formula (1) will take note of models with significantly differing core di- this fact with some accuracy. According to mensions Rc. In the model of Kozlovskaya the data from the trajectory measurements (ref. 15), Rc = 960 km; in the model of Ring- on Mariner 5 and radar studies of the wood (ref. 16), RCfxt =1720 km and Rc.ni planet, the radius of Venus Rv — 6052.5 = 1510 km; in the model of Binder and Davis ± 2.5km (ref. 29). (ref. 17), Rr = 1250 km; in the model of The period of rotation, as determined by Anderson (ref. 18), Rc = 1350 km; and in Dyce et al. (ref. 30), is 245.1 ± 0.7 days; the model of Johnston and Toksoz (ref. 19), according to the data of Carpenter (ref. 31), Rc = 1300 km. The ratios of volumes of the it is ~ 243 ± 1 days. cores of the Earth and Mars which they ob- The inclination of the axis of rotation rela- tain are 45, 26, 20, 16, and 18, respectively. tive to the pole of the , according to The dimensions of the magnetically active Dyce et al. (ref. 30), is -- 3.3 ± 0.4° ; accord- area of Jupiter are determined primarily by ing to Carpenter (ref. 31), ^ 2.2°. Trajec- the depth at which hydrogen is converted to tory measurements of Mariner 5 (ref. 32)

Table 1.—Parameters of Rapidly Rotating Planets fl Ho, gammas Earth 23h 56m 04s 50.25 "/yr 23.5° 30000 Mars 24h 37m 23s 7.40 "/yr 23.2° 64 Jupiter 9h 50m 56s 2.34 "/yr 3° 400 000 496 COSMOCHEMISTRY OF THE MOON AND PLANETS

have been used to estimate the degree of gle of inclination of the axis of rotation to dynamic flattening / = 1/100 of this terres- the plane of the orbit, widely different values trial planet. This value has been confirmed have been given. According to radar data by the data from Mariner 10 (/ = 1/30 000). (ref. 30), an angle of a =* 28° is given. The precession of Venus results from the in- According to terrestrial photographs of Mer- fluence of the dipole gravitational field of the cury (ref. 35), the axis of rotation is perpen- Sun. The rate of presession of Venus, with dicular to the plane of the orbit with an / = 1/30 000, is nB = 122"/yr. We can fur- accuracy of —- 3°. It is assumed (ref. 26) ther assume that the volume of the liquid that the axis of the dipole makes an angle of core of Venus is equal to or only slightly less 80 ± 10° with the plane of the ecliptic. The than the volume of the liquid core of Earth orbital plane of Mercury makes an angle of (ref. 20). Then formula (1) gives us, with —- 7° with the plane of the ecliptic. a = 2.2° and 3.3°, a field intensity H0 < 20- Studies of the internal structure of Mer- 30 gammas in comparison with the terrestrial cury, based on determinations of mass and magnetic field. radius, lead to comparatively large dimen- Thus, proceeding from formula (1), there sions for the core of Mercury. Plageman (ref. is no need to postulate any difference in the 36) indicates Rc = 2112 km, and Kozlovskaya internal structure of Venus from that of (ref. 37) indicates Rc = 1730 km. Earth, except on the basis of the fact that With the existing uncertainties of most it has no significant magnetic field (refs. 33 parameters, there is reason to take certain and 34). arbitrary parameters and estimate from them, the order of magnitude of a quite un- Mercury known quantity, the dynamic flattening. Let us take Hn = 250 gammas, R, = 1800 km, For Mercury, data on most of the param- a = 20°. Then, comparing with the Earth, eters included in formula (1) are quite in- from formula (1) we can estimate the rate definite. Nevertheless, there is reason to use of precession of Mercury: QMerc = 190"/yr. (1) to estimate the limits of dynamic flatten- The dynamic flattening of Mercury /Mere can ing of the planet Mercury. The flight of be estimated by comparing the rate of preces- Mariner 10 produced new experimental data sion of Mars and Mercury, since both planets on the parameters of rotation and, we can precess under the influence of the gravita- hope, in the near future it will be possible to tional field of the Sun. compare results and take into consideration /„ = ^^ . T<$ .flMer c . /j / Merc r>3 rrt „ * JO the new experimental results. K Merc *• Merc «O For the value of the dipole magnetic field, Then where Ho = 250 gammas, / rc = 5.7 • 10~5 assuming that it is intrinsic to the planet Me 5 (ref. 26), there are at present two values HH = 380 gammas, /Merc = 8.3 • 1Q- available: 227 gammas and 380 gammas. The These values are somewhat greater than first of these values corresponds to the theo- those for the dynamic flattening of Venus. retical dipole, somewhat inclined to the axis Thus, even with the current uncertainty of rotation and displaced from the center by for a number of parameters of Mercury, for- 0.47 km, to achieve best agreement with the mula (1) indicates that the fact of the exis- experimental data. The second value is pro- tence of a magnetic field of Mercury (ref. duced from the data of the gasdynamic model 26) is not too surprising. of solar wind flow around the planet. The shock front was intersected in this experi- Correspondence of the Dynamo ment at an angle of ~ 110° with respect to Model the direction to the Sun. It is now generally accepted that Mercury 1. In the theory of the terrestrial dynamo, rotates with a period of 58 days. For the an- the inclination of the axis of the dipole is not MAGNETIC DYNAMO OF THE PLANETS 497 considered to be a chance phenomenon, but 9. VERHOOGEN, J., Phys. Earth Planet. Interiors, rather a fact directly related to the mec- No. 7, 1973, pp. 47-58. 10. TUREKIAN, K. K., AND S. P. CLARK, Earth hanism of generation of the field (ref. 12). Planet. Sci. Letters, No. 6, 1969, p. 346. Paleomagnetologists, on the other hand, 11. VINOGRADOV, A. P., International Geochemical consider the current orientation of the mag- Congress, 1971 Collection, Vol. 1. netic dipole to be a brief (in the geological 12. BRAGINSKIY, S. I., Geomagnetizm i Aeronomiya, scale) deviation from the position of sym- Vol. 4, No. 4, 1964. 13. DOLGINOV, A. Z., Astron. Zhurn. Vol. 51, No. 2, metry (ref. 38). 1974. The magnetic dipoles of Jupiter, Mars, 14. LlLLEY, F. E., Pros. Ron. Soc., A. 316, No. 525, and—apparently—Mercury, are inclined in 1970, p. 153. relation to their axes of rotation. Further- 15. KOZLOVSKAYA, S. V., Izvestiya AN SSSR, Fizika more, for Jupiter, as for Earth, the dipole is Zemli, No. 7, 1972. 16. RINGWOOD, A. E., AND S. P. CLARK, Nature, Vol. eccentric. The eccentricity of the dipole ap- 234, No. 5324, 1971. pears to be particularly great for Mercury. 17. BINDER, A. B., AND DAVIS, Phys. Earth Planet. 2. Mercury, Earth, Mars, and Jupiter Interiors, Vol. 7, 1973, p. 477. have direct rotation. However, the polarity of 18. ANDERSON, D. L., Comments on Earth Sci. the magnetic fields of Mars and Jupiter is Geophys., Vol. 1, No. 5, 1971. 19. JOHNSTON, D. H., T. R. McGETCHiN, AND M. F. the reverse of the polarity of the current TOKSOZ, Preprint, April, 197.4. fields of Earth and Mercury. In kinematic 20. ZHARKOV, V. N., V. P. TRUBITSIN AND L. V. models of the terrestrial dynamo, the sign of SAMSONENKO, Fizika Zemli i Planet, 1971, p. the field is not directly related to the direc- 383. tion of rotation, and the possibility of field 21. LORELL, J., et al., Icarus, Vol. 28, No. 2, 1973, p. 304. reversals is explained by the action of insta- 22. KOVALEVSKY, J., Surface and Interiors of Planets bilities in the mechanism of generation. The and Satellites, 1970. possible relationship of terrestrial field re- 23. DOLGINOV, SH. SH., YE. G. EROSHENKO, AND L. N. versals with precession of its axis has been ZHUZGOV, DAN SSSR, Vol. 207, No. 6, 1972. indicated by (ref. 39). DAN, 1974. (In press). 24. DOLGINOV, SH. SH., YE. G. EROSHENKO, AND Experiments in the third decade of the L. N. ZHUZGOV, Kosmich. Issled., 1967. space age should clarify the degree of regu- 25. SMITH, E. J., L. J. DAVIS, D. E. JONES, D. S. COL- larity of the connections revealed between BURN, P. J. COLEMAN, P. DYAL, AND C. P. the planets' parameters of rotation, their SONETT, Science, Vol. 183, 1974, p. 305. structure, and their magnetic fields. 26. NESS, N. F., K. W. BEHANNON, R. P. LEPPING, Y. C. WHANG, AND K. H. SCHATTEN, Science, 1974. In press. 27. NESS, N. F., K. W. BEHANNON, R. P. LEPPING, References Y. C. WHANG, AND K. H. SCHATTEN, Science, Vol. 183, 1974, p. 1301. 1. BRAGINSKIY, S. I., Izvestiya AN SSSR, Fizika 28. KERN, J. W., AND E. H. VESTINE, Space Science Zemli, No. 10, 1972. Reviews, 1963, p. 136. 2. Origin of the Geomagnetic Field in World Mag- 29. ANDERSON, J. D., J. Atmosph. Sci., Vol. 25, No. netic Survey 1957-1969. JACA Bulletin, No. 6, 1968. 28, 1971. 30. DYCE, R. B., G. H. PETTENGILL, AND I. SHAPIRO, 3. UREY, H. C., The Planets. Yale Univ. Press, Astron. Zhurn. Vol. 72, No. 3, 1967, p. 371. 1952. 31. CARPENTER, R. L., Astron. Zhurn. Vol. 71, No. 9, 4. VERHOOGEN, J., Geophys. J. R.A.S., Vol. 4, No. 1966. 276, 1961. 32. ANDERSON, J. D., L. EFRON, AND G. E. PEASE, '5. BRAGINSKIY, S. I., Geomagnetizm I Aeronomiya, Astron. Zhurn., Vol. 73, No. 10, Part 2, 1968. Vol. 4, No. 5, 1964. 33. HENDERVARI, P., Observatory 86, No. 953, 1966. 6. MALKUS, W. V. R., J. Geophys. Res., No. 68, 34. RUNCORN, S. K., Geophys, J. Roy. Astron. Soc., 1963, p. 2871. Vol. 15, Nos. 1 and 2, 1968. 7. MALKUS, W. V. R., Science, Vol. 160, No. 3825, 35. MURRAY, J. B., A. DOLLFUS, AND B. SMITH, 1968. Icarus, Vol. 17, No. 3, 1972. 8. KENNEDY, G. G., AND G. H. HIGGINS, The Moon, 36. PLACEMAN, S., J. Geophys. Res., Vol. 70, No. 4, Vol. 7, Nos. 1 and 2, 1973. 1965. 498 COSMOCHEMISTRY OF THE MOON AND PLANETS

37. KOZLOVSKAYA, S. V., Physics of the Moon and 38. STACY, Fizika Zemli, 1972. Planets, International Symposium, October 39. KALININ, Yu. D., IFSO Preprint, 2 F, 1972. 15-22, 1968, Kiev, 1972.