Water undersaturated mantle plume volcanism on present‐day Mars

Based on meteorite evidence, the present‐day Martian mantle has a combined abundance of up to a few hundred ppm of H₂O, Cl, and F, which lowers the solidus and enhances the magma production rate. Adiabatic decompression melting in upwelling mantle plumes is the best explanation for young (last 200 Myr) volcanism on Mars. We explore water undersaturated mantle plume volcanism using a finite element mantle convection model coupled to a model of hydrous peridotite melting. Relative to a dry mantle, the reduction in solidus temperature due to water increases the magma production rate by a factor of 1.3–1.7 at 50 ppm water and by a factor of 1.9–3.2 at 200 ppm water. Mantle water also decreases the viscosity and increases the vigor of convection, which indirectly increases the magma production rate by thinning the thermal boundary layer and increasing the flow velocity. At conditions relevant to Mars, these indirect effects can cause an order of magnitude increase in the magma production rate. Using geologic and geophysical observations of the Late Amazonian magma production rate and geochemical observations of melt fractions in shergottite meteorites, present‐day Mars is constrained to have a core–mantle boundary temperature of ~1750 to 1800 °C and a volume‐averaged thermal Rayleigh number of 2 × 10⁶ to 10⁷, indicating that moderately vigorous mantle convection has persisted to the present day. Melting occurs at depths of 2.5–6 GPa and is controlled by the Rayleigh number at the low pressure end and by the mantle water concentration at high pressure.     

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.     

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.     

Crustal recycling, mantle dehydration, and the thermal evolution of Mars

We have reinvestigated the coupled thermal and crustal evolution of Mars taking new laboratory data concerning the flow behavior of iron-rich olivine into account. The low mantle viscosities associated with the relatively higher iron content of the martian mantle as well as the observed high concentrations of heat producing elements in a crust with a reduced thermal conductivity were found to promote phases of crustal recycling in many models. As crustal recycling is incompatible with an early separation of geochemical reservoirs, models were required to show no episodes of crustal recycling. Furthermore, admissible models were required to reproduce the martian crust formation history, to allow for the formation of partial melt under present day mantle conditions and to reproduce the measured concentrations of potassium and thorium on the martian surface. Taking dehydration stiffening of the mantle viscosity by the extraction of water from the mantle into account, we found that admissible models have low initial upper mantle temperatures around 1650 K, preferably a primordial crustal thickness of 30 km, and an initially wet mantle rheology. The crust formation process on Mars would then be driven by the extraction of a primordial crust after core formation, cooling the mantle to temperatures close to the peridotite solidus. According to this scenario, the second stage of global crust formation took place over a more extended period of time, waning at around 3500 Myr b.p., and was driven by heat produced by the decay of radioactive elements. Present-day volcanism would then be driven by mantle plumes originating at the core-mantle boundary under regions of locally thickened, thermally insulating crust. Water extraction from the mantle was found to be relatively efficient and close to 40% of the total inventory was lost from the mantle in most models. Assuming an initial mantle water content of 100 ppm and that 10% of the extracted water is supplied to the surface, this amount is equivalent to a 14 m thick global surface layer, suggesting that volcanic outgassing of H₂O could have significantly influenced the early martian climate and increased the planet’s habitability.     

Mars and its magnetic field— Some new aspects

This paper is the continuation of an earlier one on planetary evolution of Earth and Mars (Franck and Orgzall, 1987). Here models for the history of the Martian magnetic field are presented. It is concluded, that Mars once had a substantial magnetic field that was switched off approximately 1 billion years ago on the one hand due to missing gravitational energy supply— the core reached its eutectic composition earlier— and on the other hand due to dynamic inhibitions.     

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.     

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.     

Melting in the martian mantle: Shergottite formation and implications for present‐day mantle convection on Mars

