Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.

Similar documents
Ground displacement in a fault zone in the presence of asperities

Introduction to Displacement Modeling

3D Finite Element Modeling of fault-slip triggering caused by porepressure

Derivation of Table 2 in Okada (1992)

Shear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with Rigid Surface

EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE MECHANISMS SHOW MOTION

BEYOND TRAVELTIMES AND EARTHQUAKE LOCATION What else can seismograms tell us about the nature of earthquakes on faults?

! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU

1. classic definition = study of deformed rocks in the upper crust

SUPPLEMENTARY INFORMATION

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL

Interactive 3D Sketchupbook

The Frictional Regime

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications

Earthquake and Volcano Clustering at Mono Basin (California)

Earthquakes and Seismotectonics Chapter 5

Resolving Stress Components and Earthquake Triggering

Basics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole

Mechanics of Earthquakes and Faulting

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment

Stress, Strain, Mohr s Circle

Earthquake and Volcano Deformation

both an analytical approach and the pole method, determine: (a) the direction of the

DATA REPOSITORY SUPPLEMENTARY MATERIAL. We analyzed focal mechanism solutions for aftershocks of the Loma Prieta

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES

Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement

Supporting Information for Inferring field-scale properties of a fractured aquifer from ground surface deformation during a well test

Simulation of earthquake rupture process and strong ground motion

ANALYSIS OF STRAINS CONCEPT OF STRAIN

Synthetic Seismicity Models of Multiple Interacting Faults

Effect of buttress on reduction of rock slope sliding along geological boundary

Elastic models of deformation in nature: why shouldn t we use the present day fault geometry?

Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress

Case Study 1: 2014 Chiang Rai Sequence

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

DEFORMATION PATTERN IN ELASTIC CRUST

Stress and Strain. Stress is a force per unit area. Strain is a change in size or shape in response to stress

Mathematical Modelling of a Fault Slip Induced by Water Injection

Data Repository Hampel et al., page 1/5

Learning goals - January 16, Describe the geometry of a fault (1) orientation of the plane (strike and dip) (2) slip vector

Seismic Source Mechanism

A circular tunnel in a Mohr-Coulomb medium with an overlying fault

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Provided by Tasa Graphic Arts, Inc. for An Introduction to Structural Methods DVD-ROM

EXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION

Exercise: concepts from chapter 8

Tectonics. Lecture 12 Earthquake Faulting GNH7/GG09/GEOL4002 EARTHQUAKE SEISMOLOGY AND EARTHQUAKE HAZARD

Prof. B V S Viswanadham, Department of Civil Engineering, IIT Bombay

Basic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008

A METHOD FOR DETERMINING ASPERITY PARAMETERS PRODUCING SPECIFIC MAXIMUM GROUND MOTION

Application of Fault Response Modelling Fault Response Modelling theory

GG303 Lecture 17 10/25/09 1 MOHR CIRCLE FOR TRACTIONS

Finite element modelling of fault stress triggering due to hydraulic fracturing

Onto what planes should Coulomb stress perturbations be resolved?

Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller

Exercise: infinitesimal strain during the great San Francisco earthquake of 1906

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

[5] Stress and Strain

to: Interseismic strain accumulation and the earthquake potential on the southern San

Crustal Deformation Earth - Chapter Pearson Education, Inc.

Using An Introduction to Structural Methods - An Interactive CD-ROM - In and Out of the Classroom

Analytic and Numeric Tests of Fourier Deformation Model (Copyright 2003, Bridget R. Smith and David T. Sandwell)

UNIVERSITY OF PRETORIA Department of Geology STRUCTURAL GEOLOGY -GLY 254 SEMESTER EXAM

Elastic rebound theory

Advanced Structural Analysis EGF Section Properties and Bending

San Andreas Movie Can It Happen?

Two Neighbouring Strike Slip Faults and Their Interaction

SUPPLEMENTAL INFORMATION

STRAIN. Normal Strain: The elongation or contractions of a line segment per unit length is referred to as normal strain denoted by Greek symbol.

Mechanics of Earthquakes and Faulting

Culminating Review for Vectors

Jack Loveless Department of Geosciences Smith College

This paper was prepared for presentation at the Unconventional Resources Technology Conference held in San Antonio, Texas, USA, 1-3 August 2016.

FRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao

Mechanics of Earthquakes and Faulting

Using information about wave amplitudes to learn about the earthquake size.

Force and Stress. Processes in Structural Geology & Tectonics. Ben van der Pluijm. WW Norton+Authors, unless noted otherwise 1/9/ :35 PM

Globale Seismizität I

INFLUENCE OF LOCAL PERTURBATION ON REGIONAL STRESS AND ITS IMPACT ON THE DESIGN OF MAJOR UNDERGROUND STRUCTURE IN HYDROELECTRIC PROJECT

STUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Improvement in the Fault Boundary Conditions for a Staggered Grid Finite-difference Method

Seismic Monitoring of Mine Environments

A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements

Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient

Lecture 8. Stress Strain in Multi-dimension

GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise

Introduction Faults blind attitude strike dip

AVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE

Three Dimensional Simulations of Tsunami Generation and Propagation

RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE

Summary so far. Geological structures Earthquakes and their mechanisms Continuous versus block-like behavior Link with dynamics?

DEFORMATION PATTERN IN ELASTIC CRUST

Mechanics of Earthquakes and Faulting

IZMIT BAY BRIDGE SOUTH APPROACH VIADUCT: SEISMIC DESIGN NEXT TO THE NORTH ANATOLIAN FAULT

Asish Karmakar 1, Sanjay Sen 2 1 (Corresponding author, Assistant Teacher, Udairampur Pallisree Sikshayatan (H.S.), Udairampur, P.O.

Earthquakes and Faulting

Orthogonal jointing during coeval igneous degassing and normal faulting, Yucca Mountain, Nevada

MEASUREMENT OF HYDRAULICALLY ACTIVATED SUBSURFACE FRACTURE SYSTEM IN GEOTHERMAL RESERVOIR BY USING ACOUSTIC EMISSION MULTIPLET-CLUSTERING ANALYSIS

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets

Transcription:

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much lower than an elastic model predicts, but strains and the principal axes of strain are more or less correct (1). The form of the strain tensor can consequently be used together with a simple Coulomb-like failure criterion (2) to define the orientations and mechanism of predicted faulting. In this work, the shear zone or fault is modeled as a rectangular dislocation in an elastic half-space (3). A Young's Modulus of 1.0 and a Poisson's ratio of 0.5 are adopted. With these values the results can then be interpreted in terms of strain, with the Coulomb strain being numerically equivalent to the Coulomb stress. The calculations are carried out using an effective friction of 0.4. Because the model uses a non-zero friction, the failure planes associated with Coulomb strain are not orthogonal and a standard seismic focal mechanism representation of the predicted fault planes is not correct. The convention adopted to denote the faults uses a simplified stereonet representation of the fault geometry and kinematics (Figs. 1d. and 2e). For each fault mechanism a blue or red line shows the intersection of the fault with the plane of the model. The rake for each fault is indicated by a red or blue dot showing the intersection of the slip vector with a hemisphere behind the modeled plane. To aid in discriminating relative motions, blue lines represent apparent clockwise (or right-lateral) motion, while red lines indicate apparent anticlockwise (or left-lateral) motion. For all of the models, Coulomb strain is plotted as a function of master fault slip. This allows the relative activity of different parts of the slip-partitioned system to be estimated. The figures use simplified representations of the master faults consistent with the known form of the fault system. Three segments (SA1 to 3) are used to model the San Andreas fault (Fig. 2b) and four (HA1 to 4) to model the Haiyuan fault (Fig. 2d). The locations of the faults within the figures are shown by green lines. Strains are calculated for 1 meter of horizontal slip in both cases. The direction of this slip is taken to be parallel to the mean orientation of the central San Andreas (Fig. 2b, segment SA1) and to the mean orientation of the Haiyuan west of the modeled area (Fig. 2d, segment

HA1). The slip is resolved onto the other segments to give appropriate strike-slip and dip-slip components. For both models the master faults extend from 30 km to 1000 km (infinite) depth. The coordinates of the segment ends for the San Andreas model are:- Segment SA1 dip 90 : 54.8300 N, -158.345 E to 34.9734 N -119.487 E; Segment SA2 dip 80 N: 34.9734 N, -119.487 E to 33.8646 N, -116.272 E; Segment SA3 dip 90 : 33.8646 N, -116.272 E to -13.5670 N, -82.2496 E. Segments 1 and 3 are pure strikeslip. The coordinates of the fault ends for the Haiyuan model are:- Segment HA1 dip 90 : 37.7293 N, 80.3217 E to 38.2828 N, 99.5088 E; Segment HA2 dip 70 S: 38.2828 N, 99.5088 E to 37.3872 N, 102.152 E; Segment HA3 dip 80 S: 37.3872 N, 102.152 E to 37.4510 N, 105.725 E; Segment HA4 (Gulang fault) dip 80 S: 37.3872 N, 102.152 E to 36.4163 N, 125.298 E. For segments HA3 and HA4 the horizontal component of slip is 70% on HA3, 30% on HA4. The conversion from geographic to Cartesian co-ordinates is carried out using a stereographic projection. For the San Andreas model the pole is at 34 N, -108 E. For the Haiyuan model it is at 38.0 N, 101.0 E.

DipofMainFault=30 Cross-Section Map-View Reverse Component = 1.0 m Strike-Slip Component = 0.0 m Reverse Component = 1.0 m Strike-Slip Component = 0.2 m Reverse Component = 1.0 m Strike-Slip Component = 0.4 m Reverse Component = 1.0 m Strike-Slip Component = 0.6 m Reverse Component = 1.0 m Strike-Slip Component = 0.8 m Reverse Component = 1.0 m Strike-Slip Component = 1.0 m

Reverse Component = 0.8 m Strike-Slip Component = 1.0 m Reverse Component = 0.6 m Strike-Slip Component = 1.0 m Reverse Component = 0.4 m Strike-Slip Component = 1.0 m Reverse Component = 0.2 m Strike-Slip Component = 1.0 m Reverse Component = 0.0 m Strike-Slip Component = 1.0 m -1x10-5 strain 1x10-5 Figure S1. Slip partitioning resulting from strain on a 30 dipping oblique fault. Left side of figure shows strain distribution and optimal failure planes in a cross-section perpendicular to the fault. Right side shows strain distribution and optimal failure planes in map view. Blue mechanisms indicate clockwise relative displacement in cross-section and right-lateral motion in map view. Red symbols indicate counter-clockwise motion in cross-section and left-lateral motion in map view. See Fig. 1d for explanation of symbols in crosssection and Fig. 2e for explanation of symbols in map view.

Dip of Main Fault = 45 Cross-section Map View Reverse Component = 1.0 m Strike-Slip Component = 0.0 m Reverse Component = 1.0 m Strike-Slip Component = 0.2 m Reverse Component = 1.0 m Strike-Slip Component = 0.4 m Reverse Component = 1.0 m Strike-Slip Component = 0.6 m Reverse Component = 1.0 m Strike-Slip Component = 0.8 m Reverse Component = 1.0 m Strike-Slip Component = 1.0 m

Reverse Component = 0.8 m Strike-Slip Component = 1.0 m Reverse Component = 0.6 m Strike-Slip Component = 1.0 m Reverse Component = 0.4 m Strike-Slip Component = 1.0 m Reverse Component = 0.2 m Strike-Slip Component = 1.0 m Reverse Component = 0.0 m Strike-Slip Component = 1.0 m -1x10-5 strain 1x10-5 Figure S2. Slip partitioning resulting from strain on a 45 dipping oblique fault. Left side of figure shows strain distribution and optimal failure planes in a cross-section perpendicular to the fault. Right side shows strain distribution and optimal failure planes in map view. Blue mechanisms indicate clockwise relative displacement in cross-section and right-lateral motion in map view. Red symbols indicate counter-clockwise motion in cross-section and left-lateral motion in map view. See Fig. 1d for explanation of symbols in crosssection and Fig. 2e for explanation of symbols in map view.

Dip of Main Fault = 60 Cross-Section Map View Reverse Component = 1.0 m Strike-Slip Component = 0.0 m Reverse Component = 1.0 m Strike-Slip Component = 0.2 m Reverse Component = 1.0 m Strike-Slip Component = 0.4 m Reverse Component = 1.0 m Strike-Slip Component = 0.6 m Reverse Component = 1.0 m Strike-Slip Component = 0.8 m Reverse Component = 1.0 m Strike-Slip Component = 1.0 m

Reverse Component = 0.8 m Strike-Slip Component = 1.0 m Reverse Component = 0.6 m Strike-Slip Component = 1.0 m Reverse Component = 0.4 m Strike-Slip Component = 1.0 m Reverse Component = 0.2 m Strike-Slip Component = 1.0 m Reverse Component = 0.0 m Strike-Slip Component = 1.0 m -1x10-5 strain 1x10-5 Figure S3. Slip partitioning resulting from strain on a 60 dipping oblique fault. Left side of figure shows strain distribution and optimal failure planes in a cross-section perpendicular to the fault. Right side shows strain distribution and optimal failure planes in map view. Blue mechanisms indicate clockwise relative displacement in cross-section and right-lateral motion in map view. Red symbols indicate counter-clockwise motion in cross-section and left-lateral motion in map view. See Fig. 1d for explanation of symbols in crosssection and Fig. 2e for explanation of symbols in map view.

Dip of Main Fault = 75 Cross-Section Map View Reverse Component = 1.0 m Strike-Slip Component = 0.0m Reverse Component = 1.0 m Strike-Slip Component = 0.2 m Reverse Component = 1.0 m Strike-Slip Component = 0.4 m Reverse Component = 1.0 m Strike-Slip Component = 0.6 m Reverse Component = 1.0 m Strike-Slip Component = 0.8 m Reverse Component = 1.0 m Strike-Slip Component = 1.0 m

Reverse Component = 0.8 m Strike-Slip Component = 1.0 m Reverse Component = 0.6 m Strike-Slip Component = 1.0 m Reverse Component = 0.4 m Strike-Slip Component = 1.0 m Reverse Component = 0.2 m Strike-Slip Component = 1.0 m Reverse Component = 0.0 m Strike-Slip Component = 1.0 m -1x10-5 1x10-5 strain Figure S4. Slip partitioning resulting from strain on a 75 dipping oblique fault. Left side of figure shows strain distribution and optimal failure planes in a cross-section perpendicular to the fault. Right side shows strain distribution and optimal failure planes in map view. Blue mechanisms indicate clockwise relative displacement in cross-section and right-lateral motion in map view. Red symbols indicate counter-clockwise motion in cross-section and left-lateral motion in map view. See Fig. 1d for explanation of symbols in crosssection and Fig. 2e for explanation of symbols in map view.

References 1. A. Hubert-Ferrari et al., Geophys. J. Int., 153, 111 (2003). 2. G.C.P.King, M. Cocco. Advances in Geophysics, 44, 1 (2001). 3. Y. Okada, Y. Bull. seism. Soc. Amer. 75,1135 (1985).