Search

Mantle plumes provide valuable information about whole-mantle convection: they originate at the core-mantle boundary, cross Earth's mantle and interact with the lithosphere. For instance, it has been proposed that the mobility/stability of plumes depends on plume intrinsic properties, on how slabs interact with the basal boundary layer, on mantle flow, or on their proximity to mid-ocean ridges. Here, we use 3D-spherical models of mantle convection generating self-consistent plate-like behaviour to investigate the mechanisms linking tectonics and mantle convection to plume dynamics. Our models produce fully-dynamic mantle plumes that rise vertically with deflection <10°and present excess temperatures, rising speeds, buoyancy and heat fluxes comparable to observations. In the absence of plate tectonics, plumes are stable and their lifetime exceeds hundreds of million years. With plate tectonics, plumes are more mobile, and we identify four physical mechanisms controlling their stability. 1/ Fixed plumes are located at saddle points of basal mantle flow. 2/ Plumes moving at speeds between 0.5-1 cm yr−1 are slowly entrained by passive mantle flow. 3/ Fast plume motions between 2-5 cm yr−1 lasting several tens of million years are caused by slab push. 4/ Plumes occasionally drift at speeds >5 cm yr−1 over <10 Myr through plume merging. We do not observe systematic anchoring of plumes to mid-oceanic ridges. Independent of the presence of a dense basal layer, plate-like regimes decrease the lifetime of plumes compared to stagnant-lid models. Plume age, temperature excess or buoyancy flux are not diagnostic of plume lateral speed. The fraction of plumes moving by less than 0.5 cm yr−1 is >25%, which suggests that fixed hotspot reference frames can be defined from carefully selected hotspot tracks.

The relationships between plate motions and basal mantle structure remain poorly understood, with some models implying that the basal mantle structure has remained stable over time, while others suggest that it could be shaped by the aggregation and dispersal of supercontinents. Here we investigate the evolution of mantle flow driven by end-member plate tectonic models over 1 Gyr. We implement a tectonic scenario in which supercontinent reassembly occurs by introversion, and consider three distinct references frames that result in different net lithospheric rotation. Our flow models predict a dominant degree-2 mantle structure most of the time. We analyze the relationship between imposed tectonic velocities and deep mantle flow, and find that at spherical harmonic degree 2, the maxima of lower mantle radial flow and temperature follow the motion path of the maxima of surface divergence. It may take 160–240 Myr for lower mantle structure to reflect plate motion changes when the lower mantle is reorganized by slabs sinking onto basal thermochemical structures, and/or when slabs stagnate in the transition zone before sinking to the lower mantle. Basal thermochemical structures move at less than 0.6°/Myr in our models, with a temporal average of 0.16°/Myr when there is no net lithospheric rotation, and between 0.20 and 0.23°/Myr when net lithospheric rotation exists and is induced in the lower mantle. Our results suggest that basal thermochemical structures are not stationary, but rather linked to global plate motions and plate boundary reconfigurations, reflecting the dynamic nature of the coevolving plate- mantle system.

© 2020 International Association of Sedimentologists and European Association of Geoscientists and Engineers and John Wiley & Sons Ltd During the Paleozoic, sedimentary basins developed within Gondwana without evolving to diverging plate boundaries. Such intracontinental basins present long subsidence histories with multiple phases of accelerated subsidence that are not always easily explained by far-field tectonic forces, and may be driven by processes other than rifting and thermal subsidence. Here we investigate the subsidence of Paleozoic Australian intracontinental basins by comparing one-dimensional backstripped tectonic subsidence histories from the western Australian Canning and Southern Carnarvon Basins and the central Australian Cooper Basin to forward subsidence models for pure shear lithospheric thinning. We make the hypothesis that differences between observed and model subsidence may be explained by mantle-flow driven topography, in addition to tectonic forces. To test this hypothesis, we compute dynamic topography from the first geodynamic models of mantle flow spanning the entire Phanerozoic Eon, and we analyse the relationship between dynamic topography and anomalous basin subsidence to dynamic topography and mantle flow. Although reconstructions of mantle flow in deep geological times are uncertain, our results suggest that long-wavelength dynamic topography could explain aspects of the complex tectonic histories intracontinental basins. In the presented reconstruction of mantle flow, topographic rebound following the sinking of a Cambrian aged slab resulted in a minor phase of dynamic uplift in the Cooper Basin in middle Permian times. Throughout Carboniferous-Triassic times Australia was positioned above a mantle upwelling driven by a hot structure at the base of the mantle. Structural uplift in the Canning and Southern Carnarvon basins during the Triassic-Jurassic interval was augmented by dynamic uplift produced by that large-scale upwelling, a