Abstract— Radiometric age dating of the shergottite meteorites and cratering studies of lava flows in Tharsis and Elysium both demonstrate that volcanic activity has occurred on Mars in the geologically recent past. This implies that adiabatic decompression melting and upwelling convective flow in the mantle remains important on Mars at present. I present a series of numerical simulations of mantle convection and magma generation on Mars. These models test the effects of the total radioactive heating budget and of the partitioning of radioactivity between crust and mantle on the production of magma. In these models, melting is restricted to the heads of hot mantle plumes that rise from the core‐mantle boundary, consistent with the spatially localized distribution of recent volcanism on Mars. For magma production to occur on present‐day Mars, the minimum average radioactive heating rate in the martian mantle is 1.6 times 10^?12 W/kg, which corresponds to 39% of the Wanke and Dreibus (1994) radioactivity abundance. If the mantle heating rate is lower than this, the mean mantle temperature is low, and the mantle plumes experience large amounts of cooling as they rise from the base of the mantle to the surface and are, thus, unable to melt. Models with mantle radioactive heating rates of 1.8 to 2.1 times 10 ^?12 W/kg can satisfy both the present‐day volcanic resurfacing rate on Mars and the typical melt fraction observed in the shergottites. This corresponds to 43–50% of the Wanke and Dreibus radioactivity remaining in the mantle, which is geochemically reasonable for a 50 km thick crust formed by about 10% partial melting. Plausible changes to either the assumed solidus temperature or to the assumed core‐mantle boundary temperature would require a larger amount of mantle radioactivity to permit present‐day magmatism. These heating rates are slightly higher than inferred for the nakhlite source region and significantly higher than inferred from depleted shergottites such as QUE 94201. The geophysical estimate of mantle radioactivity inferred here is a global average value, while values inferred from the martian meteorites are for particular points in the martian mantle. Evidently, the martian mantle has several isotopically distinct compositions, possibly including a radioactively enriched source that has not yet been sampled by the martian meteorites. The minimum mantle heating rate corresponds to a minimum thermal Rayleigh number of 2 times 10⁶, implying that mantle convection remains moderately vigorous on present‐day Mars. The basic convective pattern on Mars appears to have been stable for most of martian history, which has prevented the mantle flow from destroying the isotopic heterogeneity.     

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     

The influence of mantle melting on the evolution of Mars

We present a parameterized convection model of Mars by incorporating a new heat-flow scaling law for stagnant-lid convection, to better understand how the evolution of Mars may be affected by mantle melting. Melting in the mantle during convection leads to the formation of a compositionally buoyant lithosphere, which may also be intrinsically more viscous by dehydration. The consequences of these melting effects on the evolution of terrestrial planets have not been explored before. The temporal evolution of crust and lithospheric mantle is modeled in a self-consistent manner considering mantle melting, convective instability, and the rewetting of dehydrated lithosphere from below by hydrogen diffusion. Though the effect of compositional buoyancy turns out to be minimal, the introduction of viscosity contrast between wet and dry mantle can considerably slow mantle cooling and sometimes lead to non-monotonic core cooling. Furthermore, with or without dehydration stiffening, our model predicts that the martian mantle must have been degassed more extensively (>80%) than previously suggested (<10%); the loss of such a large amount of water from the mantle to surface has significant implications about the role of water in the early surface and climate evolution of Mars.     

Internal structure and properties of Mars

Theoretical physical models of the Martia interior are presented in the light of new and revised data and constraints. These models include thermal evolution, densities, and seismic wave velocities. The interior of Mars appears to be Earth-like in many respects. Although thermal models indicate that Mars has passed its peak of evolution it may still have an asthenosphere and may be moderately active tectonically. Mars has an Fe-FeS core with a radius of and may be moderately active tectonically. Mars has an Fe-FeS core with a radius of 1500–2000 km. The mantle is enriched in FeO with an olivine composition of about Fo₇₅. Theoretically determined seismic wave velocities are relatively well constrained in the mantle with upper-mantle P_n velocities ranging from 7.64 to 7.80 km/sec. However, there are wide variations in V_P in the core dependent on composition. The shadow zone due to the core is larger than the Earth's.     

The thermal evolution of Mars as constrained by paleo-heat flows

Lithospheric strength can be used to estimate the heat flow at the time when a given region was deformed, allowing us to constrain the thermal evolution of a planetary body. In this sense, the high (>300 km) effective elastic thickness of the lithosphere deduced from the very limited deflection caused by the north polar cap of Mars indicates a low surface heat flow for this region at the present time, a finding difficult to reconcile with thermal history models. This has started a debate on the current heat flow of Mars and the implications for the thermal evolution of the planet. Here we perform refined estimates of paleo-heat flow for 22 martian regions of different periods and geological context, derived from the effective elastic thickness of the lithosphere or from faulting depth beneath large thrust faults, by considering regional radioactive element abundances and realistic thermal conductivities for the crust and mantle lithosphere. For the calculations based on the effective elastic thickness of the lithosphere we also consider the respective contributions of crust and mantle lithosphere to the total lithospheric strength. The obtained surface heat flows are in general lower than the equivalent radioactive heat production of Mars at the corresponding times, suggesting a limited contribution from secular cooling to the heat flow during the majority of the history of Mars. This is contrary to the predictions from the majority of thermal history models, but is consistent with evidence suggesting a currently fluid core, limited secular contraction for Mars, and recent extensive volcanism. Moreover, the interior of Mars could even have been heating up during part of the thermal history of the planet.     

