Boundary-layer mantle flow under the Dead Sea transform fault inferred from seismic anisotropy

Lithospheric-scale transform faults play an important role in the dynamics of global plate motion. Near-surface deformation fields for such faults are relatively well documented by satellite geodesy, strain measurements and earthquake source studies, and deeper crustal structure has been imaged by seismic profiling. Relatively little is known, however, about deformation taking place in the subcrustal lithosphere--that is, the width and depth of the region associated with the deformation, the transition between deformed and undeformed lithosphere and the interaction between lithospheric and asthenospheric mantle flow at the plate boundary. Here we present evidence for a narrow, approximately 20-km-wide, subcrustal anisotropic zone of fault-parallel mineral alignment beneath the Dead Sea transform, obtained from an inversion of shear-wave splitting observations along a dense receiver profile. The geometry of this zone and the contrast between distinct anisotropic domains suggest subhorizontal mantle flow within a vertical boundary layer that extends through the entire lithosphere and accommodates the transform motion between the African and Arabian plates within this relatively narrow zone.     

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.     

Resistance to mantle flow inferred from the electromagnetic strike of the Australian upper mantle

Seismic anisotropy is thought to result from the strain-induced lattice-preferred orientation of mantle minerals, especially olivine, owing to shear waves propagating faster along the a-axis of olivine crystals than along the other axes. This anisotropy results in birefringence, or 'shear-wave splitting', which has been investigated in numerous studies. Although olivine is also anisotropic with respect to electrical conductivity (with the a-axis being most conductive), few studies of the electrical anisotropy of the upper mantle have been undertaken, and these have been limited to relatively shallow depths in the lithospheric upper mantle. Theoretical models of mantle flow have been used to infer that, for progressive simple shear imparted by the motion of an overriding tectonic plate, the a-axes of olivine crystals should align themselves parallel to the direction of plate motion. Here, however, we show that a significant discrepancy exists between the electromagnetic strike of the mantle below Australia and the direction of present-day absolute plate motion. We infer from this discrepancy that the a-axes of olivine crystals are not aligned with the direction of the present-day plate motion of Australia, indicating resistance to deformation of the mantle by plate motion.     

Natural examples of olivine lattice preferred orientation patterns with a flow-normal a-axis maximum

Tectonic plate motion is thought to cause solid-state plastic flow within the underlying upper mantle and accordingly lead to the development of a lattice preferred orientation of the constituent olivine crystals. The mechanical anisotropy that results from such preferred orientation typically produces a direction of maximum seismic wave velocity parallel to the plate motion direction. This has been explained by the existence of an olivine preferred orientation with an 'a-axis' maximum parallel to the induced mantle flow direction. In subduction zones, however, the olivine a axes have been inferred to be arranged roughly perpendicular to plate motion, which has usually been ascribed to localized complex mantle flow patterns. Recent experimental work suggests an alternative explanation: under conditions of high water activity, a 'B-type' olivine preferred orientation may form, with the a-axis maximum perpendicular to the flow direction. Natural examples of such B-type preferred orientation are, however, almost entirely unknown. Here we document widespread B-type olivine preferred orientation patterns from a subduction-type metamorphic belt in southwest Japan and show that these patterns developed in the presence of water. Our discovery implies that mantle flow above subduction zones may be much simpler than has generally been thought.     

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.     

Magma-assisted rifting in Ethiopia

The rifting of continents and evolution of ocean basins is a fundamental component of plate tectonics, yet the process of continental break-up remains controversial. Plate driving forces have been estimated to be as much as an order of magnitude smaller than those required to rupture thick continental lithosphere. However, Buck has proposed that lithospheric heating by mantle upwelling and related magma production could promote lithospheric rupture at much lower stresses. Such models of mechanical versus magma-assisted extension can be tested, because they predict different temporal and spatial patterns of crustal and upper-mantle structure. Changes in plate deformation produce strain-enhanced crystal alignment and increased melt production within the upper mantle, both of which can cause seismic anisotropy. The Northern Ethiopian Rift is an ideal place to test break-up models because it formed in cratonic lithosphere with minor far-field plate stresses. Here we present evidence of seismic anisotropy in the upper mantle of this rift zone using observations of shear-wave splitting. Our observations, together with recent geological data, indicate a strong component of melt-induced anisotropy with only minor crustal stretching, supporting the magma-assisted rifting model in this area of initially cold, thick continental lithosphere.     

