Chiou and Youngs PEER-NGA Empirical Ground Motion Model for ...

Chiou and Youngs PEER-NGA Empirical Ground Motion Model for ... Chiou and Youngs PEER-NGA Empirical Ground Motion Model for ...

peer.berkeley.edu
from peer.berkeley.edu More from this publisher
08.04.2013 Views

MODEL FORMULATION All ground motion models, empirical or theoretical/numerical, attempt to capture three basic aspects of ground motion estimation, the effect of the amount of energy radiated at the source, the effect of attenuation of seismic waves along the path due to geometric spreading and energy absorption, and the effect of local modification of the seismic waves by the nearsurface materials. The empirical ground motion model developed in this study attempts to capture these elements using (relatively) simple algebraic expressions to represent the average behavior observed in the empirical strong motion data. The form of these expressions was guided by trends in the data, simple seismological models, past experience, and examination of the results of 1-D and 3-D ground motion simulations conducted as part of the PEER-NGA project. Seismic Source Scaling Effect of Earthquake Size: The simplest and most commonly used measure for correlating the amount of energy released in an earthquake with the resulting amplitudes of ground motions is the earthquake magnitude, defined in the PEER-NGA database in terms of moment magnitude M. Most empirical ground motion models have used a polynomial function for scaling the log of ground motions y with magnitude of the form n ln( y ) ∝ C2 ( T ) × M + C3 ( T ) × ( mC − M) (1) where coefficients C3 and possibly C2 vary with spectral period T, the exponent n is typically in the range of 2 to 3, and the value of mC is independent of T. We have introduced an alternative form for magnitude scaling that reflects the results obtained using seismological models for the form of the earthquake source spectrum. Figure 7 shows response spectra computed using the stochastic ground motion model (e.g. Boore 2003) with two source spectra models, a Brune (1970, 1971) source with a stress parameter of 70 bars and the Atkinson and Silva (1997) empirical source spectra model for California earthquakes. Both forms of source spectra display similar trends in magnitude scaling. For high frequency motions (> ~ 10 Hz), the magnitude scaling of ln(y) is approximately linear at a relatively flat slope in the magnitude range of 5 ≤ M ≤ 8.5. The slope of the curves reflects the scaling of the source spectrum above the corner frequency. For very low frequency motions (< ~0.3 Hz) the magnitude scaling in the range of 5 ≤ M ≤ 7 is again approximately linear, but at a much larger slope than for high frequency motions, reflecting scaling of the source spectrum below the corner frequency. At intermediate spectral frequencies, there is a transition in the scaling of ln(y) with magnitude from low-frequency scaling at low magnitudes to high-frequency scaling at large magnitudes. This transition occurs over the magnitude range where the corner frequency of the source spectrum is near the spectral frequency of the ground motion. As the spectral frequency of the ground motions decreases, the magnitude range for this transition shifts to larger magnitudes, reflecting the decrease in corner frequency with increasing magnitude. The shape of the magnitude scaling curves shown in Figure 7 is modeled by the expression: C&Y2006 Page 13

