Low strength of deep San Andreas fault gouge from SAFOD core

The San Andreas fault accommodates 28-34 mm?yr(-1) of right lateral motion of the Pacific crustal plate northwestward past the North American plate. In California, the fault is composed of two distinct locked segments that have produced great earthquakes in historical times, separated by a 150-km-long creeping zone. The San Andreas Fault Observatory at Depth (SAFOD) is a scientific borehole located northwest of Parkfield, California, near the southern end of the creeping zone. Core was recovered from across the actively deforming San Andreas fault at a vertical depth of 2.7 km (ref. 1). Here we report laboratory strength measurements of these fault core materials at in situ conditions, demonstrating that at this locality and this depth the San Andreas fault is profoundly weak (coefficient of friction, 0.15) owing to the presence of the smectite clay mineral saponite, which is one of the weakest phyllosilicates known. This Mg-rich clay is the low-temperature product of metasomatic reactions between the quartzofeldspathic wall rocks and serpentinite blocks in the fault. These findings provide strong evidence that deformation of the mechanically unusual creeping portions of the San Andreas fault system is controlled by the presence of weak minerals rather than by high fluid pressure or other proposed mechanisms. The combination of these measurements of fault core strength with borehole observations yields a self-consistent picture of the stress state of the San Andreas fault at the SAFOD site, in which the fault is intrinsically weak in an otherwise strong crust.     

Talc-bearing serpentinite and the creeping section of the San Andreas fault

The section of the San Andreas fault located between Cholame Valley and San Juan Bautista in central California creeps at a rate as high as 28 mm yr(-1) (ref. 1), and it is also the segment that yields the best evidence for being a weak fault embedded in a strong crust. Serpentinized ultramafic rocks have been associated with creeping faults in central and northern California, and serpentinite is commonly invoked as the cause of the creep and the low strength of this section of the San Andreas fault. However, the frictional strengths of serpentine minerals are too high to satisfy the limitations on fault strength, and these minerals also have the potential for unstable slip under some conditions. Here we report the discovery of talc in cuttings of serpentinite collected from the probable active trace of the San Andreas fault that was intersected during drilling of the San Andreas Fault Observatory at Depth (SAFOD) main hole in 2005. We infer that the talc is forming as a result of the reaction of serpentine minerals with silica-saturated hydrothermal fluids that migrate up the fault zone, and the talc commonly occurs in sheared serpentinite. This discovery is significant, as the frictional strength of talc at elevated temperatures is sufficiently low to meet the constraints on the shear strength of the fault, and its inherently stable sliding behaviour is consistent with fault creep. Talc may therefore provide the connection between serpentinite and creep in the San Andreas fault, if shear at depth can become localized along a talc-rich principal-slip surface within serpentinite entrained in the fault zone.     

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.     

Migration of seismic scatterers associated with the 1993 Parkfield aseismic transient event

The time-varying deformation field within a fault zone, particularly at depths where earthquakes occur, is important for understanding fault behaviour and its relation to earthquake occurrence. But detection of this temporal variation has been extremely difficult, although laboratory studies have long suggested that certain structural changes, such as the properties of crustal fractures, should be seismically detectable. Here we present evidence that such structural changes are indeed observable. In particular, we find a systematic temporal variation in the seismograms of repeat microearthquakes that occurred on the Parkfield segment of the San Andreas fault over the decade 1987-97. Our analysis reveals a change of the order of 10 m in the location of scatterers which plausibly lie within the fault zone at a depth of approximately 3 km. The motion of the scatterers is coincident, in space and time, with the onset of a well documented aseismic transient (deformation event). We speculate that this structural change is the result of a stress-induced redistribution of fluids in fluid-filled fractures caused by the transient event.     

Implications for prediction and hazard assessment from the 2004 Parkfield earthquake

Obtaining high-quality measurements close to a large earthquake is not easy: one has to be in the right place at the right time with the right instruments. Such a convergence happened, for the first time, when the 28 September 2004 Parkfield, California, earthquake occurred on the San Andreas fault in the middle of a dense network of instruments designed to record it. The resulting data reveal aspects of the earthquake process never before seen. Here we show what these data, when combined with data from earlier Parkfield earthquakes, tell us about earthquake physics and earthquake prediction. The 2004 Parkfield earthquake, with its lack of obvious precursors, demonstrates that reliable short-term earthquake prediction still is not achievable. To reduce the societal impact of earthquakes now, we should focus on developing the next generation of models that can provide better predictions of the strength and location of damaging ground shaking.     

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.     

