The role of fluids in lower-crustal earthquakes near continental rifts

The occurrence of earthquakes in the lower crust near continental rifts has long been puzzling, as the lower crust is generally thought to be too hot for brittle failure to occur. Such anomalous events have usually been explained in terms of the lower crust being cooler than normal. But if the lower crust is indeed cold enough to produce earthquakes, then the uppermost mantle beneath it should also be cold enough, and yet uppermost mantle earthquakes are not observed. Numerous lower-crustal earthquakes occur near the southwestern termination of the Taupo Volcanic Zone (TVZ), an active continental rift in New Zealand. Here we present three-dimensional tomographic imaging of seismic velocities and seismic attenuation in this region using data from a dense seismograph deployment. We find that crustal earthquakes accurately relocated with our three-dimensional seismic velocity model form a continuous band along the rift, deepening from mostly less than 10 km in the central TVZ to depths of 30-40 km in the lower crust, 30 km southwest of the termination of the volcanic zone. These earthquakes often occur in swarms, suggesting fluid movement in critically loaded fault zones. Seismic velocities within the band are also consistent with the presence of fluids, and the deepening seismicity parallels the boundary between high seismic attenuation (interpreted as partial melt) within the central TVZ and low seismic attenuation in the crust to the southwest. This linking of upper and lower-crustal seismicity and crustal structure allows us to propose a common explanation for all the seismicity, involving the weakening of faults on the periphery of an otherwise dry, mafic crust by hot fluids, including those exsolved from underlying melt. Such fluids may generally be an important driver of lower-crustal seismicity near continental rifts.     

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.     

Whole-mantle convection and the transition-zone water filter

Because of their distinct chemical signatures, ocean-island and mid-ocean-ridge basalts are traditionally inferred to arise from separate, isolated reservoirs in the Earth's mantle. Such mantle reservoir models, however, typically satisfy geochemical constraints, but not geophysical observations. Here we propose an alternative hypothesis that, rather than being divided into isolated reservoirs, the mantle is filtered at the 410-km-deep discontinuity. We propose that, as the ascending ambient mantle (forced up by the downward flux of subducting slabs) rises out of the high-water-solubility transition zone (between the 660 km and 410 km discontinuities) into the low-solubility upper mantle above 410 km, it undergoes dehydration-induced partial melting that filters out incompatible elements. The filtered, dry and depleted solid phase continues to rise to become the source material for mid-ocean-ridge basalts. The wet, enriched melt residue may be denser than the surrounding solid and accordingly trapped at the 410 km boundary until slab entrainment returns it to the deeper mantle. The filter could be suppressed for both mantle plumes (which therefore generate wetter and more enriched ocean-island basalts) as well as the hotter Archaean mantle (thereby allowing for early production of enriched continental crust). We propose that the transition-zone water-filter model can explain many geochemical observations while avoiding the major pitfalls of invoking isolated mantle reservoirs.     

Low electrical resistivity associated with plunging of the Nazca flat slab beneath Argentina

Beneath much of the Andes, oceanic lithosphere descends eastward into the mantle at an angle of about 30 degrees (ref. 1). A partially molten region is thought to form in a wedge between this descending slab and the overlying continental lithosphere as volatiles given off by the slab lower the melting temperature of mantle material. This wedge is the ultimate source for magma erupted at the active volcanoes that characterize the Andean margin. But between 28 degrees and 33 degrees S the subducted Nazca plate appears to be anomalously buoyant, as it levels out at about 100 km depth and extends nearly horizontally under the continent. Above this 'flat slab', volcanic activity in the main Andean Cordillera terminated about 9 million years ago as the flattening slab presumably squeezed out the mantle wedge. But it is unknown where slab volatiles go once this happens, and why the flat slab finally rolls over to descend steeply into the mantle 600 km further eastward. Here we present results from a magnetotelluric profile in central Argentina, from which we infer enhanced electrical conductivity along the eastern side of the plunging slab, indicative of the presence of partial melt. This conductivity structure may imply that partial melting occurs to at least 250 km and perhaps to more than 400 km depth, or that melt is supplied from the 410 km discontinuity, consistent with the transition-zone 'water-filter' model of Bercovici and Karato.     

