Geodesy constraints on the interior structure and composition of Mars

Knowledge of the interior structure of Mars is of fundamental importance to the understanding of its past and present state as well as its future evolution. The most prominent interior structure properties are the state of the core, solid or liquid, its radius, and its composition in terms of light elements, the thickness of the mantle, its composition, the presence of a lower mantle, and the density of the crust. In the absence of seismic sounding only geodesy data allow reliably constraining the deep interior of Mars. Those data are the mass, moment of inertia, and tides. They are related to Mars’ composition, to its internal mass distribution, and to its deformational response to principally the tidal forcing of the Sun. Here we use the most recent estimates of the moment of inertia and tidal Love number k₂ in order to infer knowledge about the interior structure of the Mars.We have built precise models of the interior structure of Mars that are parameterized by the crust density and thickness, the volume fractions of upper mantle mineral phases, the bulk mantle iron concentration, and the size and the sulfur concentration of the core. From the bulk mantle iron concentration and from the volume fractions of the upper mantle mineral phases, the depth dependent mineralogy is deduced by using experimentally determined phase diagrams. The thermoelastic properties at each depth inside the mantle are calculated by using equations of state. Since it is difficult to determine the temperature inside the mantle of Mars we here use two end-member temperature profiles that have been deduced from studies dedicated to the thermal evolution of Mars. We calculate the pressure and temperature dependent thermoelastic properties of the core constituents by using equations state and recent data about reference thermoelastic properties of liquid iron, liquid iron-sulfur, and solid iron. To determine the size of a possible inner core we use recent data on the melting temperature of iron-sulfur.Within our model assumptions the geodesy data imply that Mars has no solid inner core and that the liquid core contains a large fraction of sulfur. The absence of a solid inner is in agreement with the absence of a global magnetic field. We estimate the radius of the core to be 1794 ± 65 km and its core sulfur concentration to be 16 ± 2 wt%. We also show that it is possible for Mars to have a thin layer of perovskite at the bottom of the mantle if it has a hot mantle temperature. Moreover a chondritic Fe/Si ratio is shown to be consistent with the geodesy data, although significantly different value are also possible. Our results demonstrate that geodesy data alone, even if a mantle temperature is assumed, can almost not constrain the mineralogy of the mantle and the crust. In order to obtain stronger constraints on the mantle mineralogy bulk properties, like a fixed Fe/Si ratio, have to be assumed.     

The influence of Mercury’s inner core on its physical libration

As Mercury orbits the Sun, gravitational torques on its equatorial elliptical shape give rise to a planetary libration. The amplitude of Mercury’s libration, as determined from Earth-based radar speckle pattern observations, suggests that only the mantle participates in the motion. This indicates a decoupling between the core and the mantle, and therefore that the outermost part of the core must be fluid. If a solid inner core is present at the center of Mercury, the equatorial elliptical shape of the latter may become misaligned with that of Mercury’s mantle, leading to an internal gravitational torque between the two. If this torque is large, it may participate in the dynamics of Mercury’s libration. The goal of this work is to determine whether Mercury’s observed librations can be used to place constraints on the properties of its inner core. We present a comparison between predicted and observed librations for a range of interior models of Mercury, with various inner core sizes and fluid core densities. We show that a marginally better fit to observations can be achieved for interior models that have an inner core radius larger than 400 km. However, the improvement in fit is small, and it is not possible to draw robust conclusions on the size of Mercury’s inner core on the basis of existing libration data. Nevertheless, our study demonstrates that the influence of the inner core on the libration of Mercury could be detected with a decade worth of accurate observations.     

Viscosity of the Martian mantle and its initial temperature: Constraints from crust formation history and the evolution of the magnetic field

