Sellami
Sellami
Sellami
Seismische Gefährdungsanalyse
-------------------------------------------
The general theme of this course is earthquake risk. The concept of risk includes hazard
and vulnerability. The first part has dealt with earthquakes, where when they occur, how
big they are and why they happen. The second part is about the effects.
Hazard assessment is to evaluate, for a certain place, how frequent and how strong
earthquake will be felt, in order to take measure to reduce the possible damages. In other
terms, it is to qualify and quantify the level of ground motion in a site due to the earthquake.
Seismic hazard maps depict the levels of chosen ground motions that likely will, or will not,
be exceeded in specified exposure times. The ground motion can be the intensity of the
earthquake, displacement, velocity or acceleration of the seismic wave at the site.
The hazard can be estimated using deterministic or probabilistic methods. The probabilistic
method (Cornell 1968), broadly applied, will be first described. Alternative methods will be
presented in the discussion paragraph.
The three major elements of the probabilistic method are: 1) the characterisation of seismic
sources; 2) the characterisation of attenuation of ground motion; and 3) the actual
calculation of probabilities.
Probabilistic method basic principles
The different steps of the probabilistic (Cornell) method are outlined on the sketch (figure
1.) below (Ruettener 95):
Figure 1.
1) the characterisation of seismic sources is usually achieved by
• compilation of an earthquake catalogue (a)
• delineation of the seismic sources (b)
• magnitudes-frequency distribution (c)
2) the characterisation of attenuation of ground motion is described by attenuation
functions (d)
3) computation of the probability analysis (c)
The first element of seismic hazard assessment, the characterisation of seismic sources,
involves obtaining ‘robust’ answers to three questions, which have been addressed in the
first part of this course:
• Where do earthquakes occur?
• How often do earthquakes occur?
• How big can we expect these earthquakes to be?
Seismicity catalogues are the fundamental data base used to determine where, how often,
and how big earthquakes are likely to be. However and seismicity statistics are based on
geologically short catalogues. For magnitude above 6, the completeness is less than 1000
year in Switzerland (one of the best case) and less than 200 years in California). Therefore
other deformation data are examined. The results from seismic monitoring, the historic
record, geodetic monitoring, and the geologic record are combined to characterise seismic
sources. These data, when available, are used to interpret seismic source zones. Because
many interpretations of the input data are possible, large uncertainties are associated with
source characterisation.
The identification of the seismogemic sources in the region is very important. In theory
earthquake sources are faults. In most of the places, the earthquake distribution does not
coincide with known fault visible at the earth surface. So in practice, the shape of a source
can be a fault, but they are surfaces (area zones) when active fault can not be recognised,
which is the most common case (for example in Switzerland), The Figure 2 (Rutenner 1995)
shows the historical seismicicity map (time span 1300-1994 and earthquake intensity V and
above) together with the geographic distribution of the seismic sources (Sägesser and
Mayer-Rosa, 1978). The shape of the source zone depends strongly on the earthquake
distribution ‘per extenso’ on the catalogue. Their design is subject to a part of subjectivity or
expert judgement (Schenk 1996). This is illustrated, for example, by the source zoning of
the Ibero-maghreb region (Jimenez et al 1999) shown on Figure 3.
Figure 2. Seismicicity (1300-1994) and sources zones (Ruettener 1995)
where n(M) is the number of events per year having magnitudes greater than M. a and b are
constants defined by regression analysis. The slope of the magnitude-frequency Gutenberg-
Richter defines the "b value" parameter.
Maximum magnitudes
Mupp and Mlow are the upper and lower bound magnitude on the source. Mlow, is the
magnitude below which no engineering-significant damage is expected and Mupp represents
the maximum expected magnitude. The maximum magnitude is related to the tectonic
setting, geometry, and type of the seismic source. Although no standard method exists for
assigning a maximum magnitude to a given fault, empirical correlation are used based on
the length of rupture of the fault, the total length of the fault trace or the area of the fault
rupture zone. In most cases, faults cannot be clearly recognised, maximum magnitude are
than deduced either from the earthquake catalogue, from the recurrence rate (extrapolating
the Gutemberg-Richter relationship) or from paleoseismicity studies.
The ground motion at a site, for example Peak Ground Acceleration depends on the
earthquake source, the seismic wave propagation and the site response. Earthquake source
signifies the earthquake magnitude, the depth and the focal mechanism, the propagation
depends mainly on the distance to the site. The site response deals with the local geology
(site classification); it is the subject of microzonation.
The basic functional (logarithmic) form for ground motion attenuation relationship is
defined as (Reiter 1990)
These relationship are for a specific site classification (hard rock, soft rock, etc.). Hazard
values calculated for rock/stiff soil sites (the most common site classifications) are lower
than hazard values calculated for soil sites. A comprehensive review and application for
Switzerland can be found in Smit (1996).
Example of different attenuation functions used in the different countries of the Ibero-
maghreb region are shown on table 1 and figure 4 (Jimenez et al. 1999).
Table 1
Country Attenuation law
For the calculation of the seismic hazard map of the whole region, only one attenuation law
was considered Joyner and Boore (1981). This relationship is represented on the figure 4
with one standard deviation in log PGA. It could roughly represent an average of the
different laws used in the region.
To adjust the parameters of the attenuation function for intensities (ex. Portugal and Spain in
table 1), macroseismic data are needed. For example the isoseismal map of the 1356 Basel
earthquake (after Mayer-Rosa and Cadiot, 1979) shown on figure 5.
Figure 5. Isoseismal map of the 1356 Basel event (Mayer-Rosa and Cadiot, 1979)
The figure 6 depict the observed intensity distribution of an event (13 March 1964 in
Central Switzerland) together with the attenuation law, for subalpine region
(Ruettener1995). It shows a significant scattering of data around the attenuation function.
This scattering depends not only on physical effect like the directivity of the radiated energy
or the local geology but also on the quality of intensity data.
Calculation of probabilities
The third element of hazard assessment, the actual calculation of expected ground motion
values, involves determining an annual frequency of exceedance of the ground motion
parameter of interest, then summing over the time period of interest.
Probability of Exceedance
The annual mean number of events in which GM exceeds the specified ground motion level
Z is calculated by summing up the incremental contributions of the N sources, taking into
account:
• the annual mean rate of recurrence of earthquakes of magnitude Mi on each source
• given the occurrence of an earthquake of magnitude Mi, the distance from the rupture
surface to the site,
• given the occurrence of an earthquake of magnitude Mi at a distance of rj, the probability
that the value of ground motion at the site exceeds a specified level.
The probability (ex 10%) of exceeding a specified level of ground motion (such as a certain
value of peak acceleration) at a selected site within the time interval of interest (such as 50
years) is calculated by:
1/ combining the three probability functions,
2/ integrating over all possible earthquake, location, magnitude for a source
3/ integrating over all the sources
The probability of exceedance, in a specific time interval, that ground motion amplitude a*
is exceeded can be expressed as follows (McGuire, 1993)
P[A > a * in time t ] / t = ∑ v i ∫ ∫ G A m,r (a *)f m (m) f r (r m)dmdr
i
where ν i is the mean annual rate of occurence in source i, G is the probability that an
earthquake given m and r exceeds ground motion a* at a specific site. f m(m) is the
probability density function for magnitude, f r (r m) and is the probability distribution
function for distance.
Probability terms
Figure 7.a and b. Hazard outputs of the historical method for Brig (Ruettener 1995)
Return period = -T / ln(1-P(Zz))
Probability of non-
Probability of exccedance in approximate average return
exceedance in a given
a given exposure time period in years
exposure time
Incorporation of Uncertainties
The uncertainties of the basic input data must be taken into account. Uncertainties are
introduced either by lack of data or/and lack of knowledge. There are random uncertainties
which could be incorporated in the hazard curve calculation. There are also systematic or
modelling uncertainties for example for the choice of maximum magnitude, the correct
ground motion model. These uncertainties are taken into account by developing alternative
strategies and models in the interpretation of those input data for which significant
uncertainties are known to exist. For example, multiple source zone models may be defined.
Hazard calculations from each model are then combined using various schemes that produce
a weighted mean (or median) hazard value. It is the logic tree analysis (Frankel 1995).
The probabilistc method, based on a poissonian model, does not depict a possible variation
of the seismicity in time, because of the hypothesis of stationarity of the model. Some
research is developing in this direction.
Seismic hazard assessment in low seismicity areas is much more subject to large errors than
in areas with high earthquake activity. This is specifically the case if the time span of the
available data catalogue is considerably smaller than the mean return interval of large
events, for which the hazard data has to be calculated.
Alternative techniques of the probabilistic method are mainly to avoid the delimitation of
source zone. The historical method for example, applied in Switzerland by Ruettener (1995)
estimates the ground motion at a site from each event of the historical catalogue. The
probability distribution of the historic occurrences of earthquake is calculated from the
earthquake rate.
The deterministic approach evaluates the maximum expected ground motion at a site,
resulting from the strongest potential earthquake at the nearer possible distance.
Deterministic approaches are often used to evaluate the hazard for a selected site. They are
in general more conservative. But they do not take into account uncertainties nor an estimate
of frequency of occurrence which is needed by decision makers for planning purposes.
The probabilistic method allows to take uncertainties into account and is easily applied.
However, as the quality of the output depends strongly on the quality of the input
parameters (earthquake catalogues, strong motion relationships), and it might be that there is
little transparency on the quality and the integration of the data. On the other hand the
output reflects the state of the data and is subject to improvement.
Probabilistic approach can be applied to mapping the hazard for different probabilities and
exposure time and for different area sizes (local / regional (Figure 8) / global).
Global seismic hazard assessment maps: The variations in each element of the seismic
hazard assessment lead to differences in the estimated hazard along the national borders.
Some research programs (DACH GSHAP, SESAME) have been launched in order to
homogenise this hazard (Figure 9).
Figure 8. Earthquake hazard map of the Ibero-Maghreb region. PGA [m/s2] with
90% probability of non-exceedance in 50 years.
Figure 9. Horizontal peak ground acceleration seismic hazard map representing stiff
site conditions for an exceedance or occurrence rate of 10% within 50 years
(Gruenthal et al. 1999).
Links
Application projects:
• USGS national seismic hazard mapping
http://gldage.cr.usgs.gov/eq/
• Global Seismic Hazard Program
http://seismo.ethz.ch/GSHAP/
- Catalogue / Katalog / catalogue: record of historical earthquake. It can be used for the
time span of completness for different magnitudes.
- Exposure time / - / - : The period of time (period of interest ex 50 years) that a structure
or a facility is exposed to hazard. Linked to the design lifetime of the structure.
- Seismic activity/ Aktivität/ activité sismique: (of a source) it is described with the
frequency/magnitude relationship ( Häufigkeitverteilung) characterise by the b-value.
- Source / Quelle / source: the source of energy release causing an earthquake. For
computing purpose, they can be areas, faults or points the source is characterise by its
activity and the maximal an minimal magnitudes