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

ecorded within 70 km of earthquake ruptures, the region of primary concern for application of the model. We also believe that the model properly represents the attenuation of ground motion at larger distances in California. The issue we have raised points out the need to systematically collect and process all recordings from large earthquakes if there is interest in correctly modeling the attenuation of ground motions at large distances. Equation (19) provides a relationship for γ for pga. The relationships for other spectral periods was constructed by using the period-dependence of γ computed from analysis of the broadband data from the three small southern California earthquakes and scaling the values by the relative difference in γ at pga, that is: γ ( pga) California γ ( T ) California = γ ( T ) Anza− YL−BBC × (20) γ ( pga) Anza−YL−BBC Figure 30 compares the model developed for γ with values obtained from direct regression of the PEER-NGA data without imposing any distance truncation and without using truncated regression. The variation of γ with spectral period T is similar. Figure 30: Comparison of the model developed for γ (Equations 19 and 20) versus values of γ obtained from regression without truncation (black points) of the PEER-NGA data set. MODELING STEPS The parameters of the ground motion model were developed through an iterative process that involved performing regressions for the entire spectral period range with some parts of the model fixed, examining the trend of remaining parameters with spectral period, developing smoothing models for these parameters with period, and then repeating the analysis to examine the variation of the remaining parameters. The major steps in the development of the parameters are summarized below. All analyses were performed using the PEER-NGA C&Y2006 Page 39

data truncated at a maximum distance of RRUP = 70 km. In addition, the magnitude scaling parameters c2 and c3 were fixed at 1.06 and 3.45 for all analyses and in the final model. As discussed in the model development section, c3 = 3.45 represents scaling of low-frequency amplitudes proportional to seismic moment. Parameter c2 was fixed at 1.06, the value derived from the Atkinson and Silva (2006) seismic source model. All previous applications of the general function form used by Ross Sadigh have yielded far-source high-frequency magnitude scaling coefficients that are consistent with this value. Preliminary analysis of the PEER-NGA data also produced values near 1.06. Parameter c4 was examined by analyses of PEER-NGA data and TriNet data. It was concluded that a range of values for c4 would provide satisfactory fits to the data with adjustments to the rate of attenuation occurring through changes of parameters c5 and c6. The change in attenuation rate with distance was set to occur at cRB = 50 km. A value near 50 km was consistently found in all preliminary analyses. After extensive exploratory analysis of the data, we concluded that the PEER-NGA dataset does not sample a sufficiently wide range of motion and Vs30 to allow the simultaneous estimation of all 3 nonlinear soil parameters (φ2, φ3,, and φ4) We thus decided to fix φ4 through other supporting information. Both our preliminary analysis of the PEER-NGA data and our analysis of the simulated soil amplification factors from Silva (2004) suggested that φ4 is about 0.1g in the case of pga. The empirical study of Choi and Stewart (2005) also indicate 0.1g is a viable value for φ4. We thus fix pga’s φ4 at 0.1g. To fix φ4 of other periods, we anchor the rock spectral shape of M6.5 at 10 km (Silva et a. Add referrecne) at 0.1g pga. The remaining two nonlinear parameters will adjust to the fixed value of φ4 and they are more than sufficient to capture the nonlinear soil behavior existing in the PEER-NGA dataset. In the preliminary analyses it was found that the data exhibited a statistically significant dependence on source depth parameterized as the depth to top of rupture, ZTOR and that aftershocks showed a stronger dependence on depths than main shocks. Also, the data indicated that the style of faulting effects were weaker for aftershocks than for main shocks. Therefore, the analyses were performed with separate depth dependence and style of faulting effects for main shocks and aftershocks. First Phase: Step 1: The first phase of the analysis used only data from neutral sites. Sites with the hanging wall/foot wall flags set to ‘fw’ or ‘hw’ in the PEER-NGA database were removed. This simplification was performed to examine other portions of the model without interaction with the hanging wall function fHW. In particular, regression with the full model showed high correlation between estimates of the hanging wall parameters and the parameters for style of faulting and soil nonlinearity. Using a simple model for the reference motion (without style of faulting, depth of rupture, and variable magnitude scaling) the parameters of the non-linear soil model (φ2 and φ3) were obtained from fits to the NGA data. At longer periods when the soil response becomes nearly linear, estimates of the nonlinear soil parameter φ3 fluctuated wildly and had very large errors of estimation. Therefore parameter φ3 for periods longer than 1.0 seconds were fixed to produce a smooth transition to linear behavior at spectral periods near 1.0 sec and greater. C&Y2006 Page 40

