Hotspots, polar wander, Mesozoic convection and the geoid

The supercontinent cycle and Earth's long-term climate

Earth's long-term climate has been profoundly influenced by the episodic assembly and breakup of supercontinents at intervals of ca. 500 m.y. This reflects the cycle's impact on global sea level and atmospheric CO2 (and other greenhouse gases), the levels of which have fluctuated in response to variations in input from volcanism and removal (as carbonate) by the chemical weathering of silicate minerals. Supercontinent amalgamation tends to coincide with climatic cooling due to drawdown of atmospheric CO2 through enhanced weathering of the orogens of supercontinent assembly and a thermally uplifted supercontinent. Conversely, breakup tends to coincide with increased atmospheric CO2 and global warming as the dispersing continental fragments cool and subside, and weathering decreases as sea level rises. Supercontinents may also influence global climate through their causal connection to mantle plumes and large igneous provinces (LIPs) linked to their breakup. LIPs may amplify the warming trend of breakup by releasing greenhouse gases or may cause cooling and glaciation through sulfate aerosol release and drawdown of CO2 through the chemical weathering of LIP basalts. Hence, Earth's long-term climatic trends likely reflect the cycle's influence on sea level, as evidenced by Pangea, whereas its influence on LIP volcanism may have orchestrated between Earth's various climatic states.

Thermochronologic constraints on the origin of the Great Unconformity

The Great Unconformity involves a common gap of hundreds of millions to billions of years in the geologic record. The cause of this missing time has long eluded explanation, but recently two opposing hypotheses claim either a glacial or a plate tectonic origin in the Neoproterozoic. We provide thermochronologic evidence of rock cooling and multiple kilometers of exhumation in the Cryogenian Period in support of a glacial origin for erosion contributing to the composite basement nonconformity found across the North American interior. The broad synchronicity of this cooling signal at the continental scale can only be readily explained by glacial denudation.

The origin of the phenomenon known as the Great Unconformity has been a fundamental yet unresolved problem in the geosciences for over a century. Recent hypotheses advocate either global continental exhumation averaging 3 to 5 km during Cryogenian (717 to 635 Ma) snowball Earth glaciations or, alternatively, diachronous episodic exhumation throughout the Neoproterozoic (1,000 to 540 Ma) due to plate tectonic reorganization from supercontinent assembly and breakup. To test these hypotheses, the temporal patterns of Neoproterozoic thermal histories were evaluated for four North American locations using previously published medium- to low-temperature thermochronology and geologic information. We present inverse time–temperature simulations within a Bayesian modeling framework that record a consistent signal of relatively rapid, high-magnitude cooling of ∼120 to 200 °C interpreted as erosional exhumation of upper crustal basement during the Cryogenian. These models imply widespread, synchronous cooling consistent with at least ∼3 to 5 km of unroofing during snowball Earth glaciations, but also demonstrate that plate tectonic drivers, with the potential to cause both exhumation and burial, may have significantly influenced the thermal history in regions that were undergoing deformation concomitant with glaciation. In the cratonic interior, however, glaciation remains the only plausible mechanism that satisfies the required timing, magnitude, and broad spatial pattern of continental erosion revealed by our thermochronological inversions. To obtain a full picture of the extent and synchroneity of such erosional exhumation, studies on stable cratonic crust below the Great Unconformity must be repeated on all continents.

The Juvenile Hafnium Isotope Signal as a Record of Supercontinent Cycles

Hf isotope ratios measured in igneous zircon are controlled by magmatic source, which may be linked to tectonic setting. Over the 200–500 Myr periodicity of the supercontinent cycle - the principal geological phenomenon controlling prevailing global tectonic style - juvenile Hf signals, i.e. most radiogenic, are typically measured in zircon from granites formed in arc settings (crustal growth), and evolved zircon Hf signals in granites formed in continent-collision settings (crustal reworking). Interrogations of Hf datasets for excursions related to Earth events commonly use the median value, however this may be equivocal due to magma mixing. The most juvenile part of the Hf signal is less influenced by crustal in-mixing, and arguably a more sensitive archive of Earth’s geodynamic state. We analyze the global Hf dataset for this juvenile signal, statistically correlating supercontinent amalgamation intervals with evolved Hf episodes, and breakup leading to re-assembly with juvenile Hf episodes. The juvenile Hf signal is more sensitive to Pangaea and Rodinia assembly, its amplitude increasing with successive cycles to a maximum with Gondwana assembly which may reflect enhanced subduction-erosion. We demonstrate that the juvenile Hf signal carries important information on prevailing global magmatic style, and thus tectonic processes.

Thermal erosion of cratonic lithosphere as a potential trigger for mass-extinction