Destruction geodynamics of the North China craton and its Paleoproterozoic plate tectonics

Much attention has been paid in the last two decades to the physical and chemical processes as well as temporal-spatial variations of the lithospheric mantle beneath the North China Craton. In order to provide insights into the geodynamics of this variation, it is necessary to thoroughly study the state and structure of the lithospheric crust and mantle of the North China Craton and its adjacent regions as an integrated unit. Based on the velocity structure of the crust and upper mantle constrained from seismological studies, this paper presents various available geophysical results regarding the lithosphere thickness, the nature of crust-mantle boundary, the upper mantle structure and deformation characteristics as well as their tectonic features and evolution systematics. Combined with the obtained data from petrology and geochemistry, a mantle flow model is proposed for the tectonic evolution of the North China Craton during the Mesozoic-Cenozoic. We suggest that subduction of the Pacific plate made the mantle underneath the eastern Asian continent unstable and able to flow faster. Such a regional mantle flow system would cause an elevation of melt/fluid content in the upper mantle of the North China Craton and the lithospheric softening, which, subsequently resulted in destruction of the North China Craton in different ways of delamination and thermal erosion in Yanshan, Taihang Mountains and the Tan-Lu Fault zone. Multiple lines of evidence recorded in the crust of the North China Craton, such as the amalgamation of the Archean eastern and western blocks, the subduction of Paleo-oceanic crust and Paleo-continental residue, indicate that the Earth in the Paleoproterozoic had already evolved into the plate tectonic system similar to the present plate tectonics.     

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.     

Peridotite-melt interaction: A key point for the destruction of cratonic lithospheric mantle

This paper presents an overview of recent studies dealing with different ages of mantle peridotitic xenoliths and xenocrysts from the North China Craton, with aim to provide new ideas for further study on the destruction of the North China Craton. Re-Os isotopic studies suggest that the lithospheric mantle of the North China Craton is of Archean age prior to its thinning. The key reason why such a low density and highly refractory Archean lithospheric mantle would be thinned is changes in composition, thermal regime, and physical properties of the lithospheric mantle due to interaction of peridotites with melts of different origins. Inward subduction of circum craton plates and collision with the North China Craton provided not only the driving force for the destruction of the craton, but also continuous melts derived from partial melting of subducted continental or oceanic crustal materials that resulted in the compositional change of the lithospheric mantle. Regional thermal anomaly at ca. 120 Ma led to the melting of highly modified lithospheric mantle. At the same time or subsequently lithospheric extension and asthenospheric upwelling further reinforced the melting and thinning of the lithospheric mantle. Therefore, the destruction and thinning of the North China Craton is a combined result of per- idotite-melt interaction （addition of volatile）, enhanced regional thermal anomaly （temperature increase） and lithospheric extension （decompression）. Such a complex geological process finally produced a ＂mixed＂ lithospheric mantle of highly chemical heterogeneity during the Mesozoic and Cenozoic. It also resulted in significant difference in the composition of mantle peridotitic xenoliths between different regions and times.     

Arc-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge or eroded fore-arc complexes but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (63-190 mm yr(-1)) and are comparable to the magnitude of subducting Cocos plate motion (approximately 85 mm yr(-1)). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.     

Uppermost mantle structure of the North China Craton: Constraints from interstation Pn travel time difference tomography

