Change In Plate Motion Research Articles

Coupled Evolution of Plate Tectonics and Basal Mantle Structure

AbstractThe 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.

Eocene intra-plate shortening responsible for the rise of a faunal pathway in the northeastern Caribbean realm.

Intriguing latest Eocene land-faunal dispersals between South America and the Greater Antilles (northern Caribbean) has inspired the hypothesis of the GAARlandia (Greater Antilles Aves Ridge) land bridge. This landbridge, however, should have crossed the Caribbean oceanic plate, and the geological evolution of its rise and demise, or its geodynamic forcing, remain unknown. Here we present the results of a land-sea survey from the northeast Caribbean plate, combined with chronostratigraphic data, revealing a regional episode of mid to late Eocene, trench-normal, E-W shortening and crustal thickening by ∼25%. This shortening led to a regional late Eocene-early Oligocene hiatus in the sedimentary record revealing the location of an emerged land (the Greater Antilles-Northern Lesser Antilles, or GrANoLA, landmass), consistent with the GAARlandia hypothesis. Subsequent submergence is explained by combined trench-parallel extension and thermal relaxation following a shift of arc magmatism, expressed by a regional early Miocene transgression. We tentatively link the NE Caribbean intra-plate shortening to a well-known absolute and relative North American and Caribbean plate motion change, which may provide focus for the search of the remaining connection between 'GrANoLA' land and South America, through the Aves Ridge or Lesser Antilles island arc. Our study highlights the how regional geodynamic evolution may have driven paleogeographic change that is still reflected in current biology.

Young and Old Granulites: A Volatile Connection

Granulite facies metamorphism of the lower crust has decreased in scale since the Late Archean, but many of its definitive features have persisted: (1) Punctuated, sometimes relatively short-lived, episodes of high-grade metamorphism. These are recorded, in favorably simple cases, by discrete growth rims on zircons. (2) A consistent age gap of a few to several tens of millions of years between juvenile magmatism (crustal accretion) and high-temperature metamorphism. The secondary thermal pulse is an event distinct from primary crustal accretion. (3) Involvement of mineralizing pore fluids of lowered H2O activity, that is, with high CO2 and saline concentrations. Very high oxidation states of some granulites implicate sulfur as an important fluid component. (4) Transcurrent faulting as a conspicuous feature of synmetamorphic deformation. This gives rise to characteristic transposed foliation and lineation. (5) Emplacement of coeval postorogenic K-rich granites at midcrust levels. These features can be rationalized by concepts of modern plate tectonics. High-angle plate collision is succeeded by orogen-parallel transport. This change of plate motion necessarily detaches the underthrust portion of the lithosphere, liberating asthenospheric melts and/or fluids in a postorogenic resurgence. A generation of volatile-rich mafic magmas invades the continental margin; high CO2 and halogen contents cause outgassing and freezing of the magmas at depth. Liberated volatiles effect granulite facies metamorphism by leaching H2O and lithophile elements, importantly K, and transporting these components and heat upward. Extensive melting of the lower crust is inhibited by the low H2O activity of saline-carbonic pore fluids at high pressure. Melting of orthogneiss and supracrustal rocks occurs at midcrust levels by increase of H2O activity as pressure on alkali chloride solutions falls below 0.6–0.5 GPa. The foregoing hypothesis is an alternative to the classical view that granite results from fluid-absent partial melting of, and extraction from, the lower crust, leaving granulites.

A Late Cretaceous‐Eocene Geomagnetic Polarity Timescale (MQSD20) That Steadies Spreading Rates on Multiple Mid‐Ocean Ridge Flanks

