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.     

Impact demagnetization of the martian crust

The lack of magnetic anomalies within the giant martian impact basins, Hellas, Argyre, and Isidis suggests that the impacts demagnetized the crust. Our analysis of the magnetic anomaly intensity shows that the interior parts of the basins are completely demagnetized, while the outer parts and surroundings are partially demagnetized. We investigate the shock pressure and impact heating resulting from the impacts. The crust has been completely demagnetized within ∼0.8 basin radius by a combination of thermal and shock effects, and the surroundings have been partially demagnetized by shock to a distance of at least 1.4 radii. We also investigate magnetic signatures of intermediate-size craters. From the pressures generated by both the large and intermediate-sized impacts, we conclude that the remanent magnetization is carried at least in part by high coercivity rocks. Since the crust beneath the basins does not appear to have been remagnetized as it cooled following the impacts, we conclude that the martian core dynamo was inactive or very weak for at least 100 Myr following the Hellas impact.     

Lithospheric drift on early Mars: Evidence in the magnetic field

The crustal magnetic anomalies on Mars may represent hot spot tracks resulting from lithospheric drift on ancient Mars. As evidence, an analysis of lineation patterns derived from the ΔBr magnetic map is presented. The ΔBr map, largely free of external magnetic field effects, allows excellent detail of the magnetic anomaly pattern, particularly in areas of Mars where the field is relatively weak. Using cluster analysis, we show that the elongated anomalies in the martian magnetic field form concentric small circles (parallels of latitude) about two distinct north pole locations. If these pole locations represent ancient spin axes, then tidal force on the early lithosphere by former satellites in retrograde orbits may have pulled the lithosphere in an east-west direction over hot mantle plumes. With an active martian core dynamo, this may have resulted in the observed magnetic anomaly pattern of concentric small circles. As further evidence, we observe that, of the 15 martian giant impact basins that were possibly formed while the core dynamo was active, seven lie along the equators of our two proposed paleopoles. We also find that four other re-magnetized giant impact basins lie along a great circle about the mean magnetic paleopole of Mars. These 11 impact basins, likely the result of fallen retrograde satellite fragments, indicate that Mars once had moons large enough to cause tidal drag on the early martian lithosphere. The results of this study suggest that the magnetic signatures of this tidal interaction have been preserved to the present day.     

Fluid dynamics of local martian magma oceans

The evolution of a melt region produced by a large impact during Mars formation is addressed. While some impact induced melt is redistributed during crater excavation, sufficiently large impacts (much larger than basin forming impacts) generate an intact melt region which is retained beneath the excavation zone, i.e., a local magma ocean. Local magma ocean evolution depends on the effective rheology controlling large scale deformation of the solid part of the planet, mechanism of crystallization, and melt region size. Within the uncertainties of various parameters, two scenarios are possible. For sufficiently weak rheology or large melt region size, evolution is characterized by rapid extrusion and formation of a global magma ocean. For sufficiently strong rheology or small melt region size, in situ crystallization to a partially molten solid state occurs prior to isostatic adjustment. Subsequent to in situ crystallization, local magma ocean evolution depends on melt region size and efficiency of lateral redistribution compared to bulk conductive cooling. For large melt regions, lateral spreading occurs via plastic deformation and results in an asymmetric, global, partial melt layer. For small melt region size, viscous spreading viscous can result in bulk cooling below the solidus prior to formation of a global layer. A hypothesis for the origin of the hemispherical crustal dichotomy and Tharsis rise is suggested. The dichotomy is associated with a global partial melt layer produced by evolution of a large, local magma ocean. After dichotomy formation, evolution of a second, smaller, local magma ocean is related to Tharsis development.     

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.     

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.     

Concentric crater fill in Utopia Planitia: History and interaction between glacial “brain terrain” and periglacial mantle processes