The uppermost mantle is the key area for exchange of heat flux and material convection between the crust and lithospheric mantle. Spatial variations of lithospheric thinning and dynamic processes in the North China Craton could inevitably induce the velocity heterogeneity in the uppermost mantle. In this study, we used Pn arrivals from permanent seismic stations in North China and surrounding regions to construct a tomographic image of the North China Craton. The tomographic method with Pn travel time difference data were used to study the velocity variations in the uppermost mantle. Pn velocities in the uppermost mantle varied significantly in the Eastern, Central and Western blocks of the North China Craton. This suggests that the lithosphere beneath different blocks of the North China Craton have experienced distinct tectonic evolutions and dynamic processes since the Paleozoic. The current uppermost mantle has been imprinted by these tectonic and dynamic processes. Fast Pn velocities are prominent beneath the Bohai Bay Basin in the Eastern Block of the North China Craton, suggesting residuals of the Archean lithospheric mantle. Beneath the Tanlu Fault Zone and Bohai Sea, slow Pn velocities are present in the uppermost mantle, which can be attributed to significant lithospheric thinning and asthenospheric upwelling. The newly formed lithospheric mantle beneath Yanshan Mountain may be the dominant reason for the existence of slow Pn velocities in this region. Conversely, the ancient lower crust and lithospheric mantle already have been delaminated. In the Central Block, significant slow Pn velocities are present in Taihangshan Mountain, which also extends northward to the Yinchuan-Hetao Rift on the northern margin of the Ordos Block and Yinshan Orogen. This characteristic probably is a result of hot asthenospheric upwelling along the active tectonic boundary on the margin of the Western Block. The protracted thermal erosion and underplating of hot asthenospheric upwelling may induce lithospheric thinning and significant slow velocities in the uppermost mantle. Fast velocities beneath the Western Block suggest that the thick, cold and refractory Archean lithospheric keel of craton still is retained without apparent destruction.     

Global azimuthal seismic anisotropy and the unique plate-motion deformation of Australia

Differences in the thickness of the high-velocity lid underlying continents as imaged by seismic tomography, have fuelled a long debate on the origin of the 'roots' of continents. Some of these differences may be reconciled by observations of radial anisotropy between 250 and 300 km depth, with horizontally polarized shear waves travelling faster than vertically polarized ones. This azimuthally averaged anisotropy could arise from present-day deformation at the base of the plate, as has been found for shallower depths beneath ocean basins. Such deformation would also produce significant azimuthal variation, owing to the preferred alignment of highly anisotropic minerals. Here we report global observations of surface-wave azimuthal anisotropy, which indicate that only the continental portion of the Australian plate displays significant azimuthal anisotropy and strong correlation with present-day plate motion in the depth range 175-300 km. Beneath other continents, azimuthal anisotropy is only weakly correlated with plate motion and its depth location is similar to that found beneath oceans. We infer that the fast-moving Australian plate contains the only continental region with a sufficiently large deformation at its base to be transformed into azimuthal anisotropy. Simple shear leading to anisotropy with a plunging axis of symmetry may explain the smaller azimuthal anisotropy beneath other continents.     

Preservation of ancient and fertile lithospheric mantle beneath the southwestern United States

Stable continental regions, free from tectonic activity, are generally found only within ancient cratons-the centres of continents which formed in the Archaean era, 4.0-2.5 Gyr ago. But in the Cordilleran mountain belt of western North America some younger (middle Proterozoic) regions have remained stable, whereas some older (late Archaean) regions have been tectonically disturbed, suggesting that age alone does not determine lithospheric strength and crustal stability. Here we report rhenium-osmium isotope and mineral compositions of peridotite xenoliths from two regions of the Cordilleran mountain belt. We found that the younger, undeformed Colorado plateau is underlain by lithospheric mantle that is 'depleted' (deficient in minerals extracted by partial melting of the rock), whereas the older (Archaean), yet deformed, southern Basin and Range province is underlain by 'fertile' lithospheric mantle (not depleted by melt extraction). We suggest that the apparent relationship between composition and lithospheric strength, inferred from different degrees of crustal deformation, occurs because depleted mantle is intrinsically less dense than fertile mantle (due to iron having been lost when melt was extracted from the rock). This allows the depleted mantle to form a thicker thermal boundary layer between the deep convecting mantle and the crust, thus reducing tectonic activity at the surface. The inference that not all Archaean crust developed a strong and thick thermal boundary layer leads to the possibility that such ancient crust may have been overlooked because of its intensive reworking or lost from the geological record owing to preferential recycling.     