AbstractMagnetic anomalies over mid‐ocean ridge flanks record the history of geomagnetic field reversals, and the width of magnetized crustal blocks can be combined with absolute dates to generate a Geomagnetic Polarity Timescale (GPTS). We update here the current GPTS for the Late Cretaceous‐Eocene (chrons C33–C13, ~84–33 Ma) by extending to several spreading centers the analysis that originally assumed smoothly varying spreading rates in the South Atlantic. We assembled magnetic anomaly tracks from the southern Pacific (23 ship tracks), the northern Pacific (35), the southern Atlantic (45), and the Indian Ocean (51). Tracks were projected onto plate tectonic flow lines, and distances to magnetic polarity block boundaries were estimated by fitting measured magnetic anomalies with a Monte Carlo algorithm that iteratively changed block model distances and anomaly skewness angles. Distance data from each track were then assembled in summary sets of block model distances over 13 ridge flank regions. We obtained a final MQSD20 GPTS with another Monte Carlo algorithm that iteratively perturbed ages of polarity chron boundaries to minimize the variability of spreading rates over all ridge flanks and fit an up‐to‐date set of radioisotopic dates. The MQSD20 GPTS highlights a major plate motion change at 50–45 Ma, when spreading rates decreased in the Indian Ocean as India collided with Eurasia while spreading rates increased in the South Atlantic and Northern Pacific and the Hawaii‐Emperor seamount chain changed its orientation.

The impact of oblique inheritance and changes in relative plate motion on the development of rift-transform systems

In transform margins, oblique structural inheritance and plate motion vector changes have a direct impact on the margin's morphology and duration of transform activity. We investigate the effect of these two factors using numerical modelling. To simulate oblique inheritance in continental lithosphere, we model an initial rift-transform-rift configuration oriented at a range of angles (−45o to +45o) with respect to the extension direction. In a second suite of calculations, we first extend a rift-transform-rift system orthogonally and then vary the extension direction to simulate rotation of the far-field stress directions, and hence the relative plate motion. We found that transpressional deformation on the transform results in increased duration of fault activity and diffuse, longer transform zones at higher angles, while the opposite is true for the transtensional case. These observations are in good agreement with natural examples such as the Ungava Transform Zone, the Gulf of California and the Gulf of Aden, indicating that relative plate rotation plays an important role in the structural evolution of transform margins. Finally, we present a metric that links current transform margin morphology to past evolution.

Does pulsed Tibetan deformation correlate with Indian plate motion changes?

Models that aim to explain the causes of the significant Indian plate motion acceleration around 70 Ma, and the subsequent deceleration around 52 Ma predict different scenarios regarding crustal shortening of the Tibetan Plateau, which can be tested by precisely determining the timing of regional shortening events in Tibet. Here we attempt to determine this timing by presenting a high-resolution magnetostratigraphy of a ∼3.5 km thick sedimentary sequence in the syn-contractional Gonjo Basin, east-central Tibet. We successfully isolated the primary remanence as confirmed by positive fold and reversal tests. Correlation to the geomagnetic polarity time scale reveals a 69–41.5 Ma age for the Gonjo Basin sedimentary succession. Average sedimentation rates indicate two episodes of enhanced sediment accumulation rate at 69–64 Ma and 52–48 Ma, which coincide with periods of vertical axis rotation recorded in the basin fill. This coincidence suggests a tectonic cause, which given regional structures we interpret as shortening pulses. Our results are similar to those from basins elsewhere in southern, central and northern Tibet, suggesting plateau-wide, synchronous shortening pulses at ∼69–64 Ma and ∼52–48 Ma. These pulses are synchronous with major acceleration and deceleration of India-Asia convergence rate, suggesting that both the acceleration and deceleration of India-Asia convergence may be associated with enhanced crustal deformation in Tibet, which we use to evaluate previous dynamic models explaining the Indian plate motion changes and India-Asia collision processes.

Continental-scale geographic change across Zealandia during Paleogene subduction initiation

Abstract Data from International Ocean Discovery Program (IODP) Expedition 371 reveal vertical movements of 1–3 km in northern Zealandia during early Cenozoic subduction initiation in the western Pacific Ocean. Lord Howe Rise rose from deep (∼1 km) water to sea level and subsided back, with peak uplift at 50 Ma in the north and between 41 and 32 Ma in the south. The New Caledonia Trough subsided 2–3 km between 55 and 45 Ma. We suggest these elevation changes resulted from crust delamination and mantle flow that led to slab formation. We propose a “subduction resurrection” model in which (1) a subduction rupture event activated lithospheric-scale faults across a broad region during less than ∼5 m.y., and (2) tectonic forces evolved over a further 4–8 m.y. as subducted slabs grew in size and drove plate-motion change. Such a subduction rupture event may have involved nucleation and lateral propagation of slip-weakening rupture along an interconnected set of preexisting weaknesses adjacent to density anomalies.