利用接收函数方法研究大盈江断裂两侧S波速度结构

 利用研究区(24.2°~25.2°N,97.5°~98.5°E)内大盈江断裂两侧5个流动数字地震台站记录到的宽频带远震P波波形数据进行接收函数反演,得到台站下方0~100km深度范围内地壳、上地幔S波速度细结构.结果表明:研究区内,以大盈江断裂为界,其西北侧Moho面深度约为38km;东南侧Moho面深度为40~42km.断裂两侧地壳、上地幔S波速度结构存在显著差异,东南侧台站下方地壳和上地幔均存在大范围低速区;西北侧台站下方地壳内存在低速层,而上地幔中无明显低速层.研究区内的S波速度结构存在明显的横向非均匀性.     

Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system

The San Andreas fault in California is a mature continental transform fault that accommodates a significant fraction of motion between the North American and Pacific plates. The two most recent great earthquakes on this fault ruptured its northern and central sections in 1906 and 1857, respectively. The southern section of the fault, however, has not produced a great earthquake in historic times (for at least 250 years). Assuming the average slip rate of a few centimetres per year, typical of the rest of the San Andreas fault, the minimum amount of slip deficit accrued on the southern section is of the order of 7-10 metres, comparable to the maximum co-seismic offset ever documented on the fault. Here I present high-resolution measurements of interseismic deformation across the southern San Andreas fault system using a well-populated catalogue of space-borne synthetic aperture radar data. The data reveal a nearly equal partitioning of deformation between the southern San Andreas and San Jacinto faults, with a pronounced asymmetry in strain accumulation with respect to the geologically mapped fault traces. The observed strain rates confirm that the southern section of the San Andreas fault may be approaching the end of the interseismic phase of the earthquake cycle.     

Evidence of crustal ＇channel flow＇ in the eastern margin of Tibetan Plateau from MT measurements

Magnetotelluric （MT） survey has been carried out in the eastern margin of the Tibetan Plateau and its neighboring Shimian-Leshan area, Sichuan Province. Analysis of this MT data reveals that the electric structure of the Tibetan Plateau differ much from that of the Sichuan block. In general, the electric resistivity of crust beneath the Sichuan block in the east is larger than that of the eastern margin of the Tibetan Plateau in the west. The crust of the plateau is divided into upper, middle, and lower layers. The middle crust is a low resistivity layer with minimum down to 3-10Ωm about 10-15 km thick. It presumably contains partial melt and/or salt-bearing fluids with low viscosity, prone to deform and flow, producing a ＂channel flow＂ under the southeastward squeeze of the eastern Tibetan Plateau. This low-resistivity layer makes the upper crust decoupled mechanically from the lower crust. In the brittle upper crust, faults are dominated by left-lateral strike-slip and thrust motions, leading to surface rising and shallow earthquakes. The low-resistivity layer also cut the Xianshuihe-Anninghe fault zone into two sections vertically. In this region, the thicknesses of upper, middle, and lower crust vary laterally, producing a transitional zone in the eastern margin of the Tibetan Plateau characterized by thicker crust and higher elevation in the west and thinner crust and lower elevation in the east.     

Slip on 'weak' faults by the rotation of regional stress in the fracture damage zone

Slip on unfavourably oriented faults with respect to a remotely applied stress is well documented and implies that faults such as the San Andreas fault and low-angle normal faults are weak when compared to laboratory-measured frictional strength. If high pore pressure within fault zones is the cause of such weakness, then stress reorientation within or close to a fault is necessary to allow sufficient fault weakening without the occurrence of hydrofracture. From field observations of a major tectonic fault, and using laboratory experiments and numerical modelling, here we show that stress rotation occurs within the fractured damage zone surrounding faults. In particular, we find that stress rotation is considerable for unfavourably oriented 'weak' faults. In the 'weak' fault case, the damage-induced change in elastic properties provides the necessary stress rotation to allow high pore pressure faulting without inducing hydrofracture.     

Numerical simulation of GPS observed clockwise rotation around the eastern Himalayan syntax in the Tibetan Plateau