Global sea level change can be inferred from sequence stratigraphic and continental flooding data. These methods reconstruct sea level from peri-cratonic and cratonic basins that are assumed to be tectonically stable and sometimes called reference districts, and from spatio-temporal correlations across basins. However, it has been understood that long-wavelength (typically hundreds of km) and low-amplitude (km) vertical displacements of the Earth's surface due to mantle flow, namely dynamic topography, can occur in the absence of crustal deformation. Dynamic topography can drive marine inundation or regional emergence of continents and must be taken into consideration for eustasy estimates. Our analysis indicates that the long-term trend in global-scale maximum flooding over the late Paleozoic generally correlates with global sea level curves. The first-order flooding history of North America correlates with some estimates of eustasy. The Paleozoic inundation of South America does not follow long-term sea level variations. The flooding lows during the Early Carboniferous and high during the Late Carboniferous are at odds with estimates of eustasy and can be explained by dynamic uplift and subsidence, respectively. Our dynamic topography models indicate that the Yangtze Platform of South China experienced significant dynamic subsidence during the transition from Permian to Triassic largely due to proto-Pacific subduction and its northward motion to collide with North China. The reference districts - Western New York, Oklahoma and Kansas, and West Texas in North America - were to some degree affected by dynamic uplift and subsidence associated with long-lived Panthalassa subduction zones, closure of the Rheic Ocean and large-scale upwelling above the African deep-mantle structure during late Paleozoic times. This indicates that some published global sea level curves may include non-eustatic signals such as dynamic uplift or subsidence. The interpretation of stratigrap

The structure of the lithosphere is key to reconciling the dynamic topography predicted by mantle convection models with residual topography derived from observations, suggest analyses of both models and data.

Dynamic topography due to mantle flow contributes to shaping Earth's evolving landscapes by affecting sediment routing, which has rarely been explored in source-to-sink contexts. Here we design a generic model to investigate the impact of dynamic topography on both landscape evolution and stratigraphic formations. An imposed wave of dynamic topography propagates laterally under a fixed continent, exerting transient surface uplift and subsidence. We find that a migrating dynamic topography can induce significant drainage reorganizations and affect sediment routing from source to sink. Variations in sediment supply driven by the lateral migrating dynamic topography contribute to the formation of diachronous unconformities along the margin. The predicted sediment flux histories are then put into perspective with the Cretaceous sedimentary records along the Southern African margins. Finally, we demonstrate that correlating offshore depositional hiatuses and unconformities has the potential to constrain the spatiotemporal evolution of past dynamic topography events.

The absolute motion of tectonic plates since Pangea can be derived from observations of hotspot trails, paleomagnetism, or seismic tomography. However, fitting observations is typically carried out in isolation without consideration for the fit to unused data or whether the resulting plate motions are geodynamically plausible. Through the joint evaluation of global hotspot track observations (for times <80 >Ma), first‐order estimates of net lithospheric rotation (NLR), and parameter estimation for paleo-trench migration (TM), we present a suite of geodynamically consistent, data‐optimized global absolute reference frames from 220 Ma to the present. Each absolute plate motion (APM) model was evaluated against six published APM models, together incorporating the full range of primary data constraints. Model performance for published and new models was quantified through a standard statistical analyses using three key diagnostic global metrics: root‐mean square plate velocities, NLR characteristics, and TM behavior. Additionally, models were assessed for consistency with published global paleomagnetic data and for ages motion, track geometry, and time dependence. Optimized APM models demonstrated significantly improved global fit with geological and geophysical observations while performing consistently with geodynamic constraints. Critically, APM models derived by limiting average rates of NLR to ~0.05°/Myr and absolute TM velocities to ~27‐mm/year fit geological observations including hotspot tracks. This suggests that this range of NLR and TM estimates may be appropriate for Earth over the last 220 Myr, providing a key step toward the practical integration of numerical geodynamics into plate tectonic reconstructions.