Eocene to Miocene Subduction Initiation Recorded in Stratigraphy of Reinga Basin, Northwest New Zealand

AbstractReinga Basin, northwest of New Zealand, contains deformed Cretaceous and Cenozoic strata, flat unconformities interpreted as sea level‐modulated erosion surfaces, and volcanics. Marine conditions since the Late Cretaceous resulted in deposition of thick sequences of fossil‐rich sediments. We present a seismic stratigraphic analysis that is tied to dredged rock samples and International Ocean Discovery Program Site U1508. The Cenozoic tectonic evolution of Reinga Basin comprises four phases. (1) Folding was accompanied by uplift from lower bathyal water depths at 56–43 Ma along West Norfolk Ridge, where wave ravinement surfaces were cut. (2) Basin‐wide compression and reverse faulting at 39–34 Ma exposed Eocene and older strata along the northeastern basin margin and farther south in western New Zealand. (3) Emplacement of South Maria Allochthon at 34–28 Ma along the northeast margin postdated thick‐skinned reverse faulting and folding, was broadly synchronous with the onset of basin volcanism, and immediately predated wave ravinement of nearby ridges. Normal faulting disrupted the western basin margin at this time. (4) During the final phase of the Reinga Basin deformation, South Maria Ridge subsided ∼1,000–2,000 m from the middle shelf to bathyal depths after 23–19 Ma, and the deformation migrated southeastwards. The Northland Allochthon was emplaced at 23–20 Ma, and Northland volcanism with confirmed arc character occurred at 23–12 Ma. The early Eocene onset of tectonic activity recorded by deformation in Reinga Basin is shown to be older than previously recognized and broadly synchronous with other events related to subduction initiation and plate motion change elsewhere in the western Pacific.

Geochemical and geochronological records of tectonic changes along a flat-slab arc-transform junction: Circa 30 Ma to ca. 19 Ma Sonya Creek volcanic field, Wrangell Arc, Alaska

AbstractThe Sonya Creek volcanic field (SCVF) contains the oldest in situ volcanic products in the ca. 30 Ma–modern Wrangell Arc (WA) in south-central Alaska, which commenced due to Yakutat microplate subduction initiation. The WA occurs within a transition zone between Aleutian subduction to the west and dextral strike-slip tectonics along the Queen Charlotte–Fairweather and Denali–Duke River fault systems to the east. New 40Ar/39Ar geochronology of bedrock shows that SCVF magmatism occurred from ca. 30–19 Ma. New field mapping, physical volcanology, and major- and trace-element geochemistry, coupled with the 40Ar/39Ar ages and prior reconnaissance work, allows for the reconstruction of SCVF magmatic evolution. Initial SCVF magmatism that commenced at ca. 30 Ma records hydrous, subduction-related, calc-alkaline magmatism and also an adakite-like component that we interpret to represent slab-edge melting of the Yakutat slab. A minor westward shift of volcanism within the SCVF at ca. 25 Ma was accompanied by continued subduction-related magmatism without the adakite-like component (i.e., mantle-wedge melting), represented by ca. 25–20 Ma basaltic-andesite to dacite domes and associated diorites. These eruptions were coeval with another westward shift to anhydrous, transitional-tholeiitic, basaltic-andesite to rhyolite lavas and tuffs of the ca. 23–19 Ma Sonya Creek shield volcano; we attribute these eruptions to intra-arc extension. SCVF activity was also marked by a small southward shift in volcanism at ca. 21 Ma, characterized by hydrous calc-alkaline lavas. SCVF geochemical compositions closely overlap those from the <13 Ma WA, and no alkaline lavas that characterize the ca. 18–10 Ma eastern Wrangell volcanic belt exposed in Yukon Territory are observed. Calc-alkaline, transitional-tholeiitic, and adakite-like SCVF volcanism from ca. 30–19 Ma reflects subduction of oceanic lithosphere of the Yakutat microplate beneath North America. We suggest that the increase in magmatic flux and adakitic eruptions at ca. 25 Ma, align with a recently documented change in Pacific plate direction and velocity at this time and regional deformation events in southern Alaska. By ca. 18 Ma, SCVF activity ceased, and the locus of WA magmatism shifted to the south and east. The change in relative plate motions would be expected to transfer stress to strike-slip faults above the inboard margin of the subducting Yakutat slab, a scenario consistent with increased transtensional-related melting recorded by the ca. 23–19 Ma transitional-tholeiitic Sonya Creek shield volcano between the Denali and Totschunda faults. Moreover, we infer the Totschunda fault accommodated more than ∼85 km of horizontal offset since ca. 18 Ma, based on reconstructing the initial alignment of the early WA (i.e., 30–18 Ma SCVF) and temporally and chemically similar intrusions that crop out to the west on the opposite side of the Totschunda fault. Our results from the SCVF quantify spatial-temporal changes in deformation and magmatism that may typify arc-transform junctions over similar time scales (>10 m.y.).

