Paul A. Bedrosian

Along strike variations in the electrical structure of the San Andreas Fault at Parkfield, California

Magnetotelluric exploration has been used to image along strike variations in the electrical resistivity structure of the San Andreas Fault at Parkfield, California. A low resistivity wedge extending to a depth of several kilometers is continuous over a horizontal distance of 8 km. The base of the wedge is coincident with the shallowest microearthquakes. A change in the electrical and fluid connection of the San Andreas Fault with a low resistivity zone in the Franciscan formation is observed along the Parkfield segment.

Correlation between deep fluids, tremor and creep along the central San Andreas fault

The seismicity pattern along the San Andreas fault near Parkfield and Cholame, California, varies distinctly over a length of only fifty kilometres. Within the brittle crust, the presence of frictionally weak minerals, fault-weakening high fluid pressures and chemical weakening are considered possible causes of an anomalously weak fault northwest of Parkfield. Non-volcanic tremor from lower-crustal and upper-mantle depths is most pronounced about thirty kilometres southeast of Parkfield and is thought to be associated with high pore-fluid pressures at depth. Here we present geophysical evidence of fluids migrating into the creeping section of the San Andreas fault that seem to originate in the region of the uppermost mantle that also stimulates tremor, and evidence that along-strike variations in tremor activity and amplitude are related to strength variations in the lower crust and upper mantle. Interconnected fluids can explain a deep zone of anomalously low electrical resistivity that has been imaged by magnetotelluric data southwest of the Parkfield–Cholame segment. Near Cholame, where fluids seem to be trapped below a high-resistivity cap, tremor concentrates adjacent to the inferred fluids within a mechanically strong zone of high resistivity. By contrast, subvertical zones of low resistivity breach the entire crust near the drill hole of the San Andreas Fault Observatory at Depth, northwest of Parkfield, and imply pathways for deep fluids into the eastern fault block, coincident with a mechanically weak crust and the lower tremor amplitudes in the lower crust. Fluid influx to the fault system is consistent with hypotheses of fault-weakening high fluid pressures in the brittle crust.

Joint inversion for Vp, Vs, and Vp/Vs at SAFOD, Parkfield, California

We refined the three-dimensional (3-D) Vp, Vs and Vp/Vs models around the San Andreas Fault Observatory at Depth (SAFOD) site using a new double-difference (DD) seismic tomography code (tomoDDPS) that simultaneously solves for earthquake locations and all three velocity models using both absolute and differential P, S, and S-P times. This new method is able to provide a more robust Vp/Vs model than that from the original DD tomography code (tomoDD), obtained simply by dividing Vp by Vs. For the new inversion, waveform cross-correlation times for earthquakes from 2001 to 2002 were also used, in addition to arrival times from earthquakes and explosions in the region. The Vp values extracted from the model along the SAFOD trajectory match well with the borehole log data, providing in situ confirmation of our results. Similar to previous tomographic studies, the 3-D structure around Parkfield is dominated by the velocity contrast across the San Andreas Fault (SAF). In both the Vp and Vs models, there is a clear low-velocity zone as deep as 7 km along the SAF trace, compatible with the findings from fault zone guided waves. There is a high Vp/Vs anomaly zone on the southwest side of the SAF trace that is about 1–2 km wide and extends as deep as 4 km, which is interpreted to be due to fluids and fractures in the package of sedimentary rocks abutting the Salinian basement rock to the southwest. The relocated earthquakes align beneath the northeast edge of this high Vp/Vs zone. We carried out a 2-D correlation analysis for an existing resistivity model and the corresponding profiles through our model, yielding a classification that distinguishes several major lithologies.

Anatomy of the Dead Sea Transform from lithospheric to microscopic scale

Fault zones are the locations where motion of tectonic plates, often associated with earthquakes, is accommodated. Despite a rapid increase in the understanding of faults in the last decades, our knowledge of their geometry, petrophysical properties, and controlling processes remains incomplete. The central questions addressed here in our study of the Dead Sea Transform (DST) in the Middle East are as follows: (1) What are the structure and kinematics of a large fault zone? (2) What controls its structure and kinematics? (3) How does the DST compare to other plate boundary fault zones? The DST has accommodated a total of 105 km of left-lateral transform motion between the African and Arabian plates since early Miocene (similar to 20 Ma). The DST segment between the Dead Sea and the Red Sea, called the Arava/Araba Fault (AF), is studied here using a multidisciplinary and multiscale approach from the mu m to the plate tectonic scale. We observe that under the DST a narrow, subvertical zone cuts through crust and lithosphere. First, from west to east the crustal thickness increases smoothly from 26 to 39 km, and a subhorizontal lower crustal reflector is detected east of the AF. Second, several faults exist in the upper crust in a 40 km wide zone centered on the AF, but none have kilometer-size zones of decreased seismic velocities or zones of high electrical conductivities in the upper crust expected for large damage zones. Third, the AF is the main branch of the DST system, even though it has accommodated only a part (up to 60 km) of the overall 105 km of sinistral plate motion. Fourth, the AF acts as a barrier to fluids to a depth of 4 km, and the lithology changes abruptly across it. Fifth, in the top few hundred meters of the AF a locally transpressional regime is observed in a 100-300 m wide zone of deformed and displaced material, bordered by subparallel faults forming a positive flower structure. Other segments of the AF have a transtensional character with small pull-aparts along them. The damage zones of the individual faults are only 5-20 m wide at this depth range. Sixth, two areas on the AF show mesoscale to microscale faulting and veining in limestone sequences with faulting depths between 2 and 5 km. Seventh, fluids in the AF are carried downward into the fault zone. Only a minor fraction of fluids is derived from ascending hydrothermal fluids. However, we found that on the kilometer scale the AF does not act as an important fluid conduit. Most of these findings are corroborated using thermomechanical modeling where shear deformation in the upper crust is localized in one or two major faults; at larger depth, shear deformation occurs in a 20-40 km wide zone with a mechanically weak decoupling zone extending subvertically through the entire lithosphere.

