UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions
|
|
- Henry Tyler
- 6 years ago
- Views:
Transcription
1 UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions Bruce E. Shaw Lamont Doherty Earth Observatory, Columbia University Statement of the Problem In UCERF2 [WGCEP, 2008] Magnitude-Area scaling relations contributed one of the main sources of uncertainty in the final hazard estimates. In this task, we seek ways of finding additional constraints or alternative methods of using scaling which will reduce these uncertainties. Summary of Recommendations A thorough study was made of the two main uses of magnitude-area relations: one, to figure out the sizes of events, and two, to figure out the rates of events through moment or slip-rate balancing. The accompanying paper [Shaw, 2011b] discusses the issues and results more comprehensively. Here, we summarize the recommendations for UCERF3. With respect to moment or slip-rate balancing, qualitative consistency was found between the implied slip-length from magnitude-area scaling, and surface slip-length observations. However, some quantitative difference, of order 30%, remained, with the implied slip from magnitude-area being higher. This could be due to the deep coseismic slip being mapped onto shallower seismogenic layers, so that magnitude-area relations are overestimating average slip. Alternatively, it could be due to the surface slip measurements underestimating slip, perhaps due to distributed plastic deformation in unconsolidated layers or other such mechanisms. This quantitative discrepancy remains an unresolved epistemic uncertainty. The qualitative agreement, however, does give useful information in distinguishing different scaling laws. The accompanying paper discusses this more thoroughly. Given the epistemic uncertainties, we suggest two main branches. The difference in the two branches relates to how the slip-rate balancing is done. In all cases, we suggest continuing to use magnitude-area relations to estimate the sizes of events. Whatever the epistemic uncertainties associated with the depth dependence of slip, magnitude-area relations give good size estimates, and seismogenic area is a good proxy, and self-consistent when understood as such, for these relations. Three magnitude-area scaling relations are considered as most viable for estimating sizes of events: Ellsworth-B [WGCEP, 2003], Hanks-Bakun [Hanks and Bakun, 2002, 2008], and Shaw09 [Shaw, 2009]. The Shaw09 relation generalizes the Hanks and Bakun relation from two to three regimes, adding a third regime for the largest events with slip saturating with width W scaling. The relations are Ellsworth-B: M = log A +.2 (1) This has the virtue of being extremely simple, being a one parameter fit, and being a good empirical fit to the large earthquake (M > 6.5) data it was developed to fit. Hanks and Bakun [2002] developed a two regime saling: { log A A 537 km 2 ; M = 3 log A A>537 km2. This has the virtue of doing a better job fitting the data at smaller magnitudes. Shaw [2009] developed a generalization of the Hanks and Bakun model, which extends the bilinear Hanks-Bakun two-regime scaling to a three-regime scaling which has a third asymptotic regime valid for very long ruptures L W whereby S approaches W scaling asymptotically. This is done at the price of one additional (2) 1
2 scaling parameter (the scale at which the transition to the third regime occurs). The Shaw [2009] scaling relation is parameterized as follows. M = log 10 A + 2 A 3 log max(1, ) H const. (3) A (1 + max(1, ))/2 H 2 β Here H is the seismogenic depth, and β is a fitting parameter which gives a crossover scale length to the asymptotic W scaling. In the limit of β it reduces to the Hanks-Bakun scaling relation, and in this sense is a one parameter extension of it. Note also that in the limit as L it scales asymptotically as M 2 3 log 10A which gives S W (using W H), as the Knopoff [1958] solution suggests. Fits to the Hanks-Bakun data give best fitting parameters of values H = 16.1 km and β =6.6. Based on the whole magnitude range in the Hanks and Bakun database for magnitude-area of strike-slip events, Shaw09 slightly outperforms Hanks and Bakun on an AIC test, though they perform comparably, and both outperform Ellsworth-B. Based on M > 6.5, all three scaling relations perform comparably. This is based on misfits in magnitude-log area space, so all events regardless of size are weighted equally. Figure 1 shows the fit of the three scaling relations to the Hanks-Bakun data. Table 1 summarizes these results of the fits. Looking in linear weighted space, as is relevant for moment or slip-rate balancing, the picture changes substantially. There, the Hanks-Bakun relation is seen to be a poor fit of the largest events, with Shaw09 slightly outperforming Ellsworth-B, and both very much outperforming Hanks and Bakun. Figure 2a shows these results, plotting the magnitude-area data transformed into slip-length data, and plotted on a linear scale relevant for slip-rate balancing. Table 2 summarizes these results of the fits. Weights for sizes of events branches For UCERF3 purposes, all three magnitude-area scaling relations perform sufficiently for estimating sizes of events. Shaw09 performed well on both the log-weighting (across all magnitudes) and linear weighting (weighted towards largest magnitudes), while Hanks and Bakun did not perform well on the linear weighting and Ellsworth-B did not perform as well on the log weighting (though questions have been raised about the area estimates for the smaller magnitude events). Given these results, viable branch weights could be replacing the UCERF2 Hanks and Bakun weight with the Shaw09 relation which generalizes it. Alternatively, some weight on the Hanks and Bakun branch could be retained for continuity. Weights for Moment or Slip-Rate Balancing Two main branches are suggested. One branch continues to use magnitude-area scaling to estimate average slip, as was done in UCERF2. For this branch, we suggest using the Ellsworth-B and Shaw09 relations, since the Hanks and Bakun relation performed poorly on the linear weighting relevant to rate balancing. (See Table 2). For the second branch, slip rate would be balanced using slip-length scaling relations. Linear scaling with L out to the largest events was found to be a poor fit to the data. The best fits were for a constant stress drop model, which transitions from asymptotic L scaling to asymptotic W scaling, and for a square root L 1/2 model. The equations for these two models are as follows. For the square root scaling consistent with the implied Ellsworth-B scaling, S A 1/2 = C (LW ) 1/2 () 2
3 with C a fitting constant. For the constant stress drop model [Shaw, 2011b] consistent with the Shaw09 scaling S = σ 1 (5) µ 7 3L + 1 2W with µ the shear modulus and σ the stress drop. This relationship comes from a more careful treatment of how slip would be expected to transition from circular ruptures just beginning to break the surface to long rectangular ruptures. Best fit values of the constants C and σ were found for fits to the Wesnousky [2008] data in the associated paper Shaw [2011b], with C = ,and assuming µ =30GP a abestfit σ =3.91 MPa. These constants are expected to change slightly when fit to the updated UCERF3 large event dataset. The constant stress drop model somewhat outperformed the L 1/2 scaling, but they both performed well in fitting the data. Figure 2b shows the fits of these slip-length functions to surface slip data. Table 3 summarizes the results of these fits. Based on the fits, and a desire to span epistemic uncertainty, we suggest equal weighting to these two models. W scaling For use with variable downdip width W, the dependence on W has been included in the slip-length scaling relations, Equations () and (5). (See paper for details). Note that the scalings imply different W dependence. Looking at the largest events to try to distinguish between the candidate scalings, the sparse noisy data favored W scaling over W 1/2 scaling (See Figure 6 in paper). New datasets Data sets examined thus far are the Hanks and Bakun [2008] magnitude-area dataset and the Wesnousky [2008] surface slip dataset. The updated dataset for large earthquakes being compiled for UCERF3 would be very useful in further refining weightings for scaling relations. The basic framework is not expected to change, but the relative performance in fitting the data might, and this could impact recommended branch weights. Epistemic Uncertainty Currently, to fits for the two types of slip-length relations, those derived from implied magnitude-area scaling laws, and those derived from fits to surface slip data, span the range of epistemic uncertainty. This epistimic undertainty is illustrated in Figure 3, which shows superposed from Figure 2a the implied magnitude area scaling laws with dashed lines, and from Figure 2b the surface slip scaling laws with solid lines. Newer datasets are expected to raise somewhat the amplitude of the surface slip measurement scaling laws. Cybershake Cybershake [Graves et al., 2010] is somewhat outside the scope of UCERF3, but it is, nevertheless an interested user, and potentially can develop into a useful constraint. Some of the branches in UCERF3 will be more compatible with cybershake than others. Since cybershake itself has many underlying assumptions (such as the degree of coherence in the ruptures), sensitivities in cybershake need to be explored further before it is used as a constraint on UCERF relations. Nevertheless, it is a goal to have compatibility developed, and thus exploring the implications of the various relations in this context is an important task. An unresolved epistemic uncertainty in the depth dependence of slip concerns how much coseismic slip may be occurring below the seismogenic layer. Dynamic models suggest a significant fraction of 3
4 slip may be occurring below the seismogenic layer up to of order 1/3 the total moment [Shaw and Wesnousky, 2008]. This slip would be occurring mainly as long period motion, with a dearth of high frequencies [Shaw and Wesnousky, 2008]. An exponential decay or a linear decay with depth below the seismogenic layer are both reasonable parameterizations of what is seen in the dynamic models [Shaw and Wesnousky, 2008]. The scale length of the decay depends on the degree of velocity strengthening; the degree of strengthening is not well constrained, and thus the scale length remains an epistemic uncertainty. The distributions of slip, and rupture velocities, are important factors which warrant further study in cybershake. Correlations of slip at the surface have gaussian, or even faster than gaussian falloff in the tails [Shaw, 2011a], so cybershake simulations positing extreme slip values in the distributions (e.g. 50m slip in great strike-slip earthquake patches) should be reexamined. Given the implied high stress drops at moderate magnitudes in the Ellsworth-B relations, and the implied high stress drops at large magnitudes in the Hanks-Bakun relations, it would be interesting to see if constant stress drop models can be seen to be compatible with cybershake. Is, for example, the Shaw09 magnitude-area scaling relation, combined with a 30% subseismogenic high-frequency deficient moment, compatible with more high variance rupture velocity low variance slip patch versions of cybershake? References Graves, R., et al., Cybershake: A physics-based seismic hazard model for Southern California, Pageoph, Hanks, T. C., and W. H. Bakun, A blinear source-scaling model for M-log A observations of continental earthquakes, Bull. Seismol. Soc. Am., 92, 181, Hanks, T. C., and W. H. Bakun, M-log A observations of recent large earthquakes, Bull. Seismol. Soc. Am., 98, 90, Knopoff, L., Energy release in earthquakes, Geophys. J. R. A. S., 1,, Shaw, B. E., Constant stress drop from small to great earthquakes in magnitude-area scaling, Bull. Seismol. Soc. Am., 99, 871, Shaw, B. E., Surface slip gradients of large earthquakes, Bull. Seismol. Soc. Am., 101, 792, 2011a. Shaw, B. E., An alternative pathway for seismic hazard estimates based on slip-length scaling, submitted, 2011b. Shaw, B. E., and S. G. Wesnousky, Slip-length scaling in large earthquakes: The role of deep penetrating slip below the seismogenic layer, Bull. Seismol. Soc. Am., 98, 1633, Wesnousky, S. G., Displacement and geometrical characteristics of earthquake surface ruptures: Issues and implications for seismic-hazard analysis and the process of earthquake rupture, Bull. Seismol. Soc. Am., 98, 1609, WGCEP, Earthquake probabilities in the San Francisco Bay Region: 2002 to 2031, U.S. Geol. Surv. Open File Rep., 03-21, WGCEP, The Uniform California Earthquake Rupture Forecast, version 2 (ucerf2), U.S. Geol. Surv. Open File Rep., , 2008.
5 Tables of Fits of Scalings to Data Magnitude-Area Scaling Table 1: Magnitude-Area Scaling Fits Scaling Best fits Parameters Name Eqn Std. dev. AIC No. Best fit values S09 (3) C 0 =3.98 W =16.1 km β =6.6 HB (2) C 0 =3.98 A c = 537 km 2 EB (1) C 0 =.20 The standard deviation in the above table measures the difference in magnitudes of the data from the predicted curve. The data comes from the Hanks-Bakun database [Hanks and Bakun, 2008], which is reproduced in the electronic supplement S1. The number of free parameters is given by No., with the best fitted parameter values following that. The parameter C 0 is the constant in the magnitudearea scaling relations. Ranking in this table, and those that follow, is by AIC, which rewards reduced standard deviation but penalizes additional parameters. Slip-Length Scaling from Implied Magnitude Area Data Table 2: Implied Slip-Length from Magnitude-Area Scaling Best fits Parameters Name Eqn Std. dev. AIC No. Values S09 (3) C 0 =3.98 Wβ = km EB (1) C 0 =.20 HB (2) C 0 =3.98 A c = 537 km 2 The standard deviation in the above table measures the difference in slip of the data from the predicted curve. The data comes from the Hanks-Bakun database [Hanks and Bakun, 2008], which is transformed from magnitude-area to slip-length as follows. To get slip we convert magnitude to moment and divide by area and modulus to get slip. To get length we divide area by width, assuming width is the lesser of the length compared with the downdip seismogenic width. The table shows a comparison of the transformed data with scaling relations derived from transformed magnitude area to implied slip-length. scaling relations by a constant factor ρ. 5
6 Slip-Length Scaling from Surface Slip Data Table 3: Slip-Length from Surface Slip Data Scaling Best fits Parameters Name Eqn Std. dev. AIC No. Best fit values S11 (5) σ =3.91 MPa L 1/2 () C = L C 5 = km 1 The standard deviation in the above table measures the difference in slip of the data from the predicted curve. The data is derived from the Wesnousky database for surface slip in large Strike-Slip events [Wesnousky, 2008]. The table shows a comparison of the surface slip data with scaling relations proposed for surface slip-length scaling. The best overall fit to the data, in terms of AIC, is S11, the constant stress drop model. 6
7 9 8 M log A [km 2 ] 10 Figure 1: Magnitude area relations for large strike-slip events. Red circles denote magnitude and area of events from Hanks and Bakun [2008] database. Solid black line is linear Ellsworth-B [WGCEP, 2003] magnitude-area relation, Equation (1). Dashed green line is Hanks and Bakun [2002] bilinear relation, Equation (2). Blue line is new generalized scaling relation, Equation (3) [Shaw, 2009]. 7
8 S[m] 6 S[m] (a) L[km] (b) L[km] Figure 2: Fits of slip-length scaling. (a) Slip-length scaling implied by magnitude-area data. Data (light blue circles) derived from Hanks and Bakun [2008] mag-area data by rescaling as follows. Magnitude is converted to moment, then divided by area and modulus to get slip. Area is converted to length by dividing area by width assuming seismogenic depth H=15 km. Different color lines represent different magnitude-area scaling laws rescaled the same way the data has been. The curves are the implied slip-length for: Ellsworth-B [WGCEP, 2003] (red), Hanks-Bakun [Hanks and Bakun, 2008] (green), and Shaw [2009] (black). (b) Fits of surface sliplength scaling laws to surface slip observations of average slip versus length. Color indicates focal mechanism: strike-slip (blue), normal (cyan), thrust (magenta). In order of best fit to least good fit: black line shows Equation (5) S11 scaling; red line shows Equation () L 1/2 scaling. Note change in vertical axis scale in (a) versus (b), reflecting overall systematic difference, and fundamental epistemic uncertainty, in average slip estimates from magnitude-area versus surface slip observations. Some of this difference reflects different events in the two different databases, but some of it, a difference of around 30%, remains when we restrict to looking at the same events S[m] L[km] Figure 3: Slip length scaling relations for large strike-slip events. Dashed lines are implied slip-length scaling relations derived from Magnitude-Area scaling relations. Solid lines are slip-length scaling relations derived from surface slip observations. Note that the dashed lines are systematically above the solid lines. Dashed cyan line is implied slip-length from Ellsworth-B [WGCEP, 2003] magnitude-area relation, Equation (1). Dashed green line is implied slip-length from Hanks and Bakun [2002] bilinear relation, Equation (2). Dashed blue line is implied slip-length from [Shaw, 2009], Equation (3). Solid red line is L 1/2 scaling, Equation (). Solid black line is constant stress drop scaling [Shaw, 2011b], Equation (5). The dashed blue line best fits the slip-length data derived from the magnitude-area data. The solid black line best fits the surface slip-length data. 8
UCERF3 Task R2- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions
UCERF3 Task R- Evaluate Magnitude-Scaling Relationships and Depth of Rupture: Proposed Solutions Bruce E. Shaw Lamont Doherty Earth Observatory, Columbia University Statement of the Problem In UCERF Magnitude-Area
More informationProbabilities for Jumping Fault Segment Stepovers
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:1.129/, Probabilities for Jumping Fault Segment Stepovers Bruce E. Shaw Lamont-Doherty Earth Observatory, Columbia University, New York James H. Dieterich
More informationMagnitude-Area Scaling of Strike-Slip Earthquakes. Paul Somerville, URS
Magnitude-Area Scaling of Strike-Slip Earthquakes Paul Somerville, URS Scaling Models of Large Strike-slip Earthquakes L Model Scaling (Hanks & Bakun, 2002) Displacement grows with L for L > > Wmax M
More informationGLOBAL SOURCE PARAMETERS OF FINITE FAULT MODEL FOR STRONG GROUND MOTION SIMULATIONS OR PREDICTIONS
13 th orld Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2743 GLOBAL SOURCE PARAMETERS OF FINITE FAULT MODEL FOR STRONG GROUND MOTION SIMULATIONS OR PREDICTIONS
More informationGeophysical Research Letters
Geophysical Research Letters 16 JANUARY 2007 Volume 34 Number 1 American Geophysical Union Toward a more objective determination of seismic hazards Andaman Islands deformation associated with the 2004
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 informationAppendix O: Gridded Seismicity Sources
Appendix O: Gridded Seismicity Sources Peter M. Powers U.S. Geological Survey Introduction The Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) is a forecast of earthquakes that fall
More information6 Source Characterization
6 Source Characterization Source characterization describes the rate at which earthquakes of a given magnitude, and dimensions (length and width) occur at a given location. For each seismic source, the
More informationSlip-Length Scaling in Large Earthquakes: The Role of Deep-Penetrating Slip below the Seismogenic Layer
Bulletin of the Seismological Society of America, Vol. 98, No. 4, pp. 633 64, August 2008, doi: 0.785/0200709 Slip-Length Scaling in Large Earthquakes: The Role of Deep-Penetrating Slip below the Seismogenic
More informationDepth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24301, doi:10.1029/2008gl036249, 2008 Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
More informationSurface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes
Surface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes Brad Aagaard May 21, 2009 Surface Rupture from Kinematic Rupture Models Objective Holistic approach to calculating surface
More informationGround displacement in a fault zone in the presence of asperities
BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 2, pp. 95-110; JUNE 2000 Ground displacement in a fault zone in the presence of asperities S. SANTINI (1),A.PIOMBO (2) and M. DRAGONI (2) (1) Istituto
More informationPeter Shearer 1, Robin Matoza 1, Cecily Wolfe 2, Guoqing Lin 3, & Paul Okubo 4
Characterizing fault zones and volcanic conduits at Kilauea and Mauna Loa volcanoes by large-scale mapping of earthquake stress drops and high precision relocations Peter Shearer 1, Robin Matoza 1, Cecily
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 informationGPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise
GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise Strain Analysis Introduction Name: The earthquake cycle can be viewed as a process of slow strain accumulation
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 18, 16 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Earthquake Magnitude and Moment Brune Stress Drop Seismic Spectra & Earthquake Scaling laws Scaling and
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
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 informationSouthern California ground motion envelopes over ranges of magnitudes, distances, and site conditions
55 Chapter 3 Average properties of Southern California ground motion envelopes over ranges of magnitudes, distances, and site conditions In this chapter, I use the envelope attenuation relationships derived
More informationA Bilinear Source-Scaling Model for M log A Observations of Continental Earthquakes
Bulletin of the Seismological Society of America, Vol. 9, No. 5, pp. 1841 1846, June A Bilinear Source-Scaling Model for M log A Observations of Continental Earthquakes by Thomas C. Hanks and William H.
More informationWAACY Magnitude PSF Model (Wooddell, Abrahamson, Acevedo-Cabrera, and Youngs) Norm Abrahamson DCPP SSC workshop #3 Mar 25, 2014
WAACY Magnitude PSF Model (Wooddell, Abrahamson, Acevedo-Cabrera, and Youngs) Norm Abrahamson DCPP SSC workshop #3 Mar 25, 2014 UCERF3 Allowing larger magnitudes by linking rupture segments Grand inversion
More informationIntroduction to Displacement Modeling
Introduction to Displacement Modeling Introduction Deformation on the Earth surface informs us about processes and material properties below surface Observation tools: GPS (static-dynamic) InSAR (static)
More informationThe Length to Which an Earthquake will go to Rupture. University of Nevada, Reno 89557
The Length to Which an Earthquake will go to Rupture Steven G. Wesnousky 1 and Glenn P. Biasi 2 1 Center of Neotectonic Studies and 2 Nevada Seismological Laboratory University of Nevada, Reno 89557 Abstract
More informationKinematics of the Southern California Fault System Constrained by GPS Measurements
Title Page Kinematics of the Southern California Fault System Constrained by GPS Measurements Brendan Meade and Bradford Hager Three basic questions Large historical earthquakes One basic question How
More informationSurface Rupture Effects on Earthquake Moment-Area Scaling Relations
Pure Appl. Geophys. 07 Springer International Publishing DOI 0.007/s0004-07-467-4 Pure and Applied Geophysics Surface Rupture Effects on Earthquake Moment-Area Scaling Relations YINGDI LUO, JEAN-PAUL AMPUERO,
More informationSimulation-based Seismic Hazard Analysis Using CyberShake
Simulation-based Seismic Hazard Analysis Using CyberShake SCEC CyberShake Collaboration: Robert Graves, Scott Callaghan, Feng Wang, Thomas H. Jordan, Philip Maechling, Kim Olsen, Kevin Milner, En-Jui Lee,
More informationPredicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L19305, doi:10.1029/2008gl035270, 2008 Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake D. A. Schmidt
More informationDIRECT HAZARD ANALYSIS OF INELASTIC RESPONSE SPECTRA
DIRECT HAZARD ANALYSIS OF INELASTIC RESPONSE SPECTRA ABSTRACT Y. Bozorgnia, M. Hachem, and K.W. Campbell Associate Director, PEER, University of California, Berkeley, California, USA Senior Associate,
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 informationArthur Frankel, William Stephenson, David Carver, Jack Odum, Robert Williams, and Susan Rhea U.S. Geological Survey
Probabilistic Seismic Hazard Maps for Seattle: 3D Sedimentary Basin Effects, Nonlinear Site Response, and Uncertainties from Random Velocity Variations Arthur Frankel, William Stephenson, David Carver,
More informationDCPP Seismic FAQ s Geosciences Department 08/04/2011 GM1) What magnitude earthquake is DCPP designed for?
GM1) What magnitude earthquake is DCPP designed for? The new design ground motions for DCPP were developed after the discovery of the Hosgri fault. In 1977, the largest magnitude of the Hosgri fault was
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 informationCHARACTERIZING EARTHQUAKE SLIP MODELS FOR THE PREDICTION OF STRONG GROUND MOTION
CHARACTERIZING EARTHQUAKE SLIP MODELS FOR THE PREDICTION OF STRONG GROUND MOTION P G SOMERVILLE 1, K IRIKURA, N ABRAHAMSON 3, S SAWADA 4, T KAGAWA 5 And Y TATSUMI 6 SUMMARY Over the past fifteen years,
More informationEarthquake stress drop estimates: What are they telling us?
Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,
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 informationRotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications
Bulletin of the Seismological Society of America, Vol. 85, No. 5, pp. 1513-1517, October 1995 Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications by
More informationEleventh U.S. National Conference on Earthquake Engineering Integrating Science, Engineering & Policy June 25-29, 2018 Los Angeles, California
Eleventh U.S. National Conference on Earthquake Engineering Integrating Science, Engineering & Policy June 25-29, 2018 Los Angeles, California Site-Specific MCE R Response Spectra for Los Angeles Region
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 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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationUsing information about wave amplitudes to learn about the earthquake size.
Earthquake Magnitudes and Moments Using information about wave amplitudes to learn about the earthquake size. Need to correct for decrease with distance M = log(a/t) + F(h,Δ) + C A is the amplitude of
More informationGraves and Pitarka Method
Based on Sept 26 meeting, is the output of the BBP consistent with the expectations from your method? Is there anything in the simulations that stand out as problematic? The results shown at the Sept 26
More informationJohn Shaw CRUSTAL DEFORMATION MODEL: OBLIQUE CONVERGENCE IN THE INNER CALIFORNIA BORDERLANDS ACCOMMODATED BY ACTIVE STRIKE-SLIP AND REVERSE FAULTS
CRUSTAL DEFORMATION MODEL: OBLIQUE CONVERGENCE IN THE INNER CALIFORNIA BORDERLANDS ACCOMMODATED BY ACTIVE STRIKE-SLIP AND REVERSE FAULTS August 2013 SONGS SSC SSHAC Workshop # 2 (August 12-14, 2013) Crustal
More informationHybrid Empirical Ground-Motion Prediction Equations for Eastern North America Using NGA Models and Updated Seismological Parameters
Hybrid Empirical Ground-Motion Prediction Equations for Eastern North America Using NGA Models and Updated Seismological Parameters by Shahram Pezeshk, 1 Arash Zandieh, 1 and Behrooz Tavakoli 2 1 Department
More informationHybrid Empirical Ground-Motion Prediction Equations for Eastern North America Using NGA Models and Updated Seismological Parameters
Bulletin of the Seismological Society of America, Vol. 101, No. 4, pp. 1859 1870, August 2011, doi: 10.1785/0120100144 Hybrid Empirical Ground-Motion Prediction Equations for Eastern North America Using
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 information4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration
Title 4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration Authors Bloomer, D., Ikon Science Asia Pacific Reynolds, S., Ikon Science Asia Pacific Pavlova, M., Origin
More informationRUPTURE MODELS AND GROUND MOTION FOR SHAKEOUT AND OTHER SOUTHERN SAN ANDREAS FAULT SCENARIOS
RUPTURE MODELS AND GROUND MOTION FOR SHAKEOUT AND OTHER SOUTHERN SAN ANDREAS FAULT SCENARIOS ABSTRACT : L.A. Dalguer 1, S.M. Day 2, K. Olsen 3, and V.M. Cruz-Atienza 4 1 Institute of Geophysics, ETH-Zurich,
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/4/3/eaao4915/dc1 Supplementary Materials for Global variations of large megathrust earthquake rupture characteristics This PDF file includes: Lingling Ye, Hiroo
More informationEarthquake Stress Drops in Southern California
Earthquake Stress Drops in Southern California Peter Shearer IGPP/SIO/U.C. San Diego September 11, 2009 Earthquake Research Institute Lots of data for big earthquakes (rupture dimensions, slip history,
More informationSynthetic Seismicity Models of Multiple Interacting Faults
Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).
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 informationComment on Systematic survey of high-resolution b-value imaging along Californian faults: inference on asperities.
Comment on Systematic survey of high-resolution b-value imaging along Californian faults: inference on asperities Yavor Kamer 1, 2 1 Swiss Seismological Service, ETH Zürich, Switzerland 2 Chair of Entrepreneurial
More informationSource parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop
Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for
More informationOn the validity of time-predictable model for earthquake generation in north-east India
Proc. Indian Acad. Sci. (Earth Planet. Sci.), Vol. 101, No. 4, December 1992, pp. 361-368. 9 Printed in India. On the validity of time-predictable model for earthquake generation in north-east India V
More informationEstimating Earthquake-Rupture Rates on a Fault or Fault System
Bulletin of the Seismological Society of America, Vol. 101, No. 1, pp. 79 92, February 2011, doi: 10.1785/0120100004 Estimating Earthquake-Rupture Rates on a Fault or Fault System by Edward H. Field and
More informationSan Francisco Bay Area Earthquake Simulations: A step toward a Standard Physical Earthquake Model
San Francisco Bay Area Earthquake Simulations: A step toward a Standard Physical Earthquake Model Steven N. Ward Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA,
More informationFault Scaling Relationships Depend on the
Manuscript Click here to download Manuscript RL+SR_2MW_r8.pdf Fault Scaling Relationships Depend on the Average Geological Slip Rate John G. Anderson, Glenn P. Biasi, Steven G. Wesnousky Abstract This
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 16, 9 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Energy Balance of dynamic rupture Crack tip stress field Frictional Rupture Fronts Meet in the lab (522
More informationI. Locations of Earthquakes. Announcements. Earthquakes Ch. 5. video Northridge, California earthquake, lecture on Chapter 5 Earthquakes!
51-100-21 Environmental Geology Summer 2006 Tuesday & Thursday 6-9:20 p.m. Dr. Beyer Earthquakes Ch. 5 I. Locations of Earthquakes II. Earthquake Processes III. Effects of Earthquakes IV. Earthquake Risk
More informationChallenges of Applying Ground Motion Simulation to Earthquake Engineering
Challenges of Applying Ground Motion Simulation to Earthquake Engineering Methodology of simulating ground motions from crustal earthquake and mega-thrust subduction earthquakes: application to the 2016
More informationThe UCERF3 Grand Inversion: Solving for the Long-Term Rate of Ruptures in a Fault System
Bulletin of the Seismological Society of America, Vol. 104, No. 3, pp. 1181 1204, June 2014, doi: 10.1785/0120130180 The UCERF3 Grand Inversion: Solving for the Long-Term Rate of Ruptures in a Fault System
More informationPotency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0
Geophys. J. Int. (2002) 148, F1 F5 FAST TRACK PAPER Potency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0 Yehuda Ben-Zion 1, and Lupei Zhu 2 1 Department of Earth
More informationMaterials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
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
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More informationREAL-TIME TSUNAMI INUNDATION FORECAST STUDY IN CHIMBOTE CITY, PERU
REAL-TIME TSUNAMI INUNDATION FORECAST STUDY IN CHIMBOTE CITY, PERU Nabilt Moggiano Supervisor: Kenji SATAKE MEE16720 ABSTRACT For rapid forecast of tsunami inundation during a tsunamigenic event, we constructed
More informationVerification of the asperity model using seismogenic fault materials Abstract
Verification of the asperity model using seismogenic fault materials Takehiro Hirose*, Wataru Tanikawa and Weiren Lin Kochi Institute for Core Sample Research/JAMSTEC, JAPAN * Corresponding author: hiroset@jamstec.go.jp
More informationCAMPBELL-BOZORGNIA NEXT GENERATION ATTENUATION (NGA) RELATIONS FOR PGA, PGV AND SPECTRAL ACCELERATION: A PROGRESS REPORT
Proceedings of the 8 th U.S. National Conference on Earthquake Engineering April 18-22, 2006, San Francisco, California, USA Paper No. 906 CAMPBELL-BOZORGNIA NEXT GENERATION ATTENUATION (NGA) RELATIONS
More informationModifications to Existing Ground-Motion Prediction Equations in Light of New Data
Bulletin of the Seismological Society of America, Vol. 101, No. 3, pp. 1121 1135, June 2011, doi: 10.1785/0120100270 Modifications to Existing Ground-Motion Prediction Equations in Light of New Data by
More informationA GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR
A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR Annemarie CHRISTOPHERSEN 1 And Euan G C SMITH 2 SUMMARY This paper considers the distribution of aftershocks in space, abundance, magnitude and time. Investigations
More informationM 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2006jb004415, 2006 M 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model Tom Parsons 1 Received
More informationACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT
ACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT Edward H FIELD 1 And SCEC PHASE III WORKING GROUP 2 SUMMARY Probabilistic seismic hazard analysis
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 informationApplication of Phase Matched Filtering on Surface Waves for Regional Moment Tensor Analysis Andrea Chiang a and G. Eli Baker b
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 Application of Phase Matched Filtering on Surface Waves for Regional Moment Tensor Analysis Andrea Chiang a and G. Eli
More informationTwo ways to think about the dynamics of earthquake ruptures
Two ways to think about the dynamics of earthquake ruptures (1) In terms of friction (2) In terms of fracture mechanics Scholz describes conditions for rupture propagation (i.e. instability) via energy
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 informationby George E. Hilley, J Ramón Arrowsmith, and Elizabeth Stone Introduction
Bulletin of the Seismological Society of America, 91, 3, pp. 427 440, June 2001 Inferring Segment Strength Contrasts and Boundaries along Low-Friction Faults Using Surface Offset Data, with an Example
More informationFocused Observation of the San Andreas/Calaveras Fault intersection in the region of San Juan Bautista, California
Focused Observation of the San Andreas/Calaveras Fault intersection in the region of San Juan Bautista, California Clifford Thurber, Dept. of Geology and Geophysics, Univ. Wisconsin-Madison Mick Gladwin,
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 informationJournal of Geophysical Research (Solid Earth) Supporting Information for
Journal of Geophysical Research (Solid Earth) Supporting Information for Postseismic Relocking of the Subduction Megathrust Following the 2007 Pisco, Peru earthquake D.Remy (a), H.Perfettini (b), N.Cotte
More informationRate and State-Dependent Friction in Earthquake Simulation
Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity
More informationLecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller
Lecture 20: Slow Slip Events and Stress Transfer GEOS 655 Tectonic Geodesy Jeff Freymueller Slow Slip Events From Kristine Larson What is a Slow Slip Event? Slip on a fault, like in an earthquake, BUT
More informationPower-law distribution of fault slip-rates in southern California
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L23307, doi:10.1029/2007gl031454, 2007 Power-law distribution of fault slip-rates in southern California Brendan J. Meade 1 Received 31
More informationRelevance Vector Machines for Earthquake Response Spectra
2012 2011 American American Transactions Transactions on on Engineering Engineering & Applied Applied Sciences Sciences. American Transactions on Engineering & Applied Sciences http://tuengr.com/ateas
More informationCOULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES
COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES Madlazim Physics Department, Faculty Mathematics and Sciences of Surabaya State University (UNESA) Jl. Ketintang, Surabaya 60231, Indonesia. e-mail:
More informationRupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from : 1. Source Parameter Scaling Relationships
Journal of Geophysical Research Solid Earth Supporting Information for Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from 1990-2015: 1. Source Parameter Scaling Relationships
More informationRapid Earthquake Rupture Duration Estimates from Teleseismic Energy Rates, with
1 2 Rapid Earthquake Rupture Duration Estimates from Teleseismic Energy Rates, with Application to Real-Time Warning 3 Jaime Andres Convers 1 and Andrew V. Newman 1 4 5 1. School of Earth and Atmospheric
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 informationCoulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment
Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Abstract D. Weatherley University of Queensland Coulomb stress change analysis has been applied in
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi: 10.1038/ngeo739 Supplementary Information to variability and distributed deformation in the Marmara Sea fault system Tobias Hergert 1 and Oliver Heidbach 1,* 1 Geophysical
More informationThe 3 rd SCEC CSM workshop
The 3 rd SCEC CSM workshop Welcome on behalf of the organizers Jeanne Hardebeck Brad Aagaard David Sandwell Bruce Shaw John Shaw Thorsten Becker Thanks for playing! SCEC Community Stress Model (CSM) Community
More informationScaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence
Earth Planets Space, 53, 943 948, 2001 Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence Z. L. Wu Institute of Geophysics,
More informationThe ISC-GEM Global Instrumental Reference Earthquake Catalogue ( )
The ISC-GEM Global Instrumental Reference Earthquake Catalogue (1900-2009) Dmitry Storchak, Domenico Di Giacomo, István Bondár, James Harris International Seismological Centre (ISC), Thatcham, United Kingdom
More information7 Ground Motion Models
7 Ground Motion Models 7.1 Introduction Ground motion equations are often called attenution relations but they describe much more than just the attenutation of the ground motion; they describe the probability
More informationThe Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) Project Plan
UCERF3_Project_Plan_v55.doc 1 The Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) Project Plan by the Working Group on California Earthquake Probabilities (WGCEP) Notes on this version
More informationSONGS SSC. Tom Freeman GeoPentech PRELIMINARY RESULTS
SONGS SSC Tom Freeman GeoPentech PRELIMINARY RESULTS Focused Questions Summarize the tectonic setting What is the seismogenic thickness? Are you including deep ruptures in the upper mantle (~30 km)? Do
More informationExternal Grant Award Number 04HQGR0058 IMPROVED THREE-DIMENSIONAL VELOCITY MODELS AND EARTHQUAKE LOCATIONS FOR CALIFORNIA
External Grant Award Number 04HQGR0058 IMPROVED THREE-DIMENSIONAL VELOCITY MODELS AND EARTHQUAKE LOCATIONS FOR CALIFORNIA Clifford H. Thurber University of Wisconsin-Madison 1215 W. Dayton St. Madison,
More informationKnowledge of in-slab earthquakes needed to improve seismic hazard estimates for southwestern British Columbia
USGS OPEN FILE REPORT #: Intraslab Earthquakes 1 Knowledge of in-slab earthquakes needed to improve seismic hazard estimates for southwestern British Columbia John Adams and Stephen Halchuk Geological
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 informationSource Duration Scales with Magnitude Differently For Earthquakes on the San Andreas Fault and on Secondary Faults in Parkfield, CA
Source Duration Scales with Magnitude Differently For Earthquakes on the San Andreas Fault and on Secondary Faults in Parkfield, CA Rebecca M. Harrington University of California, Los Angeles Dept. of
More information