From Global Position System (GPS) measurements, there is a clockwise rotation around the eastern Himalayan syntax in the Tibetan Plateau. This phenomenon is difficult to be interpreted by simple two-dimensional modeling from a geodynamic point of view. Because of the extremely thick crust and the lower crust with relatively high temperature in the Tibetan Plateau, the lithospheric rheology in Tibet and surrounding areas present a complex structure. In general, the tectonic structure of the Tibetan Plateau consists of brittle upper crust, ductile lower crust, high viscosity lithospheric upper mantle, and low viscosity asthenosphere, the same as the case in many other continental regions. However, the lower crust in the Tibetan Plateau is much more ductile with a lower viscosity than those of its surroundings at the same depth, and the effective viscosity is low along the collision fault zone. In this study, we construct a three-dimensional Maxwell visco-elastic model in spherical coordinate system, and simulate the deformation process of the Tibetan Plateau driven by a continuous push from the Indian plate. The results show that the existence of the soft lower crust under the plateau makes the entire plateau uplift as a whole, and the Himalayas and the eastern Himalayan syntax uplift faster. Since the lower crust of surrounding blocks is harder except in the southeastern corner where the high-temperature material is much softer and forms an exit channel for material transfer, after the whole plateau reaches a certain height, the lower crustal and upper mantle material begins to move eastward or southeastward and drag the upper crust to behave same way. Thus, from the macroscopic point of view, a relative rigid motion of the plateau with a clockwise rotation around the eastern Himalayan syntax is developed. Supported by Knowledge Innovation Project of the Chinese Academy of Sciences (Grant No. KZCX2-YW-123) and National Natural Science Foundation of China (Grant Nos. 40774048 and 90814014)     

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.     

郯庐断裂带中生代火山活动与深部过程

郯庐断裂带中生代火山活动主要发生在早白垩世 ,形成了一套以橄榄安粗岩系为主的火山岩。火山岩岩石化学特征、地球化学特征与断裂带外明显不同 ,以富碱富钾、富轻稀土元素为特征 ,具有较高的 ISr值和明显偏低的εNd值 ,表明火山活动不仅与地壳混染有关 ,而且与深部地幔以及断裂带的演化有密切的关系。研究显示该区岩浆起源于地幔底侵、壳幔相互作用下的壳幔过渡带 ,也反映此时的断裂带切入到壳幔边界。     

Iron-silica interaction at extreme conditions and the electrically conducting layer at the base of Earth's mantle

The boundary between the Earth's metallic core and its silicate mantle is characterized by strong lateral heterogeneity and sharp changes in density, seismic wave velocities, electrical conductivity and chemical composition. To investigate the composition and properties of the lowermost mantle, an understanding of the chemical reactions that take place between liquid iron and the complex Mg-Fe-Si-Al-oxides of the Earth's lower mantle is first required. Here we present a study of the interaction between iron and silica (SiO2) in electrically and laser-heated diamond anvil cells. In a multianvil apparatus at pressures up to 140 GPa and temperatures over 3,800 K we simulate conditions down to the core-mantle boundary. At high temperature and pressures below 40 GPa, iron and silica react to form iron oxide and an iron-silicon alloy, with up to 5 wt% silicon. At pressures of 85-140 GPa, however, iron and SiO2 do not react and iron-silicon alloys dissociate into almost pure iron and a CsCl-structured (B2) FeSi compound. Our experiments suggest that a metallic silicon-rich B2 phase, produced at the core-mantle boundary (owing to reactions between iron and silicate), could accumulate at the boundary between the mantle and core and explain the anomalously high electrical conductivity of this region.     

Elastic wave velocity in rocks from Dabieshan and its constraints for lithospheric composition and crust-mantle recycling

P-and S-wave velocities in eclogites and granulites from the Dabie ultrahigh pressures (UHP) meta morphic belt, China, were measured at room temperature under the hydrostatic pressures up to 1.0 GPa. The ultrahigh pressure eclogites had the highest densities (3.3 ～ 3.6 g.cm-3) , high velocities and the lowest anisotropy (1.4%c ～ 2.6 % ) . The lowest densities (2.8 ～ 3.1 g. cm-3 ) and the highest Poisson' s ratios (0.28 ～ 0.29) were found in gran ulites, whereas the strongest anisotropies (6. 1 % ～8.4% ) were found in the high-pressure (HP) eclogites. Compari son of the velocities in rocks with that observed in the deep seismic sounding profile crossing Dabieshan suggests that e clogites might exist in the lower crust of Dabieshan, but the quantity might be small. The upper mantle has very similar velocities as the UHP eclogites and serpentinizated/water-bearing dunite. The formation of eclogite represents the crust mantle recycling processes. Crustal material is delaminated and sinks into the mantle by way of eclogite, whereas only a small part of the eclogite could return to the crust.     

The effect of energy feedbacks on continental strength