The rheology of the Martian mantle and the planet's initial temperature is constrained with thermal evolution models that include crust growth and test the conditions for magnetic field generation in the core. As observations we use the present-day average crustal thickness of 50-120 km as estimated from the Mars Global Surveyor gravity and topography data, the evidence for the crust being produced mostly early, with a rate declining from the Noachian to the Hesperian, and the evidence for an early magnetic field that likely existed for less than a billion years. We use the fact that the rate of crust growth is a function of temperature, which must be above the solidus in the sub-lithosphere mantle, and the mantle convection speed because the latter determines the rate at which melt can be replenished. The convection speed is a strong function of viscosity which, in turn, is a strong function of temperature and also of the water content of the mantle. We use a viscosity parameterization with a reference viscosity evaluated at 1600 K the value of which can be characteristic of either a dry or a wet mantle. We further consider the Fe-FeS phase diagram for the core and compare the core liquidus estimated for a sulphur content of 14% as suggested by the SNC meteorite compositions with the core temperatures calculated for our cooling models. Two data sets of the Fe-FeS eutectic temperature have been used that differ by about 200 K [Böhler, R., 1996. Fe-FeS eutectic temperatures at 620 kbar. Phys. Earth Planet. Inter. 96, 181-186; Fei, Y., Bertka, C.M., Finger, L.W., 1997. High-pressure iron-sulphur compound, Fe₃S₂, and melting relations in the Fe-FeS system. Science 275, 1621-1623] at Martian core-mantle boundary pressure and in the eutectic composition by 5 wt%. The differences in eutectic temperature and composition translate into a difference of about 400 K in liquidus temperature for 14 wt% sulphur.We find it premature to rule out specific mantle rheologies on the basis of the presently available crustal thickness and crust growth evidence. Rather a trade-off exists between the initial mantle temperature and the reference viscosity. Both a wet mantle rheology with a reference viscosity less than 10²⁰ Pas and a dry mantle rheology with a reference viscosity of 10²¹ Pas or more can be acceptable if initial mantle temperatures between roughly 1700 and 2000 K are allowed. To explain the magnetic field history, the differences in liquidus temperatures matter. For a liquidus temperature of about 1900 K at the Martian core-mantle boundary as calculated from the Böhler et al. eutectic, a dry mantle rheology can best explain the lack of a present-day dynamo. For a liquidus temperature of about 1500 K at the core-mantle boundary as calculated from the Fei et al. eutectic all models are consistent with the observed lack of dynamo action. The reason lies with the fact that at 14 wt% S the Martian core would be close to the eutectic composition if the Fei et al. data are correct. As inner core growth is unlikely for an almost eutectic core, the early field would have been generated by a thermally driven dynamo. Together with the measured strength of the Martian crustal magnetization this would prove the feasibility of a strong thermally driven dynamo.     

Implications from Galileo Observations on the Interior Structure and Chemistry of the Galilean Satellites

Data from the recent gravity measurements by the Galileo mission are used to construct wide ranges of interior structure and composition models for the Galilean satellites of Jupiter. These models show that mantle densities of Io and Europa are consistent with an olivine-dominated mineralogy with the ratios of Mg to Fe components depending on mantle temperature for Io and on ice shell thickness for Europa. The mantle density and composition depend relatively little on core composition. The size of the core is largely determined by the core's composition with core radius increasing with the concentration of a light component such as sulfur. For Io, the range of possible core sizes is between 38 and 53% of the satellite's radius. For Europa, there is also a substantial effect of the thickness of the ice layer which is varied between 120 and 170 km on the core size. Core sizes are between 10 and 45% of Europa's radius. The core size of Ganymede ranges between one-quarter and one-third of the surface radius depending on its sulfur content and the thickness of the ice shell. A subset of the Ganymede models is consistent with an olivine-dominated mantle mineralogy. The thickness of the silicate mantle above the core varies between 900 and 1100 km. The outermost ice shell is about 900 km in thickness and is further subdivided by pressure-induced phase transitions into ice I, ice III, ice V, and ice VI layers. Callisto should be differentiated, albeit incompletely. It is proposed that this satellite was never molten at a large scale but differentiated through the convective gradual unmixing of the ice and the metal/rock component. Bulk iron-to-silicon ratios Fe/Si calculated for the inner pair of satellites, Io and Europa, are less than the CI carbonaceous chondrite value of 1.7±0.1, whereas ratios for the outer pair, Ganymede and Callisto, cover a broad range above the chondritic value. Although the ratios are uncertain, in particular for Ganymede and Callisto, the values are sufficiently distinct to suggest a difference in composition between these two pairs of satellites. This may indicate a difference in iron-silicon fractionation during the formation of both classes of satellites in the protojovian nebula.     

Magnetism and thermal evolution of the terrestrial planets