Variations of the magnetic field near Mars caused by magnetic crustal anomalies

The possible avenues for photoelectron transport were determined during southern hemisphere winter at Mars by using a mapping analysis of the theoretical magnetic field. Magnetic field line tracing was performed by superposing two magnetic field models: (1) magnetic field derived from a three-dimensional (3D) self-consistent quasi-neutral hybrid model which does not contain the Martian crustal magnetic anomalies and (2) a 3D map of the magnetic field associated with the magnetic anomalies based on Mars Global Surveyor magnetic field measurements. It was found that magnetic field lines connected to the nightside of the planet are mainly channeled within the optical shadow of the magnetotail whereas magnetic field lines connected to the dayside of the planet are observed to form the remainder of the magnetosphere. The simulation suggests that the crustal anomalies create “a magnetic shield” by decreasing the region near Mars which is magnetically connected to the Martian magnetosphere. The rotation of Mars causes periodic changes in magnetic connectivity, but not to qualitative changes in the overall magnetic field draping around Mars.     

Divergent evolution of Earth and Venus: Influence of degassing,tectonics, and magnetic fields

Knowledge of the earliest evolution of Earth and Venus is extremely limited, but it is obvious from their dramatic contrasts today that at some point in their evolution conditions on the two planets diverged. In this paper we develop a geophysical systems box model that simulates the flux of carbon through the mantle, atmosphere, ocean, and seafloor, and the degassing of water from the mantle. Volatile fluxes, including loss to space, are functions of local volatile concentration, degassing efficiency, tectonic plate speed, and magnetic field intensity. Numerical results are presented that demonstrate the equilibration to a steady state carbon cycle, where carbon and water are distributed among mantle, atmosphere, ocean, and crustal reservoirs, similar to present-day Earth. These stable models reach steady state after several hundred million years by maintaining a negative feedback between atmospheric temperature, carbon dioxide weathering, and surface tectonics. At the orbit of Venus, an otherwise similar model evolves to a runaway greenhouse with all volatiles in the atmosphere. The influence of magnetic field intensity on atmospheric escape is demonstrated in Venus models where either a strong magnetic field helps the atmosphere to retain about 60 bars of water vapor after 4.5 Gyr, or the lack of a magnetic field allows for the loss of all atmospheric water to space in about 1 Gyr. The relative influences of plate speed and degassing rate on the weathering rate and greenhouse stability are demonstrated, and a stable to runaway regime diagram is presented. In conclusion, we propose that a stable climate-tectonic-carbon cycle is part of a larger coupled geophysical system where a moderate surface climate provides a stabilizing feedback for maintaining surface tectonics, the thermal cooling of the deep interior, magnetic field generation, and the shielding of the atmosphere over billion year time scales.     

Modeling of major martian magnetic anomalies: Further evidence for polar reorientations during the Noachian

Maps of the vector components of the Mars crustal magnetic field are constructed at the mapping altitude (360 to 410 km) using a selected set of data obtained with the Mars Global Surveyor magnetometer during 2780 orbits of the planet in 1999. Forward modeling calculations are then applied to six relatively strong and isolated, dominantly dipolar, magnetic anomalies for the primary purpose of estimating bulk directions of magnetization. Assuming that the magnetizing field was a (dipolar) core dynamo field centered in the planet, paleomagnetic pole positions are calculated for the six primary source bodies together with that for a seventh anomaly analyzed earlier. In agreement with several previous studies, it is found that six of the seven pole positions are clustered in what is now the northern lowlands in a region centered northwest of Olympus Mons (mean pole position: 34°±10° N, 202°±58° E). Assuming that the dynamo dipole moment vector was approximately parallel to the rotation axis, the modeling results therefore suggest a major reorientation of Mars relative to its rotation axis after magnetization was acquired. Such a reorientation may have been stimulated by internal mass redistributions associated with the formation of the northern lowlands and Tharsis, for example. A comparison of the mean paleo (magnetic) equator to the global distribution of crustal fields shows that magnetic anomalies tend to occur at low paleolatitudes. The same appears to be true for the Noachian-aged valley networks, which exhibit a broad spatial correlation with the magnetic anomalies. A possible interpretation is that the formation of magnetic anomalies and the valley networks was favored in the tropics where melting of water ice and snow was a stronger source of both surface valley erosion and groundwater recharge during the earliest history of the planet. This would be consistent with models in which hydrothermal alteration of crustal rocks played a role in producing the unusually strong martian magnetic anomalies.     

