Flament, N. (2016). A rapid burst in hotspot motion through the interaction of tectonics and deep mantle flow. Nature, 533 (7602), 239-242.

Even though it is well accepted that the Earth's surface topography has been affected by mantle-convection induced dynamic topography, its magnitude and time-dependence remain controversial. The dynamic influence to topographic change along continental margins is particularly difficult to unravel, because their stratigraphic record is dominated by tectonic subsidence caused by rifting. We follow a three-fold approach to estimate dynamic topographic change along passive margins based on a set of seven global mantle convection models. We first demonstrate that a geodynamic forward model that includes adiabatic and viscous heating in addition to internal heating from radiogenic sources, and a mantle viscosity profile with a gradual increase in viscosity below the mantle transition zone, provides a greatly improved match to the spectral range of residual topography end-members as compared with previous models at very long wavelengths (spherical degrees 2-3). We then combine global sea level estimates with predicted surface dynamic topography to evaluate the match between predicted continental flooding patterns and published paleo-coastlines by comparing predicted versus geologically reconstructed land fractions and spatial overlaps of flooded regions for individual continents since 140 Ma. Modelled versus geologically reconstructed land fractions match within 10% for most models, and the spatial overlaps of inundated regions are mostly between 85% and 100% for the Cenozoic, dropping to about 75-100% in the Cretaceous. Regions that have been strongly affected by mantle plumes are generally not captured well in our models, as plumes are suppressed in most of them, and our models with dynamically evolving plumes do not replicate the location and timing of observed plume products. We categorise the evolution of modelled dynamic topography in both continental interiors and along passive margins using cluster analysis to investigate how clusters of similar dynamic topography time series are distributed spatially. A subdivision of four clusters is found to best reveal end-members of dynamic topography evolution along passive margins and their hinterlands, differentiating topographic stability, longterm pronounced subsidence, initial stability over a dynamic high followed by moderate subsidence and regions that are relatively proximal to subduction zones with varied dynamic topography histories. Along passive continental margins the most commonly observed process is a gradual motion from dynamic highs towards lows during the fragmentation of Pangea, reflecting the location of many passive margins now over slabs sinking in the lower mantle. Our best-fit model results in up to 500 (± 150) m of total dynamic subsidence of continental interiors while along passive margins the maximum predicted dynamic topographic change over 140 million years is about 350 (± 150) m of subsidence. Models with plumes exhibit clusters of transient passive margin uplift of about 200 ± 200 m, but are mainly characterised by long-term subside...

In global convection models constrained by plume motions and subduction history over the last 230 Myr, plumes emerge preferentially from the edges of thermochemical structures that resemble present-day large low shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean. It has been argued that large igneous provinces (LIPs) erupting since 200 Ma may originate from plumes that emerged from the edges of the LLSVPs and numerical models have been devised to validate this hypothesis. Although qualitative assessments that are broadly in agreement with this hypothesis have been derived from numerical models, a quantitative assessment has been lacking. We present a novel plume detection scheme and derive Monte Carlo-based statistical correlations of model plume eruption sites and reconstructed LIP eruption sites. We show that models with a chemically anomalous lower mantle are highly correlated to reconstructed LIP eruption sites, whereas the confidence level obtained for a model with purely thermal plumes falls just short of 95%. A network of embayments separated by steep ridges form in the deep lower mantle in models with a chemically anomalous lower mantle. Plumes become anchored to the peaks of the chemical ridges and the network of ridges acts as a floating anchor, adjusting to slab push forces through time. The network of ridges imposes a characteristic separation between conduits that can extend into the interior of the thermochemical structures. This may explain the observed clustering of reconstructed LIP eruption sites that mostly but not exclusively occur around the present-day LLSVPs.

No abstract

Abstract. The CO2 liberated along subduction zones through intrusive/extrusive magmatic activity and the resulting active and diffuse outgassing influences global atmospheric CO2. However, when melts derived from subduction zones intersect buried carbonate platforms, decarbonation reactions may cause the contribution to atmospheric CO2 to be far greater than segments of the active margin that lacks buried carbon-rich rocks and carbonate platforms. This study investigates the contribution of carbonate-intersecting subduction zones (CISZs) to palaeo-atmospheric CO2 levels over the past 410 million years by integrating a plate motion and plate boundary evolution model with carbonate platform development through time. Our model of carbonate platform development has the potential to capture a broader range of degassing mechanisms than approaches that only account for continental arcs. Continuous and cross-wavelet analyses as well as wavelet coherence are used to evaluate trends between the evolving lengths of carbonate-intersecting subduction zones, non-carbonate-intersecting subduction zones and global subduction zones, and are examined for periodic, linked behaviour with the proxy CO2 record between 410 Ma and the present. Wavelet analysis reveals significant linked periodic behaviour between 60 and 40 Ma, when CISZ lengths are relatively high and are correlated with peaks in palaeo-atmospheric CO2, characterised by a 32–48 Myr periodicity and a ∼ 8–12 Myr lag of CO2 peaks following CISZ length peaks. The linked behaviour suggests that the relative abundance of CISZs played a role in affecting global climate during the Palaeogene. In the 200–100 Ma period, peaks in CISZ lengths align with peaks in palaeo-atmospheric CO2, but CISZ lengths alone cannot be determined as the cause of a warmer Cretaceous–Jurassic climate. Nevertheless, across the majority of the Phanerozoic, feedback mechanisms between the geosphere, atmosphere and biosphere likely played dominant roles in modulating climate. Our modelled subduction zone lengths and carbonate-intersecting subduction zone lengths approximate magmatic activity through time, and can be used as input into fully coupled models of CO2 flux between deep and shallow carbon reservoirs.