Of the terrestrial planets, Earth and probably Mercury possess substantial intrinsic magnetic fields generated by core dynamos, while Venus and Mars apparently lack such fields. Thermal histories are calculated for these planets and are found to admit several possible present states, including those which suggest simple explanations for the observations; whule the cores of Earth and Mercury are continuing to freeze, the cores of Venus and Mars may still be completely liquid. The models assume whole mantle convection, which is parameterized by a simple Nusselt-Rayleigh number relation and dictates the rate at which heat escapes from the core. It is found that completely fluid cores, devoid of intrinsic heat sources, are not likely to sustain thermal convection for the age of the solar system but cool to a subadiabatic, conductive state that can not maintain a dynamo. Planets which nucleate an inner core continue to sustain a dynamo because of the gravitational energy release and chemically driven convection that accompany inner core growth. The absence of a significant inner core can arise in Venus because of its slightly higher temperature and lower central pressure relative to Earth, while a Martian core avoids the onset of freezing if the abundance of sulfur in the core is ?15% by mass. All of the models presented assume that (I) core dynamos are driven by thermal and/or chemical convection; (ii) radiogenic heat production is confined to the mantle; (iii) mantle and core cool from initially hot states which are at the solidus and superliquidus, respectively; and (iv) any inner core excludes the light alloying material (sulfur or oxygen) which then mixes uniformly upward through the outer core. The models include realistic pressure and composition-dependent freezing curves for the core, and material parameters are chosen so that the correct present-day values of heat outflow, upper mantle temperature and viscosity, and inner core radius are obtained for the earth. It is found that Venus and Mars may have once had dynamos maintained by thermal convection alone. Earth may have had a completely fluid core and a dynamo maintained by thermal convection for the first 2 to 3 by, but an inner core nucleates and the dynamo energetics are subsequently dominated by gravitational energy release. Complete freezing of the Mercurian core is prohibited if it contains even a small amount of sulfur, and a dynamo can be maintained by chemical convection in a thin, fluid shell.     

Size and compositional constraints of Ganymede's metallic core for driving an active dynamo

Ganymede has an intrinsic magnetic field which is generally considered to originate from a self-excited dynamo in the metallic core. Driving of the dynamo depends critically on the satellite's thermal state and internal structure. However, the inferred structure based on gravity data alone has a large uncertainty, and this makes the possibility of dynamo activity unclear; variations in core size and composition significantly change the heat capacity and alter the cooling history of the core. The main objectives of this study is to explore the structural conditions for a currently active dynamo in Ganymede using numerical simulations of the thermal history, and to evaluate under which conditions Ganymede can maintain the dynamo activity at present. We have investigated the satellite's thermal history using various core sizes and compositions satisfying the mean density and moment of inertia of Ganymede, and evaluate the temperature and heat flux at the core-mantle boundary (CMB). Based on the following two conditions, we evaluate the possibility of dynamo activity, thereby reducing the uncertainty of the previously inferred interior structure. The first condition is that the temperature at the CMB must exceed the melting point of a metallic core, and the second is that the heat flux through the CMB must exceed the adiabatic temperature gradient. The mantle temperature starts to increase because of the decay of long-lived radiogenic elements in the rocky mantle. After a few Gyr, radiogenic elements are exhausted and temperature starts to decrease. As the rocky mantle cools, the heat flux at the CMB steadily increases. If the temperature and heat flux at the CMB satisfy these conditions simultaneously, we consider the case as capable of driving a dynamo. Finally, we identify the Dynamo Regime, which is the specific range of internal structures capable of driving the dynamo, based on the results of simulations with various structures. If Ganymede's self-sustained magnetic field were maintained by thermal convection, the satellite's metallic core would be relatively large and, in comparison to other terrestrial-type planetary cores, strongly enriched in sulfur. The dynamo activity and the generation of the magnetic field of Ganymede should start from a much later stage, possibly close to the present.     

Thermal evolutions of the terrestrial planets

