Mantle Dynamics in Super-Earths: Post-Perovskite Rheology and Self-Regulation of Viscosity

Simple scalings suggest that super-Earths are more likely than an equivalent Earth-sized planet to be undergoing plate tectonics. Generally, viscosity and thermal conductivity increase with pressure while thermal expansivity decreases, resulting in lower convective vigor in the deep mantle. According to conventional thinking, this might result in no convection in a super-Earth's deep mantle. Here we evaluate this. First, we here extend the density functional theory (DFT) calculations of post-perovskite activation enthalpy of to a pressure of 1 TPa. The activation volume for diffusion creep becomes very low at very high pressure, but nevertheless for the largest super-Earths the viscosity along an adiabat may approach 1030 Pa s in the deep mantle. Second, we use these calculated values in numerical simulations of mantle convection and lithosphere dynamics of planets with up to ten Earth masses. The models assume a compressible mantle including depth-dependence of material properties and plastic yielding induced plate tectonics. Results confirm the likelihood of plate tectonics and show a novel self-regulation of deep mantle temperature. The deep mantle is not adiabatic; instead internal heating raises the temperature until the viscosity is low enough to facilitate convective loss of the radiogenic heat, which results in a super-adiabatic temperature profile and a viscosity increase with depth of no more than ~3 orders of magnitude, regardless of the viscosity increase that is calculated for an adiabat. Convection in large super-Earths is characterised by large upwellings and small, time-dependent downwellings. If a super-Earth was extremely hot/molten after its formation, it is thus likely that even after billions of years its deep interior is still extremely hot and possibly substantially molten with a "super basal magma ocean" - a larger version of (Labrosse et al., 2007).Comment: 25 pages, 5 figure

The dynamical control of subduction parameters on surface topography

he long-wavelength surface deflection of Earth's outermost rocky shell is mainly controlled by large-scale dynamic processes like isostasy or mantle flow. The largest topographic amplitudes are therefore observed at plate boundaries due to the presence of large thermal heterogeneities and strong tectonic forces. Distinct vertical surface deflections are particularly apparent at convergent plate boundaries mostly due to the convergence and asymmetric sinking of the plates. Having a mantle convection model with a free surface that is able to reproduce both realistic single-sided subduction and long-wavelength surface topography self-consistently, we are now able to better investigate this interaction. We separate the topographic signal into distinct features and quantify the individual topographic contribution of several controlling subduction parameters. Results are diagnosed by splitting the topographic signal into isostatic and residual components, and by considering various physical aspects like viscous dissipation during plate bending. Performing several systematic suites of experiments, we are then able to quantify the topographic impact of the buoyancy, rheology, and geometry of the subduction-zone system to each and every topographic feature at a subduction zone and to provide corresponding scaling laws. We identify slab dip and, slightly less importantly, slab buoyancy as the major agents controlling surface topography at subduction zones on Earth. Only the island-arc high and the back-arc depression extent are mainly controlled by plate strength. Overall, his modeling study sets the basis to better constrain deep-seated mantle structures and their physical properties via the observed surface topography on present-day Earth and back through time

Quantifying the correlation between mobile continents and elevated temperatures in the subcontinental mantle

Continents influence the mantle's convective wavelength and the heat flow escaping from the planet's surface. Over the last few decades, many numerical and analytical studies have contributed to the debate about whether the continents can warm up the subcontinental mantle or not and if they do, then to what extent? However, a consensus regarding the exact nature and magnitude of this correlation between continents and elevated temperatures in the subcontinental mantle remains to be achieved. By conducting a systematic parameter study using 2‐D global mantle convection simulations with mobile continents, we provide qualitative and quantitative observations on the nature of this correlation. In our incompressible and compressible convection models, we observe the general processes of downwellings bringing cold material into the mantle along continental margins and a subsequent buildup of warm thermal anomalies underneath the continents. We compute the amplitude and degree of this correlation using spectral decomposition of the temperature and composition fields. The dominant degree of correlation evolves with time and changes with continental configuration. Using simple empirical fits, we observe that this correlation decreases with increasing core temperature, number of continents, internal heating, or decreasing reference viscosity. We also report simple regressions of the time dependence of this correlation. Additionally, we show that decompression melting as a result of a mantle upwelling or small‐scale sublithospheric convection leads to voluminous volcanism. The emplacement of this dense basalt‐eclogite material breaks the continents apart and destroys the correlation

Linking continental drift, plate tectonics and the thermal state of the Earth's mantle