Earth's large-scale topography and lower-mantle structure are linked to past tectonic motions and mantle flow, making it possible to gain insights in the properties of the solid Earth from time-dependent global convection models driven by tectonic reconstructions. Recent work suggests that the amplitude of residual topography, obtained by subtracting models of isostatic topography from total topography, may be up to ~1000 m at spherical harmonic degree two (wavelength ~13 300 km) and > 200 m at spherical harmonic degree 12 (wavelength ~3000 km). The amplitude of dynamic topography and the structure of the lower mantle predicted by time-dependent forward mantle flow models both depend on the physical assumption and on model parameters. Here, I investigate the consequence of using the Boussinesq, extended Boussinesq or truncated anelastic liquid approximation (TALA) in time-dependent flow models for predicting present-day mantle structure and dynamic topography; I characterize the sensitivity of the spectral distribution of dynamic topography amplitude to the boundary conditions and model set-up for the computation of dynamic topography; and I investigate the sensitivity of model results to parameters including the depth- and temperature dependence of viscosity, the model initial age and the density of the basal layer. Extended Boussinesq and TALA models are preferred to Boussinesq models that overpredict the volume of lower-mantle slabs and the amplitude of long-wavelength dynamic topography. The correlation between dynamic and residual topography models for spherical harmonic degrees 1-3 generally ranges between 0.4 and 0.5. The flow models better predict the geographical location of large low shear velocity provinces than that of lower-mantle slabs, which cover smaller areas in map view. I show that preserving shallow lateral viscosity variations in the computation of dynamic topography increases the amplitude of dynamic topography for wavelengths > 6000 km. Parame

Whether the latitudinal distribution of climate-sensitive lithologies is stable through greenhouse and icehouse regimes remains unclear. Previous studies suggest that the palaeolatitudinal distribution of palaeoclimate indicators, including coals, evaporites, reefs and carbonates, has remained broadly similar since the Permian period, leading to the conclusion that atmospheric and oceanic circulation control their distribution rather than the latitudinal temperature gradient. Here we revisit a global-scale compilation of lithologic indicators of climate, including coals, evaporites and glacial deposits, back to the Devonian period. We test the sensitivity of their latitudinal distributions to the uneven distribution of continental areas through time and to global tectonic models, correct the latitudinal distributions of lithologies for sampling- and continental area-bias, and use statistical methods to fit these distributions with probability density functions and estimate their high-density latitudinal ranges with 50% and 95% confidence intervals. The results suggest that the palaeolatitudinal distributions of lithologies have changed through deep geological time, notably a pronounced poleward shift in the distribution of coals at the beginning of the Permian. The distribution of evaporites indicates a clearly bimodal distribution over the past ~400 Ma, except for Early Devonian, Early Carboniferous, the earliest Permian and Middle and Late Jurassic times. We discuss how the patterns indicated by these lithologies change through time in response to plate motion, orography, evolution and greenhouse/icehouse conditions. This study highlights that combining tectonic reconstructions with a comprehensive lithologic database and novel data analysis approaches provide insights into the nature and causes of shifting climatic zones through deep time.

Some changes in the topography of eastern China since Late Jurassic times cannot be well explained by lithospheric deformation. Here we analyze global mantle flow models to investigate how mantle‐driven long‐wavelength topography may have contributed to shaping the surface topography of eastern China. Paleodrainage directions suggest that a southward tilted topography once existed in eastern North China in the latest Jurassic Period, which is different from that at present day (southeastward tilting). Our model dynamic topography reveals a southward tilting topography between 160 and 150 Ma, followed by southeastward tilting and rapid subsidence, which is compatible with paleodrainage directions and tectonic subsidence of the Ordos Basin. The Cretaceous anomalous subsidence of the Songliao and North Yellow Sea basins, as well as the Cenozoic anomalous subsidence of the East China Sea Shelf Basin can also be explained by dynamic topography. An apatite fission track study in the Taihang Mountains reveals four stages of evolution: Late Jurassic fast unroofing, Cretaceous slow unroofing, early Cenozoic fast unroofing and late Cenozoic slow unroofing. We propose that mantle flow influenced this surface unroofing because the model predicts Late Jurassic dynamic uplift, Cretaceous dynamic subsidence, early Cenozoic dynamic uplift and late Cenozoic dynamic subsidence. Apatite fission track data from northern South China are also in reasonable agreement with predicted dynamic topography between 80 and 30 Ma. The spatial and temporal agreement between geological observations and model dynamic topography indicates that mantle flow has had a significant influence in shaping the surface topography of eastern China.