The thermal evolution of the Moon, Mercury, Mars, Venus and hypothetical minor planets is calculated theoretically, taking into account conduction, solid-state convection, and differentiation. An assortment of geological, geochemical, and geophysical data is used to constrain both the present day temperatures and thermal histories of the planets' interiors. Such data imply that the planets were heated during or shortly after formation and that all the terrestrial planets started their differentiations early in their history. Initial temperatures and core formation play the most important roles in the early differentiation. The size of the planet is the primary factor in determining its present day thermal state. A planetary body with radius less than 1000 km is unlikely to reach melting given heat source concentrations similar to terrestrial values and in the absence of intensive early heating such as short half-life radioactive heating and inductive heating.Studies of individual planets are constrained by varying amounts of data. Most data exist for the Earth and Moon. The Moon is a differentiated body with a crust, a thick solid mantle and an interior region which may be partially molten. It is presently cooling rapidly and is relatively inactive tectonically.Mercury most likely has a large core. Thermal calculations indicate it may have a 500 km thick solid lithosphere, and the core may be partially molten if it contains some heat sources. If this is not the case, the planet's interior temperatures are everywhere below the melting curve for iron. The thermal evolution is dominated by core separation and the high conductivity of iron which makes up the bulk of Mercury.Mars, intermediate in size among the terrestrial planets, is assumed to have differentiated an Fe–FeS core. Differentiation and formation of an early crust is evident from Mariner and Viking observations. Theoretical models suggest that melting and differentiation of the mantle silicates has occurred at least up until 1 billion years ago. Present day temperature profiles indicate a relatively thick (250 km) lithosphere with a possible asthenosphere below. The core is molten.Venus is characterized as a planet similar to the Earth in many respects. Core formation probably occurred during the first billion years after the formation. Present day temperatures indicate a partially molten upper mantle overlain by a 100 km thick lithosphere and a molten Fe–Ni core. If temperature models are good indicators, we can expect that today, Venus has tectonic processes similar to the Earth's.Paper dedicated to Professor Hannes Alfvén on the occasion of his 70th birthday, 30 May 1978.     

Mars

Mars is the fourth planet out from the sun. It is a terrestrial planet with a density suggesting a composition roughly similar to that of the Earth. Its orbital period is 687 days, its orbital eccentricity is 0.093 and its rotational period is about 24 hours. Mars has two small moons of asteroidal shapes and sizes (about 11 and 6 km mean radius), the bigger of which, Phobos, orbits with decreasing semimajor orbit axis. The decrease of the orbit is caused by the dissipation of tidal energy in the Martian mantle. The other satellite, Deimos, orbits close to the synchronous position where the rotation period of a planet equals the orbital period of its satellite and has hardly evolved with time. Mars has a tenous atmosphere composed mostly of CO with strong winds and with large scale aeolian transport of surface material during dust storms and in sublimation-condensation cycles between the polar caps. The planet has a small magnetic field, probably not generated by dynamo action in the core but possibly due to remnant magnetization of crustal rock acquired earlier from a stronger magnetic field generated by a now dead core dynamo. A dynamo powered by thermal power alone would have ceased a few billions of years ago as the core cooled to an extent that it became stably stratified. Mars' topography and its gravity field are dominated by the Tharsis bulge, a huge dome of volcanic origin. Tharsis was the major center of volcanic activity, a second center is Elysium about 100° in longitude away. The Tharsis bulge is a major contributor to the non-hydrostaticity of the planet's figure. The moment of inertia factor together with the mass and the radius presently is the most useful constraint for geophysical models of the Martian interior. It has recently been determined by Doppler range measurements to the Mars Pathfinder Lander to be (Folkner et al. 1997). In addition, models of the interior structure use the chemistry of the SNC meteorites which are widely believed to have originated on Mars. According to the models, Mars is a differentiated planet with a 100 to 200 km thick basaltic crust, a metallic core with a radius of approximately half the planetary radius, and a silicate mantle. Mantle dynamics is essential in forming the elements of the surface tectonics. Models of mantle convection find that the pressure-induced phase transformations of -olivine to -spinel, -spinel to -spinel, and -spinel to perovskite play major roles in the evolution of mantle flow fields and mantle temperature. It is not very likely that the -spinel to perovskite transition is present in Mars today, but a few 100 km thick layer of perovskite may have been present in the lower mantle immediately above the core-mantle boundary early in the Martian history when mantle temperatures were hotter than today. The phase transitions act to reduce the number of upwellings to a few major plumes which is consistent with the bipolar distribution of volcanic centers of Mars. The phase transitions also cause a partial layering of the lower mantle which keeps the lower mantle and the core from extensive cooling over the past aeons. A relatively hot, fluid core is the most widely accepted explanation for the present lack of a self-generated magnetic field. Growth of an inner core which requires sub-liquidus temperatures in the core would have provided an efficient mechanism to power a dynamo up to the present day. Received 10 May 1997     

Tides and nutation of the Earth: II. Realistic models of the liquid core

