scholarly journals How well do surface slip measurements track slip at depth in large strike-slip earthquakes? The importance of fault structural maturity in controlling on-fault slip versus off-fault surface deformation

2014 ◽  
Vol 388 ◽  
pp. 38-47 ◽  
Author(s):  
James F. Dolan ◽  
Ben D. Haravitch

2021 ◽  
Author(s):  
Sarah Visage ◽  
Pauline Souloumiac ◽  
Nadaya Cubas ◽  
Bertrand Maillot ◽  
Arthur Delorme ◽  
...  

<p>During large strike-slip earthquakes, the displacement at the ground surface, only partially measured, is often under-estimated in comparison with the amount of slip inferred at depth. The resulting concept of shallow slip deficit is challenged by the precise measurements of surface deformation of on- and off-fault deformation by space imaging techniques, showing that a significant amount of deformation might be accommodated through distributed damage in a zone several hundred meters to kilometers wide around the fault. In this study, analogue modeling is used to quantify the distribution of on/off-fault surface deformation along strike-slip faults over the long term and to understand how it relates to the deep structure of the fault.</p><p>To do so, we used a 1.5 m x 1.34 m PVC box, and studied the deformation of a homogeneous sand pack deposited above a straight basal fault, with sand thicknesses varying from 2 to 8 cm. During strike-slip fault experiments, the first structures to appear are the Riedel shears (R-shears) followed by the synthetic shears (S-shears). These structures eventually coalesce to form an anastomosed fault zone, made of a succession of segments separated by geometrical complexities of variable size. Optical imagery is used, at every stage of the strike-slip fault formation, to (1) describe the 3D surface displacement and (2) precisely quantify on/off-fault deformation. </p><p>At the initiation of the fault before the formation of the Riedels, a zone of diffuse deformation is highlighted by a positive divergence of the displacement. This diffuse zone is also characterized by a vertical deformation that forms a bulge.</p><p>When the displacement Ux parallel to the basal fault has a gradient dUx/dy >= 0.1, we consider that it is "on-fault" deformation, and it is "off-fault, when that gradient is between 0.02 and 0.1.</p><p>At the Riedel shear stage, we find 40% of off-fault deformation over a unique Riedel fault and about 60% if deformation is distributed over two Riedels.</p><p>Once the strike-slip fault is formed, the ratio drops between 0 to 5 % of off-fault deformation over a fault segment, but the ratio increases to 20% along geometrical complexities.</p><p>Moreover, we also show that off-fault deformation around the early Riedel structures partly control the long-lived segmentation and morphology of the strike-slip fault.</p><p>Experimental results are then compared to observations and measurements of near-field and far-field deformation obtained along the 2013 Mw 7.7 Balochistan earthquake by Vallage et al. (2015) and Gold et al. (2015). Azimuthal displacements measured in a relay zone  (Vallage et al. 2015) are consistent with those observed along our experimental relay zones. Although our experiments were only run with sand, we found a similar distribution of the deformation at the surface. These observations suggest that the distribution of the surface deformation of strike-slip fault earthquakes is inherent to the fault structure, possibly inherited from the Riedel shear stage, and not induced by earthquakes dynamics.</p>



2020 ◽  
Author(s):  
Rishabh Dutta ◽  
Sigurjón Jónsson ◽  
Hannes Vasyura-Bathke

<p>Earthquake fault ruptures are typically complex and can consist of en echelon segments, have bends, large step-overs, and be curved or warped at different spatial scales. Although surface fault ruptures can be mapped using a variety of geological and geophysical techniques, the subsurface topology of faults is challenging to estimate. One of the main options is to use geodetic data (InSAR and GPS) of coseismic surface deformation to estimate the subsurface earthquake fault geometry along with the distributed slip. The general practice is to assume a planar fault surface and estimate the strike and dip of a simple rectangular fault prior to the spatially-variable slip estimation. Using such simplistic fault geometry during source fault estimations of large earthquakes rarely captures all the crustal deformation details seen in the data and can cause biased estimation results of the fault slip. Here, we show how complex non-planar fault geometry can be estimated simultaneously with spatially-variable slip from geodetic data in a Bayesian framework, where our non-planar fault geometry parametrization approach allows for various undulations of the fault surface in both the along-strike and down-dip directions.</p><p>We exemplify this approach through synthetic tests considering a checkerboard-like slip pattern on a listric non-planar fault. The results show that fault slip can be over-estimated by about 50-100% when using pre-assumed planar fault geometry. In contrast, both the non-planar fault geometry and spatially-variable slip are better retrieved when using our estimation approach. We then apply this modeling approach to the 2011 M<sub>W</sub>9.1 megathrust Tohoku-Oki (Japan) earthquake. Here we use prior information like the location of the trench and earthquake hypocenters during the Bayesian estimation to reduce the extent of the model space. The resulting fault geometry shows variations in fault dip in both the along-strike and down-dip directions that compare well with Hayes’ slab1.0 model of the subduction interface. The estimated fault slip is also comparable to the results that pre-defined the fault geometry using the slab1.0 model. In the future, the proposed method could lead to more realistic source models of major earthquakes, aided by improving computational resources and spatial resolution of geodetic data.</p>



