Episodic zircon age spectra mimic fluctuations in subduction.

Decades of geochronological work have shown the temporal distribution of zircon ages to be episodic on billion-year timescales and seemingly coincident with the lifecycle of supercontinents, but the physical processes behind this episodicity remain contentious. The dominant, end-member models of fluctuating magmatic productivity versus selective preservation of zircon during times of continental assembly have important and very different implications for long-term, global-scale phenomena, including the history of crustal growth, the initiation and evolution of plate tectonics, and the tempo of mantle outgassing over billions of years. Consideration of this episodicity has largely focused on the Precambrian, but here we analyze a large collection of Phanerozoic zircon ages in the context of global, full-plate tectonic models that extend back to the mid-Paleozoic. We scrutinize two long-lived and relatively simple active margins, and show that along both, a relationship between the regional subduction flux and zircon age distribution is evident. In both cases, zircon age peaks correspond to intervals of high subduction flux with a ~10-30 Ma time lag (zircons trailing subduction), illuminating a possibly intrinsic delay in the subduction-related magmatic system. We also show that subduction fluxes provide a stronger correlation to zircon age distributions than subduction lengths do, implying that convergence rates play a significant role in regulating the volume of melting in subduction-related magmatic systems, and thus crustal growth.

Intraoceanic subduction spanned the Pacific in the Late Cretaceous-Paleocene.

The notorious ~60° bend separating the Hawaiian and Emperor chains marked a prominent change in the motion of the Pacific plate at ~47 Ma (million years ago), but the origin of that change remains an outstanding controversy that bears on the nature of major plate reorganizations. Lesser known but equally significant is a conundrum posed by the pre-bend (~80 to 47 Ma) motion of the Pacific plate, which, according to conventional plate models, was directed toward a fast-spreading ridge, in contradiction to tectonic forcing expectations. Using constraints provided by seismic tomography, paleomagnetism, and continental margin geology, we demonstrate that two intraoceanic subduction zones spanned the width of the North Pacific Ocean in Late Cretaceous through Paleocene time, and we present a simple plate tectonic model that explains how those intraoceanic subduction zones shaped the ~80 to 47 Ma kinematic history of the Pacific realm and drove a major plate reorganization.

Pacific plate motion change caused the Hawaiian-Emperor Bend.

A conspicuous 60° bend of the Hawaiian-Emperor Chain in the north-western Pacific Ocean has variously been interpreted as the result of an abrupt Pacific plate motion change in the Eocene (∼47 Ma), a rapid southward drift of the Hawaiian hotspot before the formation of the bend, or a combination of these two causes. Palaeomagnetic data from the Emperor Seamounts prove ambiguous for constraining the Hawaiian hotspot drift, but mantle flow modelling suggests that the hotspot drifted 4-9° south between 80 and 47 Ma. Here we demonstrate that southward hotspot drift cannot be a sole or dominant mechanism for formation of the Hawaiian-Emperor Bend (HEB). While southward hotspot drift has resulted in more northerly positions of the Emperor Seamounts as they are observed today, formation of the HEB cannot be explained without invoking a prominent change in the direction of Pacific plate motion around 47 Ma.

Archaean zircons in Miocene oceanic hotspot rocks establish ancient continental crust beneath Mauritius.

A fragment of continental crust has been postulated to underlie the young plume-related lavas of the Indian Ocean island of Mauritius based on the recovery of Proterozoic zircons from basaltic beach sands. Here we document the first U-Pb zircon ages recovered directly from 5.7 Ma Mauritian trachytic rocks. We identified concordant Archaean xenocrystic zircons ranging in age between 2.5 and 3.0 Ga within a trachyte plug that crosscuts Older Series plume-related basalts of Mauritius. Our results demonstrate the existence of ancient continental crust beneath Mauritius; based on the entire spectrum of U-Pb ages for old Mauritian zircons, we demonstrate that this ancient crust is of central-east Madagascar affinity, which is presently located ∼700 km west of Mauritius. This makes possible a detailed reconstruction of Mauritius and other Mauritian continental fragments, which once formed part of the ancient nucleus of Madagascar and southern India.

Global correlation of lower mantle structure and past subduction.

Advances in global seismic tomography have increasingly motivated identification of subducted lithosphere in Earth's deep mantle, creating novel opportunities to link plate tectonics and mantle evolution. Chief among those is the quest for a robust subduction reference frame, wherein the mantle assemblage of subducted lithosphere is used to reconstruct past surface tectonics in an absolute framework anchored in the deep Earth. However, the associations heretofore drawn between lower mantle structure and past subduction have been qualitative and conflicting, so the very assumption of a correlation has yet to be quantitatively corroborated. Here we show that a significant, time-depth progressive correlation can be drawn between reconstructed subduction zones of the last 130 Myr and positive S wave velocity anomalies at 600-2300 km depth, but that further correlation between greater times and depths is not presently demonstrable. This correlation suggests that lower mantle slab sinking rates average between 1.1 and 1.9 cm yr-1.

Continental crust beneath southeast Iceland.