The temporal coincidence between large igneous provinces (LIPs) and mass extinctions has led many to pose a causal relationship between the two. However, there is still no consensus on a mechanistic model that explains how magmatism leads to the turnover of terrestrial and marine plants, invertebrates and vertebrates. Here we present a synthesis of ammonite biostratigraphy, isotopic data and high precision U-Pb zircon dates from the Triassic-Jurassic (T-J) and Pliensbachian-Toarcian (Pl-To) boundaries demonstrating that these biotic crises are both associated with rapid change from an initial cool period to greenhouse conditions. We explain these transitions as a result of changing gas species emitted during the progressive thermal erosion of cratonic lithosphere by plume activity or internal heating of the lithosphere. Our petrological model for LIP magmatism argues that initial gas emission was dominated by sulfur liberated from sulfide-bearing cratonic lithosphere before CO₂ became the dominant gas. This model offers an explanation of why LIPs erupted through oceanic lithosphere are not associated with climatic and biotic crises comparable to LIPs emitted through cratonic lithosphere.

Pangea breakup and northward drift of the Indian subcontinent reproduced by a numerical model of mantle convection

Since around 200 Ma, the most notable event in the process of the breakup of Pangea has been the high speed (up to 20 cm yr(-1)) of the northward drift of the Indian subcontinent. Our numerical simulations of 3-D spherical mantle convection approximately reproduced the process of continental drift from the breakup of Pangea at 200 Ma to the present-day continental distribution. These simulations revealed that a major factor in the northward drift of the Indian subcontinent was the large-scale cold mantle downwelling that developed spontaneously in the North Tethys Ocean, attributed to the overall shape of Pangea. The strong lateral mantle flow caused by the high-temperature anomaly beneath Pangea, due to the thermal insulation effect, enhanced the acceleration of the Indian subcontinent during the early stage of the Pangea breakup. The large-scale hot upwelling plumes from the lower mantle, initially located under Africa, might have contributed to the formation of the large-scale cold mantle downwelling in the North Tethys Ocean.

Petrological evidence for secular cooling in mantle plumes

Geological mapping and geochronological studies have shown much lower eruption rates for ocean island basalts (OIBs) in comparison with those of lavas from large igneous provinces (LIPs) such as oceanic plateaux and continental flood provinces. However, a quantitative petrological comparison has never been made between mantle source temperature and the extent of melting for OIB and LIP sources. Here we show that the MgO and FeO contents of Galapagos-related lavas and their primary magmas have decreased since the Cretaceous period. From petrological modelling, we infer that these changes reflect a cooling of the Galapagos mantle plume from a potential temperature of 1,560-1,620 degrees C in the Cretaceous to 1,500 degrees C at present. Iceland also exhibits secular cooling, in agreement with previous studies. Our work provides quantitative petrological evidence that, in general, mantle plumes for LIPs with Palaeocene-Permian ages were hotter and melted more extensively than plumes of more modern ocean islands. We interpret this to reflect episodic flow from lower-mantle domains that are lithologically and geochemically heterogeneous.

Modeling the dynamics of a free boundary on turbulent thermal convection

Based on our previous experimental study, we present a one-dimensional phenomenological model of a thermal blanket floating on the upper surface of a thermally convecting fluid. The model captures the most important interactions between the floating solid and the fluid underneath. By the thermal blanketing effect, the presence of the solid plate modifies the flow structure below; in turn, the flow exerts a viscous drag that causes the floating boundary to move. An oscillatory state and a trapped state are found in this model, which is in excellent agreement with experimental observations. The model also offers details on the transition between the states, and gives useful insights on this coupled system without the need for full-scale simulations.

The Phanerozoic Record of Global Sea-Level Change

We review Phanerozoic sea-level changes [543 million years ago (Ma) to the present] on various time scales and present a new sea-level record for the past 100 million years (My). Long-term sea level peaked at 100 +/- 50 meters during the Cretaceous, implying that ocean-crust production rates were much lower than previously inferred. Sea level mirrors oxygen isotope variations, reflecting ice-volume change on the 10(4)- to 10(6)-year scale, but a link between oxygen isotope and sea level on the 10(7)-year scale must be due to temperature changes that we attribute to tectonically controlled carbon dioxide variations. Sea-level change has influenced phytoplankton evolution, ocean chemistry, and the loci of carbonate, organic carbon, and siliciclastic sediment burial. Over the past 100 My, sea-level changes reflect global climate evolution from a time of ephemeral Antarctic ice sheets (100 to 33 Ma), through a time of large ice sheets primarily in Antarctica (33 to 2.5 Ma), to a world with large Antarctic and large, variable Northern Hemisphere ice sheets (2.5 Ma to the present).

A correlation between mid-ocean-ridge basalt chemistry and distance to continents

To fully understand the structure and dynamics of the Earth's convecting mantle, the origins of temperature variations within the mantle need to be resolved. Different hypotheses have been proposed to account for these temperature variations: for example, heat coming from the decay of radioactive elements or heat flowing out of the Earth's core. In addition, theoretical studies suggest that the thermal properties of continental masses can affect mantle convection, but quantitative data that could allow us to test these models are scarce. To address this latter problem, we have examined the chemistry of mid-ocean-ridge basalt--which reflects the temperature of the source mantle--as a function of the distance of the ridge from the closest continental margin. No correlation is observed for oceanic ridges close to subduction zones or hotspots; subduction zones probably inhibit thermal transfer between the mantle beneath continents and ocean, whereas hotspots influence the major-element chemistry of ridge basalts, which makes their interpretation with respect to mantle temperature more difficult. However, we do observe a significant correlation for mid-oceanic basalts from the Atlantic and Indian oceans. From this, we conclude that the location of continental masses relative to active ridges influences the large-scale thermal structure of the mantle and we estimate that the mantle cools by 0.05 to 0.1 degrees C per kilometre from the continental margins.

