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

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

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

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

Geochemical observations and one layer mantle convection

cited By 125Rare gas systematics are at the heart of the discrepancy between geophysical and geochemical models. Since more and more robust evidence of whole mantle convection comes from seismic tomography and geoid modeling, the interpretation of high 3He/4He in some oceanic island basalts as being primitive has to be revisited. A time dependent model with five reservoirs (bulk mantle, continental crust, atmosphere, residual deep mantle and D″) is studied for Rb/Sr and U/Pb/He systems. The dynamics of this model correspond to whole mantle convection in which subducted oceanic crust, transformed into dense assemblages, partially segregates to form a D″ layer growing with time as in Christensen and Hofmann [J. Geophys. Res. 99 (1994) 19867-19884]. A complementary cold and depleted harzburgitic lithosphere remains above D'. We assume that hotspots arise from the deep thermal boundary layer and tap, in variable proportions, material from both the residual deep mantle and D'. The difference between HIMU and Hawaiian basalts is attributed to HIMU being mostly from strongly degassed oceanic crust, though enriched in incompatibles (D'), while Hawaii is mostly from MORB source residuals that are variably degassed and depleted. We suggest that a significant part of the Earth's radioactive elements (~1/3) is trapped in the D″ layer