Although often separated as two unique subjects in science, geology and biology have been intricately intertwined since life on Earth first evolved billions of years ago.

Mantle convection shapes Earth's surface by generating dynamic topography. Observational constraints and regional convection models suggest that surface topography could be sensitive to mantle flow for wavelengths as short as 1,000 and 250 km, respectively. At these spatial scales, surface processes including sedimentation and relative sea‐level change occur on million‐year timescales. However, time‐dependent global mantle flow models do not predict small‐scale dynamic topography yet. Here we present 2‐D spherical annulus numerical models of mantle convection with large radial and lateral viscosity contrasts. We first identify the range of Rayleigh number, internal heat production rate and yield stress for which models generate plate‐like behavior, surface heat flow, surface velocities, and topography distribution comparable to Earth's. These models produce both whole‐mantle convection and small‐scale convection in the upper mantle, which results in small‐scale (km) to large‐scale (>104 km) dynamic topography, with a spectral power for intermediate scales (500 to 104 km) comparable to estimates of present‐day residual topography. Timescales of convection and the associated dynamic topography vary from five to several hundreds of millions of years. For a Rayleigh number of 107, we investigate how lithosphere yield stress variations (10-50 MPa) and the presence of deep thermochemical heterogeneities favor small‐scale (200-500 km) and intermediate‐scale (500-104 km) dynamic topography by controlling the formation of small‐scale convection and the number and distribution of subduction zones, respectively. The interplay between mantle convection and lithosphere dynamics generates a complex spatial and temporal pattern of dynamic topography consistent with constraints for Earth.

A complex history of subduction, back-arc basin formation, terrane accretion and transpressional shearing characterizes the evolution of the Caribbean and northern South American margin since Jurassic times. Quantitative plate tectonic reconstructions of the area do not include Jurassic-Cretaceous back-arc terranes of which there are both geological and geophysical observations. We developed a revised plate tectonic reconstruction based on geological observations and seismic tomography models to constrain the Jurassic-Cretaceous subduction history of eastern Panthalassa, along the western margin of the Caribbean region. This reconstruction considers the opening of a Northern Andean back-arc basin at 145 Ma, the Quebradagrande back-arc, closing at 120 Ma and followed by terrane accretion and northward translation along the South American margin starting at 100 Ma. This kinematic reconstruction is tested against two previously published tectonic reconstructions via coupling with global numerical mantle convection models using CitcomS. A comparison of modelled versus tomographically imaged mantle structure reveals that subduction outboard of the South American margin, lacking in previous tectonic models, is required to reproduce mid-mantle positive seismic anomalies imaged in P- and S-wave seismic tomography beneath South America, 500-2000 km in depth. Furthermore, we show that this subduction zone is likely produced by a back-arc basin that developed along the northern Andes during the Cretaceous via trench roll-back from 145 Ma and was closed at 100 Ma. The contemporaneous opening of the Quebradagrande back-arc basin with the Rocas Verdes back-arc basin in the southern Andes is consistent with a model that invokes return flow of mantle material behind a retreating slab and may explain why extension along the Peruvian and Chilean sections of the Andean margin did not experience full crustal break-up and back-arc opening during the late Jurassic-early Cretaceous Period

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

Traditional plate reconstruction methodologies do not allow for plate deformation to be considered. Here we present software to construct and visualize global tectonic reconstructions with deforming plates within the context of rigid plates. Both deforming and rigid plates are defined by continuously evolving polygons. The deforming regions are tessellated with triangular meshes such that either strain rate or cumulative strain can be followed. The finite strain history, crustal thickness and stretching factor of points within the deformation zones are tracked as Lagrangian points. Integrating these tools within the interactive platform GPlates enables specialized users to build and refine deforming plate models and integrate them with other models in time and space. We demonstrate the integrated platform with regional reconstructions of Cenozoic western North America, the Mesozoic South American Atlantic margin, and Cenozoic southeast Asia, embedded within global reconstructions, using different data and reconstruction strategies.