The magmatic activity (0-16 Ma) in Iceland is linked to a deep mantle plume that has been active for the past 62 My. Icelandic and northeast Atlantic basalts contain variable proportions of two enriched components, interpreted as recycled oceanic crust supplied by the plume, and subcontinental lithospheric mantle derived from the nearby continental margins. A restricted area in southeast Iceland--and especially the Öræfajökull volcano--is characterized by a unique enriched-mantle component (EM2-like) with elevated (87)Sr/(86)Sr and (207)Pb/(204)Pb. Here, we demonstrate through modeling of Sr-Nd-Pb abundances and isotope ratios that the primitive Öræfajökull melts could have assimilated 2-6% of underlying continental crust before differentiating to more evolved melts. From inversion of gravity anomaly data (crustal thickness), analysis of regional magnetic data, and plate reconstructions, we propose that continental crust beneath southeast Iceland is part of ∼350-km-long and 70-km-wide extension of the Jan Mayen Microcontinent (JMM). The extended JMM was marginal to East Greenland but detached in the Early Eocene (between 52 and 47 Mya); by the Oligocene (27 Mya), all parts of the JMM permanently became part of the Eurasian plate following a westward ridge jump in the direction of the Iceland plume.

Deep mantle structure as a reference frame for movements in and on the Earth.

Earth's residual geoid is dominated by a degree-2 mode, with elevated regions above large low shear-wave velocity provinces on the core-mantle boundary beneath Africa and the Pacific. The edges of these deep mantle bodies, when projected radially to the Earth's surface, correlate with the reconstructed positions of large igneous provinces and kimberlites since Pangea formed about 320 million years ago. Using this surface-to-core-mantle boundary correlation to locate continents in longitude and a novel iterative approach for defining a paleomagnetic reference frame corrected for true polar wander, we have developed a model for absolute plate motion back to earliest Paleozoic time (540 Ma). For the Paleozoic, we have identified six phases of slow, oscillatory true polar wander during which the Earth's axis of minimum moment of inertia was similar to that of Mesozoic times. The rates of Paleozoic true polar wander (<1°/My) are compatible with those in the Mesozoic, but absolute plate velocities are, on average, twice as high. Our reconstructions generate geologically plausible scenarios, with large igneous provinces and kimberlites sourced from the margins of the large low shear-wave velocity provinces, as in Mesozoic and Cenozoic times. This absolute kinematic model suggests that a degree-2 convection mode within the Earth's mantle may have operated throughout the entire Phanerozoic.

Plate tectonic controls on atmospheric CO2 levels since the Triassic.

Climate trends on timescales of 10s to 100s of millions of years are controlled by changes in solar luminosity, continent distribution, and atmosphere composition. Plate tectonics affect geography, but also atmosphere composition through volcanic degassing of CO2 at subduction zones and midocean ridges. So far, such degassing estimates were based on reconstructions of ocean floor production for the last 150 My and indirectly, through sea level inversion before 150 My. Here we quantitatively estimate CO2 degassing by reconstructing lithosphere subduction evolution, using recent advances in combining global plate reconstructions and present-day structure of the mantle. First, we estimate that since the Triassic (250-200 My) until the present, the total paleosubduction-zone length reached up to ∼200% of the present-day value. Comparing our subduction-zone lengths with previously reconstructed ocean-crust production rates over the past 140 My suggests average global subduction rates have been constant, ∼6 cm/y: Higher ocean-crust production is associated with longer total subduction length. We compute a strontium isotope record based on subduction-zone length, which agrees well with geological records supporting the validity of our approach: The total subduction-zone length is proportional to the summed arc and ridge volcanic CO2 production and thereby to global volcanic degassing at plate boundaries. We therefore use our degassing curve as input for the GEOCARBSULF model to estimate atmospheric CO2 levels since the Triassic. Our calculated CO2 levels for the mid Mesozoic differ from previous modeling results and are more consistent with available proxy data.

4D Arctic: A Glimpse into the Structure and Evolution of the Arctic in the Light of New Geophysical Maps, Plate Tectonics and Tomographic Models.

Knowledge about the Arctic tectonic structure has changed in the last decade as a large number of new datasets have been collected and systematized. Here, we review the most updated, publicly available Circum-Arctic digital compilations of magnetic and gravity data together with new models of the Arctic's crust. Available tomographic models have also been scrutinized and evaluated for their potential to reveal the deeper structure of the Arctic region. Although the age and opening mechanisms of the Amerasia Basin are still difficult to establish in detail, interpreted subducted slabs that reside in the High Arctic's lower mantle point to one or two episodes of subduction that consumed crust of possibly Late Cretaceous-Jurassic age. The origin of major igneous activity during the Cretaceous in the central Arctic (the Alpha-Mendeleev Ridge) and in the proximity of rifted margins (the so-called High Arctic Large Igneous Province-HALIP) is still debated. Models of global plate circuits and the connection with the deep mantle are used here to re-evaluate a possible link between Arctic volcanism and mantle plumes.

Stability of active mantle upwelling revealed by net characteristics of plate tectonics.

