The importance of temporal stress variation and dynamic disequilibrium for the initiation of plate tectonics

We use 1-D thermal history models and 3-D numerical experiments to study the impact of dynamic thermal disequilibrium and large temporal variations of normal and shear stresses on the initiation of plate tectonics. Previous models that explored plate tectonics initiation from a steady state, single plate mode of convection concluded that normal stresses govern the initiation of plate tectonics, which based on our 1-D model leads to plate yielding being more likely with increasing interior heat and planet mass for a depth-dependent Byerlee yield stress. Using 3-D spherical shell mantle convection models in an episodic regime allows us to explore larger temporal stress variations than can be addressed by considering plate failure from a steady state stagnant lid configuration. The episodic models show that an increase in convective mantle shear stress at the lithospheric base initiates plate failure, which leads with our 1-D model to plate yielding being less likely with increasing interior heat and planet mass. In this out-of-equilibrium and strongly time-dependent stress scenario, the onset of lithospheric overturn events cannot be explained by boundary layer thickening and normal stresses alone. Our results indicate that in order to understand the initiation of plate tectonics, one should consider the temporal variation of stresses and dynamic disequilibrium

Hemispheric Tectonics on super-Earth LHS 3844b

The tectonic regime of rocky planets fundamentally influences their long-term evolution and cycling of volatiles between interior and atmosphere. Earth is the only known planet with active plate tectonics, but observations of exoplanets may deliver insights into the diversity of tectonic regimes beyond the solar system. Observations of the thermal phase curve of super-Earth LHS 3844b reveal a solid surface and lack of a substantial atmosphere, with a temperature contrast between the substellar and antistellar point of around 1000 K. Here, we use these constraints on the planet's surface to constrain the interior dynamics and tectonic regimes of LHS 3844b using numerical models of interior flow. We investigate the style of interior convection by assessing how upwellings and downwellings are organized and how tectonic regimes manifest. We discover three viable convective regimes with a mobile surface: (1) spatially uniform distribution of upwellings and downwellings, (2) prominent downwelling on the dayside and upwellings on the nightside, and (3) prominent downwelling on the nightside and upwellings on the dayside. Hemispheric tectonics is observed for regimes (2) and (3) as a direct consequence of the day-to-night temperature contrast. Such a tectonic mode is absent in the present-day solar system and has never been inferred from astrophysical observations of exoplanets. Our models offer distinct predictions for volcanism and outgassing linked to the tectonic regime, which may explain secondary features in phase curves and allow future observations to constrain the diversity of super-Earth interiors.Comment: Accepted for publication in The Astrophysical Journal Letters; 9 pages, 5 figures; summary available at http://exoplanet-talks.org/talk/26

A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the ‘sticky air' method

Calculating surface topography in geodynamic models is a common numerical problem. Besides other approaches, the so-called ‘sticky air' approach has gained interest as a free-surface proxy at the top boundary. The often used free slip condition is thereby vertically extended by introducing a low density, low viscosity fluid layer. This allows the air/crust interface to behave in a similar manner to a true free surface. We present here a theoretical analysis that provides the physical conditions under which the sticky air approach is a valid approximation of a true free surface. Two cases are evaluated that characterize the evolution of topography on different timescales: (1) isostatic relaxation of a cosine perturbation and (2) topography changes above a rising plume. We quantitatively compare topographies calculated by six different numerical codes (using finite difference and finite element techniques) using three different topography calculation methods: (i) direct calculation of topography from normal stress, (ii) body-fitting methods allowing for meshing the topography and (iii) Lagrangian tracking of the topography on an Eulerian grid. It is found that the sticky air approach works well as long as the term (ηst/ηch)/(hst/L)3 is sufficiently small, where ηst and hst are the viscosity and thickness of the sticky air layer, and ηch and L are the characteristic viscosity and length scale of the model, respectively. Spurious lateral fluctuations of topography, as observed in some marker-based sticky air approaches, may effectively be damped by an anisotropic distribution of markers with a higher number of markers per element in the vertical than in the horizontal directio