On the thermal history of a moon of fission origin

Model calculations show that the thermal history of a Moon which originated by fission from the proto-Earth is the same as that for the Moon as it is currently understood. In particular, a fissioned Moon currently has a small percent of partial melt or at least near solidus temperatures below depths of 800 km in accord with the seismic data which show that the deep interior of the Moon has a very lowQ. The models have moderate (20–50%) degrees of partial melting in the upper mantle (depths < 300 or 200 km) in the period between 3 to 4 × 10⁹ years ago and, therefore, can account for the mare filling epoch. Finally the heat flow of the models is 18 ergs cm^–2 s^–1 which is close to the average of 19 ergs cm^–2 s^–1 derived from the Apollo heat flow experiments. These findings add further support for the fission origin of the Moon.     

Impact-induced mantle dynamics on Mars

At least 20 impact basins with diameters ranging from 1000 to 3380 km have been identified on Mars, with five exceeding 2500 km. The coincidental timing of the end of the sequence of impacts and the disappearance of the global magnetic field has led to investigations of impact heating crippling an early core dynamo. The rate of core cooling (and thus dynamo activity) is limited by that of the overlying mantle. Thus, the pre-existing thermal state of the mantle controls the extent to which a sequence of impacts may affect dynamo activity. Here, we examine the effects of the initial thermal structure of the core and mantle, and the location of an impact with respect to the pre-existing convective structure on the mantle dynamics and surface heat flux.We find that the impacts that formed the five largest basins dominate the impact-driven effects on mantle dynamics. A single impact of this size can alter the entire flow field of the mantle. Such an impact promotes the formation of an upwelling beneath the impact site, resulting in long-lived single-plume convection. The interval between the largest impacts is shorter than the initial recovery time for a single impact. Hence, the change in convective pattern due to each impact sets up a long term change in the global heat flow. These long-term changes are cumulative, and multiple impacts have a synergistic effect.     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

Theoretical models for Mars and their seismic properties

Theoretical seismic properties of the planet Mars are investigated on the basis of the various models which have been proposed for the internal composition of the planet. The latest interpretation of gravity field data (Reasenberg, 1977), assuming a lower value of the moment of inertia, would require a less dense mantle and a larger core than previous models. If Mars is chondritic in composition, the most reasonable models are an incompletely differentiated H-chondrite or a mixture of H-chondrites and carbonaceous chondrites. Seismic profiles, travel times, and free oscillation periods are computed for various models, with the aim of establishing which seismic data is crucial for deciding among the alternatives. A detailed discussion is given of the seismic properties which could—in principle—help answer the following two questions: Is Mars' core liquid or solid? Does Mars have a partially molten asthenosphere in its upper mantle?     

Dayside induced magnetic field in the ionosphere of Mars

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) onboard the Mars Express spacecraft has occasionally displayed surprising features. One such feature is the occurrence of a series of broadband, low-frequency echoes at equally spaced delay times after the sounder transmitter pulse. The interval between the echoes has been shown to be at the cyclotron period of electrons orbiting in the local magnetic field. The electrons are believed to be accelerated by the large voltages applied to the antenna by the sounder transmitter. Measurements of the period of these “electron cyclotron echoes” provide a simple technique for determining the magnitude of the magnetic field near the spacecraft. These measurements are particularly useful because Mars Express carries no magnetometer, so this is the only method available for measuring the magnetic field magnitude. Using this technique, results are presented showing the large scale structure of the draped field inside the magnetic pile-up boundary. The magnitude of the draped field is shown to vary from about 40 nT at a solar zenith angle of about 25°, to about 25 nT at a solar zenith angle of 90°. The results compare favorably with similar results from the Mars Global Surveyor spacecraft. A fitting technique is developed to derive the vector direction and magnitude of the draped magnetic field in cases where the spacecraft passes through regions with significant variation in the crustal field. The magnetic field directions are consistent with current knowledge of the draping geometry of the magnetic field around Mars.