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.     

Extinct isotope heterogeneities in the mantles of Earth and Mars: Implications for mantle stirring rates

Heterogeneities in terrestrial samples for ¹⁸²W/¹⁸³W and ¹⁴²Nd/¹⁴⁴Nd are only preserved in Hadean and Archean rocks while heterogeneities in ¹²⁹Xe/¹³⁰Xe and ¹³⁶Xe/¹³⁰Xe persist to very young mantle‐derived rocks. In contrast, meteorites from Mars show that the Martian mantle preserves heterogeneities in ¹⁸²W/¹⁸³W and ¹⁴²Nd/¹⁴⁴Nd up to the present. As a consequence of the probable “deep magma ocean” core formation process, we assume that the Earth and Mars both had a very early two‐mantle‐reservoir structure with different initial extinct nuclide isotopic compositions (different ¹⁸²W/¹⁸³W, ¹⁴²Nd/¹⁴⁴Nd, ¹²⁹Xe/¹³⁰Xe, ¹³⁶Xe/¹³⁰Xe ratios). Based on this assumption, we developed a simple stochastic model to trace the evolution of a mantle with two initially distinct layers for the extinct isotopes and its development into a heterogeneous mantle by convective mixing and stretching of these two layers. Using the extinct isotope system ¹⁸²Hf‐¹⁸²W, we find that the mantles of Earth and Mars exhibit substantially different mixing or stirring rates. This is consistent with Mars having cooled faster than the Earth due to its smaller size, resulting in less efficient mantle mixing for Mars. Moreover, the mantle stirring rate obtained for Earth using ¹⁸²Hf‐¹⁸²W is consistent with the mantle stirring rate of ~500 Myr constrained by the long‐lived isotope system, ⁸⁷Rb‐⁸⁷Sr and ¹⁴⁷Sm‐¹⁴³Nd. The apparent absence of ¹⁸²W/¹⁸³W isotopic heterogeneity in modern terrestrial rocks is attributed to very active mantle stirring which reduced the ¹⁸²W/¹⁸³W isotopic heterogeneity to a relatively small scale (~83 m for a mantle stirring rate of 500 Myr) compared to the common sampling scale of terrestrial basalts (~30 or 100 km). Our results also support the “deep magma ocean” core formation model as being applicable to both Mars and Earth.     

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.     

Constraints on the structure of the martian interior determined from the chemical and isotopic systematics of SNC meteorites

Abstract— The crystallization ages of martian (SNC) meteorites give evidence that martian volcanism has continued until recent times‐perhaps until the present. These meteorites also indicate that the mantle source regions of this volcanism are modestly to extremely depleted by terrestrial standards. These 2 observations produce a conundrum. How is it that such depleted source regions have produced basaltic magma for such a long time? This contribution attempts to quantify the radiogenic heat production in 2 distinct martian mantle source regions: those of the shergottites and nakhlites. Compared to the depleted upper mantle of the Earth (MORB), the nakhlite source region is depleted by about a factor of 2, and the shergottite source region is depleted by a factor of 6. According to current geophysical models, the nakhlite source contains the minimum amount of radioactive heat production to sustain whole‐mantle convection and basalt generation over geologic time. A corollary of this conclusion is that the shergottite source contains much too little radioactivity to produce recent (<200 Ma) basalts. A model martian interior with a deep nakhlite mantle that is insulated by a shallow shergottite mantle may allow basalt production from both source regions if the divide between the nakhlite‐shergottite mantles acts as a thermal boundary layer. Similarities between lunar and martian isotopic reservoirs indicate that the Moon and Mars may have experienced similar styles of differentiation.     

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.     

A new view of martian evolution

Conventional evolutionary models for Mars adopt a dry mantle solidus. Taking into account the condensation conditions in the preplanetary nebula in the accretion zone of Mars, it can be concluded that large amounts of water or hydrated silicates have condensed in those regions. Therefore, water influences significantly the melting behaviour and the viscosity of the silicatic material. A model for the calculation of the thermal history of a planet is constructed. On this basis, and use of water — saturated solidus — it is possible to derive that the core is not liquid, as given in models employing a dry mantle solidus, but solid to a large extent, which prevents the operation of a large-scale dynamo and explains in that way the lack of a magnetic field. With these assumptions one can construct a possible evolutionary scheme that covers early crust differentiation, a hot thermal past and the missing magnetic field at present.     

Asteroid Differentiation: Pyroclastic Volcanism to Magma Oceans