The classical strength profile of continents is derived from a quasi-static view of their rheological response to stress--one that does not consider dynamic interactions between brittle and ductile layers. Such interactions result in complexities of failure in the brittle-ductile transition and the need to couple energy to understand strain localization. Here we investigate continental deformation by solving the fully coupled energy, momentum and continuum equations. We show that this approach produces unexpected feedback processes, leading to a significantly weaker dynamic strength evolution. In our model, stress localization focused on the brittle-ductile transition leads to the spontaneous development of mid-crustal detachment faults immediately above the strongest crustal layer. We also find that an additional decoupling layer forms between the lower crust and mantle. Our results explain the development of decoupling layers that are observed to accommodate hundreds of kilometres of horizontal motions during continental deformation.     

Lithospheric thermal-rheological structures of the continental margin in the northern South China Sea

Thermal structures of three deep seismic profiles in the continental margin in the northern South China Sea are calculated, their "thermal" lithospheric thicknesses are evaluated based on the basalt dry solidus, and their rheological structures are evaluated with linear frictional failure criterion and power-law creep equation. "Thermal" lithosphere is about 90 km in thickness in shelf area, and thins toward the slope, lowers to 60-65 km in the lower slope, ocean crust and Xisha Trough. In the mid-west of the studied area, the lithospheric rheological structure in shelf area and Xisha Islands is of four layers: brittle, ductile, brittle and ductile. Because of uprising of heat mantle and thinning of crust and lithosphere in Xisha Trough, the bottom of the upper brittle layer is only buried at 16 km. In the eastern area, the bottom of the upper brittle layer in the north is buried at 20 km or so, while in lower slope and ocean crust, the rheological structure is of two layers of brittle and ductile, and crust and uppermost mantle form one whole brittle layer whose bottom is buried at 30-32 km. Analyses show that the characteristics of rheological structure accord with the seismic result observed. The character of rheological stratification implies that before the extension of the continent margin, there likely was a ductile layer in mid-lower crust. The influence of the existence of ductile layer to the evolution of the continent margin and the different extensions of ductile layer and brittle layer should not be overlooked. Its thickness, depth and extent in influencing continent margin's extension and evolution should be well evaluated in building a dynamic model for the area.     

Crustal P-wave velocity structure in Lower Yangtze region： Reinterpretation of Fuliji-Fengxian deep seismic sounding profile

A deep seismic sounding profile in this paper, from Fuliji in Anhui Province to Fengxian of Shanghai City, is located at eastern China (Fig. 1). The field work was jointly accomplished by the Chinese Geological and Mineral Bureau, the China Seismological …     

Seismic detection of folded, subducted lithosphere at the core-mantle boundary

Seismic tomography has been used to infer that some descending slabs of oceanic lithosphere plunge deep into the Earth's lower mantle. The fate of these slabs has remained unresolved, but it has been postulated that their ultimate destination is the lowermost few hundred kilometres of the mantle, known as the D' region. Relatively cold slab material may account for high seismic velocities imaged in D' beneath areas of long-lived plate subduction, and for reflections from a seismic velocity discontinuity just above the anomalously high wave speed regions. The D' discontinuity itself is probably the result of a phase change in relatively low-temperature magnesium silicate perovskite. Here, we present images of the D' region beneath the Cocos plate using Kirchhoff migration of horizontally polarized shear waves, and find a 100-km vertical step occurring over less than 100 km laterally in an otherwise flat D' shear velocity discontinuity. Folding and piling of a cold slab that has reached the core-mantle boundary, as observed in numerical and experimental models, can account for the step by a 100-km elevation of the post-perovskite phase boundary due to a 700 degrees C lateral temperature reduction in the folded slab. We detect localized low velocities at the edge of the slab material, which may result from upwellings caused by the slab laterally displacing a thin hot thermal boundary layer.     

Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120-180 km depth

Fluids and melts liberated from subducting oceanic crust recycle lithophile elements back into the mantle wedge, facilitate melting and ultimately lead to prolific subduction-zone arc volcanism. The nature and composition of the mobile phases generated in the subducting slab at high pressures have, however, remained largely unknown. Here we report direct LA-ICPMS measurements of the composition of fluids and melts equilibrated with a basaltic eclogite at pressures equivalent to depths in the Earth of 120-180 km and temperatures of 700-1,200 degrees C. The resultant liquid/mineral partition coefficients constrain the recycling rates of key elements. The dichotomy of dehydration versus melting at 120 km depth is expressed through contrasting behaviour of many trace elements (U/Th, Sr, Ba, Be and the light rare-earth elements). At pressures equivalent to 180 km depth, however, a supercritical liquid with melt-like solubilities for the investigated trace elements is observed, even at low temperatures. This mobilizes most of the key trace elements (except the heavy rare-earth elements, Y and Sc) and thus limits fluid-phase transfer of geochemical signatures in subduction zones to pressures less than 6 GPa.