Nicolas Flament

A review of observations and models of dynamic topography

The topography of Earth is primarily controlled by lateral differences in the density structure of the crust and lithosphere. In addition to this isostatic topography, flow in the mantle induces deformation of its surface leading to dynamic topography. This transient deformation evolves over tens of millions of years, occurs at long wavelength, and is relatively small ( 5000 km), we show that there is good agreement between published residual topography fields, including the one described here, and present-day dynamic topography predicted from mantle flow models, including a new one. Residual and predicted fields show peak-to-peak amplitudes of roughly ±2 km and a dominant degree two pattern with high values for the Pacific Ocean, southern Africa, and the North Atlantic and low values for South America, western North America, and Eurasia. The flooding of continental interiors has long been known to require both larger amplitudes and to be temporally phase-shifted compared with inferred eustatic changes. Such long-wavelength inferred vertical motions have been attributed to dynamic topography. An important consequence of dynamic topography is that long-term global sea-level change cannot be estimated at a single passive margin. As a case study, we compare the results of three published models and of our model to the subsidence history of well COST-B2 offshore New Jersey. The <400 ± 45 m amount of anomalous subsidence of this well since 85 Ma is best explained by models that predict dynamic subsidence of the New Jersey margin during that period. Explicitly including the lithosphere in future global mantle flow models should not only facilitate such comparisons between model results and data, but also further constrain the nature of the coupling between the mantle and the lithosphere.

Ridge subduction sparked reorganization of the Pacific plate-mantle system 60-50 million years ago

A reorganization centered on the Pacific plate occurred ~53–47 million years ago. A “top-down” plate tectonic mechanism, complete subduction of the Izanagi plate, as opposed to a “bottom-up” mantle flow mechanism, has been proposed as the main driver. Verification based on marine geophysical observations is impossible as most ocean crust recording this event has been subducted. Using a forward modeling approach, which assimilates surface plate velocities and shallow thermal structure of slabs into mantle flow models, we show that complete Izanagi plate subduction and margin-wide slab detachment induced a major change in sub-Pacific mantle flow, from dominantly southward before 60 Ma to north-northeastward after 50 Ma. Our results agree with onshore geology, mantle tomography, and the inferred motion of the Hawaiian hot spot and are consistent with a plate tectonic process driving the rapid plate-mantle reorganization in the Pacific hemisphere between 60 and 50 Ma. This reorganization is reflected in tectonic changes in the Pacific and surrounding ocean basins.

Insights on the kinematics of the India‐Eurasia collision from global geodynamic models

The Eocene India-Eurasia collision is a first order tectonic event whose nature and chronology remains controversial. We test two end-member collision scenarios using coupled global plate motion-subduction models. The first, conventional model, invokes a continental collision soon after ∼60 Ma between a maximum extent Greater India and an Andean-style Eurasian margin. The alternative scenario involves a collision between a minimum extent Greater India and a NeoTethyan back-arc at ∼60 Ma that is subsequently subducted along southern Lhasa at an Andean-style margin, culminating with continent-continent contact at ∼40 Ma. Our numerical models suggest the conventional scenario does not adequately reproduce mantle structure related to Tethyan convergence. The alternative scenario better reproduces the discrete slab volumes and their lateral and vertical distribution in the mantle, and is also supported by the distribution of ophiolites indicative of Tethyan intraoceanic subduction, magmatic gaps along southern Lhasa and a two-stage slowdown of India. Our models show a strong component of southward mantle return flow for the Tethyan region, suggesting that the common assumption of near-vertical slab sinking is an oversimplification with significant consequences for interpretations of seismic tomography in the context of subduction reference frames.

The GPlates Portal: Cloud-Based Interactive 3D Visualization of Global Geophysical and Geological Data in a Web Browser.

The pace of scientific discovery is being transformed by the availability of ‘big data’ and open access, open source software tools. These innovations open up new avenues for how scientists communicate and share data and ideas with each other and with the general public. Here, we describe our efforts to bring to life our studies of the Earth system, both at present day and through deep geological time. The GPlates Portal (portal.gplates.org) is a gateway to a series of virtual globes based on the Cesium Javascript library. The portal allows fast interactive visualization of global geophysical and geological data sets, draped over digital terrain models. The globes use WebGL for hardware-accelerated graphics and are cross-platform and cross-browser compatible with complete camera control. The globes include a visualization of a high-resolution global digital elevation model and the vertical gradient of the global gravity field, highlighting small-scale seafloor fabric such as abyssal hills, fracture zones and seamounts in unprecedented detail. The portal also features globes portraying seafloor geology and a global data set of marine magnetic anomaly identifications. The portal is specifically designed to visualize models of the Earth through geological time. These space-time globes include tectonic reconstructions of the Earth’s gravity and magnetic fields, and several models of long-wavelength surface dynamic topography through time, including the interactive plotting of vertical motion histories at selected locations. The globes put the on-the-fly visualization of massive data sets at the fingertips of end-users to stimulate teaching and learning and novel avenues of inquiry.

Spreading continents kick-started plate tectonics

Stresses acting on cold, thick and negatively buoyant oceanic lithosphere 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 buoyant, 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.

Lower crustal flow kept Archean continental flood basalts at sea level