Abstract— Asteroid differentiation was driven by a complex array of magmatic processes. This paper summarizes theoretical and somewhat speculative research on the physics of these processes. Partial melts in asteroids migrate rapidly, taking < 10⁶ years to reach surface regions. On relatively small (<100 km) asteroids with sufficient volatiles in partial melts (<3000 ppm), explosive volcanism accelerated melts to greater than escape velocity, explaining the apparent lack of basaltic components on the parent asteroids of some differentiated meteorites. Partial melting products include the melts (some eucrites, angrites), residues (lodranites, ureilites), and unfractionated residues (acapulcoites). The high liquidus temperatures of magmatic iron meteorites, the existence of pallasites with only olivine, and the fact that enstatite achondrites formed from ultramafic magmas argue for the existence of magma oceans on some asteroids. Asteroidal magma oceans would have been turbulently convective. This would have prevented crystals nucleated at the upper cooling surface (the only place for crystal nucleation in a low-pressure body) from settling until the magma became choked with crystals. After turbulent convection slowed, crystals and magma would have segregated, leaving a body stratified from center to surface as follows: a metallic core, a small pallasite zone, a dunite region, a feldspathic pyroxenite, and basaltic intrusions and lava flows (if the basaltic components had not been lost by explosive volcanism). The pallasite and dunite zones probably formed from coarse (0.5–1 cm) residual olivine left after formation of the magma ocean at >50% partial melting of the silicate assemblage. Iron cores crystallized dendritically from the outside to the inside. The rapid melt migration rate of silicate melts suggests that ²⁶Al could not be responsible for forming asteroidal magma oceans because it would leave the interior before a sufficient amount of melting occurred. Other heat sources are more likely candidates. Our analysis suggests that if Earth-forming planetesimals had differentiated they were either small (<100 km) and poor in volatiles (<1000 ppm) or they were rich in volatiles and large enough (>300 km) to retain the products of pyroclastic eruptions; if these conditions were not met, Earth would not have a basaltic component.     

The heat flow during the formation of ribbon terrains on Venus

Ribbons are regularly spaced, between 2 and 6 km, troughs that exist on venusian tesserae, which are mainly located in, and characterize to, venusian crustal plateaus. Independent of the geological or temporal relations with other features, regularly and similarly spaced ribbons on several tesserae strongly suggest a thermal control on the thickness of the deformed layer. This can be used to constraint the heat flow at the time of ribbon formation, which holds important implications for the viability of the hypotheses that address the origin and evolution of crustal plateaus. For a brittle–ductile transition ∼1–3 km deep (as proposed from ribbon spacing), realistic strain rates, and a present-day surface temperature of 740 K, the implied heat flow is very high, 130–780 mW m⁻². If Venus has experienced higher surface temperatures due to climate forcing by massive volcanism, then the heat flow could be greatly reduced. For surface temperatures of 850–900 K the heat flow is 190–560, 60–230 and 20–130 mW m⁻² for brittle–ductile transition depths of 1, 2 and 3 km, respectively. Heat flow values around 80–100 mW m⁻² are reasonable for venusian hotspots, based on terrestrial analogs, but hardly consistent with coldspot settings. High surface temperatures are also required to maintain the crustal solidus deeper than a few kilometers during the formation of ribbon terrains. For the obtained heat flows, a solidus deeper than ∼30 km (the likely mean value for the crustal thickness) is difficult to achieve. This suggests that a substantial proportion of the crust beneath crustal plateaus was emplaced subsequently to the time when ribbon terrains were formed. Alternatively, at that time a magma reservoir inside the crust could have existed.     

Mars Hesperian Magmatism as Revealed by Syrtis Major and the Circum-Hellas Volcanic Province

Review of morphologic, morphometric and compositional data from Mars suggests that volcanism in the early Hesperian Syrtis Major edifice was predominantly ultramafic, in contrast to the abundant basaltic volcanism of the Hesperian to Amazonian Tharsis and Elysium provinces. Comparisons of edifice characteristics between Syrtis Major and the large, circum-Hellas Noachian to Hesperian volcanoes suggest that these structures may also be formed by ultramafic volcanic activity. The data suggest that a global scale magma compositional change occurred on Mars during the late Hesperian. The occurrence of widespread ultramafic volcanism suggests that a high degree of partial melting in a relatively hot mantle characterized Mars?? early thermal history, conditions that may be analogous to those that prevailed in the Archean Earth.     

The composition and thickness of the crust of Mars estimated from rare earth elements and neodymium-isotopic compositions of Martian meteorites

