SI: Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California
|
|
- Sheena Thompson
- 5 years ago
- Views:
Transcription
1 SI: Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California T. H. W. Goebel 1*, S. M. Hosseini 2, F. Cappa 3,4, E. Hauksson 1, J. P. Ampuero 1, F. Aminzadeh 2, and J. B. Saleeby 1 1 Division of Geological and Planetary Sciences, CA Institute of Technology, Pasadena, CA, USA. 2 Department of Chemical Engineering and Material Sciences, University of Southern CA, Los Angeles, CA, USA. 3 Geoazur Laboratory, University of Nice Sophia-Antipolis, Nice, France. 4 Institut Universitaire de France, Paris, France * now at University of California, Santa Cruz Contents 1 Tectonic setting at the southern end of the Central Valley 4 2 Geologic structures involved in rupture and fluid migration processes 6 3 FMD, b-value and confidence interval 10 4 Earthquake cluster relocation 11 5 Additional earthquake detection using the template matching method 13 6 Moment tensor inversion and moment magnitudes 20 1
2 7 Hydrogeological model of fluid injection induced fluid pressure perturbations Intro: Crustal permeability structure close to injection well WD Theoretical background Sensitivity analysis of permeability and fault zone width
3 Additional Supporting Information (Files uploaded separately) File Name data.zip WF_new.zip modelresults.zip WWF_2D_v2.gif WWF_3D_v1.gif Description Seismicity catalogs and injection rates in well WD01, WD04 & WD05 Waveform plots of newly detected events and templates Results and initial conditions of pore-pressure modeling 2D animation of WWF swarm epicenter locations 3D animation of WWF swarm hypocenter locations 3
4 1 2 1 Tectonic setting at the southern end of the Central Valley In order to elucidate the complex structure in the area of the injection well and to discuss the plausible fluid pathways that could induce seismicity, we present a simplified fault map of the Tejon area (Figure 1), and generalized cross sections that show fault depth relations and fault inter-relationships (Figure 2). Figure 2 shows the position of the injection well above the approximate surface traces of buried, Early-Miocene, normal faults. A comparison of Figure 2 with the seismicity map in the main text (Figure 1a) shows a band of background seismicity along one of the principal buried normal faults, which along with geomorphic observations along the Tehachapi foothills where the fault is observed at basement levels, indicates Holocene fault-remobilization 20 (J.B. Saleeby, unpub. mapping). Recognizing the dominate north-south directed thrust or reverse sense to the M w 4.6 mainshock, we show in Figure 2 a south-north cross section that crosses the mainshock focal point (synthesized from 10,11,15 ; and J.B. Saleeby, unpub. mapping). The focal point of the M w 4.6 mainshock is deep within the footwall domains of the Pleito and Wheeler Ridge north-directed thrust faults, but lies within the depth projected White Wolf fault zone. On the south-north cross section we also denote the focal point of the north-directed reverse dominated 1952, M w 7.3 Kern County earthquake, which also lies along the depth projected White Wolf fault, at 19 km depth. These relations lead us to hypothesize that the M w 4.6 mainshock is associated with the White Wolf fault. 4
5 119.2 W W W W W W W 35.2 N 35.1 N 35.0 N 34.9 N 34.8 N Figure 1: Generalized fault map for Tejon area showing approximate surface traces of an active blind thrust and several reverse faults, buried early Miocene normal faults, and areas where faults are exposed in adjacent exhumed basement (sources: 6,14,20,21, J.B. Saleeby, unpub. mapping. Final drafting by Z.A. Saleeby.) 23 The principal basement structure of the White Wolf fault formed in the Late 24 Cretaceous in continuity with the early-phases of the Breckenridge-Kern Canyon 25 system to the north 6. Tens of kilometers of spatially varying oblique slip on the in- 26 tegrated system led to severe basement damage zones, numerous subsidiary fault 27 breaks, and rapid differential basement exhumation of up to 15 km 31. The typi- 28 cal width of the ductile to brittle basement damage zone, where exposed, is shown 29 on the south-north cross section in Figure 2. From Early-Miocene to Early-Pliocene 30 the White Wolf segment of the system functioned primarily as a north-down nor- 31 mal fault 10,11,33, resulting in its profound north wall sediment accommodation 32 space (Figure 2, south-north cross section). With the late Pliocene-Quaternary 5
6 33 34 growth of the Tehachapi-San Emigdio fold thrust belt the White Wolf zone has inverted along much of its length to a north-directed reverse/thrust zone 8, Geologic structures involved in rupture and fluid mi- gration processes Several plausible tectonic structures exist that could be associated with the M w 4.6 mainshock such as: 1) Rupture processes along a Riedel shear within the principal basement damage zone; 2) slip near the intersection of a hypothetical blind thrust fault zone, suggested by Davis and Lagoe, (1988) and Gordon & Gerke (2009), and approximated in depth as the hypothetical blind thrust on the Figure 2 south-north cross section; or 3) slip along a zone of stress concentration at the basement surface footwall cut-off of the White Wolf fault (Figure 2 south-north cross section). Potential fluid migration paths from the injection wells into the White Wolf fault zone can be visualized by projecting the structure of the Figure 2 1 west-east cross section northwestward onto the fault zone. Most notable are remobilized Early-Miocene normal faults. Fault damage zone permeability and hanging wall brecciation provide potential direct fluid pathways. Our geologic model is supplemented by detailed reservoir models based on standard petroleum industry data such as seismic profiles and well log data, available from the CA Division of Oil, Gas and Geothermal Resources (DOGGR) 3,4. The 1 In Figure 2, we included an uncertainty ellipse for the injection position shown on the west-east cross section. The uncertainty is derived from: 1) Geological uncertainty in subsurface stratigraphic position resulting from having only a mud log, and no supporting electric log nor core data for the injection well; 2) Rapid south to north facies changes in the Neogene section from the basin margin to interior 10,14,20 ; 3) Northwards projection of the structure across relatively small offset NE SW trending normal faults that are transverse to the principal NW SE trending normal faults that the projection follows; and 4) The effect of the 1500 m long NNW trending horizontal well segment for which fluids are incrementally injected over the last 520 m. 6
7 S 1 Tehachapi - San Emigdio foothills Pleito flt. Wheeler Ridge W-E section San Joaquin Basin N 1 km MSL -1? MSL W 1 Wheeler Ridge flt. Quaternary alluvium Base of Pleistocene Tulare Fm. Base of up. Pliocene San Joaquin Fm. Base of up. Miocene Santa Margarita Fm. Base of up. Miocene Monterey Fm. Base of lw. Miocene Tunis-Tecuya volcanics (Tv)? Basement Pre-Cenozoic White Wolf fault damage zone? 1952 M 7.3 reverse event at 16 km depth Hypothetical Blind Thrust M 4.7 White Wolf fault S-N section E injection well projected ~6 km 1 km northward along normal fault strike normal growth fault sedimentation? approx. depth of Tv isochronous surface km mi km MSL MSL -1-2 uncertainty on injection position N - plunging turbiditic Holocene sand lenses in Monterey remobilization Lw. Miocene normal growth faults Tunis volcanics km mi -6-6 km Figure 2: Simplified south to north, and west to east geologic cross sections for area of injection well and hypothesized induced seismicity. Locations given in Figure 1. Selected stratigraphic horizons are shown as structural form lines. Sources given in Supplementary text. Final drafting by Z.A. Saleeby. 7
8 industry data highlights many fracture and fault zones located between injection site and the White Wolf fault (Figure 3). The WNW orientation of these faults is in agreement with a linear feature in the seismicity catalog highlighting that at least one of these faults is seismically active. This fault is referred to as Tejon fault in the main text. The seismicity catalog indicates that the Tejon fault deepens toward the White Wolf fault which is also seen within some of the interpreted well-logs that extent to 3.7 km depth 4. The closest fault zones identified in the DOGGR data is located less than 0.5 km from injection site WD05. Thus, these faults likely influenced the nearby permeability structure and fluid flow path ways of injected wastewater. 8
9 M4.7 o W e hit W au lf F lt Tejon, North Pleito Fault Wh ult lf Fa o ite W Fault and Fracture Zones 0 1 WD05 Figure 3: Mapped faults and fractures within the Tejon, North oilfield. The top plot shows the approximate boundaries of the oil field (blue region) and earthquake locations. Seismicity likely associated with the Tejon fault is highlighted in black and background seismicity in gray. The bottom plot shows the more detailed oilfield boundaries (gray shaded area) based on a DOGGR technical report 3. Faults and fracture zones were identified based on industry data such as seismic profiles, geological mapping and well logs (Modified after California Division of Oil, Gas & Geothermal Resources, Technical Report 3 ). The displayed seismicity was predominantly recorded before the 2005 swarm and highlights the approximately linear trend of the Tejon fault which connects the injection site and largest magnitude epicenter. 9
10 62 3 FMD, b-value and confidence interval Frequency-magnitude distributions (FMD) of earthquakes within the study area generally followed a power-law of the form: 65 log N = a b M (1) where b is the power-law exponent, N and M are the event number and magnitude and a is a constant related to the overall seismic productivity within the area 17. To compute b-values, we used the maximum likelihood approach 1, b = 1 M M c log(e) (2) where M c is the magnitude of completeness, corrected for bin-size 39 and M is the mean magnitude of events above M c. The lower power-law bound, M c, was determined by minimizing the Kolmogorov-Smirnoff distance between the observed and modeled power-law distributions 7. We compared the observed FMD with expectations for FMDs with a b-value of unity using Monte Carlo simulations 13. The corresponding confidence limits as function of magnitude are reported as the 2.5 and 97.5 percentiles of 1,000 randomly sampled FMDs with the same parameters (a, b, M c ) as the observed distribution. In the main manuscript, we show that b-values show significant variations that coincide with peak injection rates in well WD05. b-values are close to unity before and after peak injection but reach values as low as b=0.6 during the highest injection rates in In addition, we investigate long-term temporal variations in b-value within a 10 km radius of well WD05 between 1995 and 2011 (Figure 4). b-values show strong variations between 0.6 and 1.4. We observed a systematic 10
11 decrease in b starting in 2004 which is when one of the high-rate injection wells became active. This first order correlation between injection rates and temporal b-value variations may indicate a connection between pressure changes and the tendency of earthquakes to grow to larger sizes resulting in lower b-values. However, variations in b-value should be interpreted carefully because of the inherent uncertainties in estimating power-law exponents over a limited number of decades and from seismicity records with changing magnitude of completeness over time. Moreover, a precise correlation between b-value variations and other time-series, e.g. injection rates, is complicated by different temporal resolutions. Thus, the most robust features within our data are generally low b-values during the time of peak injection rates Earthquake cluster relocation In addition to the standard single event locations within the Southern California Seismic Network (SCSN) earthquake catalog, we use a state-of-the-art waveformrelocated catalog 19. The relocated catalog is based on a detailed, 3D velocity model and precise, travel-time differences determined by waveform cross-correlations of earthquakes within clusters 26,36. The relocated seismicity catalog highlights that most of the earthquakes that occurred during the time of peak injection in WD05 occurred at depths between the injection site and the M w 4.6 hypocenter (Figure 5). The seismic activity that deepened toward the White-Wolf fault was likely a key component in maintaining a high-permeability pathway that transmitted pressure changes to large depths. 11
12 WD05 Figure 4: Temporal b-value variations within a 10 km radius around well WD05 and injection rates in wells WD01, WD04 and WD05. 12
13 0 5 Depth [km] Number of Earthquakes Figure 5: Focal depth distributions of events within a 10 km radius of well WD05 between 04/2005 to 09/2005. The red star highlights the focal depth of the M w 4.6 hypocenter Additional earthquake detection using the template matching method We extend the available earthquake catalog by cross-correlating event template waveforms with the continuous waveform record. The waveform templates are created from the earthquake record within the SCSN catalog. We use event waveforms as templates that were recorded at the closest station with high signal-tonoise ratio, i.e. station MPI. For this study, we selected 14 template event waveforms with magnitudes between 0.6 and 1.8, which occurred within a ten kilometer radius of well WD05 and were recorded and produced high-quality template waveforms on all three components at station MPI. This station was installed in Our implementation of the template-matching method largely follows Meng et al., (2013) 30 and consists of 4 main steps: Template selection and processing: We down-sample all waveforms to 20 Hz to match the continuous, broadband waveform record which is available af- ter June 2005 and bandpass-filter all waveforms between 5 to 10 Hz. This fre- 13
14 quency band is generally above the largest, coherent noise sources, providing the most stable results for event detection. After bandpass-filtering, the template waveforms are cut from 1 s before to 3 s after P-arrival on the vertical component. Waveforms on the radial components are selected between 2 s before and after the S-phase arrivals accounting for larger uncertainties in phase picks Cross-correlation: The template waveforms are then cross-correlated with the continuous waveform record between June 30th 2005 until the mainshock occurrence on September 22nd Cross-correlation coefficients (CC) between templates and continuous waveforms are computed for a 4 sec sliding window that is advanced by one sample steps for all three traces Event detection: An event detection is declared if the individual and stacked cross-correlation coefficients exceed a certain threshold. For the present application, we chose a conservative CC threshold to minimize false detection. We tested different CC thresholds and find that a factor of 7 times the standard deviation provides the best results Earthquake-catalog creation: For the newly-detected events, we use the times of highest correlation coefficients as origin times and adopt the locations and magnitudes of the template events for the new event We test our implementation of the template matching method by first evaluating its ability to self-detect events within the earthquake catalog that were recorded between June 2001 and October Expectedly, the corresponding correlationcoefficients (CC) across all traces are close to unity for these events (Figure 6). We minimize false-detection by requiring relatively high CC values that are at least 7 times larger than the standard deviation of noise correlations. In addition, 14
15 Figure 6: Example of waveform cross-correlations for a test-run of the template matching algorithm to self-detect a catalog event. Top: Template (blue) and continuous seismograms (gray) for the vertical, and the two horizontal components. CC values are shown in upper left corner of each subplot. Bottom Left: Time series of different CC values (black), detection threshold (red line) and event detection (red dot). Bottom Right: Histogram of CC values (black), detection threshold (red line) and event detection (red dot) we test if S-P times and event durations are in approximate agreement between template and newly detected events for each of the 21 new detection. Lastly, we check if the first motion polarities agree between templates and new detection. Figure 8 to 10 show three examples of successfully detected events including waveform records and CC values. The waveform plots for the other detected events are part of the Online Supplement. Seven of the 14 template waveforms were connected to additional seismic event detection resulting in a total of 21 additional event detection between the start of injection in well WD05 and the occurrence the mainshock (Figure 7). 15
16 Newly Detected Events Template Loca,on MPI Figure 7: Map-view of template (blue dots and waveforms) and newly-detected (red dots and gray waveforms) events based on cross-correlating waveforms at station MPI (red triangle). The White-Wolf swarm is highlighted by red, orange, yellow and green dots analogous to Figure 1a in the main manuscript.the injection site is highlighted by a blue triangle and the largest magnitude swarm event by a red star. 16
17 Figure 8: Same as Figure 6 but now showing an example of waveforms and CC values of a newly detected event Template matching is generally superior to simple phase picking algorithms, which rely on a notable amplitude or frequency difference between signal and noise, allowing for event detections even in high-noise environments 29,30,37. However, the limited availability of template events due to only one nearby station with continuous waveforms records that was recording over several years before and after the White Wolf swarm in 2005, prevents a systematic lowering of the completeness threshold within the area. Nevertheless, the template matching and catalog extension allow us to test if part of the Tejon fault was active during the high-rate injection period in WD05 and provides information about possibly systematic migration patterns (see main manuscript). 17
18 Figure 9: Same as Figure 6 but now showing an example of waveforms and CC values of a newly detected event. 18
19 Figure 10: Same as Figure 6 but now showing an example of waveforms and CC values of a newly detected event. 19
20 165 6 Moment tensor inversion and moment magnitudes We test the reliability of magnitude estimates in the upper tail of the FMD through determining moment magnitudes of events with M L 4.0. We find that the local magnitudes of M L = 4.5, 4.7 and 4.3 within the SCSN earthquake catalog exceed moment magnitude estimates of M w = 4.1, 4.6 and 4.0. However, three events above M w 4 are in good agreement with the expected number of M 4 events for the determined power-law parameters (b=0.6, a=2.8) and further support a significant deviation from a b-value of unity. Moment magnitudes were computed by fitting long-period waveforms in a 10 to 20 sec window immediately after P and S-phase arrivals using the Cut-And- Paste method for pure, double-couple source tensors 40,41. The synthetic far-field displacement for a double-couple source can be described as follows: 177 s(t) = M 0 3 i=1 A i (φ θ, δ, λ)g i (t), (3) here, i = 1,2,3 correspond to three fundamental faults, i.e. vertical strike slip, vertical dip slip and 45 dip slip, G i s are the Green s functions, φ is the station azimuth, θ, δ, λ are the strike, dip and rake of the source, and M 0 is the seismic moment. The parameters describing the source tensor can be determined by minimizing the misfit between observed and synthetic seismograms using a grid-search method within the limited parameter space. The CAP method fits both body and surface wave phases allowing for small time shifts (here 2 to 8 sec for P and S phases) between different phases, which account for unmodeled velocity structures and resulting uncertainties in Green s function estimates. We use a maximum fitting bandwidth of 0.02 to 0.2 Hz with a slightly narrower band for the smaller magnitude events between 0.08 to 0.2 Hz. We perform another 20
21 grid-search to solve to the focal depth of the synthetic source that minimizes the waveform misfits. For the focal depth inversion, we require reduction in waveform misfits of at least 10% and use the depth within the relocated catalog otherwise. For the largest magnitude event with M w 4.6, we find a focal depth of 9 km which is in good agreement with the relocated catalog depth. The following Figures 11 to 13 show inversion results and waveform misfits for the three events with M w
22 M4.4 M w 4.6 Mainshock, September 22 nd, 2015 M4.5 M4.5 M4.5 M4.5 M4.5 M4.5 M4.5 M4.6 M4.6 M4.6 strike = 89 dip = 30 rake = 72 sec Distance Sta@on Azimuth Time Shi= Correla@on Coefficient Figure 11: Focal mechanism (upper left), results of focal depth grid search including moment magnitude (upper right) and waveforms misfits (bottom) for the largest magnitude event that was part of the 2005 White Wolf swarm. The best-fitting focal mechanism is highlighted in red, observed waveforms are shown in black and synthetic waveforms in red. 22
23 M w 4.1 Forehock, September 22 nd, 2015 M4.1 M4.1 M4.1 M4.1 M4.1 M4.1 M4.1 M4.2 M4.2 M4.0M4.0 M4.0 M4.1 M4.0 M4.1 strike = 110 dip = 32 rake = 90 sec Distance Sta?on Azimuth Time Shi; Correla?on Coefficient Figure 12: Same as Figure 11 for an earthquake above magnitude 4 that occurred 44 sec before the M w 4.6 event. 23
24 M w 4.0 A)ershock, September 22 nd, 2015 M3.9 M3.9 M3.9 M4.0 M4.0 M4.0 M4.0 M4.0 M3.9 M3.9 M3.9 M4.0 M4.0 M4.0 M4.0 strike = 89 dip = 40 rake = 72 sec Distance StaEon Azimuth Time Shi) CorrelaEon Coefficient Figure 13: Same as Figure 11 for an earthquake above magnitude 4 that occurred 6 min after the M w 4.6 event. 24
25 Hydrogeological model of fluid injection induced fluid pressure perturbations Intro: Crustal permeability structure close to injection well WD05 We assess injection induced fluid pressure changes as a function of distance and time from injection operations in well WD05 by modeling pressure diffusion in three dimensions. The 3-D numerical diffusion model is based on the most complete available datasets within the upper 2 3 km of sedimentary basins and includes seismicity records, geologic mapping results, and industry data (i.e. welllogs, stratigraphic columns and interpreted reservoir structure). Below 3 km depth little data is available except for seismicity catalogs. The model includes three principal stratigraphic zones, in addition to the high-permeability portion of the Tejon fault which is implemented as a vertical zone of elevated permeability. These three zones are: 1) The injection zone, i.e. a m thick turbiditic sand lens in the Monterey formation with a lateral extent of up to 1.5 km, labeled as "Transition zone" in the industry data 4 ; 2) the crystalline (gneissic) basement complex, and 3) the Monterey formation (see 10,11,15 ). Permeability is generally high within the sand lenses (i.e. 1 D) and very low ( 10 4 md) above injection depth which is one of the requirements for the selection of an injection site. Similarly, permeability is low ( 10 4 md) within the basement and Monterey formations outside of the injection zone (permeabilities are reported in milli-darcy, 1 md m 2 ). High permeability within the injection zone is supported by publicly available well-log data, i.e., low Gamma-ray and high resistivity values 4. 25
26 Theoretical background In our model, injection induced pore-pressure diffusion is described by the uncoupled diffusivity equation for slightly compressible fluids for which changes in pressures do not affect the elastic deformation of the rock matrix 28,35 : p k x x 2 + k 2 p y y 2 + k 2 p z z 2 = φµc p t, (4) where k x, k y and k z are permeability in x, y and z directions, p is pore-pressure, φ is porosity, µ is fluid viscosity and c is rock matrix compressibility. To present the diffusivity equation in the presence of sources and sinks, we derive the mathematical description of single-phase fluid flow in porous media. We start with the continuity equation, Darcy s law for flow in 3D and the equation of state. The continuity equation in Cartesian coordinates with added source terms can be written as 25,28 : 231 x (ρq x) + y (ρq y) + z (ρq z) +... N i=1 [ρq wδ(x x w )δ(y y w )δ(z z w )] (5) = t (φρ), where q x, q y and q z are volumetric rates of flow per unit cross-sectional area, ρ is fluid density, x w, y w and z w are the injection locations, Q w is the flow rate per unit volume, N is the total number of injection wells in the domain, δ(x) is the Dirac delta function and φ is porosity. Darcy s law for flow in the x, y and z directions 26
27 236 can be written as 28 : q x = k xρ µ p x q y = k yρ p µ y q z = k ( ) zρ p µ z + ρg (6) (7) (8) where k x, k y and k z are permeability in x, y and z directions, p is pore-pressure, µ is fluid viscosity and ρ is fluid density. Using Darcy s flow law and Eq. 5, we can determine diffusion in the presence of sources and gravitational effects: ( ) ( ) ( ) kx ρ p x µ x + y ky ρ p µ y + z kz ρ µ ( p z + ρg)... (9) i=1 N [ρq wδ(x x w )δ(y y w )δ(z z w )] = t (φρ), In order to study isothermal single-phase flow of fluids with small and constant permeability, we use the following expression for fluid compressibility, c: c = 1 ρ ρ p (10) After integration of Eq. 10, we obtain the equation of state for slightly com- pressible fluids: ρ = ρ 0 e c(p p0). (11) as: The combined compressibility of fluid (c) and rock matrix (c r ) can be expressed c c = c + c r, (12) 254 The rock or formation compressibility can be described as change in porosity as a 27
28 255 function of change in pore pressure: c r = 1 φ φ p, (13) where φ again is porosity. Using the equation of state (Eq. 11) and Eq. 9, 12 & 13, we get: 259 k x 2 p x 2 + k y 2 p y 2 + k z 2 p z 2 + 2k z cρg p z... N i=1 [ρq wδ(x x w )δ(y y w )δ(z z w )] = φµc t p t, (14) which is the diffusivity equation used in our model. This description of the diffusive process assumes that pressure gradients and compressibility are small, and porosity does not change with time. This diffusivity equation is solved numerically using the finite-difference reservoir simulator ECLIPSE 32, Numerical model set-up For the present study, we used the single-phase reservoir simulator module of ECLIPSE ( black-oil simulator ) 34. The initial conditions of the numerical model are constrained by the available industry and geologic data described in Section 7.1 and include the high-permeability part at the top of the Tejon fault, identified through a combination of geologic mapping, seismicity locations and industry data. We assume a pressure gradient of 12 MPa/km. We use an average value for the porosity, φ = 0.2, rock compressibility, c r = Pa 1 and fluid compressibility, c = Pa 1. 3,18 The water formation volume factor at a reference pressure of 30 MPa is 1.13 (reservoir m 3 / standard condition m 3 ). The water viscosity at reference pressure is 0.5 cp. The size of the injection layer of 2.1 km 2 is determined from the publicly available data-sets on the California Division of Oil, 28
29 Gas and Geothermal Resources (DOGGR) website 3. For the present study, we assign no-flux boundary conditions. We chose a large modeling domain (x=30 km, y=5 km, z=18 km) and placed the injection well at the center of this domain to minimize potential effects of boundary conditions on the pressure solutions. The model domain is discretized in the following grid sizes: dx = m, dy = 500 m and dz = 100 m. At the vicinity of injection point and high permeability pathway, we use local grid refinement down to dx = 5.78, dy = 9.55 and dz = 1.58 m. The injection well is modeled as a vertical injector with a perforation zone of 30 m that spans the entire extent of the injection zone and exact monthly injection schedule is obtained from DOGGR 4, plotted in Figure Sensitivity analysis of permeability and fault zone width Crust and fault zone permeability vary non-linearly over several orders of magnitude both laterally and with depth so that permeability is one of the modeling parameters most difficult to constrain. The permeability within fault and fracture zones is generally higher than in the fault core and protolith and can be even higher in the presence of rupture induced damage, for example, toward the up-dip end of blind faults 12,23,27,38. The interplay between fault motion and permeability increase has repeatedly been observed in controlled injection experiments and is essential in maintaining relatively high fault permeability at seismogenic depth 2,5,9,16. This is also in agreement with observations at the continental deep drilling project (KTB), Germany where zones of elevated permeability were associated with fracture zones preferably oriented to slip within the local stress field 38,42,43. These zones of elevated permeability persisted down to depth of 9.1 km 42. Similarly, we expect regions with systematically higher seismicity density toward the upper limit of the Tejon fault to exhibit higher permeability. 29
30 In our model, we vary both permeability (k) and width (w) of the upper portion of the Tejon fault over a wide range of values (i.e. k=10 6 to 1 D and w=100 to 1000 m) to determine a plausible range of pore-pressure changes at the M w 4.6 hypocenter (Table 1 & 2). The sensitivity analysis showed that differences in permeability have a stronger influence on pore-pressure variations; however both higher fault zone permeability and narrower fault zone width can lead to larger pore-pressure changes at the M w 4.6 hypocenter. For fault zone permeability above 200 md and fault width below 1000 m, we observe a pressure increase at the M w 4.6 hypocenter of more than 0.1 MPa which is sufficient to trigger earthquakes on faults close to failure 22,24. The modeled pore-pressure increase at the M w 4.6 hypocenter is lower for a 312 fault zone with homogeneous permeability, i.e. without the high-permeability pathway toward the top of the fault. For fault zones with a width between w=300 and 600 m and permeability k >400 md, we expect an increase in pore-pressure at M w 4.6 hypocenter of less than 0.01 MPa (Figure 14). This further highlights that a narrow, high permeability pressure channel is expected to have the largest seismogenic consequences for the White Wolf fault. 30
31 Table 1: Overview of model runs, initial conditions and pore pressure increase at the M w 4.6 hypocenter for a wide range of permeabilities. Model run k i k f w IZL P p [md] [md] [m] [km] [MPa] < < < < < k i - injection zone permeability; k f - fault zone permeability; w - fault zone width; IZL - injection zone lateral extent; P p - pore pressure change at the M w 4.6 hypocenter after 150 days of high-rate injection 31
32 Table 2: Overview of model runs, initial conditions and pore pressure increase at the M w 4.6 hypocenter for a narrow range of permeabilities. Model run k i k f w IZL P p [md] [md] [m] [km] [MPa] k i - injection zone permeability; k f - fault zone permeability; w - fault zone width; IZL - injection zone lateral extent; P p - pore pressure change at the M w 4.6 hypocenter after 150 days of injection 32
33 P [MPa] p Figure 14: Expected pore-pressure increase at distances similar to M w 4.6 hypocentral distance for pressure diffusion along a fault with homogeneous permeability structure. The resulting pressure perturbations are significantly smaller than in the case of a highpermeability pathway at the top of the Tejon fault. 33
34 318 References [1] Aki, K. (1965), Maximum likelihood estimate of b in the formula log N = a - bm and its confidence limits, Bull. Earthquake Res. Inst., Tokyo Univ., 43, [2] Bourouis, S., and P. Bernard (2007), Evidence for coupled seismic and aseismic fault slip during water injection in the geothermal site of Soultz (France), and implications for seismogenic transients, Geophys. J. Int., 169(2), [3] CA Department of Conservation (1998), California oil and gas fields - Central California, DOGGR technical reports, 1, [4] CA Department of Conservation (2012), Division of Oil, Gas and Geothermal resources: Production history in California, (ftp://ftp.consrv.ca.gov/pub/oil/annual_reports, accessed Sept. 26, 2014) [5] Chabora, E., et al. (2012), Hydraulic stimulation of well 27 15, Desert Peak geothermal field, Nevada, USA, in Proceedings of thirty seventh workshop on geothermal reservoir engineering, Stanford University, Stanford, vol [6] Chapman, A., J. Saleeby, D. Wood, A. Piasecki, S. Kidder, M. Ducea, and K. Farley (2012), Regional displacement analysis and palinspastic restoration of dispersed crustal fragments in the southern sierra nevada, california, Geosphere, 8(2) [7] Clauset, A., C. R. Shalizi, and M. E. J. Newmann (2009), Power-law distribu- tions in empirical data, SIAM review, 51(4), [8] Clinton, J. F., E. Hauksson, and K. Solanki (2006), An evaluation of the SCSN 34
35 moment tensor solutions: Robustness of the m w magnitude scale, style of faulting, and automation of the method, Bull. Seis. Soc. Am., 96(5), [9] Cornet, F., J. Helm, H. Poitrenaud, and A. Etchecopar (1997), Seismic and aseismic slips induced by large-scale fluid injections, Pure appl. geophys., 150, [10] Davis, T. (1983), Late Cenozoic Structure and Tectonic History of the Western Big Bend of the San Andreas Fault and Adjacent San Emigdio Mountains, unpublished Ph.D. Dissertation, University of California, Santa Barbara, p. 578 pp [11] Davis, T. L., and M. Lagoe (1988), A structural interpretation of major tectonic events affecting the western and southern margins of the San Joaquin Valley, California, in Studies of the Geology of the San Joaquin Basin, Pacific Section, Society of Economic Paleontologists and Mineralogists, Pacific Section, Los Angeles, California, 60, 60 65, ed. S.A. Graham [12] Faulkner, D. R., C. A. L. Jackson, R. J. Lunn, R. W. Schlische, Z. K. Shipton, C. A. J. Wibberley, and M. O. Withjack (2010), A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones, J. Struct. Geol., 32(11), , doi: /j.jsg [13] Felzer, K. R., T. W. Becker, R. E. Abercrombie, G. Ekström, and J. R. Rice (2002), Triggering of the 1999 Mw 7.1 Hector Mine earthquake by aftershocks of the 1992 Mw 7.3 Landers earthquake, J. Geophys. Res., 107(B9), ESE [14] Goodman, E. D., and P. E. Malin (1992), Evolution of the Southern San Joaquin Basin and mid-tertiary transitional tectonics, central California, Tectonics, 11(3),
36 [15] Gordon, S., and H. Gerke (2009), Controls on petroleum occurrence and exploration prospectiveness in the southern San Joaquin Basin, California, Contributions to the Geology of the San Joaquin Basin, California, Pacific Section, American Association of Petroleum Geologists, Miscellaneous Publication, 48, [16] Guglielmi, Y., F. Cappa, J.-P. Avouac, P. Henry, and D. Elsworth (2015), Seis- micity triggered by fluid injection induced aseismic slip, Science, 348(6240), [17] Gutenberg, B., and C. F. Richter (1944), Frequency of earthquakes in Califor- nia, Bull. Seismol. Soc. Am., 34, [18] Hall, H. N. (1953), Compressibility of reservoir rocks, Journal of Petroleum Tech- nology, 5(01), [19] Hauksson, E., W. Yang, and P. M. Shearer (2012), Waveform relocated earth- quake catalog for Southern California (1981 to june 2011), Bull. Seism. Soc. Am., 102(5), [20] Hirst, B. (1986), Tectonic development of the Tejon and adjacent areas, Kern County, California, Pacific Section of AAPG, Guidebook, pp [21] Hirst, B. (1986), Early Miocene tectonism and associated turbidite deposys- tems of the Tejon area, Kern County, California, Pacific Section, Society of Eco- nomic Paleontologists and Mineralogists, 60, [22] Hornbach, M. J., et al. (2015), Causal factors for seismicity near Azle, Texas, Nature communications, [23] Ingebritsen, S., and C. Manning (2010), Permeability of the continental crust: 36
37 dynamic variations inferred from seismicity and metamorphism, Geofluids, 10(1-2), [24] Keranen, K., M. Weingarten, G. Abers, B. Bekins, and S. Ge (2014), Sharp in- crease in central Oklahoma seismicity since 2008 induced by massive wastew- ater injection, Science, 345(6195), [25] Lee, J., J. B. Rollins, and J. P. Spivey (2003), Pressure transient testing, vol. 9, Henry L. Doherty Memorial Fund of Aime Society of Petroleum [26] Lin, G., P. M. Shearer, and E. Hauksson (2007), Applying a three-dimensional velocity model, waveform cross correlation, and cluster analysis to locate southern California seismicity from 1981 to 2005, J. Geophys. Res., 112(B12), doi: /2007jb [27] Manga, M., I. Beresnev, E. E. Brodsky, J. E. Elkhoury, D. Elsworth, S. Ingebritsen, D. C. Mays, and C.-Y. Wang (2012), Changes in permeability caused by transient stresses: Field observations, experiments, and mechanisms, Rev. Geophys., 50(2) [28] Matthews, C., and D. G. Russel (1967), Pressure buildup and flow tests in wells, Society of Petroleum Engineers Monograph Series, 1, 163 pp [29] Meng, X., X. Yu, Z. Peng, and B. Hong (2012), Detecting earthquakes around Salton Sea following the 2010 m w 7. 2 El Mayor-Cucapah earthquake using GPU parallel computing, Procedia Computer Science, 9, [30] Meng, X., Z. Peng, and J. L. Hardebeck (2013), Seismicity around Park- field correlates with static shear stress changes following the 2003 m w 6.5 San Simeon earthquake, J. Geophys. Res., 118(7),
38 [31] Nadin, E. S., and J. B. Saleeby (2008), Disruption of regional primary structure of the Sierra Nevada batholith by the Kern Canyon fault system, California, Special Papers - Geol.L Soc. Am., 438, [32] Pettersen, Ø. (2006), Basics of reservoir simulation with the eclipse reservoir simulator, Bergen, Norway: University of Bergen, Department of Mathematics, Lec- ture Notes [33] Saleeby, J., Z. Saleeby, and L. Le Pourhiet (2013), Epeirogenic transients related to mantle lithosphere removal in the southern Sierra Nevada region, California: Part II: Implications of rock uplift and basin subsidence relations, Geosphere, 9(3), [34] Schlumberger (2012), Eclipse reservoir simulation software manual, [35] Shapiro, S. A., E. Huenges, and G. Borm (1997), Estimating the crust perme- ability from fluid-injection-induced seismic emission at the KTB site, Geophys. J. Int., 131, F15 F [36] Shearer, P., E. Hauksson, and G. Lin (2005), Southern California hypocenter relocation with waveform cross-correlation, part 2: Results using sourcespecific station terms and cluster analysis, Bull. Seismol. Soc. Am., 95(3), [37] Skoumal, R. J., M. R. Brudzinski, B. S. Currie, and J. Levy (2014), Optimizing multi-station earthquake template matching through reexamination of the Youngstown, Ohio, sequence, Earth Plant. Sc. Lett., doi: /j.epsl
39 [38] Townend, J., and M. D. Zoback (2000), How faulting keeps the crust strong, Geology, 28(5), [39] Utsu, T., Y. Ogata, and M. Ritsuko (1965), The centenary of Omori formula for a decay law of afterhock activity, Journal of Physics of the Earth, 43, [40] Zhao, L.-S., and D. V. Helmberger (1994), Source estimation from broadband regional seismograms, Bullet. Seismol. Soc. Am., 84(1), [41] Zhu, L., and D. V. Helmberger (1996), Advancement in source estimation tech- niques using broadband regional seismograms, Bullet. Seismol. Soc. Am., 86(5), [42] Zoback, M. D., and H.-P. Harjes (1997), Injection-induced earthquakes and crustal stress at 9 km depth at the KTB deep drilling site, Germany, J. Geophys. Res., 102(B8), [43] Zoback, M. D., and J. Townend (2001), Implications of hydrostatic pore pres- sures and high crustal strength for the deformation of intraplate lithosphere, Tectonophysics, 336(1),
Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California T. H. W. Goebel, 1* S. M. Hosseini, 2
More informationSUPPLEMENTAL INFORMATION
GSA DATA REPOSITORY 2013310 A.M. Thomas et al. MOMENT TENSOR SOLUTIONS SUPPLEMENTAL INFORMATION Earthquake records were acquired from the Northern California Earthquake Data Center. Waveforms are corrected
More informationEmpirical Green s Function Analysis of the Wells, Nevada, Earthquake Source
Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad
More informationof other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out
4. Focal Mechanism Solutions A way to investigate source properties of the 2001 sequence is to attempt finding well-constrained focal mechanism solutions to determine if they are consistent with those
More informationA comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses
A comparison of seismicity rates and fluid injection operations in Oklahoma and California: Implications for crustal stresses Thomas Goebel California Institute of Technology Seismological Laboratory,
More informationLimitations of Earthquake Triggering Models*
Limitations of Earthquake Triggering Models* Peter Shearer IGPP/SIO/U.C. San Diego September 16, 2009 Earthquake Research Institute * in Southern California Why do earthquakes cluster in time and space?
More informationCharacterization of Induced Seismicity in a Petroleum Reservoir: A Case Study
Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study Edmond Sze, M. Nafi Toksöz, and Daniel R. Burns Earth Resources Laboratory Dept. of Earth, Atmospheric and Planetary Sciences
More informationA century of oil-field operations and earthquakes in the greater Los Angeles Basin, southern California
A century of oil-field operations and earthquakes in the greater Los Angeles Basin, southern California Egill Hauksson 1, Thomas Göbel 1, Jean-Paul Ampuero 1, and Elizabeth Cochran 2 Downloaded 06/03/15
More informationOverview of the Seismic Source Characterization for the Palo Verde Nuclear Generating Station
Overview of the Seismic Source Characterization for the Palo Verde Nuclear Generating Station Scott Lindvall SSC TI Team Lead Palo Verde SSC SSHAC Level 3 Project Tuesday, March 19, 2013 1 Questions from
More informationAutomatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and
More informationAn Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events
An Investigation on the Effects of Different Stress Regimes on the Magnitude Distribution of Induced Seismic Events Afshin Amini, Erik Eberhardt Geological Engineering, University of British Columbia,
More informationKinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy)
Kinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy) Buttinelli M., Improta L., Bagh S., Chiarabba C. 1/10 The Val
More informationOn May 4, 2001, central Arkansas experienced an M=4.4 earthquake followed by a
1. Introduction On May 4, 2001, central Arkansas experienced an M=4.4 earthquake followed by a surprisingly large number of small earthquakes. We recorded about 2500 above the ambient noise level on a
More informationOverview of Seismic Source Characterization for the Diablo Canyon Power Plant
Overview of Seismic Source Characterization for the Diablo Canyon Power Plant Steve Thompson (LCI and SSC TI Team), for SWUS GMC Workshop 1, March 19, 2013 Questions from TI Team Summarize tectonic setting.
More informationto: Interseismic strain accumulation and the earthquake potential on the southern San
Supplementary material to: Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system by Yuri Fialko Methods The San Bernardino-Coachella Valley segment of the
More informationA comparison of seismicity rates and fluid-injection operations in Oklahoma and California: Implications for crustal stresses
A comparison of seismicity rates and fluid-injection operations in Oklahoma and California: Implications for crustal stresses Thomas Göbel 1 Downloaded 06/03/15 to 130.118.44.150. Redistribution subject
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationUpdated Graizer-Kalkan GMPEs (GK13) Southwestern U.S. Ground Motion Characterization SSHAC Level 3 Workshop 2 Berkeley, CA October 23, 2013
Updated Graizer-Kalkan GMPEs (GK13) Southwestern U.S. Ground Motion Characterization SSHAC Level 3 Workshop 2 Berkeley, CA October 23, 2013 PGA Model Our model is based on representation of attenuation
More informationDETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION IN A PRODUCING INDONESIAN GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 9-11, 2009 SGP-TR-187 DETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION
More informationImaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges
2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara LONG-PERIOD (3 TO 10 S) GROUND MOTIONS IN AND AROUND THE
More informationEarthquakes and Seismotectonics Chapter 5
Earthquakes and Seismotectonics Chapter 5 What Creates Earthquakes? The term Earthquake is ambiguous: Applies to general shaking of the ground and to the source of the shaking We will talk about both,
More informationMIGRATING SWARMS OF BRITTLE-FAILURE EARTHQUAKES IN THE LOWER CRUST BENEATH MAMMOTH MOUNTAIN, CALIFORNIA
MIGRATING SWARMS OF BRITTLE-FAILURE EARTHQUAKES IN THE LOWER CRUST BENEATH MAMMOTH MOUNTAIN, CALIFORNIA David Shelly and David Hill GRL, October 2011 Contents Tectonic Setting Long Valley Caldera Mammoth
More informationWhat We Know (and don t know)
What We Know (and don t know) about the M5.1 La Habra Earthquake Dr. Robert Graves U.S. Geological Survey In collaboration with: Dr. Egill Hauksson and Dr. Thomas Göbel Caltech Dr. Elizabeth Cochran U.S.
More information3D Finite Element Modeling of fault-slip triggering caused by porepressure
3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a
More informationShort Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake
Bulletin of the Seismological Society of America, Vol. 100, No. 6, pp. 3269 3277, December 2010, doi: 10.1785/0120100087 Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood,
More informationAftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007586, 2010 Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress
More informationThe Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring
The Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring Bruce R. Julian U. S. Geological Survey, Menlo Park, CA 94025 USA julian@usgs.gov Gillian R. Foulger Dept. Earth
More informationDeformation of Rocks. Orientation of Deformed Rocks
Deformation of Rocks Folds and faults are geologic structures caused by deformation. Structural geology is the study of the deformation of rocks and its effects. Fig. 7.1 Orientation of Deformed Rocks
More informationEarthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College
Earthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College What Is an Earthquake? An earthquake is the vibration of Earth, produced by the rapid release of energy.
More informationChapter 15 Structures
Chapter 15 Structures Plummer/McGeary/Carlson (c) The McGraw-Hill Companies, Inc. TECTONIC FORCES AT WORK Stress & Strain Stress Strain Compressive stress Shortening strain Tensional stress stretching
More informationEarthquakes. Building Earth s Surface, Part 2. Science 330 Summer What is an earthquake?
Earthquakes Building Earth s Surface, Part 2 Science 330 Summer 2005 What is an earthquake? An earthquake is the vibration of Earth produced by the rapid release of energy Energy released radiates in all
More informationToday: Basic regional framework. Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure
Today: Basic regional framework Western U.S. setting Eastern California Shear Zone (ECSZ) 1992 Landers EQ 1999 Hector Mine EQ Fault structure 1 2 Mojave and Southern Basin and Range - distribution of strike-slip
More informationEarthquakes. Chapter Test A. Multiple Choice. Write the letter of the correct answer on the line at the left.
Earthquakes Chapter Test A Multiple Choice Write the letter of the correct answer on the line at the left. 1. Stress that pushes a mass of rock in two opposite directions is called a. shearing. b. tension.
More informationRegional tectonic stress near the San Andreas fault in central and southern California
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S11, doi:10.1029/2003gl018918, 2004 Regional tectonic stress near the San Andreas fault in central and southern California J. Townend School of Earth Sciences,
More informationCHAPTER 1 BASIC SEISMOLOGY AND EARTHQUAKE TERMINOLGY. Earth Formation Plate Tectonics Sources of Earthquakes...
CHAPTER 1 BASIC SEISMOLOGY AND EARTHQUAKE TERMINOLGY Earth Formation... 1-2 Plate Tectonics... 1-2 Sources of Earthquakes... 1-3 Earth Faults... 1-4 Fault Creep... 1-5 California Faults... 1-6 Earthquake
More informationA mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S13, doi:10.1029/2004gl019521, 2004 A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements Jean Chéry Laboratoire Dynamique de la
More informationWhat allows seismic events to grow big?: Insights from fault roughness and b-value analysis in stick-slip experiments
What allows seismic events to grow big?: Insights from fault roughness and b-value analysis in stick-slip experiments T. H. W. Goebel 1, G. Kwiatek 2, T. W. Becker 3, G. Dresen 2, E.E. Brodsky 1 1 UC Santa
More informationAddressing the risks of induced seismicity in sub-surface energy operations
Addressing the risks of induced seismicity in sub-surface energy operations Richard Porter a, Alberto Striolo a, Haroun Mahgerefteh a, Joanna Faure Walker b a Department of Chemical Engineering, University
More informationSeismic Reflection Views of the 1994 Northridge Earthquake Hypocentral Region Using Aftershock Data and Imaging Techniques
Seismic Reflection Views of the 1994 Northridge Earthquake Hypocentral Region Using Aftershock Data and Imaging Techniques SERGIO CHÁVEZ-PÉREZ and JOHN N. LOUIE SEISMOLOGICAL LABORATORY/174, MACKAY SCHOOL
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationEarthquakes Chapter 19
Earthquakes Chapter 19 Does not contain complete lecture notes. What is an earthquake An earthquake is the vibration of Earth produced by the rapid release of energy Energy released radiates in all directions
More informationCalifornia foreshock sequences suggest aseismic triggering process
GEOPHYSICAL RESEARCH LETTERS, VOL. 0, 1, doi:10.100/grl.50, 013 California foreshock sequences suggest aseismic triggering process Xiaowei Chen 1 and Peter M. Shearer 1 Received January 013; revised 9
More informationEarthquake patterns in the Flinders Ranges - Temporary network , preliminary results
Earthquake patterns in the Flinders Ranges - Temporary network 2003-2006, preliminary results Objectives David Love 1, Phil Cummins 2, Natalie Balfour 3 1 Primary Industries and Resources South Australia
More informationstress direction are less stable during both drilling and production stages (Zhang et al., 2006). Summary
Inversion and attribute-assisted hydraulically-induced microseismic fracture prediction: A North Texas Barnett Shale case study Xavier E. Refunjol *, Katie M. Keranen, and Kurt J. Marfurt, The University
More informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationMultiscale Seismic Signature of a Small Fault Zone in a Carbonate Reservoir: Relationships Between V P Imaging, Fault Zone Architecture and Cohesion*
Multiscale Seismic Signature of a Small Fault Zone in a Carbonate Reservoir: Relationships Between V P Imaging, Fault Zone Architecture and Cohesion* Pierre Jeanne 1,2, Yves Guglielmi 2, and Frédéric Cappa
More informationFOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA
FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA Siti Norbaizura MAT SAID Supervisor: Tatsuhiko HARA MEE10505 ABSTRACT Since November 30, 2007, small local earthquakes have been observed
More informationA TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES. Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2
Printed: September 1, 2005 A TESTABLE FIVE-YEAR FORECAST OF MODERATE AND LARGE EARTHQUAKES IN SOUTHERN CALIFORNIA BASED ON SMOOTHED SEISMICITY Yan Y. Kagan 1,David D. Jackson 1, and Yufang Rong 2 1 Department
More informationEarthquakes Earth, 9th edition, Chapter 11 Key Concepts What is an earthquake? Earthquake focus and epicenter What is an earthquake?
1 2 3 4 5 6 7 8 9 10 Earthquakes Earth, 9 th edition, Chapter 11 Key Concepts Earthquake basics. "" and locating earthquakes.. Destruction resulting from earthquakes. Predicting earthquakes. Earthquakes
More informationModule 7: Plate Tectonics and Earth's Structure Topic 4 Content : Earthquakes Presentation Notes. Earthquakes
Earthquakes 1 Topic 4 Content: Earthquakes Presentation Notes Earthquakes are vibrations within the Earth produced by the rapid release of energy from rocks that break under extreme stress. Earthquakes
More informationrevised October 30, 2001 Carlos Mendoza
Earthquake Sources in the circum-caribbean Region Puerto Rico Tsunami Mitigation and Warning Program Federal Emergency Management Agency Preliminary Report: Task 3 revised October 30, 2001 Carlos Mendoza
More informationGeologic Structures. Changes in the shape and/or orientation of rocks in response to applied stress
Geologic Structures Changes in the shape and/or orientation of rocks in response to applied stress Figure 15.19 Can be as big as a breadbox Or much bigger than a breadbox Three basic types Fractures >>>
More informationEarthquakes in Barcelonnette!
Barcelonnette in the Ubaye valley : the landscape results of large deformations during the alpine orogene (40 5 Myr in this area) and the succession of Quaternary glaciations. The sedimentary rocks are
More informationToward interpretation of intermediate microseismic b-values
Toward interpretation of intermediate microseismic b-values Abdolnaser Yousefzadeh, Qi Li, Schulich School of Engineering, University of Calgary, Claudio Virues, CNOOC- Nexen, and Roberto Aguilera, Schulich
More informationAnswers: Internal Processes and Structures (Isostasy)
Answers: Internal Processes and Structures (Isostasy) 1. Analyse the adjustment of the crust to changes in loads associated with volcanism, mountain building, erosion, and glaciation by using the concept
More informationEARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE MECHANISMS SHOW MOTION
6-1 6: EARTHQUAKE FOCAL MECHANISMS AND PLATE MOTIONS Hebgen Lake, Montana 1959 Ms 7.5 1 Stein & Wysession, 2003 Owens Valley, California 1872 Mw ~7.5 EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE
More informationFULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
More informationSection 19.1: Forces Within Earth Section 19.2: Seismic Waves and Earth s Interior Section 19.3: Measuring and Locating.
CH Earthquakes Section 19.1: Forces Within Earth Section 19.2: Seismic Waves and Earth s Interior Section 19.3: Measuring and Locating Earthquakes Section 19.4: Earthquakes and Society Section 19.1 Forces
More informationSome aspects of seismic tomography
Some aspects of seismic tomography Peter Shearer IGPP/SIO/U.C. San Diego September 7, 2009 Earthquake Research Institute Part 1: Global Tomography P velocity perturbations 200 km 1000 km 2700 km MIT 2006
More informationLecture # 6. Geological Structures
1 Lecture # 6 Geological Structures ( Folds, Faults and Joints) Instructor: Dr. Attaullah Shah Department of Civil Engineering Swedish College of Engineering and Technology-Wah Cantt. 2 The wavy undulations
More informationVariations in Tremor Activity and Implications for Lower Crustal Deformation Along the Central San Andreas Fault
Variations in Tremor Activity and Implications for Lower Crustal Deformation Along the Central San Andreas Fault David R. Shelly USGS, Menlo Park Shelly and Hardebeck, GRL, 2010 Collaborators: Jeanne Hardebeck
More informationDi#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka
More informationCAP M S Wallace. Vol. 27 No. 2 Jun EARTHQUAKE RESEARCH IN CHINA M S 4. 8 CAP. 3km - - P315
27 2 207 ~ 214 2011 6 EARTHQUAKE RESEARCH IN CHINA Vol. 27 No. 2 Jun. 2011 2011 CAP 4. 8 27 2 207 ~ 214 CAP 4. 8 1 2 1 2 3 1 1 1 558 230031 2 96 230026 3 100081 CAP 2011 1 19 M S 4. 8 M W = 4. 1 I 16 74
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationModeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection
Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection S. Raziperchikolaee 1 and J. F. Miller 1 Abstract Analysis
More informationSeismological monitoring of the GRT1 hydraulic stimulation (Rittershoffen, Alsace, France)
Seismological monitoring of the GRT1 hydraulic stimulation (Rittershoffen, Alsace, France) Mohamed Boubacar 1, Olivier Lengliné 2 1. University of Quebec, INRS 2. University of Strasbourg, EOST Summary
More informationWhat Causes Rock to Deform?
Crustal Deformation Earth, Chapter 10 Chapter 10 Crustal Deformation What Causes Rock to Deform? Deformation is a general term that refers to all changes in the shape or position of a rock body in response
More informationRELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE
RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT
More informationHow will induced seismicity in Oklahoma respond to decreased saltwater injection rates?
SEISMOLOGY How will induced seismicity in Oklahoma respond to decreased saltwater injection rates? Cornelius Langenbruch* and Mark D. Zoback In response to the marked number of injection-induced earthquakes
More informationConstrained Fault Construction
Constrained Fault Construction Providing realistic interpretations of faults is critical in hydrocarbon and mineral exploration. Faults can act as conduits or barriers to subsurface fluid migration and
More informationDoes Aftershock Duration Scale With Mainshock Size?
GEOPHYSICAL RESEARCH LETTERS, VOL.???, NO., PAGES 1 16, Does Aftershock Duration Scale With Mainshock Size? A. Ziv A. Ziv, Ben-Gurion University of the Negev, Beer-Sheva, 84105, Israel. (e-mail: zival@bgu.ac.il)
More informationInduced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France. Jean Schmittbuhl 1
Induced seismicity in Enhanced Geothermal Systems (EGS) in Alsace, France Jean Schmittbuhl 1 Olivier Lengliné 1, Mohamed Boubacard 1, Lucie Lamourette 1, Lilas Vivin 1, Nicolas Cuenot 2 1 EOST-IPGS, Université
More informationHigh Resolution Imaging of Fault Zone Properties
Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California
More informationVelocity contrast along the Calaveras fault from analysis of fault zone head waves generated by repeating earthquakes
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L01303, doi:10.1029/2007gl031810, 2008 Velocity contrast along the Calaveras fault from analysis of fault zone head waves generated by
More informationSection Forces Within Earth. 8 th Grade Earth & Space Science - Class Notes
Section 19.1 - Forces Within Earth 8 th Grade Earth & Space Science - Class Notes Stress and Strain Stress - is the total force acting on crustal rocks per unit of area (cause) Strain deformation of materials
More informationMeasurements in the Creeping Section of the Central San Andreas Fault
Measurements in the Creeping Section of the Central San Andreas Fault Introduction Duncan Agnew, Andy Michael We propose the PBO instrument, with GPS and borehole strainmeters, the creeping section of
More informationStress and Strain. Stress is a force per unit area. Strain is a change in size or shape in response to stress
Geologic Structures Geologic structures are dynamically-produced patterns or arrangements of rock or sediment that result from, and give information about, forces within the Earth Produced as rocks change
More informationExtending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks
Extending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks Gisela Viegas*, ESG, Kingston, Ontario, Canada Gisela.Fernandes@esgsolutions.com and
More informationChapter 6. Conclusions. 6.1 Conclusions and perspectives
Chapter 6 Conclusions 6.1 Conclusions and perspectives In this thesis an approach is presented for the in-situ characterization of rocks in terms of the distribution of hydraulic parameters (called SBRC
More informationMicroseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties of rocks
GEOPHYSICS, VOL. 68, NO. 2 (MARCH-APRIL 2003); P. 685 689, 5 FIGS. 10.1190/1.1567239 Short Note Microseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties
More informationEarthquakes.
Earthquakes http://quake.usgs.gov/recenteqs/latestfault.htm An earthquake is a sudden motion or shaking of the Earth's crust, caused by the abrupt release of stored energy in the rocks beneath the surface.
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING EVENTS
ABSTRACT ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING EVENTS Don V. Helmberger 1, Jeroen Tromp 1, and Arthur J. Rodgers 2 California Institute of Technology 1 and Lawrence Livermore National Laboratory
More informationScience Starter. Describe in your own words what an Earthquake is and what causes it. Answer The MSL
Science Starter Describe in your own words what an Earthquake is and what causes it. Answer The MSL WHAT IS AN EARTHQUAKE AND HOW DO WE MEASURE THEM? Chapter 8, Section 8.1 & 8.2 Looking Back Deserts Wind-shaped
More informationGEOLOGY MEDIA SUITE Chapter 13
UNDERSTANDING EARTH, SIXTH EDITION GROTZINGER JORDAN GEOLOGY MEDIA SUITE Chapter 13 Earthquakes 2010 W.H. Freeman and Company Three different types of seismic waves are recorded by seismographs Key Figure
More informationCrustal Deformation. Earth Systems 3209
Crustal Deformation Earth Systems 3209 Crustal Deformation pg. 415 Refers to all changes in the original form and/or size of a rock body. May also produce changes in the location and orientation of rocks.
More informationApparent and True Dip
Apparent and True Dip Cross-bedded building stone. The contact immediately below A appears to dip gently to the right, but at B, the contact appears to dip to the left. But it's not a syncline! Both of
More informationComprehensive analysis of earthquake source spectra and swarms in the Salton Trough, California
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2011jb008263, 2011 Comprehensive analysis of earthquake source spectra and swarms in the Salton Trough, California X. Chen 1 and P. M. Shearer 1
More informationA systematic measurement of shear wave anisotropy in southern California Report for SCEC Award #15081 Submitted March 15, 2016
A systematic measurement of shear wave anisotropy in southern California Report for SCEC Award #15081 Submitted March 15, 2016 Investigators: Zhigang Peng (Georgia Tech) I. Project Overview... i A. Abstract...
More informationEarth Science, (Tarbuck/Lutgens) Chapter 10: Mountain Building
Earth Science, (Tarbuck/Lutgens) Chapter 10: Mountain Building 1) A(n) fault has little or no vertical movements of the two blocks. A) stick slip B) oblique slip C) strike slip D) dip slip 2) In a(n) fault,
More informationTectonic Seismogenic Index of Geothermal Reservoirs
Tectonic Seismogenic Index of Geothermal Reservoirs C. Dinske 1, F. Wenzel 2 and S.A. Shapiro 1 1 Freie Universität Berlin 2 KIT Karlsruhe November 27, 2012 Introduction M max reservoir location -2.0 Barnett
More informationPlate Tectonics. entirely rock both and rock
Plate Tectonics I. Tectonics A. Tectonic Forces are forces generated from within Earth causing rock to become. B. 1. The study of the origin and arrangement of Earth surface including mountain belts, continents,
More informationSeismogenic structure of 1935 Hsinchu-Taichung (M GR =7.1) earthquake, Miaoli, western Taiwan 1935 (M GR =7.1)
Seismogenic structure of 1935 Hsinchu-Taichung (M GR =7.1) earthquake, Miaoli, western Taiwan 1935 (M GR =7.1) Y.N. Nina Lin; Y.G. Chen; Y.M. Wu (Inst. of Geosciences, NTU); K.M. Yang (Exploration and
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11492 Figure S1 Short-period Seismic Energy Release Pattern Imaged by F-net. (a) Locations of broadband seismograph stations in Japanese F-net used for the 0.5-2.0 Hz P wave back-projection
More informationSendai Earthquake NE Japan March 11, Some explanatory slides Bob Stern, Dave Scholl, others updated March
Sendai Earthquake NE Japan March 11, 2011 Some explanatory slides Bob Stern, Dave Scholl, others updated March 14 2011 Earth has 11 large plates and many more smaller ones. Plates are 100-200 km thick
More informationObservational analysis of correlations between aftershock productivities and regional conditions in the context of a damage rheology model
Geophys. J. Int. (2009) 177, 481 490 doi: 10.1111/j.1365-246X.2009.04145.x FAST TRACK PAPER Observational analysis of correlations between aftershock productivities and regional conditions in the context
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO177 The Long Precursory Phase of Most Large Interplate Earthquakes Supplementary Information Supplementary Methods 1. Data and Classification We select the earthquakes
More informationSupporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system
GEOPHYSICAL RESEARCH LETTERS Supporting Information for Break of slope in earthquake-size distribution reveals creep rate along the San Andreas fault system Inessa Vorobieva, 1,2 Peter Shebalin, 1,2 Clément
More informationIntroduction Faults blind attitude strike dip
Chapter 5 Faults by G.H. Girty, Department of Geological Sciences, San Diego State University Page 1 Introduction Faults are surfaces across which Earth material has lost cohesion and across which there
More information