Faulting induced by precipitation of water at grain boundaries in hot subducting oceanic crust

Dehydration embrittlement has been proposed to explain both intermediate- and deep-focus earthquakes in subduction zones. Because such earthquakes primarily occur at shallow depths or within the core of the subducting plate, dehydration at relatively low temperatures has been emphasized. However, recent careful relocation of subduction-zone earthquakes shows that at depths of 100-250 km, earthquakes continue in the uppermost part of the slab (probably the former oceanic crust that has been converted to eclogite) where temperatures are higher. Here we show that at such pressures and temperatures, eclogite lacking hydrous phases but with significant hydroxyl incorporated as defects in pyroxene and garnet develops a faulting instability associated with precipitation of water at grain boundaries and the production of very small amounts of melt. This new faulting mechanism satisfactorily explains high-temperature earthquakes in subducting oceanic crust and could potentially be involved in much deeper earthquakes in connection with similar precipitation of water in the mantle transition zone (400-700 km depth). Of potential importance for all proposed high-pressure earthquake mechanisms is the very small amount of fluid required to trigger this instability.     

Electromagnetic detection of a 410-km-deep melt layer in the southwestern United States

A deep-seated melt or fluid layer on top of the 410-km-deep seismic discontinuity in Earth's upper mantle, as proposed in the transition-zone 'water filter' hypothesis, may have significant bearing on mantle dynamics and chemical differentiation. The geophysical detection of such a layer has, however, proved difficult. Magnetotelluric and geomagnetic depth sounding are geophysical methods sensitive to mantle melt. Here we use these methods to search for a distinct structure near 410-km depth. We calculate one-dimensional forward models of the response of electrical conductivity depth profiles, based on mineral physics studies of the effect of incorporating hydrogen in upper-mantle and transition-zone minerals. These models indicate that a melt layer at 410-km depth is consistent with regional magnetotelluric and geomagnetic depth sounding data from the southwestern United States (Tucson). The 410-km-deep melt layer in this model has a conductance of 3.0 x 10(4) S and an estimated thickness of 5-30 km. This is the only regional data set that we have examined for which such a melt layer structure was found, consistent with regional seismic studies. We infer that the hypothesized transition-zone water filter occurs regionally, but that such a layer is unlikely to be a global feature.     

Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode

Seafloor spreading centres show a regular along-axis segmentation thought to be produced by a segmented magma supply in the passively upwelling mantle. On the other hand, continental rifts are segmented by large offset normal faults, and many lack magmatism. It is unclear how, when and where the ubiquitous segmented melt zones are emplaced during the continental rupture process. Between 14 September and 4 October 2005, 163 earthquakes (magnitudes greater than 3.9) and a volcanic eruption occurred within the approximately 60-km-long Dabbahu magmatic segment of the Afar rift, a nascent seafloor spreading centre in stretched continental lithosphere. Here we present a three-dimensional deformation field for the Dabbahu rifting episode derived from satellite radar data, which shows that the entire segment ruptured, making it the largest to have occurred on land in the era of satellite geodesy. Simple elastic modelling shows that the magmatic segment opened by up to 8 m, yet seismic rupture can account for only 8 per cent of the observed deformation. Magma was injected along a dyke between depths of 2 and 9 km, corresponding to a total intrusion volume of approximately 2.5 km3. Much of the magma appears to have originated from shallow chambers beneath Dabbahu and Gabho volcanoes at the northern end of the segment, where an explosive fissural eruption occurred on 26 September 2005. Although comparable in magnitude to the ten year (1975-84) Krafla events in Iceland, seismic data suggest that most of the Dabbahu dyke intrusion occurred in less than a week. Thus, magma intrusion via dyking, rather than segmented normal faulting, maintains and probably initiated the along-axis segmentation along this sector of the Nubia-Arabia plate boundary.     