At martian mid-to-high latitudes, the surfaces of potentially ice-rich features, including concentric crater fill, lobate debris aprons, and lineated valley fill, typically display a complex texture known as “brain terrain,” due to its resemblance to the complex patterns on brain surfaces. In order to determine the structure and developmental history of concentric crater fill and overlying latitude-dependent mantle (LDM) material, “brain terrain” and polygonally-patterned LDM surfaces are analyzed using HiRISE images from four craters in Utopia Planitia containing concentric crater fill. “Brain terrain” and mantle surface textures are classified based on morphological characteristics: (1) closed-cell “brain terrain,” (2) open-cell “brain terrain,” (3) high-center mantle polygons, and (4) low-center mantle polygons. A combined glacial and thermal-contraction cracking model is proposed for the formation and modification of the “brain terrain” texture of concentric crater fill. A similar model, related to thermal contraction cracking and differential sublimation of underlying ice, is proposed for the formation and development of polygonally patterned mantle material. Both models require atmospheric deposition of ice, likely during periods of high obliquity, but do not require wet active layer processes. Crater dating of “brain terrain” and mantled surfaces suggests a transition at martian mid-latitudes from peak “glacial” conditions occurring within the past ∼10-100 My to a quiescent period followed by a cold-desert “periglacial” period during the past ∼1-2 My.     

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.     

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.     

Impact megadomes and the origin of the martian crustal dichotomy

We show that a sufficiently energetic impact can generate a melt volume which, after isostatic adjustment and differentiation, forms a spherical cap of crust with underlying depleted mantle. Depending on impact energy and initial crustal thickness, a basin may be retained or impact induced crust may be topographically elevated. Retention of a martian lowland scale impact basin at impact energies ∼3 × 10²⁸-3 × 10²⁹ J requires an initial crustal thickness greater than 10 km. Formation of impact induced crust with size comparable to the martian highlands requires a larger impact energy, ∼1-3 × 10³⁰ J, and initial crustal thickness <20 km. Furthermore, we show that the boundary of impact induced crust can be elliptical due to a spatially asymmetric impact melt volume caused by an oblique impact. We suggest the term “impact megadome” for topographically elevated, impact induced crust and propose that processes involved in megadome formation may play an important role in the origin of the martian crustal dichotomy.     

Chemical composition of the Earth after the giant impact

The giant impact hypothesis for the origin of the Moon has been widely accepted. One of the most important features of this hypothesis is that the impactor's metallic core was incorporated in the Earth after impact. If the mass of the impactor is 0.82 × 10²⁷ g, the mass of the impactor core was estimated to be 0.19 × 10²⁷ g, which is about 1/10 of present Earth's core. Liu (1982) derived the bulk composition of the Earth from CI chondrites, and concluded that the Fe content of his model appears to be low in comparison with the present Earth, which, however, can be rationalized by the addition of impactor core into the proto-Earth developed by Liu (1982). If the impactor's mantle contains 14 wt% FeO as suggested, the mass ratio of impactor/proto-Earth should not exceed 0.22. The same ratio is not likely to exceed 0.30, if a giant blowoff did not occur during impact.     

How rapidly did Mars accrete? Uncertainties in the Hf-W timing of core formation

Estimates for the martian core formation timescale based on the hafnium-tungsten (Hf-W) isotopic system have varied by almost an order of magnitude, because of uncertainties in the martian mantle Hf/W ratio. Here we argue that the Hf/W ratio is ∼4 but is uncertain by ∼25%, resulting in (instantaneous) martian core formation timescales ranging from 0 to 10 Myr; accordingly, Hf-W isotope observations currently have limited utility in distinguishing between scenarios in which Mars formed as a stranded embryo and scenarios in which Mars suffered a prolonged accretion history.     

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.     

Formation of the terrestrial planets