Large basaltic provinces as much as 15 km thick are common in Archean cratons. Many of these flood basalts erupted through continental crust but remained at sea level. Although common in the Archean record, subaqueous continental flood basalts (CFBs) are rare to absent in the post-Archean. Here we show that gravity-driven lower crustal flow may have contributed to maintaining Archean CFBs close to sea level. Our numerical experiments reveal that the characteristic time to remove the thickness anomaly associated with a CFB decreases with increasing Moho temperature (T(M)), from 500 m.y. for T(M) approximate to 320 degrees C to 1 m.y. for T(M) approximate to 900 degrees C. This strong dependency offers the opportunity to assess, from the subsidence history of CFBs, whether continental geotherms were significantly hotter in the Archean. In particular, we show that the subsidence history of the ca. 2.7 Ga upper Fortescue Group in the East Pilbara Craton, Western Australia, requires Moho temperatures >>700 degrees C. Applied to eight other unambiguous subaqueous Archean CFBs, our results indicate Moho temperatures >>650 degrees C at the time of eruption. We suggest that the decrease in the relative abundance of subaqueous CFBs over Earths history could reflect the secular cooling of the continental lithosphere due to the decrease in radiogenic heat production.

Circum‐Arctic mantle structure and long‐wavelength topography since the Jurassic

The circum-Arctic is one of the most tectonically complex regions of the world, shaped by a history of ocean basin opening and closure since the Early Jurassic. The region is characterized by contemporaneous large-scale Cenozoic exhumation extending from Alaska to the Atlantic, but its driving force is unknown. We show that the mantle flow associated with subducted slabs of the South Anuyi, Mongol-Okhotsk, and Panthalassa oceans have imparted long-wavelength deflection on overriding plates. We identify the Jurassic-Cretaceous South Anuyi slab under present-day Greenland in seismic tomography and numerical mantle flow models. Under North America, we propose the “Farallon” slab results from Andean-style ocean-continent convergence around ~30°N and from a combination of ocean-continent and intraoceanic subduction north of 50°N. We compute circum-Arctic dynamic topography through time from subduction-driven convection models and find that slabs have imparted on average <1–16 m/Myr of dynamic subsidence across the region from at least 170 Ma to ~50 Ma. With the exception of Siberia, the main phase of circum-Arctic dynamic subsidence has been followed either by slowed subsidence or by uplift of <1–6 m/Myr on average to present day. Comparing these results to geological inferences suggest that subduction-driven dynamic topography can account for rapid Middle to Late Jurassic subsidence in the Slave Craton and North Slope (respectively, <15 and 21 m/Myr, between 170 and 130 Ma) and for dynamic subsidence (<7 m/Myr, ~170–50 Ma) followed by dynamic uplift (<6 m/Myr since 50 Ma) of the Barents Sea region. Combining detailed kinematic reconstructions with geodynamic modeling and key geological observations constitutes a powerful tool to investigate the origin of vertical motion in remote regions.

Provenance of plumes in global convection models

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.

A rapid burst in hotspot motion through the interaction of tectonics and deep mantle flow

Volcanic hotspot tracks featuring linear progressions in the age of volcanism are typical surface expressions of plate tectonic movement on top of narrow plumes of hot material within Earth’s mantle. Seismic imaging reveals that these plumes can be of deep origin—probably rooted on thermochemical structures in the lower mantle. Although palaeomagnetic and radiometric age data suggest that mantle flow can advect plume conduits laterally, the flow dynamics underlying the formation of the sharp bend occurring only in the Hawaiian–Emperor hotspot track in the Pacific Ocean remains enigmatic. Here we present palaeogeographically constrained numerical models of thermochemical convection and demonstrate that flow in the deep lower mantle under the north Pacific was anomalously vigorous between 100 million years ago and 50 million years ago as a consequence of long-lasting subduction systems, unlike those in the south Pacific. These models show a sharp bend in the Hawaiian–Emperor hotspot track arising from the interplay of plume tilt and the lateral advection of plume sources. The different trajectories of the Hawaiian and Louisville hotspot tracks arise from asymmetric deformation of thermochemical structures under the Pacific between 100 million years ago and 50 million years ago. This asymmetric deformation waned just before the Hawaiian–Emperor bend developed, owing to flow in the deepest lower mantle associated with slab descent in the north and south Pacific.

Origin and evolution of the deep thermochemical structure beneath Eurasia

A unique structure in the Earths lowermost mantle, the Perm Anomaly, was recently identified beneath Eurasia. It seismologically resembles the large low-shear velocity provinces (LLSVPs) under Africa and the Pacific, but is much smaller. This challenges the current understanding of the evolution of the plate–mantle system in which plumes rise from the edges of the two LLSVPs, spatially fixed in time. New models of mantle flow over the last 230 million years reproduce the present-day structure of the lower mantle, and show a Perm-like anomaly. The anomaly formed in isolation within a closed subduction network ∼22,000 km in circumference prior to 150 million years ago before migrating ∼1,500 km westward at an average rate of 1 cm year−1, indicating a greater mobility of deep mantle structures than previously recognized. We hypothesize that the mobile Perm Anomaly could be linked to the Emeishan volcanics, in contrast to the previously proposed Siberian Traps.