Mapping the Hawaiian plume conduit with converted seismic waves

The volcanic edifice of the Hawaiian islands and seamounts, as well as the surrounding area of shallow sea floor known as the Hawaiian swell, are believed to result from the passage of the oceanic lithosphere over a mantle hotspot. Although geochemical and gravity observations indicate the existence of a mantle thermal plume beneath Hawaii, no direct seismic evidence for such a plume in the upper mantle has yet been found. Here we present an analysis of compressional-to-shear (P-to-S) converted seismic phases, recorded on seismograph stations on the Hawaiian islands, that indicate a zone of very low shear-wave velocity (< 4 km s(-1)) starting at 130-140 km depth beneath the central part of the island of Hawaii and extending deeper into the upper mantle. We also find that the upper-mantle transition zone (410-660 km depth) appears to be thinned by up to 40-50 km to the south-southwest of the island of Hawaii. We interpret these observations as localized effects of the Hawaiian plume conduit in the asthenosphere and mantle transition zone with excess temperature of approximately 300 degrees C. Large variations in the transition-zone thickness suggest a lower-mantle origin of the Hawaiian plume similar to the Iceland plume, but our results indicate a 100 degrees C higher temperature for the Hawaiian plume.     

A mechanism to thin the continental lithosphere at magma-poor margins

Where continental plates break apart, slip along multiple normal faults provides the required space for the Earth's crust to thin and subside. After initial rifting, however, the displacement on normal faults observed at the sea floor seems not to match the inferred extension. Here we show that crustal thinning can be accomplished in such extensional environments by a system of conjugate concave downward faults instead of multiple normal faults. Our model predicts that these concave faults accumulate large amounts of extension and form a very thin crust (< 10 km) by exhumation of mid-crustal and mantle material. This transitional crust is capped by sub-horizontal detachment surfaces over distances exceeding 100 km with little visible deformation. Our rift model is based on numerical experiments constrained by geological and geophysical observations from the Alpine Tethys and Iberia/Newfoundland margins. Furthermore, we suggest that the observed transition from broadly distributed and symmetric extension to localized and asymmetric rifting is directly controlled by the existence of a strong gabbroic lower crust. The presence of such lower crustal gabbros is well constrained for the Alpine Tethys system. Initial decoupling of upper crustal deformation from lower crustal and mantle deformation by progressive weakening of the middle crust is an essential requirement to reproduce the observed rift evolution. This is achieved in our models by the formation of weak ductile shear zones.     

Mantle compensation of active metamorphic core complexes at Woodlark rift in Papua New Guinea

In many highly extended rifts on the Earth, tectonic removal of the upper crust exhumes mid-crustal rocks, producing metamorphic core complexes. These structures allow the upper continental crust to accommodate tens of kilometres of extension, but it is not clear how the lower crust and underlying mantle respond. Also, despite removal of the upper crust, such core complexes remain both topographically high and in isostatic equilibrium. Because many core complexes in the western United States are underlain by a flat Moho discontinuity, it has been widely assumed that their elevation is supported by flow in the lower crust or by magmatic underplating. These processes should decouple upper-crust extension from that in the mantle. In contrast, here we present seismic observations of metamorphic core complexes of the western Woodlark rift that show the overall crust to be thinned beneath regions of greatest surface extension. These core complexes are actively being exhumed at a rate of 5-10 km Myr(-1), and the thinning of the underlying crust appears to be compensated by mantle rocks of anomalously low density, as indicated by low seismic velocities. We conclude that, at least in this case, the development of metamorphic core complexes and the accommodation of high extension is not purely a crustal phenomenon, but must involve mantle extension.     

A lithospheric instability origin for Columbia River flood basalts and Wallowa Mountains uplift in northeast Oregon