Electrical conductivity images of active and fossil fault zones

Abstract We compare recent magnetotelluric investigations of four large fault systems: (i) the actively deforming, ocean-continent interplate San Andreas Fault (SAF); (ii) the actively deforming, continent-continent interplate Dead Sea Transform (DST); (iii) the currently inactive, trench-linked intraplate West Fault (WF) in northern Chile; and (iv) the Waterberg Fault/Omaruru Lineament (WF/OL) in Namibia, a fossilized intraplate shear zone formed during early Proterozoic continental collision. These fault zones show both similarities and marked differences in their electrical subsurface structure. The central segment of the SAF is characterized by a zone of high conductivity extending to a depth of several kilometres and attributed to fluids within a highly fractured damage zone. The WF exhibits a less pronounced but similar fault-zone conductor (FZC) that can be explained by meteoric waters entering the fault zone. The DST appears different as it shows a distinct lack of a FZC and seems to act primarily as an impermeable barrier to cross-fault fluid transport. Differences in the electrical structure of these faults within the upper crust may be linked to the degree of deformation localization within the fault zone. At the DST, with no observable fault-zone conductor, strain may have been localized for a considerable time span along a narrow, metre-scale damage zone with a sustained strength difference between the shear plane and the surrounding host rock. In the case of the SAF, a positive correlation of conductance and fault activity is observed, with more active fault segments associated with wider, deeper and more conductive fault-zone anomalies. Fault-zone conductors, however, do not uniquely identify specific architectural or hydrological units of a fault. A more comprehensive whole-fault picture for the brittle crust can be developed in combination with seismicity and structural information. Giving a window into lower-crustal shear zones, the fossil WF/OL in Namibia is imaged as a subvertical, 14 km-deep, 10 km-wide zone of high and anisotropic conductivity. The present level of exhumation suggests that the WF/OL penetrated the entire crust as a relatively narrow shear zone. Contrary to the fluid-driven conductivity anomalies of active faults, the anomaly here is attributed to graphitic enrichment along former shear planes. Once created, graphite is stable over very long time spans and thus fault/shear zones may remain conductive long after activity ceases.

Segmentation of plate coupling, fate of subduction fluids, and modes of arc magmatism in Cascadia, inferred from magnetotelluric resistivity

Five magnetotelluric (MT) profiles have been acquired across the Cascadia subduction system and transformed using 2-D and 3-D nonlinear inversion to yield electrical resistivity cross sections to depths of ∼200 km. Distinct changes in plate coupling, subduction fluid evolution, and modes of arc magmatism along the length of Cascadia are clearly expressed in the resistivity structure. Relatively high resistivities under the coasts of northern and southern Cascadia correlate with elevated degrees of inferred plate locking, and suggest fluid- and sediment-deficient conditions. In contrast, the north-central Oregon coastal structure is quite conductive from the plate interface to shallow depths offshore, correlating with poor plate locking and the possible presence of subducted sediments. Low-resistivity fluidized zones develop at slab depths of 35–40 km starting ∼100 km west of the arc on all profiles, and are interpreted to represent prograde metamorphic fluid release from the subducting slab. The fluids rise to forearc Moho levels, and sometimes shallower, as the arc is approached. The zones begin close to clusters of low-frequency earthquakes, suggesting fluid controls on the transition to steady sliding. Under the northern and southern Cascadia arc segments, low upper mantle resistivities are consistent with flux melting above the slab plus possible deep convective backarc upwelling toward the arc. In central Cascadia, extensional deformation is interpreted to segregate upper mantle melts leading to underplating and low resistivities at Moho to lower crustal levels below the arc and nearby backarc. The low- to high-temperature mantle wedge transition lies slightly trenchward of the arc.