Subduction controls the distribution and fragmentation of Earth’s tectonic plates

International audienceThe theory of plate tectonics describes how the surface of the Earth is split into an organized jigsaw of seven large plates 1 of similar sizes and a population of smaller plates, whose areas follow a fractal distribution 2,3. The reconstruction of global tectonics during the past 200 My 4 suggests that this layout is probably a long-term feature of our planet, but the forces governing it are unknown. Previous studies 3,5,6 , primarily based on statistical properties of plate distributions, were unable to resolve how the size of plates is determined by lithosphere properties and/or underlying mantle convection. Here, we demonstrate that the plate layout of the Earth is produced by a dynamic feedback between mantle convection and the strength of the lithosphere. Using 3D spherical models of mantle convection with plate-like behaviour that match the plate size-frequency distribution observed for Earth, we show that subduction geometry drives the tectonic fragmentation that generates plates. The spacing between slabs controls the layout of large plates, and the stresses caused by the bending of trenches, break plates into smaller fragments. Our results explain why the fast evolution in small back-arc plates 7,8 reflects the dramatic changes in plate motions during times of major reorganizations. Our study opens the way to use convection simulations with plate-like behaviour to unravel how global tectonics and mantle convection are dynamically connected

Convection in colloidal suspensions with particle-concentration-dependent viscosity

The onset of thermal convection in a horizontal layer of a colloidal suspension is investigated in terms of a continuum model for binary-fluid mixtures where the viscosity depends on the local concentration of colloidal particles. With an increasing difference between the viscosity at the warmer and the colder boundary the threshold of convection is reduced in the range of positive values of the separation ratio psi with the onset of stationary convection as well as in the range of negative values of psi with an oscillatory Hopf bifurcation. Additionally the convection rolls are shifted downwards with respect to the center of the horizontal layer for stationary convection (psi>0) and upwards for the Hopf bifurcation (psi<0).Comment: 8 pages, 6 figures, submitted to European Physical Journal

Thin films of fluorinated 3d-metal phthalocyanines as chemical sensors of ammonia: an optical spectroscopy study

A comparative study of the sensor response toward gaseous ammonia of hexadecafluorinated 3d-metal phthalocyanine (MPcF16, MCu(II), Co(II), Zn(II), Ni(II)) thin films was performed using complementary experimental (viz., surface plasmon resonance, SPR, and IR absorption spectroscopy) along with theoretical (density functional theory calculations, DFT) techniques. SPR measurements revealed changes of both thickness and optical parameters (refraction indices and extinction coefficients) of the MPcF16 films caused by adsorption of NH3. The MPcF16 species studied exhibited the following order of sensor response: ZnPcF16>CoPcF16≥CuPcF16>NiPcF16. A good correlation was found between the DFT calculated (B3LYP/6-311++G(2df,p)) binding energies, experimentally measured shift of the selected IR bands, and the optical sensor response. Apart from this, we performed a detailed assignment of all intense..

Hydrous upwelling across the mantle transition zone beneath the Afar Triple Junction

The mechanisms that drive the upwelling of chemical heterogeneity from the lower to upper mantle (e.g., thermal versus compositional buoyancy) are key to our understanding of whole mantle con- vective processes. We address these issues through a receiver function study on new seismic data from recent deployments located on the Afar Triple Junction, a location associated with deep mantle upwelling. The detailed images of upper mantle and mantle transition zone structure illuminate features that give insights into the nature of upwelling from the deep Earth. A seismic low-velocity layer directly above the mantle transition zone, interpreted as a stable melt layer, along with a prominent 520 km discontinuity sug- gest the presence of a hydrous upwelling. A relatively uniform transition zone thickness across the region suggests a weak thermal anomaly (<100 K) may be present and that upwelling must be at least partly driven by compositional buoyancy. The results suggest that the lower mantle is a source of volatile rich, chemically distinct upwellings that influence the structure of the upper mantle, and potentially the chemis- try of surface lavas

Generic acquisition protocol for quantitative MRI of the spinal cord

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols. The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition