Characterizing the mode of growth in crustal normal faults

Earthquake ruptures along tectonically active faults nucleate predominantly at depths of 5 to 12km in the crust, so the portions of faults that slip in these events cannot be directly observed. The geometry and composition of seismogenic faults controls the nucleation, propagation and termination of the earthquake rupture process. This study aims to place constraints on the geometry and composition of seismogenic faults by examining ancient faults exhumed from the depths at which earthquakes are observed to nucleate. Faults exposed in the Sierra Nevada, California, show that the internal architecture of earthquake faults is heterogeneous at a variety of scales. Field and microstructural observations are used to describe in detail the architecture of two pseudotachylyte-bearing fault systems in the Granite Pass region of Sequoia and Kings Canyon National Park; the Granite Pass fault (GPF) and associated faults, and the Glacier Lakes fault (GLF) and faults that splay from the GLF. The GPF and sub-parallel faults are 1 to 6.7km long with left-lateral strike-slip displacements up to 80m. The GPF and GPF-parallel faults have architectures that are heterogeneous along strike. They are composed of one to four fault core strands containing cataclasites and ultracataclasites that cross-cut early localized crystal-plastic deformation. Slip surfaces developed at the edges of, within and between fault cores are defined by pseudotachylytes and cataclasites with thicknesses of ~0.01 to 20mm. Fault-related subsidiary structures are developed on either side of fault cores, and comprise damage zones with widths orthogonal to the fault of up to 30m. The GLF and splay faults have architectures that are more homogeneous along strike. These faults are composed of a tabular volume of heavily fractured and altered host rock between approximately planar fault core strands. The fault cores are centimetres wide and contain cataclasites and foliated cataclasites that are cross-cut by pseudotachylytes. Fault-related damage is limited in extent to several metres beyond the bounding fault cores. The GLF contains additional cataclasites, ultracataclasites and pseudotachylytes in a fault core strand within the tabular zone of fractured rock. Thermochronologic analyses of the host rock granodiorite, combined with previously published palaeogeobarometry and apatite fission track data, define the temperature and pressure changes associated with cooling and exhumation of the pluton. The P-T conditions prevalent during the deformation history of the GPF fault system are evaluated by relating recrystallisation mechanisms in quartz to temperature, showing that the early stages of deformation occurred at temperatures of 450 to 600oC. Dating of pseudotachylytes by the K-Ar isotopic method suggests subsequent brittle deformation took place at temperatures <350oC and pressures ≤150MPa. A model for the architecture of the GPF architecture therefore has well constrained environmental controls, and should be transferrable to faults with comparable deformation histories. Small faults (cumulative displacements <1m) in the Mount Abbot Quadrangle, 55km north of Granite Pass, have been examined to illustrate the processes associated with the earliest stages of slip in the Sierra Nevada faults. The faults have branched or straight fault traces. Pseudotachylytes in branching faults show that these faults accumulated displacement in high velocity slip events, rather than by quasi-static fault growth. Branching faults without pseudotachylytes contain chlorite breccias interpreted to have formed in response to slip along faults with elevated pore fluid pressure. Straight faults also likely underwent slip events, but contain cataclased chlorite and epidote, suggesting low fluid pressures during slip. The small faults show that fluid-rock interactions are critical to fault geometry, and that lateral structural heterogeneity is established after small finite displacements. Field and thin section observations of exhumed seismogenic faults show that fault architecture and fault rock assemblage are critical to the earthquake rupture process. The heterogeneous composition of slip surfaces in the GPF faults imply that melt lubrication cannot account for all of the dynamic slip weakening as there are no continuous pseudotachylyte generation surfaces through the fault zones. Multiple slip weakening mechanisms must have been active during single rupture events. Slip weakening mechanisms also change at a given point on the fault in response to continued deformation. Splay faults at the GLF termination suggest that structural complexity observed at the terminations of fault surface traces can also be expected at depth. The off-fault damage at the termination of the GLF will change the bulk elastic properties of the host rock and must be accounted for in models of rupture propagation beyond fault terminations, or across geometrical discontinuities. Additionally, aftershock distributions and focal mechanisms may be controlled by the geometry of structures present at fault terminations.

Echelon formation

Earthquake rupture

Source

Cite

Citations (0)

Crustal collapse and the formation of synthetic normal faults

EAEJA (2003)

Roger Buck T. Nagel L. L. Lavier

Normal fault

Source

Cite

Citations (0)

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

Arben Pitarka

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

Strong ground motion

Earthquake rupture

10.2172/1599564

Cite

Citations (0)

Scaling Law of Outer Fault Parameters of Crustal Earthquakes Caused by Normal Faults

Journal of Japan Association for Earthquake Engineering (2016)

Dianshu Ju Kazuo DAN Atsuko OANA Haruhiko Torita Kazunari Motomura

Normal fault

Scaling law

10.5610/jaee.16.3_30

Cite