The early phases of formation in the inner solar system were dominated by collisions and short-range dynamical interactions among planetesimals. But the later phases, which account for most of the differences among planets, are unsure because the dynamics are more subtle. Jupiter's influence became more important, leading to drastic clearing out of the asteroid belt and the stunting of Mars's growth. Further in, the effect of Jupiter-- both directly and indirectly, through ejection of mass in the outer solar system-- was probably to speed up the process without greatly affecting the outcome. The great variety in bulk properties of the terrestrial bodies indicate a terminal phase of great collisions, so that the outcome is the result of small-N statistics. Mercury, 65 percent iron, appears to be a residual core from a high-velocity collision. All planets appear to require a late phase of high energy impacts to erode their atmospheres: including the Earth, to remove CO₂ so that its ocean could form by condensation of water.Consistent with this model is that the largest collision, about 0.2 Earth masses, was into the proto-Earth, although the only property that appears to require it is the great lack of iron in the Moon. The other large differences between the Earth and Venus, angular momentum (spin plus satellite) and inert gas abundances, must arise from origin circumstances, but neither require nor forbid the giant impact. Venus's higher ratio of light to heavy inert gases argues for it receiving a large icy impactor, about 10^–6 Earth masses from far out, requiring some improbable dynamics to get a low enough approach velocity. Core formation in both planets probably started rather early during accretion.Some geochemical evidences argue for the Moon coming from the Earth's mantle, but are inconclusive. Large scale melting of the mantle by the giant impact would plausibly have led to stratification. But the "lock-up" at the end of turbulent mantle convection is a trade-off between rates: crystallization of constituents of small density difference versus overall freezing. Also, factors such as differences in melting temperatures and densities, melt compressibilities, and phase transitions may have had homogenizing effects in the subsequent mantle convection.     

Magnetic anomalies near Apollinaris Patera and the Medusae Fossae Formation in Lucus Planum, Mars

The nature of strong martian crustal field sources is investigated by mapping and modeling of Mars Global Surveyor magnetometer data near Apollinaris Patera, a previously proposed volcanic source, supplemented by large-scale correlative studies. Regional mapping yields evidence for positive correlations of orbital anomalies with both Apollinaris Patera and Lucus Planum, a nearby probable extrusive pyroclastic flow deposit that is mapped as part of the Medusae Fossae Formation. Iterative forward modeling of the Apollinaris Patera magnetic anomaly assuming a source model consisting of one or more uniformly magnetized near-surface disks indicates that the source is centered approximately on the construct with a scale size several times larger and comparable to that of the Apollinaris Patera free-air gravity anomaly. A significantly lower rms deviation is obtained using a two-disk model that favors a concentration of magnetization near the construct itself. Estimates for the dipole moment per unit area of the Lucus Planum source together with maximum thicknesses of ∼3 km based on topographic and radar sounding data lead to an estimated minimum magnetization intensity of ∼50 A/m within the pyroclastic deposits. Intensities of this magnitude are similar to those obtained experimentally for Fe-rich Mars analog basalts that cooled in an oxidizing (high fO₂) environment in the presence of a strong (?10 μT) surface field. Further evidence for the need for an oxidizing environment is provided by a broad spatial correlation of the locations of phyllosilicate exposures identified to date using Mars Express OMEGA data with areas containing strong crustal magnetic fields and valley networks in the Noachian-aged southern highlands. This indicates that the presence of liquid water, which is a major crustal oxidant, was an important factor in the formation of strong magnetic sources. The evidence discussed here for magnetic sources associated with relatively young volcanic units suggests that a martian dynamo existed during the late Noachian/early Hesperian, after the last major basin-forming impacts and the formation of the northern lowlands.     

Central magnetic anomalies of Nectarian-aged lunar impact basins: Probable evidence for an early core dynamo

A re-examination of all available low-altitude LP magnetometer data confirms that magnetic anomalies are present in at least four Nectarian-aged lunar basins: Moscoviense, Mendel-Rydberg, Humboldtianum, and Crisium. In three of the four cases, a single main anomaly is present near the basin center while, in the case of Crisium, anomalies are distributed in a semi-circular arc about the basin center. These distributions, together with a lack of other anomalies near the basins, indicate that the sources of the anomalies are genetically associated with the respective basin-forming events. These central basin anomalies are difficult to attribute to shock remanent magnetization of a shocked central uplift and most probably imply thermoremanent magnetization of impact melt rocks in a steady magnetizing field. Iterative forward modeling of the single strongest and most isolated anomaly, the northern Crisium anomaly, yields a paleomagnetic pole position at 81° ± 19°N, 143° ± 31°E, not far from the present rotational pole. Assuming no significant true polar wander since the Crisium impact, this position is consistent with that expected for a core dynamo magnetizing field. Further iterative forward modeling demonstrates that the remaining Crisium anomalies can be approximately simulated assuming a multiple source model with a single magnetization direction equal to that inferred for the northernmost anomaly. This result is most consistent with a steady, large-scale magnetizing field. The inferred mean magnetization intensity within the strongest basin sources is ∼1 A/m assuming a 1-km thickness for the source layer. Future low-altitude orbital and surface magnetometer measurements will more strongly constrain the depth and/or thicknesses of the sources.     