A Lower Mantle Source for Central European Volcanism

Cenozoic rifting and volcanism in Europe have been associated with either passive or active mantle upwellings. Tomographic images show a low velocity structure between 660- and 2000-kilometer depth, which we propose to represent a lower mantle upwelling under central Europe that may feed smaller upper-mantle plumes. The position of the rift zones in the foreland of the Alpine belts and the relatively weak volcanism compared to other regions with plume-associated volcanism are probably the result of the past and present subduction under southern Europe.

Long-wavelength variations in Earth’s geoid: physical models and dynamical implications

The seismic velocity anomalies resolved by seismic tomography are associated with variations in density that lead to convective flow and to dynamically maintained topography at the Earth’s surface, the core-mantle boundary (CMB), and any interior chemical boundaries that might exist. The dynamic topography resulting from a given density field is very sensitive to viscosity structure and to chemical stratification. The mass anomalies resulting from dynamic topography have a major effect on the geoid, which places strong constraints on mantle structure. Almost 90% of the observed geoid can be explained by density anomalies inferred from tomography and a model of subducted slabs, along with the resulting dynamic topography predicted for an Earth model with a low-viscosity asthenosphere ( ca . 10 20 Pa s) overlying a moderate viscosity ( ca . 10 22.5 Pa s) lower mantle. This viscosity stratification would lead to rapid mixing in the asthenosphere, with little mixing in the lower mantle. Chemically stratified models can also explain the geoid, but they predict hundreds of kilometres of dynamic topography at the 670 km discontinuity, a prediction currently unsupported by observation. A low-viscosity or chemically distinct D" layer tends to decouple CMB topography from convective circulation in the overlying mantle. Dynamic topography at the surface should result in long-term changes in eustatic sea level.

Mantle Convection with Plates and Mobile, Faulted Plate Margins

A finite-element formulation of faults has been incorporated into time-dependent models of mantle convection with realistic rheology, continents, and phase changes. Realistic tectonic plates naturally form with self-consistent coupling between plate and mantle dynamics. After the initiation of subduction, trenches rapidly roll back with subducted slabs temporarily laid out along the base of the transition zone. After the slabs have penetrated into the lower mantle, the velocity of trench migration decreases markedly. The inhibition of slab penetration into the lower mantle by the 670-kilometer phase change is greatly reduced in these models as compared to models without tectonic plates.

Relation of Major Volcanic Center Concentration on Venus to Global Tectonic Patterns

Global analysis of NASA Magellan image data indicates that a major concentration of volcanic centers covering approximately 40 percent of the surface of Venus occurs between the Beta, Atla, and Themis regiones. Associated with this enhanced concentration are geological characteristics commonly interpreted as rifting and mantle upwelling. Interconnected low plains in an annulus around this concentration are characterized by crustal shortening and infrequent volcanic centers that may represent sites of mantle return flow and net down-welling. Together, these observations suggest the existence of relatively simple, largescale patterns of mantle circulation similar to those associated with concentrations of intraplate volcanism on Earth.

A Cold Suboceanic Mantle Belt at the Earth's Equator

An exceptionally low degree of melting of the upper mantle in the equatorial part of the mid-Atlantic Ridge is indicated by the chemical composition of mantle-derived mid-ocean ridge peridotites and basalts. These data imply that mantle temperatures below the equatorial Atlantic are at least approximately 150 degrees C cooler than those below the normal mid-Atlantic Ridge, suggesting that isotherms are depressed and the mantle is downwelling in the equatorial Atlantic. An equatorial minimum of the zero-age crustal elevation of the East Pacific Rise suggests a similar situation in the Pacific. If so, an oceanic upper mantle cold equatorial belt separates hotter mantle regimes and perhaps distinct chemical and isotopic domains in the Northern and Southern hemispheres. Gravity data suggest the presence of high density material in the oceanic equatorial upper mantle, which is consistent with its inferred low temperature and undepleted composition. The equatorial distribution of cold, dense upper mantle may be ultimately an effect of the Earth's rotation.

Plate Tectonics and Hotspots: The Third Dimension

High-resolution seismic tomographic models of the upper mantle provide powerful new constraints on theories of plate tectonics and hotspots. Midocean ridges have extremely low seismic velocities to a depth of 100 kilometers. These low velocities imply partial melting. At greater depths, low-velocity and high-velocity anomalies record, respectively, previous positions of migrating ridges and trenches. Extensional, rifting, and hotspot regions have deep (> 200 kilometers) low-velocity anomalies. The upper mantle is characterized by vast domains of high temperature rather than small regions surrounding hotspots; the asthenosphere is not homogeneous or isothermal. Extensive magmatism requires a combination of hot upper mantle and suitable lithospheric conditions. High-velocity regions of the upper 200 kilometers of the mantle correlate with Archean cratons.