2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yang Gao ◽  
HuRong Duan ◽  
YongZhi Zhang ◽  
JiaYing Chen ◽  
HeTing Jian ◽  
...  

AbstractThe 2019 Ridgecrest, California seismic sequence, including an Mw6.4 foreshock and Mw7.1 mainshock, represent the largest regional seismic events within the past 20 years. To obtain accurate coseismic fault-slip distribution, we used precise positioning data of small earthquakes from January 2019 to October 2020 to determine the dip parameters of the eight fault geometry, and used the Interferometric Synthetic Aperture Radar (InSAR) data processed by Xu et al. (Seismol Res Lett 91(4):1979–1985, 2020) at UCSD to constrain inversion of the fault-slip distribution of both earthquakes. The results showed that all faults were sinistral strike-slips with minor dip-slip components, exception for dextral strike-slip fault F2. Fault-slip mainly occurred at depths of 0–12 km, with a maximum slip of 3.0 m. The F1 fault contained two slip peaks located at 2 km of fault S4 and 6 km of fault S5 depth, the latter being located directly above the Mw7.1hypocenter. Two slip peaks with maximum slip of 1.5 m located 8 and 20 km from the SW endpoint of the F2 fault were also identified, and the latter corresponds to the Mw6.4 earthquake. We also analyzed the influence of different inversion parameters on the fault slip distribution, and found that the slip momentum smoothing condition was more suitable for the inversion of the earthquakes slip distribution than the stress-drop smoothing condition.



2013 ◽  
Vol 184 (4-5) ◽  
pp. 299-306 ◽  
Author(s):  
Richard J. Lisle

AbstractThe assumption is widely made that slip on faults occurs in the direction of maximum resolved shear stress, an assumption known as the Wallace-Bott hypothesis. This assumption is used to theoretically predict slip directions from known in situ stresses, and also as the basis of palaeostress inversion from fault-slip data. This paper examines different situations in relation to the appropriateness of this assumption. Firstly, it is shown that the magnitude of the shear stress resolved within a plane is a function with a poorly defined maximum direction, so that shear stress values greater than 90% of the maximum occur within a wide angular range (± 26°) degrees. The situation of simultaneous movement on pairs of faults requires slip on each fault to be parallel to their mutual line of intersection. However, the resolved shear stresses arising from a homogeneous state of stress do not accord with such a slip arrangement except in the case of pairs of perpendicular faults. Where fault surfaces are non-planar, the directions of resolved shear stress in general give, according to the Wallace-Bott hypothesis, a set of slip directions of rigid fault blocks, which is generally kinematically incompatible. Finally, a simple model of a corrugated fault suggests that any anisotropy of the shear strength of the fault such as that arising from fault surface topography, can lead to a significant angular difference between the directions of maximum shear stress and the slip direction.These findings have relevance to the design of procedures used to estimate palaeostresses and the amount of data required for this type of analysis.



1994 ◽  
Vol 84 (3) ◽  
pp. 668-691 ◽  
Author(s):  
David J. Wald ◽  
Thomas H. Heaton

Abstract We have determined a source rupture model for the 1992 Landers earthquake (MW 7.2) compatible with multiple data sets, spanning a frequency range from zero to 0.5 Hz. Geodetic survey displacements, near-field and regional strong motions, broadband teleseismic waveforms, and surface offset measurements have been used explicitly to constrain both the spatial and temporal slip variations along the model fault surface. Our fault parameterization involves a variable-slip, multiple-segment, finite-fault model which treats the diverse data sets in a self-consistent manner, allowing them to be inverted both independently and in unison. The high-quality data available for the Landers earthquake provide an unprecedented opportunity for direct comparison of rupture models determined from independent data sets that sample both a wide frequency range and a diverse spatial station orientation with respect to the earthquake slip and radiation pattern. In all models, consistent features include the following: (1) similar overall dislocation patterns and amplitudes with seismic moments of 7 to 8 × 1026 dyne-cm (seismic potency of 2.3 to 2.7 km3); (2) very heterogeneous, unilateral strike slip distributed over a fault length of 65 km and over a width of at least 15 km, though slip is limited to shallower regions in some areas; (3) a total rupture duration of 24 sec and an average rupture velocity of 2.7 km/sec; and (4) substantial variations of slip with depth relative to measured surface offsets. The extended rupture length and duration of the Landers earthquake also allowed imaging of the propagating rupture front with better resolution than for those of prior shorter-duration, strike-slip events. Our imaging allows visualization of the rupture evolution, including local differences in slip durations and variations in rupture velocity. Rupture velocity decreases markedly at shallow depths, as well as near regions of slip transfer from one fault segment to the next, as rupture propagates northwestward along the multiply segmented fault length. The rupture front slows as it reaches the northern limit of the Johnson Valley/Landers faults where slip is transferred to the southern Homestead Valley fault; an abrupt acceleration is apparent following the transfer. This process is repeated, and is more pronounced, as slip is again passed from the northern Homestead Valley fault to the Emerson fault. Although the largest surface offsets were observed at the northern end of the rupture, our modeling indicates that substantial rupture was also relatively shallow (less than 10 km) in this region.