Flood basalts appear to form during the initiation of hotspot magmatism. The Columbia River basalts (CRB) represent the largest volume of flood basalts associated with the Yellowstone hotspot, yet their source appears to be in the vicinity of the Wallowa Mountains, about 500 km north of the projected hotspot track. These mountains are composed of a large granitic pluton intruded into a region of oceanic lithosphere affinity. The elevation of the interface between Columbia River basalts and other geological formations indicates that mild pre-eruptive subsidence took place in the Wallowa Mountains, followed by syn-eruptive uplift of several hundred metres and a long-term uplift of about 2 km. The mapped surface uplift mimics regional topography, with the Wallowa Mountains in the centre of a 'bull's eye' pattern of valleys and low-elevation mountains. Here we present the seismic velocity structure of the mantle underlying this region and erosion-corrected elevation maps of lava flows, and show that an area of reduced mantle melt content coincides with the 200-km-wide topographic uplift. We conclude that convective downwelling and detachment of a compositionally dense plutonic root can explain the timing and magnitude of Columbia River basalt magmatism, as well as the surface uplift and existence of the observed melt-depleted mantle.     

Geophysical evidence for reduced melt production on the Arctic ultraslow Gakkel mid-ocean ridge

Most models of melt generation beneath mid-ocean ridges predict significant reduction of melt production at ultraslow spreading rates (full spreading rates &<20 mm x yr(-1)) and consequently they predict thinned oceanic crust. The 1,800-km-long Arctic Gakkel mid-ocean ridge is an ideal location to test such models, as it is by far the slowest portion of the global mid-ocean-ridge spreading system, with a full spreading rate ranging from 6 to 13 mm x yr(-1) (refs 4, 5). Furthermore, in contrast to some other ridge systems, the spreading direction on the Gakkel ridge is not oblique and the rift valley is not offset by major transform faults. Here we present seismic evidence for the presence of exceptionally thin crust along the Gakkel ridge rift valley with crustal thicknesses varying between 1.9 and 3.3 km (compared to the more usual value of 7 km found on medium- to fast-spreading mid-ocean ridges). Almost 8,300 km of closely spaced aeromagnetic profiles across the rift valley show the presence of discrete volcanic centres along the ridge, which we interpret as evidence for strongly focused, three-dimensional magma supply. The traces of these eruptive centres can be followed to crustal ages of approximately 25 Myr off-axis, implying that these magma production and transport systems have been stable over this timescale.     

Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth's upper mantle

The mineral olivine dominates the composition of the Earth's upper mantle and hence controls its mechanical behaviour and seismic anisotropy. Experiments at high temperature and moderate pressure, and extensive data on naturally deformed mantle rocks, have led to the conclusion that olivine at upper-mantle conditions deforms essentially by dislocation creep with dominant [100] slip. The resulting crystal preferred orientation has been used extensively to explain the strong seismic anisotropy observed down to 250 km depth. The rapid decrease of anisotropy below this depth has been interpreted as marking the transition from dislocation to diffusion creep in the upper mantle. But new high-pressure experiments suggest that dislocation creep also dominates in the lower part of the upper mantle, but with a different slip direction. Here we show that this high-pressure dislocation creep produces crystal preferred orientations resulting in extremely low seismic anisotropy, consistent with seismological observations below 250 km depth. These results raise new questions about the mechanical state of the lower part of the upper mantle and its coupling with layers both above and below.     

Geochemical tracing of Pacific-to-Atlantic upper-mantle flow through the Drake passage

