A sequential data assimilation approach for the joint reconstruction of mantle convection and surface tectonics

International audienceWith the progress of mantle convection modelling over the last decade, it now becomes possible to solve for the dynamics of the interior flow and the surface tectonics to first order. We show here that tectonic data (like surface kinematics and seafloor age distribution) and mantle convection models with plate-like behaviour can in principle be combined to reconstruct mantle convection. We present a sequential data assimilation method, based on suboptimal schemes derived from the Kalman filter, where surface velocities and seafloor age maps are not used as boundary conditions for the flow, but as data to assimilate. Two stages (a forecast followed by an analysis) are repeated sequentially to take into account data observed at different times. Whenever observations are available, an analysis infers the most probable state of the mantle at this time, considering a prior guess (supplied by the forecast) and the new observations at hand, using the classical best linear unbiased estimate. Between two observation times, the evolution of the mantle is governed by the forward model of mantle convection. This method is applied to synthetic 2-D spherical annulus mantle cases to evaluate its efficiency. We compare the reference evolutions to the estimations obtained by data assimilation. Two parameters control the behaviour of the scheme: the time between two analyses, and the amplitude of noise in the synthetic observations. Our technique proves to be efficient in retrieving temperature field evolutions provided the time between two analyses is 10 Myr. If the amplitude of the a priori error on the observations is large (30 per cent), our method provides a better estimate of surface tectonics than the observations, taking advantage of the information within the physics of convection

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

Global warming of the mantle beneath continents back to the Archaean

Throughout its history, the Earth has experienced global magmatic events that correlate with the formationof supercontinents. This suggests that the distribution of continents at the Earth's surface is fundamental inregulating mantle temperature. Nevertheless, most large igneous provinces (LIPs) are explained in terms ofthe interaction of a hot plume with the lithosphere, even though some do not show evidence for such amechanism. The aggregation of continents impacts on the temperature and flow of the underlying mantlethrough thermal insulation and enlargement of the convection wavelength. Both processes tend to increasethe temperature below the continental lithosphere, eventually triggering melting events without theinvolvement of hot plumes. This model, called mantle global warming, has been tested using 3D numericalsimulations of mantle convection [Coltice, N., Phillips, B.R., Bertrand, H., Ricard, Y., Rey, P. (2007) Globalwarming of the mantle at the origin of flood basalts over supercontinents. Geology 35, 391–394.]. Here, weapply this model to several continental flood basalts (CFBs) ranging in age from the Mesozoic to theArchaean. Our numerical simulations show that the mantle global warming model could account for thepeculiarities of magmatic provinces that developed during the formation of Pangea and Rodinia, as well asputative Archaean supercontinents such as Kenorland and Zimvaalbara

On the impact of true polar wander on heat flux patterns at the core–mantle boundary

The heat flux across the core–mantle boundary (CMB) is a fundamental variable for Earth evolution and internal dynamics. Seismic tomography provides access to seismic heterogeneities in the lower mantle, which can be related to present-day thermal heterogeneities. Alternatively, mantle convection models can be used to either infer past CMB heat flux or to produce statistically realistic CMB heat flux patterns in self-consistent models. Mantle dynamics modifies the inertia tensor of the Earth, which implies a rotation of the Earth with respect to its spin axis, a phenomenon called true polar wander (TPW). This rotation must be taken into account to link the dynamics of the mantle to the dynamics of the core. In this study, we explore the impact of TPW on the CMB heat flux over long timescales (∼1 Gyr) using two recently published mantle convection models: one model driven by a plate reconstruction and a second that self-consistently produces a plate-like behaviour. We compute the geoid in both models to correct for TPW. In the plate-driven model, we compute a total geoid and a geoid in which lateral variations of viscosity and density are suppressed above 350 km depth. An alternative to TPW correction is used for the plate-driven model by simply repositioning the model in the original paleomagnetic reference frame of the plate reconstruction. The average TPW rates range between 0.4 and 1.8° Myr−1, but peaks up to 10° Myr−1 are observed. We find that in the plate-driven mantle convection model used in this study, the maximum inertia axis produced by the model does not show a long-term consistency with the position of the magnetic dipole inferred from paleomagnetism. TPW plays an important role in redistributing the CMB heat flux, notably at short timescales (≤10 Myr). Those rapid variations modify the latitudinal distribution of the CMB heat flux, which is known to affect the stability of the magnetic dipole in geodynamo simulations. A principal component analysis (PCA) is computed to obtain the dominant CMB heat flux pattern in the different cases. These heat flux patterns are representative of the mantle convection cases studied here and can be used as boundary conditions for geodynamo models.</p

Nickel and helium evidence for melt above the core–mantle boundary