Structure of the Altyn Tagh Fault and Daxue Shan from magnetotelluric surveys: Implications for faulting associated with the rise of the Tibetan Plateau

Magnetotelluric measurements in the foreland of the Daxue Shan have imaged low-angle (<20°) thrust faults that extend to a depth of 3 km. The surface projections of the faults at depth coincide with the mapped traces of the Chang Ma (M = 7.6, 1932) and North Hills thrusts. Minimum cumulative offsets of 7 km on the Chang Ma thrust and 14 km on the North Hills thrust are estimated from the horizontal extent of underthrust sediment. Assuming that regional thrusting began at 5–6 Ma, this corresponds to minimum convergence rates of 1.3 mm/yr and 2.5 mm/yr on the Chang Ma and North Hills thrusts, respectively. These slip rates correspond to ∼ 3 mm/yr of sinistral slip on the Altyn Tagh Fault between 96°E and 97°E. This is consistent with the 4 ± 2 mm/yr of slip estimated by geological studies. Assuming comparable slip rates and similar fault geometry in the neighboring ranges, this requires a minimum of 11 mm/yr of shortening parallel to the Altyn Tagh fault between the Dang He Nan Shan and Qilian Shan. Both the style of thrusting and rate of shortening are in agreement with geologic studies that favor a relatively high rate of slip on the Altyn Tagh fault. This, in turn, implies that lithospheric extrusion contributes significantly in accommodating the ongoing convergence between India and Asia. Farther west, the Altyn Tagh Fault is imaged on four magnetotelluric profiles as a vertical resistivity contrast extending to a depth of at least 8 km. Two strands of the North Altyn Tagh Fault are imaged east of the asperity near Subei (39°30′N, 95°E).

On the geoelectric structure of major strike-slip faults and shear zones

Magnetotelluric imaging of the San Andreas Fault has shown that seismically-active segments are characterized by a zone of low resistivity in the upper crust. Similar resistivity features are observed on other major strike-slip faults, and may have a common origin in a region of fractured rock, partially or fully saturated with groundwater. Other strike-slip faults show possible zones of reduced resistivity in the mid and lower crust that may be related to zones of ductile shear. Additional MT surveys are required to elucidate the role of fluids in controlling the seismic behaviour of major faults, both in and below the seismogenic zone. A set of synthetic inversions show that MT data is sensitive to the geoelectric structure of a shear zone at mid-crustal depths.

Joint Inversion of Seismic and Magnetotelluric Data in the Parkfield Region of California Using the Normalized Cross-Gradient Constraint

We present jointly inverted models of P-wave velocity (Vp) and electrical resistivity for a two-dimensional profile centered on the San Andreas Fault Observatory at Depth (SAFOD). Significant structural similarity between main features of the separately inverted Vp and resistivity models is exploited by carrying out a joint inversion of the two datasets using the normalized cross-gradient constraint. This constraint favors structurally similar Vp and resistivity images that adequately fit the seismic and magnetotelluric (MT) datasets. The new inversion code, tomoDDMT, merges the seismic inversion code tomoDD and the forward modeling and sensitivity kernel subroutines of the MT inversion code OCCAM2DMT. TomoDDMT is tested on a synthetic dataset and demonstrates the code’s ability to more accurately resolve features of the input synthetic structure relative to the separately inverted resistivity and velocity models. Using tomoDDMT, we are able to resolve a number of key issues raised during drilling at SAFOD. We are able to infer the distribution of several geologic units including the Salinian granitoids, the Great Valley sequence, and the Franciscan Formation. The distribution and transport of fluids at both shallow and great depths is also examined. Low values of velocity/resistivity attributed to a feature known as the Eastern Conductor (EC) can be explained in two ways: the EC is a brine-filled, high porosity region, or this region is composed largely of clay-rich shales of the Franciscan. The Eastern Wall, which lies immediately adjacent to the EC, is unlikely to be a fluid pathway into the San Andreas Fault’s seismogenic zone due to its observed higher resistivity and velocity values.

Acquisition of a Unique Onshore/Offshore Geophysical and Geochemical Dataset in the Northern Malawi (Nyasa) Rift

The Study of Extension and maGmatism in Malawi aNd Tanzania (SEGMeNT) project acquired a comprehensive suite of geophysical and geochemical datasets across the northern Malawi (Nyasa) rift in the East Africa rift system. Onshore/offshore active and passive seismic data, long‐period and wideband magnetotelluric data, continuous Global Positioning System data, and geochemical samples were acquired between 2012 and 2016. This combination of data is intended to elucidate the sedimentary, crustal, and upper‐mantle architecture of the rift, patterns of active deformation, and the origin and age of rift‐related magmatism. A unique component of our program was the acquisition of seismic data in Lake Malawi, including seismic reflection, onshore/offshore wide‐angle seismic reflection/refraction, and broadband seismic data from lake‐bottom seismometers, a towed streamer, and a large towed air‐gun source.