The Earth's convecting upper mantle can be viewed as comprising three main reservoirs, beneath the Pacific, Atlantic and Indian oceans. Because of the uneven global distribution and migration of ridges and subduction zones, the surface area of the Pacific reservoir is at present contracting at about 0.6 km2 x y(r-1), while the Atlantic and Indian reservoirs are growing at about 0.45 km2 x yr(-1) and 0.15 km2 x yr(-1), respectively. Garfunkel and others have argued that there must accordingly be net mantle flow from the Pacific to the Atlantic and Indian reservoirs (in order to maintain mass balance), and Alvarez further predicted that this flow should be restricted to the few parts of the Pacific rim (here termed 'gateways') where there are no continental roots or subduction zones that might act as barriers to shallow mantle flow. The main Pacific gateways are, according to Alvarez, the southeast Indian Ocean, the Caribbean Sea and the Drake passage. Here we report geochemical data which confirm that there has been some outflow of Pacific mantle into the Drake passage--but probably in response to regional tectonic constraints, rather than global mass-balance requirements. We also show that a mantle domain boundary, equivalent to the Australian-Antarctic discordance, must lie between the Drake passage and the east Scotia Sea.     

Variation in styles of rifting in the Gulf of California

Constraints on the structure of rifted continental margins and the magmatism resulting from such rifting can help refine our understanding of the strength of the lithosphere, the state of the underlying mantle and the transition from rifting to seafloor spreading. An important structural classification of rifts is by width, with narrow rifts thought to form as necking instabilities (where extension rates outpace thermal diffusion) and wide rifts thought to require a mechanism to inhibit localization, such as lower-crustal flow in high heat-flow settings. Observations of the magmatism that results from rifting range from volcanic margins with two to three times the magmatism predicted from melting models to non-volcanic margins with almost no rift or post-rift magmatism. Such variations in magmatic activity are commonly attributed to variations in mantle temperature. Here we describe results from the PESCADOR seismic experiment in the southern Gulf of California and present crustal-scale images across three rift segments. Over short lateral distances, we observe large differences in rifting style and magmatism--from wide rifting with minor synchronous magmatism to narrow rifting in magmatically robust segments. But many of the factors believed to control structural evolution and magmatism during rifting (extension rate, mantle potential temperature and heat flow) tend to vary over larger length scales. We conclude instead that mantle depletion, rather than low mantle temperature, accounts for the observed wide, magma-poor margins, and that mantle fertility and possibly sedimentary insulation, rather than high mantle temperature, account for the observed robust rift and post-rift magmatism.     

南海成因：岩石圈破裂与俯冲带相互作用新认识

 南海深部计划与国际大洋钻探航次取得了一系列创新进展与重大突破：1）发现南海陆缘岩石圈减薄之初未出现地幔蛇纹岩出露，且岩浆迅速出现；2）新提出南海不是“小大西洋”，而是“板缘张裂”盆地，与经典的大西洋型“板内张裂”陆缘模式不同；3）揭示南海受到俯冲带的强烈控制，提出俯冲诱发地幔上涌并影响南海岩浆活动。     

Deep crustal structure of Baiyun Sag,northern South China Sea revealed from deep seismic reflection profile

This paper discusses deep crustal architecture of the Baiyun Sag of the Pearl River Mouth Basin, northern South China Sea based on velocity analysis, time-depth conversion and seismic interpretation of the deep seismic reflection profile DSRP-2002. The profile was acquired and processed to 14 S TWT by the China National Offshore Oil Corp. (CNOOC) in 2002. It extends across the Baiyun Sag of the Pearl River Mouth Basin, from the northern continental shelf of the SCS to the deepwater province. As the first deep seismic reflection profile in the Pearl River Mouth Basin, this profile reveals seismic phases from basement down to upper most mantle. The Moho surface appears in the profile as an undulating layer of varying thickness of 1-3 km. It is not a single reflector interface, but a velocity gradient or interconversion layer. The crust thins stepwisely from the shelf to the continental slope and the abyssal plain (from north to south), and also thins under depocenters. The crustal thickness is only 7 km in the depocenter of the main Baiyun Sag, which corresponds to a Moho upwelling mirroring the basement topography. In the lower slope and the ocean-continental transition zone of the southernmost portion of the profile, three sub-parallel, NW-dipping strong reflectors found at depths around 10--21 km are interpreted as indications of a subducted Mesozoic oceanic crust. Crustal faults exist in the northern and southern boundaries of the Baiyun Sag. The intense and persistent subsidence of the Baiyun Sag might be related to the long-term activity of the crustal faults.     