Thermal evolution of the early Moon

The early thermal evolution of Moon has been numerically simulated to understand the magnitude of the impact-induced heating and the initially stored thermal energy of the accreting moonlets. The main objective of the present study was to understand the nature of processes leading to core–mantle differentiation and the production and cooling of the initial convective magma ocean. The accretion of Moon was commenced over a time scale of 100 yr after the giant impact event around 30–100 million years in the early solar system. We studied the dependence of the planetary processes on the impact scenarios, the initial average temperature of the accreting moonlets, and the size of the protomoon that accreted rapidly beyond the Roche limit within the initial 1 yr after the giant impact. The simulations indicate that the accreting moonlets should have a minimum initial averaged temperature around 1600 K. The impacts would provide additional thermal energy. The initial thermal state of the moonlets depends upon the environment prevailing within the Roche limit that experienced episodes of extensive vaporization and recondensation of silicates. The initial convective magma ocean of depth more than 1000 km is produced in the majority of simulations along with the global core–mantle differentiation in case the melt percolation of the molten metal through porous flow from bulk silicates was not the major mode of core–mantle differentiation. The possibility of shallow magma oceans cannot be ruled out in the presence of the porous flow. Our simulations indicate the core–mantle differentiation within the initial 10² to 10³ yr of the Moon accretion. The majority of the convective magma ocean cooled down for crystallization within the initial 10³ to 10⁴ yr.     

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 Manicouagan impact structure as a terrestrial analogue site for lunar and martian planetary science

The 90 km diameter, late Triassic Manicouagan impact structure of Québec, Canada, is a well-preserved, undeformed complex crater possessing an anorthositic central uplift and a 55 km diameter melt sheet. As such, it provides a valuable terrestrial analogue for impact structures developed on other planetary bodies, especially the Moon and Mars, which are currently the focus of exploration initiatives. The scientific value of Manicouagan has recently been enhanced due to the production, between 1994 and 2006, of ∼18 km of drill core from 38 holes by the mineral exploration industry. Three of these holes are in excess of 1.5 km deep, with the deepest reaching 1.8 km. Here we combine recent field work, sampling and the drill core data with previous knowledge to provide insight into processes occurring at Manicouagan and, by inference, within extraterrestrial impact structures. Four areas of comparative planetology are discussed: impact melt sheets, central uplifts, impact-generated hydrothermal regimes and footwall breccias. Human training and instrument testing opportunities are also considered. The drill core reveals that the impact melt and clast-bearing impact melts in the centre of the structure reach thicknesses of 1.4 km. The 1.1 km thick impact melt has undergone differentiation to yield a lower monzodiorite, a transitional quartz monzodiorite and an upper quartz monzonite sequence. This calls into question the previous citing of Manicouagan as an exemplar of a relatively large crater possessing an undifferentiated melt sheet, which was used as a rationale for assigning different composition lunar impact melts and clast-bearing impact melts to separate cratering events. The predominantly anorthositic central uplift at Manicouagan is comparable to certain lunar highlands material, with morphometric analogies to the King, Tycho, Pythagoras, Jackson, and Copernicus impact structures, which have similar diameters and uplift structure. Excellent exposure of the Manicouagan uplift facilitates mapping and an appraisal of its formation and collapse mechanisms, enhanced by drill core data from the centre of the structure. Recent field studies at the edge of the central island at Manicouagan, and multiple drill core sections through footwall lithologies, provide insight into allochthonous (clastic and suevitic) and autochthonous breccia formation, as well as shock effects. The hydrothermal regimes developed at Manicouagan are akin to systems proposed for Noachian (>3.5 Ga) Mars that involve alteration of impact melts via meteoritic and surface waters, with the generation of phyllosilicates, zeolites, hematite, sulfates and sulfides that can contribute to martian soil formation and sedimentation processes.