In the first part of this study (Molodensky, 2004; hereinafter, paper I), a survey of the theory of tides and nutation of the Earth was given with the inclusion of the dynamical effects of the elastic mantle and liquid core in the framework of a very simple model of a homogeneous, incompressible core. Here, effects of the inhomogeneity, compressibility, and viscosity of the liquid core are considered, along with those of electromagnetic coupling of the liquid core with the mantle and solid inner core. Errors in the present-day measurements of the amplitudes of forced nutation (of the order of 20 arc s) are well below not only the dynamical effects of the Earths liquid core but also the effects of the inelasticity of the mantle, the dynamical effects of the solid inner core, and the possible effects of electromagnetic coupling between the liquid core, solid inner core, and mantle. This opens up new avenues for astrometric studies of the mechanical properties and electrical conductivity of the mantle and solid inner core at very low frequencies. The modern data on the amplitudes and phases of the Earths forced nutation cannot be interpreted entirely unambiguously, since the following factors remain unknown: (i) the role of the core-mantle electromagnetic coupling compared to the coupling due to core viscosity near the core-mantle boundary and (ii) the frequency dependence of the creep function of the mantle at low frequencies. In large measure, the effects of electromagnetic and viscous coupling can be separated if high-precision data on the tidal variations in the gravitational force at periods of about a day are invoked and allowances are made for the differences between the effects of viscous and electromagnetic coupling on the amplitudes and phases of forced nutation and on the tidal variations in the gravitational force. Here, ranges of possible values of the creep function are constructed for periods from one hour to one day; for these calculations, values consistent with the entire set of data on the forced nutation of the Earth are assumed for the effective dynamical flattening of the liquid core-mantle and liquid core-solid inner core boundaries (these values describe the ellipticity of the boundaries of the inner core, outer core, and mantle, as well as the electromagnetic coupling between the liquid core, mantle, and solid inner core).Translated from Astronomicheskii Vestnik, Vol. 39, No. 1, 2005, pp. 61–80.Original Russian Text Copyright © 2005 by Molodensky.     

Could giant impacts cripple core dynamos of small terrestrial planets?

Large impacts not only create giant basins on terrestrial planets but also heat their interior by shock waves. We investigate the impacts that have created the largest basins existing on the planets: Utopia on Mars, Caloris on Mercury, Aitken on Moon, all formed at ∼4 Ga. We determine the impact-induced temperature increases in the interior of a planet using the “foundering” shock heating model of Watters et al. (Watters, W.A., Zuber, M.T., Hager, B.H. [2009]. J. Geophys. Res. 114, E02001. doi:10.1029/2007JE002964). The post-impact thermal evolution of the planet is investigated using 2D axi-symmetric convection in a spherical shell of temperature-dependent viscosity and thermal conductivity, and pressure-dependent thermal expansion. The impact heating creates a superheated giant plume in the upper mantle which ascends rapidly and develops a strong convection in the mantle of the sub-impact hemisphere. The upwelling of the plume rapidly sweeps up the impact-heated base of the mantle away from the core-mantle boundary and replaces it with the colder surrounding material, thus reducing the effects of the impact-heated base of the mantle on the heat flux out of core. However, direct shock heating of the core stratifies the core, suppresses the pre-existing thermal convection, and cripples a pre-existing thermally-driven core dynamo. It takes about 17, 4, and 5 Myr for the stratified cores of Mars, Mercury, and Moon to exhaust impact heat and resume global convection, possibly regenerating core dynamos.     

Highly reducing conditions during core formation on Mercury: Implications for internal structure and the origin of a magnetic field

The high average density and low surface FeO content of the planet Mercury are shown to be consistent with very low oxygen fugacity during core segregation, in the range 3-6 log units below the iron-wüstite buffer. These low oxygen fugacities, and associated high metal content, are characteristic of high-iron enstatite (EH) and Bencubbinite (CB) chondrites, raising the possibility that such materials may have been important building blocks for this planet. With this idea in mind we have explored the internal structure of a Mercury sized planet of EH or CB bulk composition. Phase equilibria in the silicate mantle have been modeled using the thermodynamic calculator p-MELTS, and these simulations suggest that orthopyroxene will be the dominant mantle phase for both EH and CB compositions, with crystalline SiO₂ being an important minor phase at all pressures. Simulations for both compositions predict a plagioclase-bearing “crust” at low pressure, significant clinopyroxene also being calculated for the CB bulk composition. Concerning the core, comparison with recent high pressure and high temperature experiments relevant to the formation of enstatite meteorites, suggest that the core of Mercury may contain several wt.% silicon, in addition to sulfur. In light of the pressure of the core-mantle boundary on Mercury (∼7 GPa) and the pressure at which the immiscibility gap in the system Fe-S-Si closes (∼15 GPa) we suggest that Mercury’s core may have a complex shell structure comprising: (i) an outer layer of Fe-S liquid, poor in Si; (ii) a middle layer of Fe-Si liquid, poor in S; and (iii) an inner core of solid metal. The distribution of heat-producing elements between mantle and core, and within a layered core have been quantified. Available data for Th and K suggest that these elements will not enter the core in significant amounts. On the other hand, for the case of U both recently published metal/silicate partitioning data, as well as observations of U distribution in enstatite chondrites, suggest that this element behaves as a chalcophile element at low oxygen fugacity. Using these new data we predict that U will be concentrated in the outer layer of the mercurian core. Heat from the decay of U could thus act to maintain this part of Mercury’s core molten, potentially contributing to the origin of Mercury’s magnetic field. This result contrasts with the Earth where the radioactive decay of U represents a negligible contribution to core heating.     