Abstract— Isotopic and trace element compositions of Martian meteorites show that early differentiation of Mars produced complementary crustal and mantle reservoirs that were sampled by later magmatic events. This paper describes a mass balance model that estimates the rare earth element (REE) content and thickness of the crust of Mars from the compositions of shergottites. The diverse REE and Nd isotopic compositions of shergottites are most easily explained by variable addition of light rare earth element (LREE)–enriched crust to basaltic magmas derived from LREE-depleted mantle source regions. Antarctic shergottites EET 79001, ALH 77005, LEW 88516, and QUE 94201 all have strongly LREE-depleted patterns and positive initial ¹⁴³Nd isotopic compositions, which is consistent with the generation of these magmas from depleted mantle sources and little or no interaction with enriched crust. In contrast, Shergotty and Zagami have negative initial ¹⁴³Nd isotopic compositions and less pronounced depletions of the LREE, which have been explained by incorporation of enriched crustal components into mantle-derived magmas (Jones, 1989; Longhi, 1991; Borg et al., 1997). The mass balance model presented here derives the REE composition of the crustal component in Shergotty by assuming it represents a mixture between a mantle-derived magma similar in composition to EET 79001A and a LREE-enriched crustal component. The amount of crust in Shergotty is constrained by mixing relations based on Nd-isotopic compositions, which allows the REE pattern of the crustal component to be calculated by mass balance. The effectiveness of this model is demonstrated by the successful recovery of important characteristics of the Earth's continental crust from terrestrial Columbia River basalts. Self-consistent results for Nd-isotopic compositions and REE abundances are obtained if Shergotty contains ～10–30% of LREE-enriched crust with >10 ppm Nd. This crustal component would have moderately enriched LREE (Sm/Nd = 0.25–0.27; ¹⁴⁷Sm/¹⁴⁴Nd = 0.15–0.17; La/Yb = 2.7–3.8), relatively unfractionated heavy rare earth elements (HREE), and no Eu anomaly. Crust with these characteristics can be produced from a primitive lherzolitic Martian mantle by modest amounts (2–8%) of partial melting, and it would have a globally averaged thickness of <45 km, which is consistent with geophysical estimates. Mars may serve as a laboratory to investigate planetary differentiation by extraction of a primary basaltic crust.     

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.     

Heterogeneous distribution of H2O in the Martian interior: Implications for the abundance of H2O in depleted and enriched mantle sources

We conducted a petrologic study of apatite within 12 Martian meteorites, including 11 shergottites and one basaltic regolith breccia. These data were combined with previously published data to gain a better understanding of the abundance and distribution of volatiles in the Martian interior. Apatites in individual Martian meteorites span a wide range of compositions, indicating they did not form by equilibrium crystallization. In fact, the intrasample variation in apatite is best described by either fractional crystallization or crustal contamination with a Cl‐rich crustal component. We determined that most Martian meteorites investigated here have been affected by crustal contamination and hence cannot be used to estimate volatile abundances of the Martian mantle. Using the subset of samples that did not exhibit crustal contamination, we determined that the enriched shergottite source has 36–73 ppm H₂O and the depleted source has 14–23 ppm H₂O. This result is consistent with other observed geochemical differences between enriched and depleted shergottites and supports the idea that there are at least two geochemically distinct reservoirs in the Martian mantle. We also estimated the H₂O, Cl, and F content of the Martian crust using known crust‐mantle distributions for incompatible lithophile elements. We determined that the bulk Martian crust has ~1410 ppm H₂O, 450 ppm Cl, and 106 ppm F, and Cl and H₂O are preferentially distributed toward the Martian surface. The estimate of crustal H₂O results in a global equivalent surface layer (GEL) of ~229 m, which can account for at least some of the surface features on Mars attributed to flowing water and may be sufficient to support the past presence of a shallow sea on Mars' surface.     

Ancient geologic events on Mars revealed by zircons and apatites from the Martian regolith breccia NWA 7034

Zircons and apatites in clasts and matrix from the Martian breccia NWA 7034 are well documented, timing ancient geologic events on Mars. Furthermore, in this study, zircon trace elemental content, apatite volatile content, and apatite volatile isotopic compositions measured in situ could constrain the evolution of those geologic events. The U‐Pb dates of zircons in basalt, basaltic andesite, trachyandesite igneous clasts, and the matrix are similar (4.4 Ga) suggesting intense volcanism on ancient Mars. However, two metamict zircon grains found in the matrix have an upper intercept date of ~4465 Ma in crystalline, whereas amorphous areas have a lower intercept date of 1634 ± 93 Ma. The younger date is consistent with the date of apatites (1530 ± 65 Ma), suggesting a metamorphic event that completely reset the U‐Pb system in both the amorphous areas of zircon and all apatites. δD values in all apatites negatively correlate with water content in a two‐endmember mixing trend. The D (δD up to 2459‰) and ³⁷Cl heavy core (3.8‰) of a large apatite grain suggest a D‐, ³⁷Cl‐rich fluid during the metamorphic event ~1.6 Ga ago, consistent with the trace elements Y, Hf and Ti and P in zircons. The fluid was also therefore P‐rich. The D‐, ³⁷Cl‐poor H₂O‐rich rim (<313‰) suggests the degassing of water from the Martian Cl‐poor interior at a later time. This D‐, ³⁷Cl‐poor Martian mantle reservoir could have derived from volcanic intrusions postdating the younger metamorphic event recorded in NWA 7034.     

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.     