Seismic reflection imaging of two megathrust shear zones in the northern Cascadia subduction zone

At convergent continental margins, the relative motion between the subducting oceanic plate and the overriding continent is usually accommodated by movement along a single, thin interface known as a megathrust. Great thrust earthquakes occur on the shallow part of this interface where the two plates are locked together. Earthquakes of lower magnitude occur within the underlying oceanic plate, and have been linked to geochemical dehydration reactions caused by the plate's descent. Here I present deep seismic reflection data from the northern Cascadia subduction zone that show that the inter-plate boundary is up to 16 km thick and comprises two megathrust shear zones that bound a >5-km-thick, approximately 110-km-wide region of imbricated crustal rocks. Earthquakes within the subducting plate occur predominantly in two geographic bands where the dip of the plate is inferred to increase as it is forced around the edges of the imbricated inter-plate boundary zone. This implies that seismicity in the subducting slab is controlled primarily by deformation in the upper part of the plate. Slip on the shallower megathrust shear zone, which may occur by aseismic slow slip, will transport crustal rocks into the upper mantle above the subducting oceanic plate and may, in part, provide an explanation for the unusually low seismic wave speeds that are observed there.     

Imaging the Indian subcontinent beneath the Himalaya

The rocks of the Indian subcontinent are last seen south of the Ganges before they plunge beneath the Himalaya and the Tibetan plateau. They are next glimpsed in seismic reflection profiles deep beneath southern Tibet, yet the surface seen there has been modified by processes within the Himalaya that have consumed parts of the upper Indian crust and converted them into Himalayan rocks. The geometry of the partly dismantled Indian plate as it passes through the Himalayan process zone has hitherto eluded imaging. Here we report seismic images both of the decollement at the base of the Himalaya and of the Moho (the boundary between crust and mantle) at the base of the Indian crust. A significant finding is that strong seismic anisotropy develops above the decollement in response to shear processes that are taken up as slip in great earthquakes at shallower depths. North of the Himalaya, the lower Indian crust is characterized by a high-velocity region consistent with the formation of eclogite, a high-density material whose presence affects the dynamics of the Tibetan plateau.     

Continental collision slowing due to viscous mantle lithosphere rather than topography

Because the inertia of tectonic plates is negligible, plate velocities result from the balance of forces acting at plate margins and along their base. Observations of past plate motion derived from marine magnetic anomalies provide evidence of how continental deformation may contribute to plate driving forces. A decrease in convergence rate at the inception of continental collision is expected because of the greater buoyancy of continental than oceanic lithosphere, but post-collisional rates are less well understood. Slowing of convergence has generally been attributed to the development of high topography that further resists convergent motion; however, the role of deforming continental mantle lithosphere on plate motions has not previously been considered. Here I show that the rate of India's penetration into Eurasia has decreased exponentially since their collision. The exponential decrease in convergence rate suggests that contractional strain across Tibet has been constant throughout the collision at a rate of 7.03?×?10(-16)?s(-1), which matches the current rate. A constant bulk strain rate of the orogen suggests that convergent motion is resisted by constant average stress (constant force) applied to a relatively uniform layer or interface at depth. This finding follows new evidence that the mantle lithosphere beneath Tibet is intact, which supports the interpretation that the long-term strain history of Tibet reflects deformation of the mantle lithosphere. Under conditions of constant stress and strength, the deforming continental lithosphere creates a type of viscous resistance that affects plate motion irrespective of how topography evolved.     

Dynamics of Mid-Palaeocene North Atlantic rifting linked with European intra-plate deformations