High ^(3)He/^(4)He ratios in some basalts have generally been interpreted as originating in an incompletely degassed lower-mantle source. This helium source may have been isolated at the core–mantle boundary region since Earth’s accretion. Alternatively, it may have taken part in whole-mantle convection and crust production over the age of the Earth; if so, it is now either a primitive refugium at the core–mantle boundary or is distributed throughout the lower mantle. Here we constrain the problem using lavas from Baffin Island, West Greenland, the Ontong Java Plateau, Isla Gorgona and Fernandina (Galapagos). Olivine phenocryst compositions show that these lavas originated from a peridotite source that was about 20 per cent higher in nickel content than in the modern mid-ocean-ridge basalt source. Where data are available, these lavas also have high ^(3)He/^(4)He. We propose that a less-degassed nickel-rich source formed by core–mantle interaction during the crystallization of a melt-rich layer or basal magma ocean, and that this source continues to be sampled by mantle plumes. The spatial distribution of this source may be constrained by nickel partitioning experiments at the pressures of the core–mantle boundary

Spreading continents kick-started plate tectonics

International audienceStresses acting on cold, thick and negatively buoyant oceanic litho- sphere are thought to be crucial to the initiation of subduction and the operation of plate tectonics, which characterizes the present- day geodynamics of the Earth. Because the Earth’s interior was hotter in the Archaean eon, the oceanic crust may have been thicker, thereby making the oceanic lithosphere more buoyant than at present, and whether subduction and plate tectonics occurred during this time is ambiguous, both in the geological record and in geodynamic models. Here we show that because the oceanic crust was thick and buoyant5, early continents may have produced intra-lithospheric gravitational stresses large enough to drive their gravitational spreading, to initiate subduction at their margins and to trigger episodes of subduction. Our model predicts the co-occurrence of deep to progressively shallower mafic volcanics and arc magmatism within continents in a self-consistent geodynamic framework, explaining the enigmatic multimodal volcanism and tectonic record of Archaean cratons. Moreover, our model predicts a petrological stratification and tectonic structure of the sub-continental lithospheric mantle, two predictions that are consistent with xenolith and seismic studies, respectively, and consistent with the existence of a mid-lithospheric seismic discontinuity. The slow gravitational collapse of early continents could have kick-started transient episodes of plate tectonics until, as the Earth’s interior cooled and oceanic lithosphere became heavier, plate tectonics became self-sustaining

Persistence of strong silica-enriched domains in the Earth's lower mantle

The composition of the lower mantle—comprising 56% of Earth’s volume—remains poorly constrained. Among the major elements, Mg/Si ratios ranging from ∼0.9–1.1, such as in rocky Solar-System building blocks (or chondrites), to ∼1.2–1.3, such as in upper-mantle rocks (or pyrolite), have been proposed. Geophysical evidence for subducted lithosphere deep in the mantle has been interpreted in terms of efficient mixing, and thus homogenous Mg/Si across most of the mantle. However, previous models did not consider the effects of variable Mg/Si on the viscosity and mixing efficiency of lower-mantle rocks. Here, we use geodynamic models to show that large-scale heterogeneity associated with a 20-fold change in viscosity, such as due to the dominance of intrinsically strong (Mg, Fe)SiO3–bridgmanite in low-Mg/Si domains, is sufficient to prevent efficient mantle mixing, even on large scales. Models predict that intrinsically strong domains stabilize mantle convection patterns, and coherently persist at depths of about 1,000–2,200 km up to the present-day, separated by relatively narrow up-/downwelling conduits of pyrolitic material. The stable manifestation of such bridgmanite-enriched ancient mantle structures (BEAMS) may reconcile the geographical fixity of deep-rooted mantle upwelling centres, and geophysical changes in seismic-tomography patterns, radial viscosity, rising plumes and sinking slabs near 1,000 km depth. Moreover, these ancient structures may provide a reservoir to host primordial geochemical signatures

A tectonically driven Ediacaran oxygenation event.

The diversification of complex animal life during the Cambrian Period (541-485.4 Ma) is thought to have been contingent on an oxygenation event sometime during ~850 to 541 Ma in the Neoproterozoic Era. Whilst abundant geochemical evidence indicates repeated intervals of ocean oxygenation during this time, the timing and magnitude of any changes in atmospheric pO₂ remain uncertain. Recent work indicates a large increase in the tectonic CO₂ degassing rate between the Neoproterozoic and Paleozoic Eras. We use a biogeochemical model to show that this increase in the total carbon and sulphur throughput of the Earth system increased the rate of organic carbon and pyrite sulphur burial and hence atmospheric pO₂. Modelled atmospheric pO₂ increases by ~50% during the Ediacaran Period (635-541 Ma), reaching ~0.25 of the present atmospheric level (PAL), broadly consistent with the estimated pO₂ > 0.1-0.25 PAL requirement of large, mobile and predatory animals during the Cambrian explosion