Earthquakes track subduction fluids from slab source to mantle wedge sink.

Subducting plates release fluids as they plunge into Earth's mantle and occasionally rupture to produce intraslab earthquakes. It is debated whether fluids and earthquakes are directly related. By combining seismic observations and geodynamic models from western Greece, and comparing across other subduction zones, we find that earthquakes effectively track the flow of fluids from their slab source at >80 km depth to their sink at shallow (<40 km) depth. Between source and sink, the fluids flow updip under a sealed plate interface, facilitating intraslab earthquakes. In some locations, the seal breaks and fluids escape through vents into the mantle wedge, thereby reducing the fluid supply and seismicity updip in the slab. The vents themselves may represent nucleation sites for larger damaging earthquakes.

Pathway from subducting slab to surface for melt and fluids beneath Mount Rainier.

Convergent margin volcanism originates with partial melting, primarily of the upper mantle, into which the subducting slab descends. Melting of this material can occur in one of two ways. The flow induced in the mantle by the slab can result in upwelling and melting through adiabatic decompression. Alternatively, fluids released from the descending slab through dehydration reactions can migrate into the hot mantle wedge, inducing melting by lowering the solidus temperature. The two mechanisms are not mutually exclusive. In either case, the buoyant melts make their way towards the surface to reside in the crust or to be extruded as lava. Here we use magnetotelluric data collected across the central state of Washington, USA, to image the complete pathway for the fluid-melt phase. By incorporating constraints from a collocated seismic study into the magnetotelluric inversion process, we obtain superior constraints on the fluids and melt in a subduction setting. Specifically, we are able to identify and connect fluid release at or near the top of the slab, migration of fluids into the overlying mantle wedge, melting in the wedge, and transport of the melt/fluid phase to a reservoir in the crust beneath Mt Rainier.

Geophysical detection of relict metasomatism from an Archean (approximately 3.5 Ga) subduction zone.

When plate tectonics started on Earth has been uncertain, and its role in the assembly of early continents is not well understood. By synthesizing coincident seismic and electrical profiles, we show that subduction processes formed the Archean Slave craton in Canada. The spatial overlap between a seismic discontinuity and a conductive anomaly at approximately 100 kilometers depth, in conjunction with the occurrence of mantle xenoliths rich in secondary minerals representative of a metasomatic front, supports cratonic assembly by subduction and accretion of lithospheric fragments. Although evidence of cratonic assembly is rarely preserved, these results suggest that plate tectonics was operating as early as Paleoarchean times, approximately 3.5 billion years ago (Ga).

A sharp lithosphere-asthenosphere boundary imaged beneath eastern North America.

Plate tectonic theory hinges on the concept of a relatively rigid lithosphere moving over a weaker asthenosphere, yet the nature of the lithosphere-asthenosphere boundary remains poorly understood. The gradient in seismic velocity that occurs at this boundary is central to constraining the physical and chemical properties that create differences in mechanical strength between the two layers. For example, if the lithosphere is simply a thermal boundary layer that is more rigid owing to colder temperatures, mantle flow models indicate that the velocity gradient at its base would occur over tens of kilometres. In contrast, if the asthenosphere is weak owing to volatile enrichment or the presence of partial melt, the lithosphere-asthenosphere boundary could occur over a much smaller depth range. Here we use converted seismic phases in eastern North America to image a very sharp seismic velocity gradient at the base of the lithosphere-a 3-11 per cent drop in shear-wave velocity over a depth range of 11 km or less at 90-110 km depth. Such a strong, sharp boundary cannot be reconciled with a purely thermal gradient, but could be explained by an asthenosphere that contains a few per cent partial melt or that is enriched in volatiles relative to the lithosphere.