Sound velocities of majorite garnet and the composition of the mantle transition region

The composition of the mantle transition region, characterized by anomalous seismic-wave velocity and density changes at depths of approximately 400 to 700 km, has remained controversial. Some have proposed that the mantle transition region has an olivine-rich 'pyrolite' composition, whereas others have inferred that it is characterized by pyroxene- and garnet-rich compositions ('piclogite'), because the sound velocities in pyrolite estimated from laboratory data are substantially higher than those seismologically observed. Although the velocities of the olivine polymorphs at these pressures (wadsleyite and ringwoodite) have been well documented, those of majorite (another significant high-pressure phase in the mantle transition region) with realistic mantle compositions have never been measured. Here we use combined in situ X-ray and ultrasonic measurements under the pressure and temperature conditions of the mantle transition region to show that majorite in a pyrolite composition has sound velocities substantially lower than those of earlier estimates, owing to strong nonlinear decreases at high temperature, particularly for shear-wave velocity. We found that pyrolite yields seismic velocities more consistent with typical seismological models than those of piclogite in the upper to middle parts of the region, except for the potentially larger velocity jumps in pyrolite relative to those observed at a depth of 410 km. In contrast, both of these compositions lead to significantly low shear-wave velocities in the lower part of the region, suggesting possible subadiabatic temperatures or the existence of a layer of harzburgite-rich material supplied by the subducted slabs stagnant at these depths.     

Decreasing overflow from the Nordic seas into the Atlantic Ocean through the Faroe Bank channel since 1950.

The overflow of cold, dense water from the Nordic seas, across the Greenland-Scotland ridge and into the Atlantic Ocean is the main source for the deep water of the North Atlantic Ocean. This flow also helps drive the inflow of warm, saline surface water into the Nordic seas. The Faroe Bank channel is the deepest path across the ridge, and the deep flow through this channel accounts for about one-third of the total overflow. Previous work has demonstrated that the overflow has become warmer and less saline over time. Here we show, using direct measurements and historical hydrographic data, that the volume flux of the Faroe Bank channel overflow has also decreased. Estimating the volume flux conservatively, we find a decrease by at least 20 per cent relative to 1950. If this reduction in deep flow from the Nordic seas is not compensated by increased flow from other sources, it implies a weakened global thermohaline circulation and reduced inflow of Atlantic water to the Nordic seas.     

Multiple seismic discontinuities near the base of the transition zone in the Earth's mantle

The seismologically defined boundary between the transition zone in the Earth's mantle (410-660 km depth) and the underlying lower mantle is generally interpreted to result from the breakdown of the gamma-spinel phase of olivine to magnesium-perovskite and magnesiowustite. Laboratory measurements of these transformations of olivine have determined that the phase boundary has a negative Clapeyron slope and does indeed occur near pressures corresponding to the base of the transition zone. But a computational study has indicated that, because of the presence of garnet minerals, multiple seismic discontinuities might exist near a depth of 660 km (ref. 4), which would alter the simple negative correlation of changes in temperature with changes in the depth of the phase boundary. In particular, garnet minerals undergo exothermic transformations near this depth, acting to complicate the phase relations and possibly effecting mantle convection processes in some regions. Here we present seismic evidence that supports the existence of such multiple transitions near a depth of 660 km beneath southern California. The observations are consistent with having been generated by garnet transformations coupling with the dissociation of the gamma-spinel phase of olivine. Temperature anomalies calculated from the imaged discontinuity depths--using Clapeyron slopes determined for the various transformations--generally match those predicted from an independent P-wave velocity model of the region.