International audienceContinents slowly drift at the top of the mantle, sometimes colliding, splitting and aggregating. The evolutions of the continent configuration, as well as oceanic plate tectonics, are surface expressions of mantle convection and closely linked to the thermal state of the mantle; however, quantitative studies are so far lacking. In the present study we use 3D spherical numerical simulations with self-consistently generated plates and compositionally and rheologically distinct continents floating at the top of the mantle in order to investigate the feedbacks between continental drift, oceanic plate tectonics and the thermal state of the Earth's mantle, by using different continent configurations ranging from one supercontinent to six small continents. With the presence of a supercontinent we find a strong time-dependence of the oceanic surface heat flow and suboceanic mantle temperature, driven by the generation of new plate boundaries. Very large oceanic plates correlate with periods of hot suboceanic mantle, while the mantle below smaller oceanic plates tends to be colder. Temperature fluctuations of subcontinental mantle are significantly smaller than in oceanic regions and are caused by a time-variable efficiency of thermal insulation of the continental convection cell. With the presence of multiple continents the temperature below individual continents is generally lower than below supercontinent and is more time-dependent, with fluctuations as large as 15% that are caused by continental assembly and dispersal. The periods featuring a hot subcontinental mantle correlate with strong clustering of the continents and periods characterized by cold subcontinental mantle, at which it can even be colder than suboceanic mantle, with a more dispersed continent configuration. Our findings with multiple continents imply that periods of partial melting and strong magmatic activity inside the continents, which may contribute to continental rifting and pronounced growth of continental crust, might be episodic processes related to the supercontinent cycle. Finally, we observe an influence of continents on the wavelength of convection: for a given strength of the lithosphere we observe longer-wavelength flow components, when continents are present. This observation is regardless of the number of continents, but most pronounced for a single supercontinent

Seafloor spreading evolution in response to continental growth

International audienc

Modelling Earth's surface topography: Decomposition of the static and dynamic components

Contrasting results on the magnitude of the dynamic component of topography motivate us to analyse the sources of uncertainties affecting long wavelength topography modelling. We obtain a range of mantle density structures from thermo-chemical interpretation of available seismic tomography models. We account for pressure, temperature and compositional effects as inferred by mineral physics to relate seismic velocity with density. Mantle density models are coupled to crustal density distributions obtained with a similar methodology. We compute isostatic topography and associated residual topography maps and perform instantaneous mantle flow modelling to calculate the dynamic topography. We explore the effects of proposed mantle 1-D viscosities and also test a 3D pressure- and temperature-dependent viscosity model. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient (r) of respectively ∼0.74 and ∼0.71, upper-lower quartile ranges of 0.86–0.65 for residual topography and 0.83–0.62 for dynamic topography maps. The amplitudes are, on the contrary, poorly constrained. For the static component, the inferred density models of lithospheric mantle give an interquartile range of isostatic topography that is always higher than 100 m, reaching 1.7 km in some locations, and averaging ∼720 m. Crustal density models satisfying the same compressional velocity structure lead to variations in isostatic topography averaging 350 m, with peaks of 1 km in thick crustal regions, and always higher than 100 m. The uncertainties on isostatic topography are strong enough to mask, if present, the contribution of mantle convection to surface topography. For the dynamic component, we obtain a peak-to-peak dynamic topography amplitude exceeding 3 km for all our mantle density and viscosity models. These extremely high values would be associated with a magnitude of geoid undulations that is not in agreement with observations. Considering chemical heterogeneities in correspondence with the lower mantle Large Low Shear wave Velocity Provinces (LLSVPs) helps to decrease the peak-to-peak amplitudes of dynamic topography and geoid, but significantly reduces the correlation between synthetic and observed geoid. The correlation coefficients between all our residual and dynamic topography maps (a total of 220 and 198, respectively) is <0.55 (average = ∼0.19). The correlation slightly improves when considering only the very long-wavelength components of the maps (average = ∼0.23). We therefore conclude that a robust determination of dynamic topography is not feasible since current uncertainties affecting crustal density, mantle density and mantle viscosity are still too large. A truly interdisciplinary approach, combining constraints from the geological record with a multi-methodological interpretation of geophysical observations, is required to tackle the challenging task of linking the surface topography to deep processes

Statistical cyclicity of the supercontinent cycle

International audienc

Using pattern recognition to infer parameters governing mantle convection

The results of mantle convection simulations are fully determined by the input parameters and boundary conditions used. These input parameters can be for initialisation, such as initial mantle temperature, or can be constant values, such as viscosity exponents. However, knowledge of Earth-like values for many input parameters are very poorly constrained, introducing large uncertainties into the simulation of mantle flow. Convection is highly non-linear, therefore linearised inversion methods cannot be used to recover past configurations over more than very short periods of time, which makes finding both initial and constant simulation input parameters very difficult. In this paper, we demonstrate a new method for making inferences about simulation input parameters from observations of the mantle temperature field after billions of years of convection. The method is fully probabilistic. We use prior sampling to construct probability density functions for convection simulation input parameters, which are represented using neural networks. Assuming smoothness, we need relatively few samples to make inferences, making this approach much more computationally tractable than other probabilistic inversion methods. As a proof of concept, we show that our method can invert the amplitude spectra of temperature fields from 2D convection simulations, to constrain yield stress, surface reference viscosity and the initial thickness of primordial material at the CMB, for our synthetic test cases. The best constrained parameter is yield stress. The reference viscosity and initial thickness of primordial material can also be inferred reasonably well after several billion years of convection