Viscous convection within the mantle is linked to tectonic plate motions and deforms Earth's surface across wide areas. Such close links between surface geology and deep mantle dynamics presumably operated throughout Earth's history, but are difficult to investigate for past times because the history of mantle flow is poorly known. Here we show that the time dependence of global-scale mantle flow can be deduced from the net behaviour of surface plate motions. In particular, we tracked the geographic locations of net convergence and divergence for harmonic degrees 1 and 2 by computing the dipole and quadrupole moments of plate motions from tectonic reconstructions extended back to the early Mesozoic era. For present-day plate motions, we find dipole convergence in eastern Asia and quadrupole divergence in both central Africa and the central Pacific. These orientations are nearly identical to the dipole and quadrupole orientations of underlying mantle flow, which indicates that these 'net characteristics' of plate motions reveal deeper flow patterns. The positions of quadrupole divergence have not moved significantly during the past 250 million years, which suggests long-term stability of mantle upwelling beneath Africa and the Pacific Ocean. These upwelling locations are positioned above two compositionally and seismologically distinct regions of the lowermost mantle, which may organize global mantle flow as they remain stationary over geologic time.

Greater India Basin hypothesis and a two-stage Cenozoic collision between India and Asia.

Cenozoic convergence between the Indian and Asian plates produced the archetypical continental collision zone comprising the Himalaya mountain belt and the Tibetan Plateau. How and where India-Asia convergence was accommodated after collision at or before 52 Ma remains a long-standing controversy. Since 52 Ma, the two plates have converged up to 3,600 ± 35 km, yet the upper crustal shortening documented from the geological record of Asia and the Himalaya is up to approximately 2,350-km less. Here we show that the discrepancy between the convergence and the shortening can be explained by subduction of highly extended continental and oceanic Indian lithosphere within the Himalaya between approximately 50 and 25 Ma. Paleomagnetic data show that this extended continental and oceanic "Greater India" promontory resulted from 2,675 ± 700 km of North-South extension between 120 and 70 Ma, accommodated between the Tibetan Himalaya and cratonic India. We suggest that the approximately 50 Ma "India"-Asia collision was a collision of a Tibetan-Himalayan microcontinent with Asia, followed by subduction of the largely oceanic Greater India Basin along a subduction zone at the location of the Greater Himalaya. The "hard" India-Asia collision with thicker and contiguous Indian continental lithosphere occurred around 25-20 Ma. This hard collision is coincident with far-field deformation in central Asia and rapid exhumation of Greater Himalaya crystalline rocks, and may be linked to intensification of the Asian monsoon system. This two-stage collision between India and Asia is also reflected in the deep mantle remnants of subduction imaged with seismic tomography.

Diamonds sampled by plumes from the core-mantle boundary.

Diamonds are formed under high pressure more than 150 kilometres deep in the Earth's mantle and are brought to the surface mainly by volcanic rocks called kimberlites. Several thousand kimberlites have been mapped on various scales, but it is the distribution of kimberlites in the very old cratons (stable areas of the continental lithosphere that are more than 2.5 billion years old and 300 kilometres thick or more) that have generated the most interest, because kimberlites from those areas are the major carriers of economically viable diamond resources. Kimberlites, which are themselves derived from depths of more than 150 kilometres, provide invaluable information on the composition of the deep subcontinental mantle lithosphere, and on melting and metasomatic processes at or near the interface with the underlying flowing mantle. Here we use plate reconstructions and tomographic images to show that the edges of the largest heterogeneities in the deepest mantle, stable for at least 200 million years and possibly for 540 million years, seem to have controlled the eruption of most Phanerozoic kimberlites. We infer that future exploration for kimberlites and their included diamonds should therefore be concentrated in continents with old cratons that once overlay these plume-generation zones at the core-mantle boundary.

Absolute plate motions and true polar wander in the absence of hotspot tracks.

The motion of continents relative to the Earth's spin axis may be due either to rotation of the entire Earth relative to its spin axis--true polar wander--or to the motion of individual plates. In order to distinguish between these over the past 320 Myr (since the formation of the Pangaea supercontinent), we present here computations of the global average of continental motion and rotation through time in a palaeomagnetic reference frame. Two components are identified: a steady northward motion and, during certain time intervals, clockwise and anticlockwise rotations, interpreted as evidence for true polar wander. We find approximately 18 degrees anticlockwise rotation about 250-220 Myr ago and the same amount of clockwise rotation about 195-145 Myr ago. In both cases the rotation axis is located at about 10-20 degrees W, 0 degrees N, near the site that became the North American-South American-African triple junction at the break-up of Pangaea. This was followed by approximately 10 degrees clockwise rotation about 145-135 Myr ago, followed again by the same amount of anticlockwise rotation about 110-100 Myr ago, with a rotation axis in both cases approximately 25-50 degrees E in the reconstructed area of North Africa and Arabia. These rotation axes mark the maxima of the degree-two non-hydrostatic geoid during those time intervals, and the fact that the overall net rotation since 320 Myr ago is nearly zero is an indication of long-term stability of the degree-two geoid and related mantle structure. We propose a new reference frame, based on palaeomagnetism, but corrected for the true polar wander identified in this study, appropriate for relating surface to deep mantle processes from 320 Myr ago until hotspot tracks can be used (about 130 Myr ago).