Did 26Al and impact‐induced heating differentiate Mercury?

Numerical models dealing with the planetary scale differentiation of Mercury are presented with the short‐lived nuclide, ²⁶Al, as the major heat source along with the impact‐induced heating during the accretion of planets. These two heat sources are considered to have caused differentiation of Mars, a planet with size comparable to Mercury. The chronological records and the thermal modeling of Mars indicate an early differentiation during the initial ~1 million years (Ma) of the formation of the solar system. We theorize that in case Mercury also accreted over an identical time scale, the two heat sources could have differentiated the planets. Although unlike Mars there is no chronological record of Mercury's differentiation, the proposed mechanism is worth investigation. We demonstrate distinct viable scenarios for a wide range of planetary compositions that could have produced the internal structure of Mercury as deduced by the MESSENGER mission, with a metallic iron (Fe‐Ni‐FeS) core of radius ~2000 km and a silicate mantle thickness of ~400 km. The initial compositions were derived from the enstatite and CB (Bencubbin) chondrites that were formed in the reducing environments of the early solar system. We have also considered distinct planetary accretion scenarios to understand their influence on thermal processing. The majority of our models would require impact‐induced mantle stripping of Mercury by hit and run mechanism with a protoplanet subsequent to its differentiation in order to produce the right size of mantle. However, this can be avoided if we increase the Fe‐Ni‐FeS contents to ~71% by weight. Finally, the models presented here can be used to understand the differentiation of Mercury‐like exoplanets and the planetary embryos of Venus and Earth.     

Long-term evolution of the spin of Mercury: I. Effect of the obliquity and core-mantle friction

The present obliquity of Mercury is very low (less than 0.1°), which led previous studies to always adopt a nearly zero obliquity during the planet’s past evolution. However, the initial orientation of Mercury’s rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long-term evolution of the spin, including the obliquity contribution.The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.     

The production of Ganymede's magnetic field

One of the great discoveries of NASA's Galileo mission was the presence of an intrinsically produced magnetic field at Ganymede. Generation of the relatively strong (750 nT) field likely requires dynamo action in Ganymede's metallic core, but how such a dynamo has been maintained into the present epoch remains uncertain. Using a one-dimensional, three layer thermal model of Ganymede, we find that magnetic field generation can only occur if the sulfur mass fraction in Ganymede's core is very low (?3%) or very high (?21%), and the silicate mantle can cool rapidly (i.e. it has a viscosity like wet olivine). However, these requirements are not necessarily compatible with cosmochemical and physical models of the satellite. We therefore investigate an alternative scenario for producing Ganymede's magnetic field in which passage through an eccentricity pumping Laplace-like resonance in Ganymede's past enables present day dynamo action in the metallic core. If sufficient tidal dissipation occurs in Ganymede's silicate mantle during resonance passage, silicate temperatures can undergo a runaway which prevents the core from cooling until the resonance passage ends. The rapid silicate and core cooling that follows resonance escape triggers dynamo action via thermal and/or compositional convection. To test the feasibility of this mechanism we couple our thermal model with an orbital evolution model to examine the effects of resonance passage on Ganymede's silicate mantle and metallic core. We find that, contrary to expectations, there are no physically plausible scenarios in which tidal heating in the silicates is sufficient to cause the thermal runaway necessary to prevent core cooling. These findings are robust to variations in the silicate rheology, tidal dissipation factor of Jupiter (QJ), structure of the ice shell, and the inclusion of partial melting in the silicate mantle. Resonance passage therefore appears unlikely to explain Ganymede's magnetic field and we must appeal to the special conditions described above to explain the presence of the field.     