Moderately and slightly siderophile element constraints on the depth and extent of melting in early Mars

The thermal history of Mars during accretion and differentiation is important for understanding some fundamental aspects of its evolution such as crust formation, mantle geochemistry, chronology, volatile loss and interior degassing, and atmospheric development. In light of data from new Martian meteorites and exploration rovers, we have made a new estimate of Martian mantle siderophile element depletions. New high pressure and temperature metal–silicate experimental partitioning data and expressions are also available. Using these new constraints, we consider the conditions under which the Martian mantle may have equilibrated with metallic liquid. The resulting conditions that best satisfy six siderophile elements—Ni, Co, W, Mo, P, and Ga—and are consistent with the solidus and liquidus of the Martian mantle phase diagram are a pressure of 14 ± 3 GPa and temperature of 2100 ± 200 K. The Martian mantle depletions of Cr and V are also consistent with metal–silicate equilibration in this pressure and temperature range if deep mantle silicate phases are also taken into account. The results are not consistent with either metal–silicate equilibrium at the surface or at the current‐day Martian core–mantle boundary. Recent measurements and modeling have concluded that deep (～17 GPa or 1350 km) mantle melting is required to explain isotopic data for Martian meteorites and the nature of differentiation into core, mantle, and crust. This is in general agreement with our estimates of the conditions of Martian core formation based on siderophile elements that result in an intermediate depth magma ocean scenario for metal–silicate equilibrium.     

Trace elements in olivine and the petrogenesis of the intermediate,olivine‐phyric shergottite NWA 10170

Olivine‐phyric shergottites represent primitive basaltic to picritic rocks, spanning a large range of Mg# and olivine abundances. As primitive olivine‐bearing magmas are commonly representative of their mantle source on Earth, understanding the petrology and evolution of olivine‐phyric shergottites is critical in our understanding of Martian mantle compositions. We present data for the olivine‐phyric shergottite Northwest Africa (NWA) 10170 to constrain the petrology with specific implications for magma plumbing‐system dynamics. The calculated oxygen fugacity and bulk‐rock REE concentrations (based on modal abundance) are consistent with a geochemically intermediate classification for NWA 10170, and overall similarity with NWA 6234. In addition, we present trace element data using laser ablation ICP‐MS for coarse‐grained olivine cores, and compare these data with terrestrial and Martian data sets. The olivines in NWA 10170 contain cores with compositions of Fo₇₇ that evolve to rims with composition of Fo₅₈, and are characterized by cores with low Ni contents (400–600 ppm). Nickel is compatible in olivine and such low Ni content for olivine cores in NWA 10170 suggests either early‐stage fractionation and loss of olivine from the magma in a staging chamber at depth, or that Martian magmas have lower Ni than terrestrial magmas. We suggest that both are true in this case. Therefore, the magma does not represent a primary mantle melt, but rather has undergone 10–15% fractionation in a staging chamber prior to extrusion/intrusion at the surface of Mars. This further implies that careful evaluation of not only the Mg# but also the trace element concentrations of olivine needs to be conducted to evaluate pristine mantle melts versus those that have fractionated olivine (±pyroxene and oxide minerals) in staging chambers.     

Constraints on mantle plumes on Venus: Implications for volatile history