Geomorphology ◽  
2021 ◽  
pp. 107965
Author(s):  
Yifei Li ◽  
Mian Liu ◽  
Huai Zhang ◽  
Yaolin Shi




2020 ◽  
Vol 110 (1) ◽  
pp. 154-165 ◽  
Author(s):  
Yuexin Li ◽  
Roland Bürgmann ◽  
Bin Zhao

ABSTRACT The Mw 6.5 Jiuzhaigou earthquake occurred on 8 August 2017 36 km west-southwest of Yongle, Sichuan, China. We use both ascending and descending Interferometric Synthetic Aperture Radar (InSAR) data from Sentinel-1 and coseismic offsets of four Global Positioning System sites to obtain the coseismic surface deformation field and invert for the fault geometry and slip distribution. Most slip of the left-lateral strike-slip earthquake occurred in the 3–10 km depth interval with a maximum slip of about 1 m and a large shallow slip deficit (SSD). An eight-month InSAR time-series analysis documents a lack of resolvable postseismic deformation, and inversions for the distribution of postseismic slip demonstrate the lack of shallow afterslip. We argue that the observations of a pronounced SSD and no early afterslip of the Jiuzhaigou earthquake are indicative of an immature fault and that all incipient young strike-slip faults likely feature a SSD. We would expect a complex rupture geometry with distributed coseismic failure in the uppermost part of the brittle crust during the fault-zone development. As faults mature, they straighten out, develop a localized fault-zone core, and the SSD diminishes. By calculating the static Coulomb stress change and nine-year viscoelastic stress change caused by the Wenchuan earthquake, we also show that the 2008 Wenchuan earthquake did not significantly affect the time of occurrence of the 2017 Jiuzhaigou earthquake.



2019 ◽  
Vol 220 (1) ◽  
pp. 79-95
Author(s):  
Lucile Bruhat

SUMMARY Most geodetic inversions of surface deformation rates consider the depth distribution of interseismic fault slip-rate to be time invariant. However, some numerical simulations show downdip penetration of dynamic rupture into regions with velocity-strengthening friction, with subsequent updip propagation of the locked-to-creeping transition. Recently, Bruhat and Segall developed a new method to characterize interseismic slip rates, that allows slip to penetrate up dip into the locked region. This simple model considered deep interseismic slip as a crack loaded at its downdip end, and provided analytical expressions for stress drop within the crack, slip and slip rate along the fault. This study extends this approach to strike-slip fault environments, and includes coupling of creep to viscoelastic flow in the lower crust and upper mantle. I use this model to investigate interseismic deformation rates along the Carrizo Plain section of the San Andreas fault. This study reviews possible models, elastic and viscoelastic, for fitting horizontal surface rates. Using this updated approach, I develop a physics-based solution for deep interseismic creep which accounts for possible slow vertical propagation, and investigate how it improves the fit of the horizontal deformation rates in the Carrizo Plain region. I found solutions for fitting the surface deformation rates that allow for reasonable estimates for earthquake rupture depth and coseismic displacement and improves the overall fit to the data. Best-fitting solutions present half-space relaxation time around 70 yr, and very low propagation speeds, less than a metre per year, suggesting a lack of creep propagation.



2013 ◽  
Vol 53 (1) ◽  
pp. 255 ◽  
Author(s):  
Ernest Swierczek ◽  
Cui Zhen-dong ◽  
Simon Holford ◽  
Guillaume Backe ◽  
Rosalind King ◽  
...  

The Rosedale Fault System (RFS) bounds the northern margin of the Gippsland Basin on the Southern Australian Margin. It comprises an anastomosing system of large, Cretaceous-age normal faults that have been variably reactivated during mid Eocene-Recent inversion. A number of large oil and gas fields are located in anticlinal traps associated with the RFS, and in the future these fields may be considered as potential storage sites for captured CO2. Given the evidence for geologically recent fault reactivation along the RFS, it is thus necessary to evaluate the potential impacts of CO2 injection on fault stability. The analysis and interpretation of 3D seismic data allowed the authors to create a detailed structural model of the western section of the RFS. Petroleum geomechanical data indicates that the in-situ stress in this region is characterised by hybrid strike-slip to reverse faulting conditions where SHmax (40.5 MPa/km) > SV (21 MPa/km) ~ Shmin (20 MPa/km). The authors performed geomechanical modelling to assess the likelihood of fault reactivation assuming that both strike-slip and reverse-stress faulting regimes exist in the study area. The authors’ results indicate that the northwest to southeast and east-northeast to west-southwest trending segments of the RFS are presently at moderate and high risks of reactivation. The authors’ results highlight the importance of fault surface geometry in influencing fault reactivation potential, and show that detailed structural models of potential storage sites must be developed to aid risk assessments before injection of CO2.



Sign in / Sign up

Export Citation Format

Share Document