Detailed global plate motion models that provide a continuous description of plate boundaries through time are an effective tool for exploring processes both at and below the Earth's surface. A new generation of numerical models of mantle dynamics pre- and post-Pangea timeframes requires global kinematic descriptions with full plate reconstructions extending into the Paleozoic (410 Ma). Current plate models that cover Paleozoic times are characterised by large plate speeds and trench migration rates because they assume that lowermost mantle structures are rigid and fixed through time. When used as a surface boundary constraint in geodynamic models, these plate reconstructions do not accurately reproduce the present-day structure of the lowermost mantle. Building upon previous work, we present a global plate motion model with continuously closing plate boundaries ranging from the early Devonian at 410 Ma to present day. We analyse the model in terms of surface kinematics and predicted lower mantle structure. The magnitude of global plate speeds has been greatly reduced in our reconstruction by modifying the evolution of the synthetic Panthalassa oceanic plates, implementing a Paleozoic reference frame independent of any geodynamic assumptions, and implementing revised models for the Paleozoic evolution of North and South China and the closure of the Rheic Ocean. Paleozoic (410-250 Ma) RMS plate speeds are on average ∼8 cm/yr, which is comparable to Mesozoic-Cenozoic rates of ∼6 cm/yr on average. Paleozoic global median values of trench migration trend from higher speeds (∼2.5 cm/yr) in the late Devonian to rates closer to 0 cm/yr at the end of the Permian (∼250 Ma), and during the Mesozoic-Cenozoic (250-0 Ma) generally cluster tightly around ∼1.1 cm/yr. Plate motions are best constrained over the past 130 Myr and calculations of global trench convergence rates over this period indicate median rates range between 3.2 cm/yr and 12.4 cm/yr with a present day median rate

The paleogeography of New Guinea indicates fluctuating periods of flooding and emergence since the Jurassic, which are inconsistent with estimates of global sea level change since the Eocene. The role of deep Earth dynamics in explaining these discrepancies has not been explored, despite the strongly time-dependent geodynamic setting within which New Guinea has evolved. We aim to investigate the role of subduction-driven mantle flow in controlling long-wavelength dynamic topography and its manifestation in the regional sedimentary record, within a tectonically complex region leading to orogeny. We couple regionally refined global plate reconstructions with forward geodynamic models to compare trends of dynamic topography with estimates of eustasy and regional paleogeography. Qualitative corroboration of modelled mantle structure with equivalent tomographic profiles allows us to ground-truth the models. We show that predicted dynamic topography correlates with the paleogeographic record of New Guinea from the Jurassic to the present. We find that subduction at the East Gondwana margin locally enhanced the high eustatic sea levels from the Early Cretaceous (∼145 Ma) to generate long-term regional flooding. During the Miocene, however, dynamic subsidence associated with subduction of the Maramuni Arc played a fundamental role in causing long-term inundation of New Guinea during a period of global sea level fall.

Northern Africa underwent widespread inundation during the Late Cretaceous. Changes in eustasy do not explain the absence of this inundation across the remainder of Africa, and the timing and location of documented tectonic deformation do not explain the large-scale paleogeographic evolution. We investigate the combined effects of vertical surface displacements predicted by a series of mantle flow models and eustasy on northern African paleoenvironmental change. We compare changes in base level computed as the difference between eustasy and long-wavelength dynamic topography arising from sources of buoyancy deeper than 350 km to the evolution of paleoshorelines derived from two interpolated global data sets since the mid-Cretaceous. We also compare the predicted mantle temperature field of these mantle flow models at present-day to several seismic tomography models. This approach reveals that dynamic subsidence, related to Africa's northward motion away from the buoyant regions overlying the African large low shear velocity province, amplified sea level rise, resulting in maximum inundation of northern Africa during the Cenomanian and Turonian. By the Cenozoic, decreased magnitudes of dynamic subsidence, reflecting the reduced drawdown effects of slab material beneath northern Africa associated with the impact of the Africa-Eurasia collision, combined with a comparatively pronounced progressive sea level fall resulted in ongoing region-wide regression along coastal regions. The temporal match between our preferred model and the paleoshoreline data sets suggests that the paleogeographic evolution of this region since the Late Cretaceous has mainly been influenced by the interplay between eustasy and long-wavelength dynamic topography arising from large-scale, subduction-driven, lower mantle convection.