Planetary core size: A seismological approach

A key parameter for understanding the geodynamics of a terrestrial planet is the size of its core. Numerical evaluation of 28 different interior structure models of Mercury, Venus, Earth, the Moon, and Mars suggests that there is an almost linear relationship between the core radius and the extent of the seismic P-wave core shadow. A scaling law is derived from a simple mantle density and velocity model that permits the interpretation of respective seismic measurements on terrestrial planetary bodies.     

Implications of large elastic thicknesses for the composition and current thermal state of Mars

The martian elastic lithosphere thickness Te has recently been constrained by modeling the geodynamical response to loading at the martian polar caps and Te was found to exceed 300 km at the north pole today. Geological evidence suggests that Mars has been volcanically active in the recent past and we have reinvestigated the martian thermal evolution, identifying models which are consistent with Te>300 km and the observed recent magmatic activity. We find that although models satisfying both constraints can be constructed, special assumptions regarding the concentration and distribution of radioactive elements, the style of mantle convection and/or the mantle's volatile content need to be made. If a dry mantle rheology is assumed, strong plumes caused by, e.g., a strongly pressure dependent mantle viscosity or endothermic phase transitions near the core-mantle boundary are required to allow for decompression melting in the heads of mantle plumes. For a wet mantle, large mantle water contents of the order of 1000 ppm are required to allow for partial mantle melting. Also, for a moderate crustal enrichment of heat producing, elements the planet's bulk composition needs to be 25 and 50% sub-chondritic for dry and wet mantle rheologies, respectively. Even then, models resulting in a globally averaged elastic thicknesses of Te>300 km are difficult to reconcile with most elastic thickness estimates available for the Hesperian and Amazonian periods. It therefore seems likely that large elastic thicknesses in excess of 300 km are not representative for the bulk of the planet and that Te possibly shows a large degree of spatial heterogeneity.     

Mercury radar speckle dynamics

Current data reveal that Mercury is a dynamic system with a core which has not yet solidified completely and is at least partially decoupled from the mantle. Radar speckle displacement experiments have demonstrated that the accuracy in spin-dynamics determination for Earth-like planets can approach 10⁻⁵. The extended analysis of space-time correlation properties of radar echoes shows that the behavior of speckles does not prevent estimation of Mercury’s instantaneous spin-vector components to accuracy of a few parts in 10⁷. This limit can be reached with more powerful radar facilities and leads to constraining the interior in more detail from effects of spin dynamics, e.g., from observation of the core-mantle interplay through high precision monitoring of the 88-day spin-variation of Mercury’s crust.     

Construction of Martian Interior Model

