Asymmetric thermal evolution of the Moon Laneuville, M.; Wieczorek, M. A.; Breuer, D. ...
Journal of geophysical research. Planets,
July 2013, Letnik:
118, Številka:
7
Journal Article
Recenzirano
Odprti dostop
The Moon possesses a clear dichotomy in geological processes between the nearside and farside hemispheres. The most pronounced expressions of this dichotomy are the strong concentration of ...radioactive heat sources on the nearside in a region known as the Procellarum KREEP Terrane (PKT) and the mare basaltic lava flows that erupted in or adjacent to this terrane. We model the thermochemical evolution of the Moon using a 3‐D spherical thermochemical convection code in order to assess the consequences of a layer enriched in heat sources below the PKT on the Moon's global evolution. We find that in addition to localizing most of the melt production on the nearside, such an enriched concentration of heat sources in the PKT crust has an influence down to the core‐mantle boundary and leaves a present‐day temperature anomaly within the nearside mantle. Moderate gravitational and topographic anomalies that are predicted in the PKT, but not observed, may be masked either by crustal thinning or gravitational anomalies from dense material in the underlying mantle. Our models also predict crystallization of an inner core for sulfur concentrations less than 6 wt %.
Key Points
The thermochemical consequences of the PKT are consistent with the observations
A thermal anomaly is present today in the mantle below the PKT
Heat sources enrichment in the PKT has an influence down to the CMB
The present‐day thermal state, interior structure, composition, and rheology of Mars can be constrained by comparing the results of thermal history calculations with geophysical, petrological, and ...geological observations. Using the largest‐to‐date set of 3‐D thermal evolution models, we find that a limited set of models can satisfy all available constraints simultaneously. These models require a core radius strictly larger than 1,800 km, a crust with an average thickness between 48.8 and 87.1 km containing more than half of the planet's bulk abundance of heat producing elements, and a dry mantle rheology. A strong pressure dependence of the viscosity leads to the formation of prominent mantle plumes producing melt underneath Tharsis up to the present time. Heat flow and core size estimates derived from the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission will increase the set of constraining data and help to confine the range of admissible models.
Plain Language Summary
We constrain the thermal state and interior structure of Mars by combining a large number of observations with thermal evolution models. Models that match the available observations require a core radius larger that half the planetary radius and a crust thicker than 48.8 km but thinner than 87.1 km on average. All best‐fit models suggest that more than half of the planet's bulk abundance of heat producing elements is located in the crust. Mantle plumes may still be active today in the interior of Mars and produce partial melt underneath the Tharsis volcanic province. Our results have important implications for the thermal evolution of Mars. Future data from the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission can be used to validate our models and further improve our understanding of the thermal evolution of Mars.
Key Points
We combine the largest‐to‐date set of 3‐D dynamical models with observations to constrain the thermal state and interior structure of Mars
Best‐fit models suggest a core radius strictly larger than 1,800 km and an average crustal thickness 48.8 km < dc < 87.1 km
Models suggest a large pressure dependence of the viscosity and a crust containing 65‐70% of the total amount of heat producing elements
A number of observations performed by the MESSENGER spacecraft can now be employed to better understand the evolution of Mercury's interior. Using recent constraints on interior structure, surface ...composition, volcanic and tectonic histories, we modeled the thermal and magmatic evolution of the planet. We ran a large set of Monte Carlo simulations based on one‐dimensional parametrized models, spanning a wide range of parameters. We complemented these simulations with selected calculations in 2‐D cylindrical and 3‐D spherical geometry, which confirmed the validity of the parametrized approach and allowed us to gain additional insight into the spatiotemporal evolution of mantle convection. Core radii of 1940 km, 2040 km, and 2140 km have been considered, and while in the first two cases several models satisfy the observational constraints, no admissible models were found for a radius of 2140 km. A typical thermal evolution scenario consists of an initial phase of mantle heating accompanied by planetary expansion and the production of a substantial amount of partial melt. The evolution subsequent to 2 Gyr is characterized by secular cooling that proceeds approximately at a constant rate and implies that planetary contraction should be ongoing today. Most of the models predict mantle convection to cease after 3–4 Gyr, indicating that Mercury may be no longer dynamically active. Finally, assuming the observed surface abundance of radiogenic elements to be representative for the entire crust, we determined bulk silicate concentrations of 35–62 ppb Th, 20–36 ppb U, and 290–515 ppm K, similar to those of other terrestrial planets.
Key Points
Models predict initial heating phase followed by cooling and contraction
Stiff rheology, thin crust and present‐day conductive mode are preferred
1D parametrized models agree well with 2D and 3D dynamic simulations
Fitting the isotopic composition of lunar rocks to a new thermal and crystallization model shows the Moon formed 4.40 to 4.45 Ga.
A giant impact onto Earth led to the formation of the Moon, resulted ...in a lunar magma ocean (LMO), and initiated the last event of core segregation on Earth. However, the timing and temporal link of these events remain uncertain. Here, we demonstrate that the low thermal conductivity of the lunar crust combined with heat extraction by partial melting of deep cumulates undergoing convection results in an LMO solidification time scale of 150 to 200 million years. Combining this result with a crystallization model of the LMO and with the ages and isotopic compositions of lunar samples indicates that the Moon formed 4.425 ± 0.025 billion years ago. This age is in remarkable agreement with the U-Pb age of Earth, demonstrating that the U-Pb age dates the final segregation of Earth’s core.
Solar radiation and geological processes over the first few million years of Earth’s history, followed soon thereafter by the origin of life, steered our planet towards an evolutionary trajectory of ...long-lived habitability that ultimately enabled the emergence of complex life. We review the most important conditions and feedbacks over the first 2 billion years of this trajectory, which perhaps represent the best analogue for other habitable worlds in the galaxy. Crucial aspects included: (1) the redox state and volatile content of Earth’s building blocks, which determined the longevity of the magma ocean and its ability to degas H
2
O and other greenhouse gases, in particular CO
2
, allowing the condensation of a water ocean; (2) the chemical properties of the resulting degassed mantle, including oxygen fugacity, which would have not only affected its physical properties and thus its ability to recycle volatiles and nutrients via plate tectonics, but also contributed to the timescale of atmospheric oxygenation; (3) the emergence of life, in particular the origin of autotrophy, biological N
2
fixation, and oxygenic photosynthesis, which accelerated sluggish abiotic processes of transferring some volatiles back into the lithosphere; (4) strong stellar UV radiation on the early Earth, which may have eroded significant amounts of atmospheric volatiles, depending on atmospheric CO
2
/N
2
ratios and thus impacted the redox state of the mantle as well as the timing of life’s origin; and (5) evidence of strong photochemical effects on Earth’s sulfur cycle, preserved in the form of mass-independent sulfur isotope fractionation, and potentially linked to fractionation in organic carbon isotopes. The early Earth presents itself as an exoplanet analogue that can be explored through the existing rock record, allowing us to identify atmospheric signatures diagnostic of biological metabolisms that may be detectable on other inhabited planets with next-generation telescopes. We conclude that investigating the development of habitable conditions on terrestrial planets, an inherently complex problem, requires multi-disciplinary collaboration and creative solutions.
The impact heat accumulated during the late stage of planetary accretion can melt a significant part or even the entire mantle of a terrestrial body, giving rise to a global magma ocean. The ...subsequent cooling of the interior causes the magma ocean to freeze from the core–mantle boundary (CMB) to the surface due to the steeper slope of the mantle adiabat compared to the slope of the solidus.
Assuming fractional crystallization of the magma ocean, dense cumulates are produced close to the surface, largely due to iron enrichment in the evolving magma ocean liquid. A gravitationally unstable mantle thus forms, which is prone to overturn. We investigate the cumulate overturn and its influence on the thermal evolution of Mars using mantle convection simulations in 2D cylindrical geometry. We present a suite of simulations using different initial conditions and a strongly temperature-dependent viscosity. We assume that all radiogenic heat sources have been enriched during the freezing-phase of the magma ocean in the uppermost 50 km and that the initial steam-atmosphere created by the degassing of the freezing magma ocean was rapidly lost, implying that the surface temperature is set to present-day values. In this case, a stagnant lid quickly forms on top of the convective interior preventing the uppermost dense cumulates to sink, even when allowing for a plastic yielding mechanism. Below this dense stagnant lid, the mantle chemical gradient settles to a stable configuration. The convection pattern is dominated by small-scale structures, which are difficult to reconcile with the large-scale volcanic features observed over Mars' surface and partial melting ceases in less than 900 Ma. Assuming that the stagnant lid can break because of additional mechanisms and allowing the uppermost dense layer to overturn, a stable density gradient is obtained, with the densest material and the entire amount of heat sources lying above the CMB. This stratification leads to a strong overheating of the lowermost mantle, whose temperature increases to values that exceed the liquidus. The iron-rich melt would most likely remain trapped in the lower part of the mantle. The upper mantle in that scenario cools rapidly and only shows partial melting during the first billion year of evolution. Therefore a fractionated global and deep magma ocean is difficult to reconcile with observations. Different scenarios assuming, for instance, a hemispherical or shallow magma ocean, or a crystallization sequence resulting in a lower density gradient than that implied by pure fractional crystallization will have to be considered.
•Numerical modeling of magma ocean overturn after fractional crystallization.•Mantle overturn strongly depends on initial density distribution and assumed rheology.•Overturn occurs most likely below a stagnant lid.•Fractionally crystallized magma ocean can hardly explain Mars' thermo-chemical history.
Context. Plate tectonics is considered a fundamental component for the habitability of the Earth. Yet whether it is a recurrent feature of terrestrial bodies orbiting other stars or unique to the ...Earth is unknown. The stagnant lid may rather be the most common tectonic expression on such bodies. Aims. To understand whether a stagnant-lid planet can be habitable, i.e. host liquid water at its surface, we model the thermal evolution of the mantle, volcanic outgassing of H2O and CO2, and resulting climate of an Earth-like planet lacking plate tectonics. Methods. We used a 1D model of parameterized convection to simulate the evolution of melt generation and the build-up of an atmosphere of H2O and CO2 over 4.5 Gyr. We then employed a 1D radiative-convective atmosphere model to calculate the global mean atmospheric temperature and the boundaries of the habitable zone (HZ). Results. The evolution of the interior is characterized by the initial production of a large amount of partial melt accompanied by a rapid outgassing of H2O and CO2. The maximal partial pressure of H2O is limited to a few tens of bars by the high solubility of water in basaltic melts. The low solubility of CO2 instead causes most of the carbon to be outgassed, with partial pressures that vary from 1 bar or less if reducing conditions are assumed for the mantle to 100–200 bar for oxidizing conditions. At 1 au, the obtained temperatures generally allow for liquid water on the surface nearly over the entire evolution. While the outer edge of the HZ is mostly influenced by the amount of outgassed CO2, the inner edge presents a more complex behaviour that is dependent on the partial pressures of both gases. Conclusions. At 1 au, the stagnant-lid planet considered would be regarded as habitable. The width of the HZ at the end of the evolution, albeit influenced by the amount of outgassed CO2, can vary in a non-monotonic way depending on the extent of the outgassed H2O reservoir. Our results suggest that stagnant-lid planets can be habitable over geological timescales and that joint modelling of interior evolution, volcanic outgassing, and accompanying climate is necessary to robustly characterize planetary habitability.
The Interior Exploration using Seismic Investigations, Geodesy and Heat Transport mission, to be launched in 2018, will perform a comprehensive geophysical investigation of Mars in situ. The Seismic ...Experiment for Interior Structure package aims to detect global and regional seismic events and in turn offer constraints on core size, crustal thickness, and core, mantle, and crustal composition. In this study, we estimate the present‐day amount and distribution of seismicity using 3‐D numerical thermal evolution models of Mars, taking into account contributions from convective stresses as well as from stresses associated with cooling and planetary contraction. Defining the seismogenic lithosphere by an isotherm and assuming two end‐member cases of 573 K and the 1073 K, we determine the seismogenic lithosphere thickness. Assuming a seismic efficiency between 0.025 and 1, this thickness is used to estimate the total annual seismic moment budget, and our models show values between 5.7 × 1016 and 3.9 × 1019 Nm.
Key Points
We compute the cumulative seismic moment due to convective and cooling stresses from 3‐D thermal evolution models of Mars
The annual seismic moment budget calculated from our models is between 5.7 × 1016 and 3.9 × 1019 Nm
A new, self‐consistent model for the spatial distribution of seismicity is derived from 3‐D thermal evolution models
The evolution of Earth's early atmosphere and the emergence of habitable conditions on our planet are intricately coupled with the development and duration of the magma ocean (MO) phase during the ...early Hadean period (4-4.5 Ga). In this paper, we study the evolution of the steam atmosphere during the MO period. We obtain the outgoing longwave radiation (OLR) using the line-by-line radiative transfer code GARLIC. Our study suggests that an atmosphere consisting of pure H2O, built as a result of outgassing, extends the MO lifetime to several million years. The thermal emission as a function of the solidification timescale of an MO is shown. We study the effect of thermal dissociation of H2O at higher temperatures by applying atmospheric chemical equilibrium, which results in the formation of H2 and O2 during the early phase of the MO. A 1%-6% reduction in the OLR is seen. We also obtain the effective height of the atmosphere by calculating the transmission spectra for the whole duration of the MO. An atmosphere of depth 100 km is seen for pure water atmospheres. The effect of thermal dissociation on the effective height of the atmosphere is also shown. Due to the difference in the absorption behavior at different altitudes, the spectral features of H2 and O2 are seen at different altitudes of the atmosphere. Therefore, these species, along with H2O, have a significant contribution to the transmission spectra and could be useful for placing observational constraints on MO exoplanets.
The atmospheric noble gas isotope and elemental bulk ratios on Venus and Earth provide important information on their origin and evolution. If the protoplanets grew to a certain mass (i.e. ...> 0.5 MEarth), they could have captured H2-dominated primordial atmospheres by accreting gas from the circumstellar disk during the formation of the Solar System, which were then quickly lost by hydrodynamic escape after the disk dissipated. In such a case, the EUV-driven hydrodynamic flow of H atoms dragged heavier elements with it at different rates, leading to changes in their initial isotope ratios. For reproducing Earth and Venus present atmospheric 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne, isotope and bulk K/U ratios we applied hydrodynamic upper atmosphere escape and Smooth Particle Hydrodynamics (SPH) impact models for the loss calculations of captured H2-dominated primordial atmospheres for various protoplanetary masses. We investigated a wide range of possible EUV evolution tracks of the young Sun and initial atmospheric compositions based on mixtures of captured nebula gas, outgassed and delivered material from ureilite, enstatite and carbonaceaous chondrites. Depending on the disk lifetime of ≈ 3–5 Myr (Bollard et al., 2017; Wang et al., 2017) and the composition of accreted material after disk dispersal, we find from the reproduction of the present atmospheric Ar, Ne, and bulk K/U ratios, that early Earth's evolution can be explained if proto-Earth had accreted masses between ≈ 0.53 − 0.58 MEarth by the time the nebula gas dissipated. If proto-Earth would have accreted a higher mass during the disk lifetime the present atmospheric Ar and Ne isotope ratios can not be reproduced with our model approach. For masses > 0.75 MEarth, Earth would have had a problem to get get rid of its primordial atmosphere. If proto-Earth accreted ≈ 0.53 − 0.58MEarth of enstatite-dominated material as suggested by Dauphas (2017) during the disk lifetime, it would have captured a tiny primordial atmosphere that was lost ≈ 3 Myr after the disk dissipated. In such a case we find that the present-day atmospheric Ar and Ne isotope ratios can be best reproduced if the post-nebula impactors contained ≈ 5% weakly depleted carbonaceous chondritic material and ≈ 95% enstatite chondrites that are strongly depleted in Ar, Ne and moderately volatile elements like potassium. If higher amounts of carbonaceous chondrites were involved in early Earth's accretion as recently suggested by Schiller et al. (2018), then the Earth's present atmospheric Ar and Ne ratios can only be reproduced if the involved carbonaceous chondritic post-nebula material was also highly depleted in these noble gases and/or had to be partially be delivered as long as the primordial atmosphere was yet escaping. As long as primordial atmospheres surround the growing protoplanets the abundance of their volatile elements is overwritten by their respective captured solar-like atmospheric abundances. Therefore the initial composition of the protoplanets at the disk dispersal time can not be identified by our method. For masses < 0.5 MEarth atmospheric escape cannot explain the present-day ratios, i.e. if Earth grew slower then these ratios have to be explained differently (Marty, 2012). If proto-Venus captured a primordial atmosphere it should have grown to masses of ≈ 0.85− 1.0 MVenus during the time until the disk dissipated and if early Venus accreted its main mass during the disk lifetime than the present atmospheric Ar and Ne isotope ratios and the observed K/U ratios on Venus surface can also be reproduced by the escape of a captured primordial atmosphere that is lost within ≤ 100 Myr, if the Sun was born as a weakly to moderately active young G star. New precise re-measurements of atmospheric noble gases are necessary by future Venus missions to better constrain the material that was involved in the planet's accretion history and possibly also the EUV activity evolution of the young Sun. In addition, measurements of other moderately volatile element and isotope ratios on the surface such as Rb/U, 64Zn/66Zn, and 39K/41K can give an insight on whether Venus accreted slow or fast, i.e. almost to its final mass within the disk lifetime.
•Measured atmospheric 36Ar/38Ar and 20Ne/22Ne noble gas isotope and bulk K/U elemental ratios on Venus and Earth can be reproduced by a combination of EUV-driven hydrodynamic escape of H atoms from an initial nebular-based H2-envelope and impacts, if•the young Sun originated as a slowly to moderately rotating early G-type star,•according to a disk lifetime of 4 Myr, proto-Venus and Earth were released from the nebula with masses of ≈0.85-1.0MVenus and ≈0.53-0.58MEarth, and small H2-envelopes,•proto-Earth accreted most of its mass (≈0.7-0.85MEarth) as long as it was surrounded by nebula gas within ≤7 Myr after the beginning of the Solar System.•a more massive solar-like H2-dominated atmosphere surrounded early Venus compared to the less massive proto-Earth, for which carbonaceous chondrites modified the evolution of the noble gas ratios to todayŠs measured atmospheric values, mainly reflected in the 20Ne/22Ne and 36Ar/22Ne ratios.•The late stage impactors occurring after the captured H2-envelope is lost have to be depleted in noble gases and moderately volatile elements like K. In addition, EarthŠs present-day atmospheric 36Ar/22Ne ratio can determine the composition of the material that was accreted after the primordial atmosphere was lost in more detail. This ratio can only be reproduced if the post-H2-envelope impactors were an admixture of 95% strongly depleted ECs and 5% CCs or ? 70% strongly depleted CCs and ? 30% depleted NCC-like material.