Reconstruction of Subduction and Back‐Arc Spreading in the NW Pacific and Aleutian Basin: Clues to Causes of Cretaceous and Eocene Plate Reorganizations

AbstractThe Eocene (~50–45 Ma) major absolute plate motion change of the Pacific plate forming the Hawaii‐Emperor bend is thought to result from inception of Pacific plate subduction along one of its modern western trenches. Subduction is suggested to have started either spontaneously, or result from subduction of the Izanagi‐Pacific mid‐ocean ridge, or from subduction polarity reversal after collision of the Olyutorsky arc that was built on the Pacific plate with NE Asia. Here we provide a detailed plate‐kinematic reconstruction of back‐arc basins and accreted terranes in the northwest Pacific region, from Japan to the Bering Sea, since the Late Cretaceous. We present a new tectonic reconstruction of the intraoceanic Olyutorsky and Kronotsky arcs, which formed above two adjacent, oppositely dipping subduction zones at ~85 Ma within the north Pacific region, during another Pacific‐wide plate reorganization. We use our reconstruction to explain the formation of the submarine Shirshov and Bowers Ridges and show that if marine magnetic anomalies reported from the Aleutian Basin represent magnetic polarity reversals, its crust most likely formed in an ~85‐ to 60‐Ma back‐arc basin behind the Olyutorsky arc. The Olyutorsky arc was then separated from the Pacific plate by a spreading ridge, so that the ~55‐ to 50‐Ma subduction polarity reversal that followed upon Olyutorsky‐NE Asia collision initiated subduction of a plate that was not the Pacific. Hence, this polarity reversal may not be a straightforward driver of the Eocene Pacific plate motion change, whose causes remain enigmatic.

How to Create New Subduction Zones: A Global Perspective

The association of deep-sea trenches—steeply angled, planar zones where earthquakes occur deep into Earth’s interior—and chains, or arcs, of active, explosive volcanoes had been recognized for 90 years prior to the development of plate tectonic theory in the 1960s. Oceanic lithosphere is created at mid-ocean ridge spreading centers and recycled into the mantle at subduction zones, where down-going lithospheric plates dynamically sustain the deep-sea trenches. Study of subduction zone initiation is a challenge because evidence of the processes involved is typically destroyed or buried by later tectonic and crust-forming events. In 2014 and 2017, the International Ocean Discovery Program (IODP) specifically targeted these processes with three back-to-back expeditions to the archetypal Izu-Bonin-Mariana (IBM) intra-oceanic arcs and one expedition to the Tonga-Kermadec (TK) system. Both subduction systems were initiated ~52 million years ago, coincident with a proposed major change of Pacific plate motion. These expeditions explored the tectonism preceding and accompanying subduction initiation and the characteristics of the earliest crust-forming magmatism. Lack of compressive uplift in the overriding plate combined with voluminous basaltic seafloor magmatism in an extensional environment indicates a large component of spontaneous subduction initiation was involved for the IBM. Conversely, a complex range of far-field uplift and depression accompanied the birth of the TK system, indicative of a more distal forcing of subduction initiation. Future scientific ocean drilling is needed to target the three-dimensional aspects of these processes at new converging margins.

