Bulletin of the Seismological Society of America, Vol. 94, No. 6, pp. 2110–2123, December 2004
Nonlinear Soil-Site Effects in Probabilistic Seismic-Hazard Analysis by Paolo Bazzurro* and C. Allin Cornell
Abstract This study presents effective probabilistic procedures for evaluating ground-motion hazard at the free-field surface of a nonlinear soil deposit located at a specific site. Ground motion at the surface, or at any depth of interest within the soil formation (e.g., at the structure foundation level), is defined here in terms either of a suite of oscillator-frequency-dependent hazard curves for spectral acceleration, s Sa(f), or of one or more spectral acceleration uniform-hazard spectra, each associated with a given mean return period. It is presumed that similar information is available for the rock-outcrop input. The effects of uncertainty in soil properties are directly included. This methodology incorporates the amplification of the local soil deposit into the framework of probabilistic seismic hazard analysis (PSHA). The soil amplification is characterized by a frequency-dependent amplification function, AF(f), where f is a s generic oscillator frequency. AF(f) is defined as the ratio ofSa(f) to the spectral s acceleration at the bedrock level,Sa(f) . The estimates of the statistics of the ampli- fication function are obtained by a limited number of nonlinear dynamic analyses of the soil column with uncertain properties, as discussed in a companion article in this issue (Bazzurro and Cornell, 2004). The hazard at the soil surface (or at any desired depth) is computed by convolving the site-specific hazard curve at the bedrock level with the probability distribution of the amplification function. The approach presented here provides more precise surface ground-motion-hazard estimates than those found by means of standard attenuation laws for generic soil conditions. The use of generic ground-motion predictive equations may in fact lead to inaccurate results especially for soft-clay-soil sites, where considerable amplifi- cation is expected at long periods, and for saturated sandy sites, where high-intensity ground shaking may cause loss of shear strength owing to liquefaction or to cyclic mobility. Both such cases are considered in this article. In addition to the proposed procedure, two alternative, easier-to-implement but approximate techniques for obtaining hazard estimates at the soil surface are also briefly discussed. One is based on running a conventional PSHA with a rock- attenuation relationship modified to include the soil response, whereas the other consists of using a simple, analytical, closed-form solution that appropriately mod- ifies the hazard results at the rock level.
Introduction
In probabilistic seismic-hazard analysis (PSHA) the ef- many references in Seismological Research Letters, 1997). fects of local soil deposits on the seismic hazard at the sur- This approach ignores virtually all site-specific information face are often treated with less rigor than this critical aspect and, therefore, produces only a broad, generic assessment of deserves. the hazard. A simple and widely used approach for noncritical fa- Because of soil nonlinearity in the soil response for se- cilities assumes that the soil conditions at the site resemble vere bedrock motions, the use of a generic soil-attenuation those at the stations in the database considered for the de- relation may, however, yield inaccurate results even when velopment of one of the many soil-site ground-motion pre- the assumption of generic soil conditions is appropriate. dictive equations available in the literature (e.g., see the Most empirical predictive equations do not explicitly incor- porate the effects of site response beyond a simple term that *Present address: AIR Worldwide, San Francisco, California. includes a macro-parameter of the soil deposit. For example,
2110 Nonlinear Soil-Site Effects in Probabilistic Seismic-Hazard Analysis 2111
Campbell (1997) uses the depth to basement rock, whereas relation to attain a site-specific soil-surface predictive equa- Boore et al. (1997) adopt the average shear-wave velocity tion to be used in PSHA; the second modifies the hazard at in the first 30 m of soil, ¢V30. These two models predict the the rock level by means of a generalization of the so-called same amplification for mild and severe bedrock shaking; that risk equation (e.g., Cornell, 1994, 1996). is, they do not include any nonlinearity in the site response. To our knowledge, the attenuation equation by Abrahamson Methodology and Silva (1997) is the only one published to date that ex- plicitly allows for the nonlinear response of local soil. To The methodology that follows is a fully probabilistic accommodate soil nonlinearity, however, this equation uses procedure to account for nonlinear soil response in PSHA. Preliminary versions were presented in Cornell and Bazzurro the expected peak bedrock ground acceleration,PGArock , which in our companion article in this issue (Bazzurro and (unpublished manuscript, 1997), Bazzurro (1998), and Baz- Cornell, 2004) has been shown to be a rather poor predictor zurro et al. (1998). This body of work has also inspired, of soil amplification at low frequencies. directly or indirectly, similar studies on this subject by Lee Another method (used almost exclusively for important et al. (1998, 1999), Lee (2000), and, more recently, Cramer facilities) consists of explicitly characterizing the soil re- (2003). Tsai (2000) also recently investigated a similar issue. sponse via a deterministic amplification function (usually the The methodology presented by the authors in their early mean or the median) computed for a generally small suite work is an integral part of the procedure prepared for the of appropriate seismograms driven through a computer U.S. Nuclear Regulatory Commission (USNRC) for nuclear- model of the soil column with best estimates of the soil pa- facility sites (e.g., see section 6 and appendix I in USNRC, rameters. The equivalent-linear SHAKE program (Schnabel 2001; and applications in USNRC, 2002) for developing haz- et al., 1972) is used for such purposes in the majority of ard-consistent spectra on soil. cases. The seismograms are either real recordings from The objective is to estimate the soil-surface-hazard s events consistent with those dominating the site hazard or curve for the spectral acceleration,Sa (f, n), at a generic os- synthetic time histories that match a predefined target- cillator frequency, f, and damping, n. This is simply a rela- s response spectrum. In most such applications the desired tionship betweenSa (f, n) and the annual mean rate of ex- hazard at the soil surface is obtained by multiplying the bed- ceedance (MRE) or alternatively versus the mean return rock hazard (in the form either of uniform hazard spectra or period (MRP). The MRP is defined as the reciprocal of the 4מ implies an MRP of 10 ן of hazard curves) by this deterministic amplification func- annual MRE (e.g., an MRE of 4 tion. This approach, which is the hybrid method in Cramer 2500 years). For small values of engineering interest the (2003) and is similar to approach 1 in U.S. Nuclear Regu- rare-event assumption holds. This assumption implies that latory Commission (2001), does not explicitly account for the likelihood of two or more events occurring in the period the amplification function record-to-record variability. It of interest is small in comparison to the likelihood that one produces surface ground-motion levels whose exceedance event happens. The rare-event assumption implies that the rates are unknown, non-uniform, inconsistent across fre- MRE of a value s is numerically equal to the probability quency, and generally nonconservative. Similarly, soil sur- that the annual maximum spectral acceleration exceeds s. face ground motions with unknown exceedance rates are in Throughout this article we make this assumption, and there- general obtained if one multiplies the bedrock hazard by the fore we are permitted to use the notation of GS(s) for this average National Earthquake Hazards Reduction Program complementary cumulative distribution function (CCDF), (NEHRP) amplification factors (e.g., see Tables 4.1.2a and understanding that it is numerically identical to the hazard 4.1.2b in Federal Emergency Management Agency, 2001). or MRE curve. In more refined studies (e.g., Cramer, 2003) the uncer- The uniform hazard spectrum (UHS) at the soil surface tainty of the site response is indeed considered but often is for any desired MRP can be obtained by interpolating hazard not properly coupled with the ground-motion variability in curves, if available, for a suitably dense range of frequencies. the PSHA computations that lead to the hazard estimates at Throughout this study the damping is assumed equal to 5% the bedrock. The result of this operation, as pointed out in of critical, and therefore n is dropped hereafter from the Silva et al. (2000), is an artificially inflated hazard at the soil notation. Note that the following procedure can be applied surface. This subject will be clarified below. to obtain the hazard at any depth (e.g., at the structure- The approach presented here, which is heavily based on foundation level), provided that the soil response is com- the findings presented in our companion article, overcomes puted accordingly. the shortcomings previously outlined. Applications are shown Consistently with our companion article (Bazzurro and for two very different offshore soil sites, one sandy and one Cornell, 2004), the effect of the soil layers on the intensity clayey. For illustration purposes, the hazard results from this of the ground motion at the surface is incorporated via a site- method are compared with those obtained by using an atten- specific, frequency-dependent amplification function, AF(f): uation equation for generic soil. We also discuss the use of s (Sa(f ס two alternative, efficient but approximate techniques: the AF(f) r , (1) first requires modifying a rock ground-motion-attenuation Sa(f) 2112 P. Bazzurro and C. A. Cornell
sr s where f is a generic oscillator frequency,Saa(f) andS (f) are for example, GZ(z) is the sought hazard curve forSa(f) , that the 5%-damped spectral-acceleration values at the soil sur- is, the annual probability of exceeding level z, and GY|X is -AF(f), conditional on a rock-level ampli ס face and at the bedrock, respectively. the CCDF of Y
The prediction of AF(f), given the knowledge of bed- tude x—and fX(x) is the probability density function (pdf) of r rock-ground-motion parameters, was studied in our compan- Sa(f). Note that fX(x) is the absolute value of the derivative ion article by running multiple recordings through different of GX(x) with respect to x, and recall that GX(x) is numeri- realizations of the soil column with uncertain properties. cally equal to MRE and fX(x) is referred to as the mean rate r That study has shown thatSa(f) is the most effective pre- density. Because in practice GZ(z) is always numerically dictor variable for estimating AF(f) (at the same frequency computed, this equation is more readily useful in discretized f) among different bedrock-ground-motion parameters such form: as magnitude, M, of the causative event; source-to-site dis- tance, R; PGA ; and spectral-acceleration values,Sr(f ) , at סՆ zz ס r a sc the initial resonant frequency, f , of the soil column. GZ(z) ͚ PY΄Η΅xpjXj(x ) ͚ GxpY|XjXj Η(x ). sc all x xxall x r jj Furthermore, results showed that once theSa(f) value (3) of a record at the bedrock is known, the additional knowl- edge of M and R, which implicitly define its average The term pX(xj) represents the probability that the rock- response spectrum shape, do not appreciably improve the input level is equal to (or better, in the neighborhood of) xj. estimation of AF(f) at the same frequency f. In other words, r This term can be approximately derived by differentiating AF(f), conditioned onSa(f) , is virtually independent on M the rock-hazard curve in discrete or numerical form. An easy and R. The interested reader is referred to our companion r way of interpreting this equation is considering that if the article for a quantification of such a statement.Sa(f) can be random variable X takes on some value x (say, 0.25 g) then said to be, therefore, a sufficient estimator of AF(f) (Luco, Z will exceed the value z (say, 0.5 g) if and only if Y takes We consider the probability of .2 ס The procedure for estimating surface hazard in this on some value y Ն z/x .(2001 case is given in the subsection that follows. Ն this event by the product P[Y z/x|xj]pX(xj), and then we Strictly speaking, this lack of M and R conditional de- sum over all possible values of xj. pendence does not hold for frequencies f below fsc, where a Y X r Assuming the lognormality of given , when all the mild dependence on M—or, more significantly, onSa(fsc) — RVs are expressed in continuous rather than discrete form was found. We will treat this latter, more complicated case the GY|X is given by: in a subsequent section. Finally, if a particular application AF f M R warrants keeping the dependence of ( )on and , be- z [(ln[mˆ (x מ΅΄sides onSrs(f) , then the procedure to estimateS (f) versus ln aa x Y|X ˆ ס MRE requires knowledge of the hazard curve at bedrock plus z GxY|XΗ U (4) results that can be extracted from standard hazard disaggre- x rlnY|X gation. The interested reader can find relevant details in Baz- -U(•) is the widely tabulated comple מ 1 ס (•) ˆzurro (1998). in which U mentary standard Gaussian CDF. Estimates of the distribu- Convolution: AF(f) Dependent on Sr(f) tion parameters of Y (i.e., the conditional median of Y, mˆ Y|X, a and the conditional standard deviation of the natural log- The proposed method for computing surface-hazard arithm of Y, r ) can be found by driving a suite of n rock- s lnY|X ס curves for Z Sa(f) convolves the site-specific rock-hazard ground-motion records through a sample of soil-column rep- r ס curves for X Sa(f), which may be exogenously provided resentations (recall that the soil properties are uncertain) and (e.g., by a site-specific PSHA analysis or by the hazard maps then regressing, for each frequency f, the values of ln Y on at the U.S. Geological Survey’s Web site (http://geohazards ln X. These regressions, based on a large suite of records, ס .cr.usgs.gov/eq/), with the probability distribution of Y were shown for two soil sites in figure 8 of our companion AF(f), an estimate of which is obtained through nonlinear article. We emphasize again, however, that in real applica- dynamic analyses of the soil. The change in notation is tions, given the relative small values of rlnY|X, only 10 soil adopted to make the following equations less cumbersome. response analyses are sufficient to achieve a level of accu- • ס in estimating the median AF(f). The large 10%ע Note that Z X Y. Based on the elementary probability racy of theory (e.g., Benjamin and Cornell, 1970, pp. 112–124), the database was used only to validate the procedure. convolution goes as follows: From the computational standpoint, implementing the convolution described by equations (2) (or 3) and (4) is an zz Ύ GxfΗ (x)dx, easy task. A word of caution, however, is in order regarding ס Ύ PY΄Η΅Ն xf(x)dx ס (G (z r ס ZXY|XX 0 x 0 x the wide range of X Sa(f) values to be considered in (2) achieving accuracy in the computation of GZ(z). In extreme cases, the exceedance of a particular soil-ground-motion where GW(w) is the CCDF of any random variable (RV), W— level equal to z may in fact occur both for a very small Nonlinear Soil-Site Effects in Probabilistic Seismic-Hazard Analysis 2113
rr bedrock ground motion significantly amplified by the nearly dependent onSaa(f) andS (fsc) and coupling it with the site- rr linear soil response or for a large bedrock ground motion specific joint hazard ofSaa(f) andS (fsc) rather than with the r greatly deamplified by the nonlinear behavior of the soil col- rock-hazard curves forSa(f) from conventional scalar PSHA. umn. Both the amplification function and the rock seismic- The site-specific joint hazard can be computed via vector- hazard curves need to span over a domain of X that is con- valued probabilistic seismic-hazard analysis (VPSHA) (Baz- siderably larger than the domain of Z for which the results zurro and Cornell, 2002). are desired. In this framework the convolution method for comput- ס s ס In the development of equation (2) we have used the ing surface-hazard curves for Z Sa(f), based on X1 סrr rare-event assumption. Relaxing this assumption implies Saa(f) and X2 S (fsc), can be described by the following permitting the occurrence of more than one event in the pe- equation: riod of interest. In this case, one would have to carefully distinguish the contributions to the random AF(f) that are z (Gx;xf (x ;x )dx dx , (5 ס (G (z the same for every event (e.g., the soil properties) from those ZYΎΎ |X12;X Η12 X 1;X 2 12 1 2 00 x1 that vary randomly from event to event (e.g., the ground- motion details). The occurrence of two or more events may or in discretized form: become important for MRE values greater than, say, 0.1, and this case is seldom of earthquake-engineering significance. z .[ ͚͚GxΗ;xp [x ;x ס (G (z Z Y|X12;X 1,j 2,iX 12;X 1,j 2,i allx1,j allx2,i x1,j rr Convolution: AF(f) Dependent onSaa(f) and S (fsc) (6) As pointed out in our companion article, in the fre- Thef (x ;x ) term in equation (5) is (numerically quency range below f the AF(f) is jointly negatively cor- X12;X 1 2 ס sc r related withSr(f) and positively correlated with M. M seems equal to) the site-specific mean rate density of ln Sa(f) ln ס a r to explain part of the variability in AF(f) because it carries X1 and ln Sa(fsc) ln X2 for a site obtained via VPSHA (Bazzurro and Cornell, 2002). In equation (6) the G information about the spectral shape and therefore about Y|X12;X r represents the CCDF of Y conditional on X1 and X2, and Sa(fsc), the bedrock spectral acceleration at the initial reso- r pX ;X [x1,j;x2,i] denotes the probability that the rock spectral- nance frequency of the soil. The higherSa(fsc) is, the more 12 significant the shift of the peak of the AF(f) function is to- acceleration values at frequencies f and fsc are in the neigh- ward lower f values, a shift that in turn tends to increase the borhood of x1,j and x2,i, respectively. As before, assuming Y X X G lognormality of , given 1 and 2, theY|X12;X in continuous value of AF(f)atf fsc. Hence, in this frequency range r form is given by: Sa(f) alone is somewhat an insufficient parameter for the estimation of AF(f) (Luco, 2001). z [( ln[mˆ (x ;x מ΅΄The consequence of this insufficiency is that the esti- ln Y|X12;X 12 mates ofmˆ rr andr that are obtained using a z x1 Uˆ סAF(f)Saa(f)lnAF(f)S (f) GxΗ;x limited dataset of real recordings are not necessarily accu- Y|X12;X 12 x1lnr Y|x12;x rate, and this may potentially lead to imprecise estimates of (7) s the hazard curves forSa(f) (Bazzurro and Cornell, 2002). The inaccuracy may arise because the record set does not where the estimates of the distribution parameters of Y, necessarily have the right, hazard-consistent spectral shape mˆ Y|X12;X (x1;x2), andrlnY|x 12;x , are obtained through the mul- ס r ס for a site as determined by the combinations of M–R events tiple regression of ln AF(f) ln Y on ln Sa(f) ln X1 and ס that nearby faults are more likely to generate. The inaccuracy r ln Sa(fsc) ln X2 (see companion article). can be reduced if the records are selected with care to ensure that the M reflects the controlling scenario for the site. Even with an infinitely large database of seismograms, however, PSHA with Soil-Specific Attenuation Equation a perfect record selection for a site may prove to be an im- The most straightforward approach for estimating possible task, because the scenario that dominates the hazard seismic-hazard curves at the surface of a soil deposit is to generally changes with hazard level. Note that this potential develop an attenuation relationship for that soil condition source of inaccuracy is a limitation of any statistical ap- and use it during conventional PSHA calculations (W. J. proach that provides estimates based on limited (and poten- Silva, personal comm., 1996). Developing a robust soil- tially biased) samples of data. specific attenuation law requires that a significant number of On the basis of the experience gained so far, we believe rock-ground-motion accelerograms (more than 100) be r that, givenSa(f) , the conditional dependence of AF(f)on propagated (and perhaps deconvolved first through a generic r Sa(fsc) in the frequency range below fsc is not strong enough rock profile; see, e.g., Kramer, 1996) through a numerical s to generate large errors in the estimate of the soil-hazard model of the sediment deposit. The computedSa(f) data curves, if neglected. Potential errors can be reduced or pos- could then be processed by using, for example, advanced sibly avoided altogether by keeping the soil response jointly regression-analysis techniques (e.g., Abrahamson and 2114 P. Bazzurro and C. A. Cornell
Youngs, 1992). These techniques account for correlations in Our companion article has shown that AF(f), condi- r the data, owing to multiple recordings of a single earthquake, tional on the spectral value at the rock level,Sa(f) , is vir- s for obtaining the desired attenuation equation forSˆa(f) , that tually independent of M and R, and other ground-motion s r is, the median ofSa(f) , given M, R, and h, where the vector parameters such asSa(fsc) or PGA. Assuming that a linear h represents source mechanism, wave travel path, and local regression in logarithmic space is appropriate to describe the r site conditions, characteristics that are commonly included dependence ofAFˆ(f) onSa(f) , at least within a limited range r in modern attenuation relations. ofSa(f) values, then: This direct method, although feasible, is rather imprac- ם ס tical. It produces, however, hazard curves at the soil surface ln AF(f) ln AFˆ(f) erlnAF(f)lnAF(f) ם r ם -that may be considered “exact.” The soil-hazard curves pro c01c ln Sa(f) erlnAF(f)lnAF(f), vided by this exact method can be used as a benchmark for (10) testing the validity of the proposed convolution approach mentioned in the two previous subsections. Such a validation where c0 and c1 are coefficients of the linear regression in was conducted in Bazzurro (1998) and will not be repeated logarithmic space of AF(f)onSr(f) ,e is a standard in this article. In contrast we explore here the use of existing a lnAF(f) normal variable, and rlnAF(f) represents the standard error of rock-attenuation laws amended with an additional term to estimation from the regression that has been found in work account for local soil conditions with uncertain properties. to date (Bazzurro and Cornell, 2004; U.S. Nuclear Regula- As pointed out earlier, this frequency-specific correction fac- tory Commission, 2001) to be not significantly functionally tor, which is simply the amplification factor, AF(f), can be r dependent onSa(f) (or M and R). obtained without driving an extensive dataset of records From these equations it follows that: through the soil column (Bazzurro and Cornell, 2004). If a linear predictive model for (log) AF(f) in terms of sr (ln AF(f םerrr ם (ln Sˆˆ(f ס (ln S (f (log)Sr(f) is appropriate, then closed-form equations can be aalnSaa(f)lnS (f) ם a derived to integrate AF(f) directly into an existing rock- erlnAF(f)lnAF(f). (11) r attenuation equation forSa(f) , transforming it into a site- specific soil-attenuation equation. If the linearity condition However from equation (10): does not hold, difficulties arise in finding a closed-form re- r ם lationship when coupling the attenuation error term for ln AFˆ(f) c01c ln Sa(f) ם r ם ס Sr(f), given M, R, and h, with the error term of the regression a c01c [ln Sˆa(f) erlnSrr(f)lnS (f)]. (12) r aa of (log) AF(f) on (log)Sa(f) . It is emphasized however, that the equations that follow can be applied also to piecewise– Substituting equation (12) into equation (11): linear models to approximate, for example, a quadratic be- (ln Sˆ(f (1 ם c) ם havior such as that shown in figure 8 in our companion ar- ln Ssr(f) c ticle (see also Fig. 2, to come). a 01 a .er םerrr(1 ם c) ם In logarithmic terms, equation (1) can be written as: 1lnSaa(f)lnS (f)lnAF(f)lnAF(f) (13) ם סsr ln Saa(f) ln S (f) ln AF(f). (8) s Equation (13) states that the medianSa(f) can be pre- r r dicted from the median ofSa(f) as: Analytical expressions for lnSa(f) in terms of M, R, and h are numerous in the literature (for a review, see Abraham- (ln S (f). (14 (1 ם c) ם ln Sˆˆsr(f) c son and Shedlock, 1997). Such attenuation functions for a a 01 a specific frequency value f have typically the form: R. K. McGuire, (personal comm., 2002) found that
r r elnAF(f) andelnSa(f) are mildly negatively correlated, mostly rr ם ס ln Sa(f) g1ln(M, R, h) erSaa(f)lnS (f) r for large values of lnSa(f) . However, if we (conservatively) r rr ם ס ln Sˆa(f) erlnSaa(f)lnS (f), (9) neglect this negative correlation, then the dispersion measure s for lnSa(f) becomes: rr where lnSˆaa(f) is the median of lnS (f) , given M, R, and h. 22 2 ם ם sr Ί The remaining variation unexplained by the nonlinear re- rlnSaa(f)1(c 1) r lnS (f)lnr AF(f). (15) gression is captured by a (standardized Gaussian) random
r variable,elnSa(f) , which can be defined as the number of (log- Equations (14) and (15) permit the assembly of an arithmic) standard deviations by which the random (loga- approximate attenuation relationship for the soil-surface- s rithmic) spectral acceleration deviates from its median value motion parameter,Sa(f) , by coupling the available attenua- r as predicted by an attenuation equation in terms of M, R, tion equation for the rock-motion parameter,Sa(f) (equation s r and h. Values ofrlnSa(f) are available for standard attenuation 14), with the site-specific regression for AF(f)onSa(f) relationships. (equation 10). Nonlinear Soil-Site Effects in Probabilistic Seismic-Hazard Analysis 2115
r It is interesting to note that because, in general, AF(f) model of ln AF(f) on lnSa(f) (equation 10). The constant r is negatively correlated with the amplitude ofSa(f) —that is, k1 is the slope (in log–log scale) of the straight-line tangent r ˆr the higher the value ofSa(f) , the lower the value of AF(f), to the hazard curve at the pointSa,z(f) (or, alternatively, the at least at frequencies f higher than 0.5 Hz (see companion slope of the secant straight line that approximates the hazard article), the coefficient c1 is negative. Hence, from equation curve in the neighborhood of the tangent point). Such a lin- 2 (15) it follows that, depending on how largerlnAF(f) is, the earization is usually an adequate approximation over an sr variability inSaa(f) may be less than the variability inS (f) . order-of-magnitude probability range around the tangency ס -This phenomenon, which is presumably due to the nonlin point. To quantify this term, note that k1 1/log a10, where earity in the soil response during severe ground shaking, has r a10 is the factor by which one must increaseSa(f) in order in fact been observed from real recordings on rock and soil to increase the hazard (or reduce the mean return period) by conditions (see Abrahamson and Sykora, 1993; Toro et al., one order of magnitude (U.S. Department of Energy, 1993). 1997). As mentioned before, accounting for the negative cor- r Hence, loosely speaking, the annual probability that relation between AF(f) andSa(f) also avoids the concern Sr(f) exceeds a given level of z is equal to the value of the pointed out by Silva et al. (2000) that artificially inflated a bedrock-hazard curve at a bedrock value x that induces a soil soil-hazard curves may be obtained owing to overestimated ground motion in exceedance of z (50% of the time) multi- uncertainty onSs(f) if one simply adds the amplification a plied by an exponential correction factor. Given the typical uncertainty to that already contained in the rock-hazard r in the order of 0.3, and of the slopes ס values of r curve without recognition of this negative correlation. lnAF(f) of the hazard curves for most sites commonly varying from An application of this method will be shown later. מ 2.5 to 6 and for c1 ranging from 0 to approximately 0.8 (values observed for the two soil sites in our companion An Analytical Estimate of the Soil Hazard article), such a correction factor is often not larger than 10. Under the same assumptions of a previous subsection, However, for some combinations of these parameter val- ס namely, that the dependence of AF(f)onM and R, given ues—for example, large negative c1 and large r rlnAF(f)— r Sa(f), is weak, the convolution in equation (2) can be re- it may become very large. We do not advise the use of this written without the explicit use of AF(f) as follows: method when the correction factor exceeds 10.