The process of continental break-up provides a large-scale experiment that can be used to test causal relations between plate tectonics and the dynamics of the Earth's deep mantle. Detailed diagnostic information on the timing and dynamics of such events, which are not resolved by plate kinematic reconstructions, can be obtained from the response of the interior of adjacent continental plates to stress changes generated by plate boundary processes. Here we demonstrate a causal relationship between North Atlantic continental rifting at approximately 62 Myr ago and an abrupt change of the intra-plate deformation style in the adjacent European continent. The rifting involved a left-lateral displacement between the North American-Greenland plate and Eurasia, which initiated the observed pause in the relative convergence of Europe and Africa. The associated stress change in the European continent was significant and explains the sudden termination of a approximately 20-Myr-long contractional intra-plate deformation within Europe, during the late Cretaceous period to the earliest Palaeocene epoch, which was replaced by low-amplitude intra-plate stress-relaxation features. The pre-rupture tectonic stress was large enough to have been responsible for precipitating continental break-up, so there is no need to invoke a thermal mantle plume as a driving mechanism. The model explains the simultaneous timing of several diverse geological events, and shows how the intra-continental stratigraphic record can reveal the timing and dynamics of stress changes, which cannot be resolved by reconstructions based only on plate kinematics.     

Active foundering of a continental arc root beneath the southern Sierra Nevada in California

Seismic data provide images of crust-mantle interactions during ongoing removal of the dense batholithic root beneath the southern Sierra Nevada mountains in California. The removal appears to have initiated between 10 and 3 Myr ago with a Rayleigh-Taylor-type instability, but with a pronounced asymmetric flow into a mantle downwelling (drip) beneath the adjacent Great Valley. A nearly horizontal shear zone accommodated the detachment of the ultramafic root from its granitoid batholith. With continuing flow into the mantle drip, viscous drag at the base of the remaining approximately 35-km-thick crust has thickened the crust by approximately 7 km in a narrow welt beneath the western flank of the range. Adjacent to the welt and at the top of the drip, a V-shaped cone of crust is being dragged down tens of kilometres into the core of the mantle drip, causing the disappearance of the Moho in the seismic images. Viscous coupling between the crust and mantle is therefore apparently driving present-day surface subsidence.     

Evolution and diversity of subduction zones controlled by slab width

Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges.     

Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling

The Colorado plateau is a large, tectonically intact, physiographic province in the southwestern North American Cordillera that stands at ～1,800-2,000?m elevation and has long been thought to be in isostatic equilibrium. The origin of these high elevations is unclear because unlike the surrounding provinces, which have undergone significant Cretaceous-Palaeogene compressional deformation followed by Neogene extensional deformation, the Colorado plateau is largely internally undeformed. Here we combine new seismic tomography and receiver function images to resolve a vertical high-seismic-velocity anomaly beneath the west-central plateau that extends more than 200?km in depth. The upper surface of this anomaly is seismically defined by a dipping interface extending from the lower crust to depths of 70-90?km. The base of the continental crust above the anomaly has a similar shape, with an elevated Moho. We interpret these seismic structures as a continuing regional, delamination-style foundering of lower crust and continental lithosphere. This implies that Pliocene (2.6-5.3?Myr ago) uplift of the plateau and the magmatism on its margins are intimately tied to continuing deep lithospheric processes. Petrologic and geochemical observations indicate that late Cretaceous-Palaeogene (～90-40?Myr ago) low-angle subduction hydrated and probably weakened much of the Proterozoic tectospheric mantle beneath the Colorado plateau. We suggest that mid-Cenozoic (～35-25?Myr ago) to Recent magmatic infiltration subsequently imparted negative compositional buoyancy to the base and sides of the Colorado plateau upper mantle, triggering downwelling. The patterns of magmatic activity suggest that previous such events have progressively removed the Colorado plateau lithosphere inward from its margins, and have driven uplift. Using Grand Canyon incision rates and Pliocene basaltic volcanism patterns, we suggest that this particular event has been active over the past ～6?Myr.