Are segmented fracture zones weak? Analytical and numerical models constrain anomalous bathymetry at the Clarion and Murray fracture zones

Between 150°W and 135°W, the Clarion Fracture Zone (FZ) coalesces from six discrete FZ traces to a single FZ over a period of ∼30 million years, coinciding with a change in plate motion that placed the associated transform fault (TF) in compression. Between 160° and 157°, the Murray FZ is also comprised of several closely-spaced (<50 km) discrete FZs and the associated TF experienced similar transpressional motion. Analysis of newly collected bathymetry data along the Clarion and Murray FZs reveals FZ structures that are inconsistent with predictions from a simple locked fault thermal subsidence model. These structures include FZ-bounded lithosphere that dips towards the old side of the FZ, whereas locked fault models predict it should dip towards the young side, and reversed fault scarp relationships where younger lithosphere lies deeper than older lithosphere. We investigate these anomalous observations using a combination of analytical and numerical models of two closely-spaced fracture zones to test how FZ strength and tectonic motions (compression or extension) affect the evolution of the intra-transform lithosphere. Our models predict that interior blocks within a segmented FZ, initially in isostatic equilibrium, can develop anomalous tilts and scarp depth relationships if the intervening FZs are relatively weak. Model results suggest that unusually low apparent coefficients of friction (<0.01) and tectonic compression are required to reproduce the magnitude of reverse tilt and scarp reversals observed in the Clarion and Murray FZs. We suggest that these low apparent friction coefficients might be the result of near lithostatic pore-fluid pressures associated with compression and alteration processes. Our results provide a tool to constrain histories of tectonic compression across FZs in places where rotation poles or other measures of plate motion may not be well constrained. Our results also imply that FZs are tectonically weak regions of hydrothermal alteration that may carry large volumes of volatiles into the mantle when subducted.

Failed rifting and fast drifting: Midcontinent Rift development, Laurentia’s rapid motion and the driver of Grenvillian orogenesis

Author(s): Swanson-Hysell, NL; Ramezani, J; Fairchild, LM; Rose, IR | Abstract: The late Mesoproterozoic was a time of large-scale tectonic activity both in the interior and on the margins of Laurentia-most notably the development of the Midcontinent Rift and the Grenvillian orogeny. Volcanism within the North American Midcontinent Rift between ca. 1109 and 1083 Ma, as well as other contemporaneous volcanism within Laurentia, has provided an opportunity to develop extensive paleomagnetic data sets spanning this time period. These data result in an apparent polar wander path (APWP) for Laurentia that goes from a high-latitude apex known as the Logan Loop into a swath known as the Keweenawan Track. A longstanding challenge of these data was the appearance of asymmetry between relatively steep reversed polarity directions from older rift rocks and relatively shallow normal polarity directions from younger rift rocks. This asymmetry was used to support an interpretation that there were large non-dipolar components to the geomagnetic field at the time. Recent data sets support the interpretation that this directional change was progressive and therefore a result of very rapid motion of Laurentia from high to low latitudes rather than a stepwise change across non-dipolar reversals. We present high-precision U-Pb dates from Midcontinent Rift volcanics that result in an improved chronostratigraphic framework for rift volcanics and unconformities that improves correlations as well as constraints on rift development. We use these dates in volcanostratigraphic context to temporally constrain a new compilation of Midcontinent Rift paleomagnetic poles. These paleomagnetic poles include new data from the North Shore Volcanic Group, Minnesota, USA and the Osler Volcanic Group, Ontario, Canada. The U-Pb dates constrain the rate of implied plate motion more precisely than has previously been possible. We apply a novel Bayesian approach to assess the rate of implied plate motion through inverting for paleomagnetic Euler poles. If the path is to be explained by a single Euler pole these inversions reveal that motion of the continent exceeded 27 cm/yr. The path is particularly well-explained by a model wherein there is continuous true polar wander in addition to rapid plate motion that changes direction and slows ca. 1096 Ma. Laurentia's movement from high to low latitudes resulted in collisional tectonics on its leading margin which could be associated with such a change in plate motion. We propose that upwelling of the Keweenawan mantle plume was associated with an avalanche of subducted slab material and associated downwelling that drove fast plate motion. This fast plate motion was followed by the Grenvillian orogeny from ca. 1090 to 980 Ma. Prolonged collisional orogenesis could have been sustained due to this strong convective cell that therefore played an integral role in the assembly of the supercontinent Rodinia.