1 ) ∝ C2 × M + ( C2 − C ) × ln M T − C ( T ) ln( y 3 n n [ 1+ exp{ C ( T ) × ( C ( ) M) } ] In Equation (2) parameter C2 is the slope of the magnitude scaling relationship above the corner frequency, C3 is the slope of the magnitude scaling relationship below the corner frequency. Parameter Cn controls the magnitude range over which the transition from C2 scaling to C3 scaling occurs. Parameter CM is the magnitude at the midpoint of this transition and its value varies with the spectral period of the ground motion parameter y. The function form of Equation (2) was tested by fitting a set of ground motion response spectra simulated using the program SMSIM (Boore, 2005) and the Atkinson and Silva (2000) stochastic ground motion. Figure 8 shows the magnitude scaling obtained from the simulated ground motions. As indicated on the figure, the functional form of Equation (2) can provide a close fit to the simulations. Figure 9 shows peak acceleration data in the distance range of 30 to 50 km and the VS30 range of 300 to 400 m/s from the PEER-NGA database and from the TriNet data (Appendix D). These values were fit using three alternative scaling relationships. The quadratic form of Equation (1) uses an exponent n of 2. The Sadigh et al. (1997) form of Equation (1) uses n equal to 2.5 and mC equal to 8.5. In the magnitude scaling defined by Equation (2), parameter C3 represents scaling of ground motion spectra below the corner frequency where the source spectra are directly proportional to seismic moment, M0. Because M ∝ ⅔M0, theoretically, C3 should equal 1.5×ln(10) or 3.45 and the fit of Equation (2) to the data shown on Figure 9 used this value of C3. Parameter C2 represents the scaling of the source spectra above the corner frequency. The source spectral model defined by Atkinson and Silva (2000) produces scaling of ln(y) ∝ 1.06M. As described in the next section, in the Sadigh et al. (1997) form of Equation (1), parameter C2 also represents the magnitude scaling of the source spectra at distances unaffected by extended-source effects and in the Sadigh et al. (1997) values of C2 were 1.0 to 1.1. The fit of the Sadigh et al. (1997) form of Equation (1) to the data on Figure 9 also produced C2 equal to 1.0. In fitting Equation (2) to the data on Figure 9, C2 was constrained to be 1.06, the value obtained from the magnitude scaling of source spectra above the corner frequency defined by Atkinson and Silva (2000). All three forms provide essentially equally good fits to the data. We prefer the scaling form of Equation (2) because we believe that it better represents our current concept of the variation in earthquake source spectra with earthquake magnitude. C&Y2006 Page 14 (2)

MODEL FORMULATION<br />

All ground motion models, empirical or theoretical/numerical, attempt to capture three basic<br />

aspects of ground motion estimation, the effect of the amount of energy radiated at the<br />

source, the effect of attenuation of seismic waves along the path due to geometric spreading<br />

<strong>and</strong> energy absorption, <strong>and</strong> the effect of local modification of the seismic waves by the nearsurface<br />

materials. The empirical ground motion model developed in this study attempts to<br />

capture these elements using (relatively) simple algebraic expressions to represent the<br />

average behavior observed in the empirical strong motion data. The <strong>for</strong>m of these<br />

expressions was guided by trends in the data, simple seismological models, past experience,<br />

<strong>and</strong> examination of the results of 1-D <strong>and</strong> 3-D ground motion simulations conducted as part<br />

of the <strong>PEER</strong>-<strong>NGA</strong> project.<br />

Seismic Source Scaling<br />

Effect of Earthquake Size: The simplest <strong>and</strong> most commonly used measure <strong>for</strong> correlating<br />

the amount of energy released in an earthquake with the resulting amplitudes of ground<br />

motions is the earthquake magnitude, defined in the <strong>PEER</strong>-<strong>NGA</strong> database in terms of<br />

moment magnitude M. Most empirical ground motion models have used a polynomial<br />

function <strong>for</strong> scaling the log of ground motions y with magnitude of the <strong>for</strong>m<br />

n<br />

ln( y ) ∝ C2<br />

( T ) × M + C3<br />

( T ) × ( mC<br />

− M)<br />

(1)<br />

where coefficients C3 <strong>and</strong> possibly C2 vary with spectral period T, the exponent n is typically<br />

in the range of 2 to 3, <strong>and</strong> the value of mC is independent of T.<br />

We have introduced an alternative <strong>for</strong>m <strong>for</strong> magnitude scaling that reflects the results<br />

obtained using seismological models <strong>for</strong> the <strong>for</strong>m of the earthquake source spectrum. Figure<br />

7 shows response spectra computed using the stochastic ground motion model (e.g. Boore<br />

2003) with two source spectra models, a Brune (1970, 1971) source with a stress parameter<br />

of 70 bars <strong>and</strong> the Atkinson <strong>and</strong> Silva (1997) empirical source spectra model <strong>for</strong> Cali<strong>for</strong>nia<br />

earthquakes. Both <strong>for</strong>ms of source spectra display similar trends in magnitude scaling. For<br />

high frequency motions (> ~ 10 Hz), the magnitude scaling of ln(y) is approximately linear at<br />

a relatively flat slope in the magnitude range of 5 ≤ M ≤ 8.5. The slope of the curves reflects<br />

the scaling of the source spectrum above the corner frequency. For very low frequency<br />

motions (< ~0.3 Hz) the magnitude scaling in the range of 5 ≤ M ≤ 7 is again approximately<br />

linear, but at a much larger slope than <strong>for</strong> high frequency motions, reflecting scaling of the<br />

source spectrum below the corner frequency.<br />

At intermediate spectral frequencies, there is a transition in the scaling of ln(y) with<br />

magnitude from low-frequency scaling at low magnitudes to high-frequency scaling at large<br />

magnitudes. This transition occurs over the magnitude range where the corner frequency of<br />

the source spectrum is near the spectral frequency of the ground motion. As the spectral<br />

frequency of the ground motions decreases, the magnitude range <strong>for</strong> this transition shifts to<br />

larger magnitudes, reflecting the decrease in corner frequency with increasing magnitude.<br />

The shape of the magnitude scaling curves shown in Figure 7 is modeled by the expression:<br />

C&Y2006 Page 13

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!