We present the results of extensive numerical modeling of the Martian interior. Yoder et al. in 2003 reported a mean moment of inertia of Mars that was somewhat smaller than the previously used value and the Love number k₂ obtained from observations of solar tides on Mars. These values of k₂ and the mean moment of inertia impose a strong new constraint on the model of the planet. The models of the Martian interior are elastic, while k₂ contains both elastic and inelastic components. We thoroughly examined the problem of partitioning the Love number k₂ into elastic and inelastic components. The information necessary to construct models of the planet (observation data, choice of a chemical model, and the cosmogonic aspect of the problem) are discussed in the introduction. The model of the planet comprises four submodels—a model of the outer porous layer, a model of the consolidated crust, a model of the silicate mantle, and a core model. We estimated the possible content of hydrogen in the core of Mars. The following parameters were varied while constructing the models: the ferric number of the mantle (Fe#) and the sulfur and hydrogen content in the core. We used experimental data concerning the pressure and temperature dependence of elastic properties of minerals and the information about the behavior of Fe(γ-Fe ), FeS, FeH, and their mixtures at high P and T. The model density, pressure, temperature, and compressional and shear velocities are given as functions of the planetary radius. The trial model M13 has the following parameters: Fe#=0.20; 14 wt % of sulfur in the core; 50 mol % of hydrogen in the core; the core mass is 20.9 wt %; the core radius is 1699 km; the pressure at the mantle-core boundary is 20.4 GPa; the crust thickness is 50 km; Fe is 25.6 wt %; the Fe/Si weight ratio is 1.58, and there is no perovskite layer. The model gives a radius of the Martian core within 1600–1820 km while ≥30 mol % of hydrogen is incorporated into the core. When the inelasticity of the Martian interior is taken into account, the Love number k₂ increases by several thousandths; therefore, the model radius of the planetary core increases as well. The prognostic value of the Chandler period of Mars is 199.5 days, including one day due to inelasticity. Finally, we calculated parameters of the equilibrium figure of Mars for the M13 model: J ₂ ⁰ = 1.82 × 10^?3, J ₄ ⁰ = ?7.79 × 10^?6, e _c-m ^D = 1/242.3 (the dynamical flattening of the core-mantle boundary).     

The evolution of the moon

The thermal evolution of the Moon as it can be defined by the available data and theoretical calculations is discussed. A wide assortment of geological, geochemical and geophysical data constrain both the present-day temperatures and the thermal history of the lunar interior. On the basis of these data, the Moon is characterized as a differentiated body with a crust, a 1000-km-thick solid mantle (lithosphere) and an interior region (core) which may be partially molten. The presence of a crust indicates extensive melting and differentiation early in the lunar history. The ages of lunar samples define the chronology of igneous activity on the lunar surface. This covers a time span of about 1.5 billion yr, from the origin to about 3.16 billion yr ago. Most theoretical models require extensive melting early in the lunar history, and the outward differentiation of radioactive heat sources.Thermal history calculations, whether based on conductive or convective computation codes define relatively narrow bounds for the present day temperatures in the lunar mantle. In the inner region of the 700 km radius, the temperature limits are wider and are between about 100 and 1600°C at the center of the Moon. This central region could have a partially or totally molten core.The lunar heat flow values (about 30 ergs/cm²s) restrict the present day average uranium abundance to 60 ± 15 ppb (averaged for the whole Moon) with typical ratios of K/U = 2000 and Th/U = 3.5. This is consistent with an achondritic bulk composition for the Moon.The Moon, because of its smaller size, evolved rapidly as compared to the Earth and Mars. The lunar interior is cooling everywhere at the present and the Moon is tectonically inactive while Mars could be and the Earth is definitely active.     

The equatorial shape and gravity field of Mercury from MESSENGER flybys 1 and 2

On 14 January and 6 October 2008 the MESSENGER spacecraft passed within 200 km of the surface of Mercury. These flybys by MESSENGER provided the first observations of Mercury from a spacecraft since the Mariner 10 flybys in 1974 and 1975. Data from the Mercury Laser Altimeter (MLA) provided new information on the equatorial shape of Mercury, and Doppler tracking of the spacecraft through the flybys provided new data on the planet’s gravity field. The MLA passes were on opposite hemispheres of the planet and span collectively ∼40% of the equatorial circumference. The mean elevation of topography observed during flyby 1, in the longitude range 0-90°E, is greater than that seen during flyby 2 in the longitude range 180-270°E, indicating an offset between centers of mass and figure having a magnitude and phase in general agreement with topography determined by Earth-based radar. Both MLA profiles are characterized by slopes of ∼0.015° downward to the east, which is consistent with a long-wavelength equatorial shape defined by a best-fitting ellipse. The Doppler tracking data show sensitivity to the gravitational structure of Mercury. The equatorial ellipticity of the gravitational field, C_2,2, is well determined and correlates with the equatorial shape. The S_2,2 coefficient is ∼0, as would be expected if Mercury’s coordinate system, defined by its rotational state, is aligned along its principal axes of inertia. The recovered value of the polar flattening of the gravitational potential, J₂, is considerably lower in magnitude than the value obtained from Mariner 10 tracking, a result that is problematic for internal structure models. This parameter is not as well constrained as the equatorial ellipticity because the flyby trajectories were nearly in the planet’s equatorial plane. The residuals from the Doppler tracking data suggest the possibility of mascons on Mercury, but flyby observations are of insufficient resolution for confident recovery. For a range of assumptions on degree of compensation and crustal and mantle densities, the allowable crustal thickness is consistent with the upper limit of about 100 km estimated from the inferred depth of faulting beneath a prominent lobate scarp, an assumed ductile flow law for crustal material, and the condition that temperature at the base of the crust does not exceed the solidus temperature. The MESSENGER value of C_2,2 has allowed an improved estimate of the ratio of the polar moment of inertia of the mantle and crust to the full polar moment (C_m/C), a refinement that strengthens the conclusion that Mercury has at present a fluid outer core.