Source and possible tectonic driver for Jurassic–Cretaceous gold deposits in the West Qinling Orogen, China

The West Qinling Orogen (WQO) in Central China Orogenic Belt contains numerous metasedimentary rock-hosted gold deposits (>2000 t Au), which mainly formed during two pulses: one previously recognized in the Late Triassic to Early Jurassic (T3–J1) and one only recently identified in the Late Jurassic to Early Cretaceous (J3–K1). Few studies have focused on the origin and geotectonic setting of the J3–K1 gold deposits.Textural relationships, LA-ICP-MS trace element and sulfur isotope compositions of pyrites in hydrothermally altered T3 dykes within the J3–K1 Daqiao deposit were used to constrain relative timing relationships between mineralization and pyrite growth in the dykes, and to characterize the source of ore fluid. These results are integrated with an overview of the regional geodynamic setting, to advance understanding of the tectonic driver for J3–K1 hydrothermal gold systems. Pyrite in breccia- and dyke-hosted gold ores at Daqiao have similar chemical and isotopic compositions and are considered to be representative of J3–K1 gold deposits in WQO. Co/Ni and sulfur isotope ratios suggest that ore fluids were derived from underlying Paleozoic Ni- and Se-rich carbonaceous sedimentary rocks. The geochemical data do not support the involvement of magmatic fluids. However, in the EQO (East Qinling Orogen), J3–K1 deposits are genetically related to magmatism. Gold mineralization in WQO is contemporaneous with magmatic deposits in the EQO and both are mainly controlled by NE- and EW-trending structures produced by changes in plate motion of the Paleo-Pacific plate as it was subducted beneath the Eurasian continent. We therefore infer that the J3–K1 structural regime facilitated the ascent of magma in the EQO and metamorphic fluids in the WQO with consequent differences in the character of contemporaneous ore deposits. If this is correct, then the far-field effects of subduction along the eastern margin of NE Asia extended 1000's of km into the continental interior.

Mid-Cenozoic Pacific plate motion change: Implications for the Northwest Hawaiian Ridge and circum-Pacific

Mid-Cenozoic Pacific plate motion change: Implications for the Northwest Hawaiian Ridge and circum-Pacific

Deciphering the origin of the Cenozoic intracontinental rifting and volcanism in eastern China using integrated evidence from the Jianghan Basin

Intracontinental rifting and low-volume volcanism are a globally common phenomenon, yet the underlying driving mechanisms and whether they can be explained through classic plate tectonic concepts, remain hotly debated. A prominent example is the Cenozoic rift and volcanic province in eastern China. Using an integration of geological, geophysical and geochemical data, we unravel the spatial and temporal variations of the rifting and volcanism in the Jianghan Basin. Both rifting and volcanism in the Jianghan Basin show two intense-to-weak cycles (65–50 Ma and 50–26 Ma, respectively) with significant enhancement in activity during the late rift phase. Moreover, rifting and depocentres progressively migrated eastward. The Jianghan basalts all share an asthenospheric origin while the source of the late phase basalts is slightly more enriched and heterogenous in Nd-Hf isotopes than that of the early phase basalts. The late phase basalts also display a smaller extent of partial melting even under a thinner lithosphere, likely indicating a significant decrease of volatile content in the mantle source. Based on regional tectonic correlations, the main stages of tectonic evolution of the Jianghan Basin and eastern China are not synchronous with changes in Pacific plate motion, while they are coincident with India-Asia collision processes. These observations lead us to propose that the asthenospheric flow driven by India-Asia collision rather than the rollback of the subducted Pacific slab has caused the widespread rifting and volcanism in eastern China. The variations of rifting and volcanism in the Jianghan Basin suggest a multiphase and eastward asthenospheric flow beneath eastern China driven by India-Asia collision, with an intense upwelling when passing through the North-South Gravity Lineament (NSGL). The much more intense rifting and volcanism during the late rift phase may indicate a much larger scale of volatile-poor asthenospheric flow than the early rift phase which could result in a more intense erosion of ancient enriched lithospheric mantle and the volatile content in the mantle source dropping sharply. This study provides an improved model based on our multidisciplinary observations for asthenospheric flow which may be an alternative driving mechanism for intracontinental rifting and low-volume volcanism in the regions where there are step changes in lithospheric thickness globally.