The analysis of Venus’ gravity field and topography suggests the presence of a small number of deep mantle plumes (～9). This study predicts the number of plumes formed at the core–mantle boundary, their characteristics, and the production of partial melt from adiabatic decompression. Numerical simulations are performed using a 3D spherical code that includes large viscosity variations and internal heating. This study investigates the effect of several parameters including the core–mantle boundary temperature, the amount of internal heating, and the mantle viscosity. The smallest number of plumes is achieved when no internal heating is present. However, scaling Earth’s radiogenic heating to Venus suggests a value of ～16 TW. Cases with internal heating produce more realistic lid thickness and partial melting, but produce either too many plumes or no plumes if a high mantle temperature precludes the formation of a hot thermal boundary layer. Mantle viscosity must be reduced to at least 10²⁰ Pa s in order to include significant internal heating and still produce hot plumes. In all cases that predict melting, melting occurs throughout the upper mantle. Only cases with high core temperature (>1700 K) produce dry melting. Over time the upper mantle may have lost significant volatiles. Depending on the water content of the lower mantle, deep plumes may contribute to present-day atmospheric water via volcanic outgassing. Assuming 50 ppm water in mantle, 10 plumes with a buoyancy flux of 500 kg/s continuously erupting for 4 myr will outgas an amount of water on the order of that in the lower atmosphere. A higher level of internal heating than achieved to date, as well as relatively low mantle viscosity, may be required to achieve simulations with ～10 plumes and a thinner lid. Alternatively, if the mantle is heating up due to the stagnant lid, the effect is equivalent to having lower rates of internal heating. A temperature increase of 110 K/byr is equivalent to ?13 TW. This value along with the internal heating of 3 TW used in this study may represent the approximate heat budget of Venus’ mantle.     

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.     

Mars’ atmospheric 40Ar: A tracer for past crustal erosion

Noble gas ⁴⁰Ar may be used as a tracer of the past evolution of volatiles in Mars’ crust, mantle and atmosphere. ⁴⁰Ar is formed by the radioactive decay of ⁴⁰K in the mantle and in the crust and is released from the mantle to the atmosphere due to volcanism and from the crust by erosion such as eolian and hydrothermal erosion. Furthermore, ⁴⁰Ar can escape from the atmosphere into space via atmospheric escape mechanisms. The evolution of the atmospheric abundance of ⁴⁰Ar thus depends on these three processes whose efficiencies vary with time.In the present study we reconsider atmospheric escape mechanism efficiencies and describe various possible scenarios of the evolution of ⁴⁰Ar with a model describing the three main reservoirs of ⁴⁰Ar, the mantle, crust and atmosphere. First, we show that atmospheric escape, which is stronger in the early evolution, does not significantly influence the present abundance of the atmospheric ⁴⁰Ar. In the early evolution the atmospheric concentration of ⁴⁰Ar is very low as the outgassing of ⁴⁰Ar from the mantle occurs relatively late in the martian evolution. Thus, the atmospheric ⁴⁰Ar concentration is essentially a tracer of Mars’ outgassing history and not of the escape processes. Second, using the results of the most recent published crustal formation models, the calculated present ⁴⁰Ar atmospheric abundance is smaller than its observed value. This discrepancy may be explained by a significant ⁴⁰Ar supply from the crust by erosion (16–30% of the ⁴⁰Ar content of the upper first 10 km of crust). The knowledge of the fraction of crustal ⁴⁰Ar outgassed to the atmosphere is an important constraint for any future global modelling of past Mars’ hydrothermal activity aiming at better characterizing the role of subsurface aqueous alteration processes in Mars climate evolution. One of the main sources of the uncertainty of these results is the present uncertainty in the measured atmospheric ⁴⁰Ar value (±20%). More precise measurements of ⁴⁰Ar and ³⁶Ar in the martian atmosphere are therefore required to better constrain the model.     

Petrology and mineralogy of volcanic glass in meteorite Northwest Africa 11801: Implications for their petrogenesis

The study of lunar magma evolution holds significant importance within the scientific community due to its relevance in understanding the Moon's thermal and geological history. However, the intricate task of unraveling the history of early volcanic activity on the Moon is hindered by the high flux of impactors, which have substantially changed the morphology of pristine volcanic constructs. In this study, we focus on a unique volcanic glass found in the lunar meteorite Northwest Africa 11801. This kind of volcanic glass is bead-like in shape and compositionally similar to the Apollo-14 and Apollo-17 very low-Ti glass. Our research approach involves conducting a comprehensive analysis of the petrology and mineralogy of the volcanic glass, coupled with multiple thermodynamic modeling techniques. Through the investigation, we aim to shed light on the petrological characteristics and evolutionary history of the glass. The results indicate that the primitive magma of the glass was created at 1398–1436°C and 8.3–11.9 kbar (166–238 km) from an olivine+orthopyroxene mantle source region. Then, the magma ascended toward the surface along a non-adiabatic path with an ascent rate of ~40 m s⁻¹ or 0.2 MPa s⁻¹. During the magma ascent, only olivine crystallized and the onset of magma eruption occurred at ~1320–1343°C. Finally, the glass cooled rapidly on the lunar surface with a cooling rate ranging between 20 and 200 K min⁻¹. Considerable evidence from petrology, mineralogy, cooling rate, and the eruption rate of the glass beads strongly supports the occurrence of ancient explosive volcanism on the Moon.