Paleogeographic reconstructions are important to understand Earth's tectonic evolution, past eustatic and regional sea level change, paleoclimate and ocean circulation, deep Earth resources and to constrain and interpret the dynamic topography predicted by mantle convection models. Global paleogeographic maps have been compiled and published, but they are generally presented as static maps with varying map projections, different time intervals represented by the maps and different plate motion models that underlie the paleogeographic reconstructions. This makes it difficult to convert the maps into a digital form and link them to alternative digital plate tectonic reconstructions. To address this limitation, we develop a workflow to restore global paleogeographic maps to their present-day coordinates and enable them to be linked to a different tectonic reconstruction. We use marine fossil collections from the Paleobiology Database to identify inconsistencies between their indicative paleoenvironments and published paleogeographic maps, and revise the locations of inferred paleo-coastlines that represent the estimated maximum transgression surfaces by resolving these inconsistencies. As a result, the consistency ratio between the paleogeography and the paleoenvironments indicated by the marine fossil collections is increased from an average of 75 % to nearly full consistency (100 %). The paleogeography in the main regions of North America, South America, Europe and Africa is significantly revised, especially in the Late Carboniferous, Middle Permian, Triassic, Jurassic, Late Cretaceous and most of the Cenozoic. The global flooded continental areas since the Early Devonian calculated from the revised paleogeography in this study are generally consistent with results derived from other paleoenvironment and paleo-lithofacies data and with the strontium isotope record in marine carbonates. We also estimate the terrestrial areal change over time associated with transferri

Tectonic plates and plate boundaries migrate substantially through time and mantle plumes are generally accepted to be mobile within the convecting mantle, but it has been proposed that large low shear velocity provinces (LLSVPs) could have been fixed and rigid for as much as 540 million years (Myr) The hypotheses of fixed and rigid LLSVPs cannot be easily tested in the absence of constraints on the past location of lowermost mantle structures. We evaluated the hypothesis of lower mantle thermochemical structure fixity with numerical experiments. As in earlier studies, we argue11 that the location of lower mantle thermochemical structures has changed through time.

Recent studies of the past eastern Australian landscape from present-day longitudinal river profiles and from mantle flow models suggest that the interaction of plate motion with mantle convection accounts for the two phases of large-scale uplift of the region since 120 Ma. We coupled the dynamic topography predicted from one of these mantle flow models to a surface process model to study the evolution of the eastern Australian landscape since the Jurassic Period. We varied the rainfall regime, erodibility, sea level variations, dynamic topography magnitude, and elastic thickness across a series of experiments. The approach accounts for erosion and sedimentation and simulates catchment dynamics. Despite the relative simplicity of our model, the results provide insights on the fundamental links between dynamic topography and continental-scale drainage evolution. Based on temporal and spatial changes in longitudinal river profiles as well as erosion and deposition maps, we show that the motion of the Australian plate over the convecting mantle has resulted in significant reorganization of the eastern Australian drainage. The model predicts that the Murray river drained eastward between 150 and ∼120 Ma, and switched to westward draining due to the tilting of the Australian plate from ∼120 Ma. First order comparisons of eight modeled river profiles and of the catchment shape of modeled Murray-Darling Basin are in agreement with present-day observations. The predicted denudation of the eastern highlands is compatible with thermochronology data and sedimentation rates along the southern Australian margin are consistent with cumulative sediment thickness.

The present-day seismic structure of the mantle under the North Atlantic Ocean indicates that the Iceland hotspot represents the surface expression of a deep mantle plume, which is thought to have erupted in the North Atlantic domain during the Palaeocene. The spatial and temporal evolution of the plume since its eruption is still highly debated, and little is known about its deep mantle history. Here, we use palaeogeographically constrained global mantle flow models to investigate the evolution of deep Earth flow beneath the North Atlantic since the Jurassic. The models show that over the last ∼ 100 Myr a remarkably stable pattern of convergent flow has prevailed in the lowermost mantle near the tip of the African Large Low-Shear Velocity Province (LLSVP), making it an ideal plume nucleation site. We extract model dynamic topography representative of a plume beneath the North Atlantic region since eruption at ∼ 60 Ma to present day and compare its evolution to available offshore geological and geophysical observations across the region. This comparison confirms that a widespread episode of Palaeocene transient uplift followed by early Eocene anomalous subsidence can be explained by the mantle-driven effects of a plume head ∼ 2500 km in diameter, arriving beneath central eastern Greenland during the Palaeocene. The location of the model plume eruption beneath eastern Greenland is compatible with several previous models. The predicted dynamic topography is within a few hundred metres of Palaeocene anomalous subsidence derived from well data. This is to be expected given the current limitations involved in modelling the evolution of Earth's mantle flow in 3-D, particularly its interactions with the base of a heterogeneous lithosphere as well as short-wavelength advective upper mantle flow, not captured in the presented global models.

A unique structure in the Earth's 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.