Sr–Pb–Nd–Hf isotopes and 40Ar/39Ar ages reveal a Hawaii–Emperor-style bend in the Rurutu hotspot

Current Pacific absolute plate motion (APM) models include two major, long-lived hotspot tracks: the ∼80 Ma Hawaii–Emperor and ∼79 Ma Louisville seamount tracks. Prior to ∼50 Ma, these two tracks show significant inter-hotspot drift, most likely due to a strong southern motion of the Hawaiian plume. With the addition of a third long-lived Pacific hotspot track, featuring a well-constrained hotspot location and a Hawaii–Emperor-style bend, we can achieve a more robust evaluation of the relationship between APM models and inter-hotspot drift. Such a candidate exists in the Rurutu hotspot track. The youngest volcanism (∼10 to 0.2 Ma) associated with the Rurutu hotspot is located in the Cook-Austral Islands and is known as the “Atiu” or “Young Rurutu” trend. Here we show with 40Ar/39Ar ages and Sr–Pb–Nd–Hf “isotopic fingerprints” that islands and seamounts ranging between ∼70 to 40 Ma fit the compositional trend and APM model-predicted locations of the Rurutu hotspot.The Rurutu hotspot has a distinctive HIMU- to FOZO-like composition, characterized by high 206Pb/204Pb (≥20.0) and low 87Sr/86Sr (≤0.704) ratios. HIMU is a rare composition, found in only two other hotspots worldwide. Uniquely, elevated Rurutu 208Pb/204Pb ratios relative to 206Pb/204Pb define a steeper regression that more closely reproduces the Northern Hemisphere Reference Line than the shallower regression defined by traditional HIMU. Rurutu also is the only HIMU-type hotspot known to continue into the West Pacific, allowing its unique composition to effectively be used as a fingerprint. The Rurutu signature can be traced through non-Samoan “interloper seamounts” in the Samoan region, to a hotspot bend in the southern Tuvalu Islands, continuing through the ∼74 to 63 Ma Gilbert Ridge and the ∼120 to 78 Ma Wake seamounts. However, precisely locating the Hawaii–Emperor style bend in Rurutu is complicated by the younger Samoan hotspot partially overlapping the Rurutu hotspot track, and Cretaceous volcanic structures predating regional Rurutu volcanism. Therefore, we present a computational method that tests a range of preexisting APM models, all generating Hawaii–Emperor style bend locations for Rurutu within error of each other. Our results confirm the presence of a third major hotspot track in the Pacific as originally proposed by W.J. Morgan (1972), with a ∼50.5–49.0 Ma bend at 178.48°E and 8.56°S.Paleomagnetic data for Hawaii indicate that changes in southern hotspot motion may explain its bend, while negligible latitudinal motion in Louisville and an unchanging distance to Rurutu (Tarduno et al., 2003; Koppers et al., 2012; Konrad et al., 2018) are consistent with subduction-driven plume motion explaining the Hawaii–Emperor bend. The presence of a defined ∼135° bend in the Rurutu hotspot track, coupled with minimal inferred latitudinal motion at that time (Konrad et al., 2018), leaves room for a small change in Pacific plate motion at ∼50 Ma.

Regional strike-slip and initial subsidence of Korea Plateau, East Sea: tectonic implications for the opening of back-arc basins