data truncated at a maximum distance of RRUP = 70 km. In addition, the magnitude scaling<br />

parameters c2 <strong>and</strong> c3 were fixed at 1.06 <strong>and</strong> 3.45 <strong>for</strong> all analyses <strong>and</strong> in the final model. As<br />

discussed in the model development section, c3 = 3.45 represents scaling of low-frequency<br />

amplitudes proportional to seismic moment. Parameter c2 was fixed at 1.06, the value<br />

derived from the Atkinson <strong>and</strong> Silva (2006) seismic source model. All previous applications<br />

of the general function <strong>for</strong>m used by Ross Sadigh have yielded far-source high-frequency<br />

magnitude scaling coefficients that are consistent with this value. Preliminary analysis of the<br />

<strong>PEER</strong>-<strong>NGA</strong> data also produced values near 1.06. Parameter c4 was examined by analyses of<br />

<strong>PEER</strong>-<strong>NGA</strong> data <strong>and</strong> TriNet data. It was concluded that a range of values <strong>for</strong> c4 would<br />

provide satisfactory fits to the data with adjustments to the rate of attenuation occurring<br />

through changes of parameters c5 <strong>and</strong> c6. The change in attenuation rate with distance was<br />

set to occur at cRB = 50 km. A value near 50 km was consistently found in all preliminary<br />

analyses.<br />

After extensive exploratory analysis of the data, we concluded that the <strong>PEER</strong>-<strong>NGA</strong> dataset<br />

does not sample a sufficiently wide range of motion <strong>and</strong> Vs30 to allow the simultaneous<br />

estimation of all 3 nonlinear soil parameters (φ2, φ3,, <strong>and</strong> φ4) We thus decided to fix φ4<br />

through other supporting in<strong>for</strong>mation. Both our preliminary analysis of the <strong>PEER</strong>-<strong>NGA</strong> data<br />

<strong>and</strong> our analysis of the simulated soil amplification factors from Silva (2004) suggested that<br />

φ4 is about 0.1g in the case of pga. The empirical study of Choi <strong>and</strong> Stewart (2005) also<br />

indicate 0.1g is a viable value <strong>for</strong> φ4. We thus fix pga’s φ4 at 0.1g. To fix φ4 of other periods,<br />

we anchor the rock spectral shape of M6.5 at 10 km (Silva et a. Add referrecne) at 0.1g pga.<br />

The remaining two nonlinear parameters will adjust to the fixed value of φ4 <strong>and</strong> they are<br />

more than sufficient to capture the nonlinear soil behavior existing in the <strong>PEER</strong>-<strong>NGA</strong><br />

dataset.<br />

In the preliminary analyses it was found that the data exhibited a statistically significant<br />

dependence on source depth parameterized as the depth to top of rupture, ZTOR <strong>and</strong> that<br />

aftershocks showed a stronger dependence on depths than main shocks. Also, the data<br />

indicated that the style of faulting effects were weaker <strong>for</strong> aftershocks than <strong>for</strong> main shocks.<br />

There<strong>for</strong>e, the analyses were per<strong>for</strong>med with separate depth dependence <strong>and</strong> style of faulting<br />

effects <strong>for</strong> main shocks <strong>and</strong> aftershocks.<br />

First Phase:<br />

Step 1: The first phase of the analysis used only data from neutral sites. Sites with the<br />

hanging wall/foot wall flags set to ‘fw’ or ‘hw’ in the <strong>PEER</strong>-<strong>NGA</strong> database were removed.<br />

This simplification was per<strong>for</strong>med to examine other portions of the model without interaction<br />

with the hanging wall function fHW. In particular, regression with the full model showed high<br />

correlation between estimates of the hanging wall parameters <strong>and</strong> the parameters <strong>for</strong> style of<br />

faulting <strong>and</strong> soil nonlinearity. Using a simple model <strong>for</strong> the reference motion (without style<br />

of faulting, depth of rupture, <strong>and</strong> variable magnitude scaling) the parameters of the non-linear<br />

soil model (φ2 <strong>and</strong> φ3) were obtained from fits to the <strong>NGA</strong> data. At longer periods when the<br />

soil response becomes nearly linear, estimates of the nonlinear soil parameter φ3 fluctuated<br />

wildly <strong>and</strong> had very large errors of estimation. There<strong>for</strong>e parameter φ3 <strong>for</strong> periods longer<br />

than 1.0 seconds were fixed to produce a smooth transition to linear behavior at spectral<br />

periods near 1.0 sec <strong>and</strong> greater.<br />

C&Y2006 Page 40

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

Saved successfully!

Ooh no, something went wrong!