This study focuses on the initial subsidence mechanism of Korea Plateau in the East Sea. New interpretation of digitized multichannel seismic profiles (vertically unexaggerated) in the southern part of the South Korea Plateau delineates that acoustic basement comprises a number of normal-fault blocks, forming intra-plateau basins and troughs. In the western part of the plateau, the normal faults are oblique to the continental margin, whereas in the eastern part, they are bounded by an orthogonal transfer fault. The entire faults are suggestive of a large-scale dextral strike-slip fault system. The dextral fault system in the plateau was accompanied by a regional strike-slip fault along the North Korean continental margin, which offsets the prominent NE-SW granite belts on land and under the sea. The regional dextral strike-slip fault system activated during the late Oligocene/early Miocene to middle Miocene. The regional dextral strike-slip along the western continental margin of the sea was coupled with the complex dextral fault on the eastern margin from Sakhalin to Honshu. The regional strike-slip deformation was due to the changes in plate motion, i.e., southward drift (late Oligocene/early Miocene) and clockwise rotation (middle Miocene) of the Japanese Arc away from the Eurasian Plate.

Decoupling of plate-asthenosphere motion caused by non-linear viscosity during slab folding in the transition zone

Although most present-day subduction zones are in trench retreat, plate reconstructions and geological observations show that individual margins experience episodes of advancing, retreating or stationary trench motion with time-variable subduction rates. However, most laboratory and numerical simulations predict steady plate velocities and sustained trench retreat unless the slab experiences folding in the transition zone. Using 2D dynamical models of subduction with a mobile trench and overriding plate, we find that rapid sinking of the slab during folding causes a reduction in asthenosphere viscosity through the non-linear rheology, which allows the overriding plate to move in the opposite direction of the asthenosphere. This decoupling of the direction of plate and asthenosphere flow allows for episodes of rapid trench advance after each slab folding event. By analyzing the interaction between slab deformation (sinking direction and speed), stress-induced changes in asthenosphere viscosity, asthenosphere flow and plate motions, we show that there are three modes of slab-flow-plate interaction: 1) coupled trench retreat during rapid vertical sinking, 2) coupled trench advance during prograde sinking of the slab, and 3) decoupled, rapid trench advance during folding with prograde motion of the shallow slab and retrograde motion of the deep slab. These results show that non-linear viscosity plays an important role in determining the force balance controlling trench motion and conversely that trench motion can be used as a constraint on the asthenosphere viscosity underlying the overriding plate. In addition, cooling by several hundreds of degrees during episodes of fast subduction could lead to a reduction in slab dehydration and fluid-induced melting in the mantle wedge. Such cold episodes would also likely lead to time-variability in the water content and related geochemical tracers in erupted lavas, as well as the amount of water being transported by slabs into the deep mantle.

Indian Ocean floor deformation induced by the Reunion plume rather than the Tibetan Plateau

The central Indian Ocean is considered the archetypal diffuse oceanic plate boundary. Data from seismic stratigraphy and deep-sea drilling indicate that the contractional deformation of the Indian Ocean lithosphere commenced at 15.4–13.9 Ma, but experienced a sharp increase at 8–7.5 Ma. This has been maintained through to the present day, with over 80% of the shortening accrued over the past 8 Myr. Here we build on previous efforts to refine the form, timing and magnitude of the regional plate-motion changes by mitigating the noise in reconstructed Indian and Capricorn plate motions relative to Somalia. Our noise-mitigated reconstructions tightly constrain the significant speed up of the Capricorn plate over the past 8 Myr and demonstrate that the history of the Indian Ocean floor deformation cannot be explained without this plate-motion change. We propose that the Capricorn plate-motion change is driven by an increase in the eastward-directed asthenospheric flow associated with the adjacent Reunion plume, and quantitatively demonstrate the viability of this hypothesis. Our inference is supported by volcanic age distributions along the Reunion hotspot track, the anomalously high residual bathymetry of the Central Indian Ridge, full-waveform seismic tomography of the underlying asthenosphere and geochemical observations from the Central Indian Ridge. These findings challenge the commonly accepted link between the deformation of the Indian Ocean floor and the Tibetan Plateau’s orogenic evolution and demonstrate that temporal variations in upwelling mantle flow can drive major tectonic events at the Earth’s surface. Deformation of the Indian Ocean floor over the past 8 million years was caused by a change in plate motions linked to flow of the Reunion mantle plume, according to analyses of forces upon plates.