Probabilistic Seismic Hazard Analysis
Probabilistic Seismic Hazard Analysis
Overview
History
1969 - Allin Cornell BSSA paper
Rapid development since that time
Probabilistic Seismic Hazard Analysis
Overview
Deterministic (DSHA)
Assumes a single “scenario”
Select a single magnitude, M Ground
Ground
Select a single distance, R motion
motion
Assume effects due to M, R parameters
parameters
Probabilistic (PSHA)
Assumes many scenarios
Consider all magnitudes Ground
Ground
Consider all distances motion
motion
Consider all effects parameters
parameters
Probabilistic Seismic Hazard Analysis
Overview
Why?
Why? Because
Because we
wedon’t
don’tknow
know when
when
earthquakes
earthquakeswill
will occur,
occur,wewe don’t
don’t
know
knowwhere
wherethey
they will
will occur,
occur, and
andwe
we
don’t
don’t know
knowhow
how big
big they
theywill
willbebe
Probabilistic (PSHA)
Assumes many scenarios
Consider all magnitudes Ground
Ground
Consider all distances motion
motion
Consider all effects parameters
parameters
Probabilistic Seismic Hazard Analysis
PSHA
PSHAcharacterizes
characterizesuncertainty
uncertainty in
in
location,
location,size,
size, frequency,
frequency,and
andeffects
effects
of
of earthquakes,
earthquakes, and and combines
combinesallall of
of
them
them toto compute
computeprobabilities
probabilitiesof
of
different
differentlevels
levelsofof ground
ground shaking
shaking
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Vertical
VerticalFaults
Faults
rjb
rrup
rseis Seismogenic
depth
rhypo
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Dipping
DippingFaults
Faults
rjb=0 rjb
rseis
rseis & rrup
rrup rhypo
rhypo
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Source-site
Source-site
distance
distancedepends
depends
on
on where
whererupture
rupture
occurs
occurs
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Source-site
Source-site
distance
distancedepends
depends
on
on where
whererupture
rupture
occurs
occurs
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
rmin
rmax fR(r)
rmin rmax r
pdf
pdffor
forsource-
source-
site
sitedistance
distance
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Measure
Measurelength
lengthof
offault,
fault,
rmax LLi,,between each pair
i between each pair
of
ofadjacent
adjacentcircles
circles
Assign
Assignweight
weightequal
equalto
to
LLi/L to each
i/L to each
corresponding
corresponding
distance
distance
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Divide
Divideentire
entirefault
faultinto
into
equal
equallength
lengthsegments
segments
rmin
Compute
Computedistance
distancefrom
from
site
siteto
tocenter
centerof
ofeach
each
rmax segment
segment
Create
Createhistogram
histogramof of
Linear source source-site
source-sitedistance.
distance.
Accuracy
Accuracyincreases
increases
with
withincreasing
increasing
number
numberof ofsegments
segments
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Divide
Dividesource
sourceinto
intoequal
equal
area
areaelements
elements
Compute
Computedistance
distancefrom
from
Areal Source center
centerof
ofeach
eachelement
element
Create
Createhistogram
histogramofof
source-site
source-sitedistance
distance
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Divide
Dividesource
sourceinto
intoequal
equal
volume
volumeelements
elements
Compute
Computedistance
distancefrom
from
center
centerof
ofeach
eachelement
element
Create
Createhistogram
histogramofof
source-site
source-sitedistance
distance
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Unequal
Unequalelement
elementareas?
areas?
Create
Createhistogram
histogramusing
using
weighting
weightingfactors
factors--weight
weight
according
accordingtotofraction
fractionof
oftotal
total
source
sourcearea
area
Probabilistic Seismic Hazard Analysis
Uncertainty in source-site distance
Quick
Quickvisualization
visualizationof
ofpdf?
pdf?
Use
Useconcentric
concentriccircle
circle
approach
approach--lets
letsyou
you“see”
“see”
basic
basicshape
shapeofofpdf
pdfquickly
quickly
Probabilistic Seismic Hazard Analysis
Characterization of maximum magnitude
Also
Alsoneed
needtotoknow
know
distribution
distributionofof magnitudes
magnitudes
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
Gutenberg-Richter
Southern California earthquake data - many faults
Counted number of earthquakes exceeding different
magnitude levels over period of many years
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
NM log NM
M M
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
log M
Mean
Meanannual
annual rate
rate
of
ofexceedance
exceedance
log M
Return
Returnperiod
period
(recurrence
(recurrenceinterval)
interval) 0.001 1000 yrs
log TR
TTRR == 11/
/MM M
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
log M
10a
b
Gutenberg-Richter
Gutenberg-Richter
Recurrence
RecurrenceLaw
Law
log TR
log
logMM== aa-- bM
bM 0 M
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
log M = a - bM
ln M = - M
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
Then
fM(m) = e-m
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
Then,
fM(m) = e-m-mo)
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
M = 12 m = 0.437 TR = 2.29
M
M>> 12
12every
everytwo
twoyears?
years?
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
Mmax M
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
Paleoseismic investigations
Show similar displacements in each earthquake
Inividual faults produce characteristic earthquakes
Characteristic earthquake occur at or near Mmax
Could be caused by geologic constraints
More research, field observations needed
Probabilistic Seismic Hazard Analysis
Distribution of earthquake magnitudes
log m
Characteristic
Seismicity data Earthquake
Recurrence
Law
Geologic data
Mmax M
Probabilistic Seismic Hazard Analysis
Predictive relationships
Standard error
log m- use to evaluate conditional probability
Y = Y*
ln Y
M = M*
R = R* log R
Mmax M
Probabilistic Seismic Hazard Analysis
Predictive relationships
ln Y
P[Y > Y*| M=M*, R=R*]
ln Y
Y = Y*
M = M*
R = R* log R
M
Probabilistic Seismic Hazard Analysis
Temporal uncertainty
Poisson process
n e
P [N n ]
n!
Poisson process
Letting = t
( t ) e t
n
P [N n ]
n!
Then
P [N 0] P [N 1] P [N 2] P [N 3] ... P [n ]
1 P [N 0 ]
1 e t
Probabilistic Seismic Hazard Analysis
Temporal uncertainty
Poisson process
P 1 e t
= 1/1000 = 0.001
P = 1 - exp[-(0.001)(100)] = 0.0952
Probabilistic Seismic Hazard Analysis
Temporal uncertainty
P = 1 - exp[-(0.001)(1000)] = 0.632
Solving for ,
ln(1 p)
t
Probabilistic Seismic Hazard Analysis
Temporal uncertainty
ln(1 0.1)
0.0021
50
U U U
P[A] = P[A B1] + P[A B2] + … + P[A B N]
B1 B2 B3 Total
Total
Probability
Probability
A Theorem
Theorem
B5
B4
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
P [Y y *] P [Y y * | X]P [ X] P [Y y * | X] f X ( X)dx
P [Y y *] P [Y y * | m, r ] f M (m ) f R (r )dmdr
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
P [Y y *] P [Y y * | m, r ] f M ( m ) f R ( r )dmdr
y * P [Y y * | m, r ] f M (m ) f R (r )dmdr
Ns
y * i P [Y y * | m, r ] f Mi (m ) f Ri (r )dmdr
i 1
NS NM NR
y * i P [Y y * | m j , r k ]P [M m j ]P [R r k ]
i 1 j 1 k 1
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
NS NM NR
y * i P [Y y * | m j , r k ]P [M m j ]P [R r k ]
i 1 j 1 k 1
All All
Allpossible
possiblemagnitudes
magnitudesareare
Allsites
sitesare
are
considered considered
considered--contribution
contributionof
of
considered All
Allpossible
possibleeffects
effectsare
are each
eachisisweighted
weightedbybyits
its
considered
considered--each
eachweighted
weighted probability
probabilityof
ofoccurrence
occurrence
by
byits
itsconditional
conditionalprobability
probability
of
ofoccurrence
occurrence
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
NM x NR possible combinations
Each produces some probability of exceeding y*
Must compute P[Y > y*|M=mj,R=rk] for all mj, rk
r NR
r1
m1 m 2 m 3 mNM
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
r1 r2 log R
r3
rN
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
“Build” hazard by:
computing conditional probability for each element
multiplying conditional probability by P[mj], P[rk], i
Repeat for each source - place values in same cells
rNR
P[Y > y*| M=m2, R=r3]
P[Y > y*| M=m2, R=r2]
P[Y > y*| M=m2, R=r1]
r1
m 1 m2 m3 mN M
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
When complete (all cells filled for all sources),
rNR
P[Y > y*| M=m2, R=r3]
P[Y > y*| M=m2, R=r2]
P[Y > y*| M=m2, R=r1]
r1
m 1 m2 m3 mN M
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
Seismic
Hazard
rNR Curve
P[Y > y*| M=m2, R=r3]
log y*
log TR
P[Y > y*| M=m2, R=r2]
P[Y > y*| M=m2, R=r1]
r1
m 1 m2 m3 y* mN M
Probabilistic Seismic Hazard Analysis
Combining uncertainties - probability computations
Seismic hazard
curve shows the
mean annual rate of
log amax
exceedance of a
log y*
log TR
log TR
particular ground
motion parameter. A
seismic hazard curve
is the ultimate result
of a PSHA.
amax y*
Probabilistic Seismic Hazard Analysis
Using seismic hazard curves
Probability of exceeding
amax = 0.30g in a 50 yr
period?
log amax
log TR
P = 1 - e-t
= 1 - exp[-(0.001)(50)]
= 0.049 = 4.9%
In a 500 yr period?
0.001
P = 0.393 = 39.3%
amax=0.30g
Probabilistic Seismic Hazard Analysis
Using seismic hazard curves
log TR
or
TR = 475 yrs
0.0021 475 yrs
Use seismic hazard
curve to find amax value
corresponding to
= 0.0021
amax=0.21g
Probabilistic Seismic Hazard Analysis
Using seismic hazard curves
Contribution of sources
each source
log TR
2 Plot seismic hazard
1 curves for each source
and total seismic
hazard curve (equal to
3 sum of source curves)
Curves may not be
parallel, may cross
amax Shows which source(s)
most important
Probabilistic Seismic Hazard Analysis
Using seismic hazard curves
Spectral accelerations
log TR
Other
Choose desired -value
Read corresponding
parameter values from
seismic hazard curves
amax, Sa
Probabilistic Seismic Hazard Analysis
Using seismic hazard curves
Spectral accelerations
log TR
Other
Choose desired -value
Read corresponding
parameter values from
seismic hazard curves
amax, Sa
Probabilistic Seismic Hazard Analysis
0.1
2%
2%in
in50
50yrs
yrs
Crustal
Peak
Peakacceleration
acceleration
0.01
amax
Intraplate
0.001
Interplate
0.0001
Probabilistic Seismic Hazard Analysis
0.1
2%
2%in
in50
50yrs
yrs
SSa(T = 3 sec)
a(T = 3 sec)
0.01 Crustal
amax Intraplate
0.001
Interplate
0.0001
Probabilistic Seismic Hazard Analysis
Uniform hazard spectrum (UHS)
Sa
Uniform
Uniform
Hazard
Hazard
Spectrum
Spectrum
T
Probabilistic Seismic Hazard Analysis
Disaggregation (De-aggregation)
Common question:
What magnitude & distance does that amax value correspond to?
Common question:
What magnitude & distance does that amax value correspond to?
USGS disaggregations
Seattle,
Seattle,WA
WA
2%
2%in in50
50yrs
yrs
(T
(TRR==2475
2475yrs)
yrs)
USGS disaggregations
Olympia,
Olympia,WA
WA
2%
2%in in50
50yrs
yrs
(T
(TRR==2475
2475yrs)
yrs)
USGS disaggregations
Olympia,
Olympia,WA
WA
2%
2%in in50
50yrs
yrs
(T
(TRR==2475
2475yrs)
yrs)
= 1.2 = 2.2
ln Y
Y = y*
r1 r2 log R
r3
rN
Probabilistic Seismic Hazard Analysis
Logic tree methods
7.0 (0.2)
G-R 7.2 (0.6)
(0.7) 7.5 (0.2)
A&S
(0.5) 7.0 (0.2)
Char. 7.2 (0.6)
(0.3) 7.5 (0.2)
Probabilistic Seismic Hazard Analysis
Logic tree methods
Attenuation Magnitude
Mmax
Model Distribution
7.0 (0.2)
Sum
Sumof
ofweighting
weightingfactors
factors G-R 7.2 (0.6)
coming
comingout
outof
ofeach
eachnode
node (0.7) 7.5 (0.2)
must
mustequal
equal1.0
1.0 BJF
(0.5) 7.0 (0.2)
Char. 7.2 (0.6)
(0.3) 7.5 (0.2)
7.0 (0.2)
G-R 7.2 (0.6)
(0.7) 7.5 (0.2)
A&S
(0.5) 7.0 (0.2)
Char. 7.2 (0.6)
(0.3) 7.5 (0.2)
Probabilistic Seismic Hazard Analysis
Logic tree methods
Attenuation Magnitude
Mmax
Model Distribution
7.0 (0.2)
G-R 7.2 (0.6)
(0.7) 7.5 (0.2)
BJF
(0.5) 7.0 (0.2)
Char. 7.2 (0.6)
(0.3) 7.5 (0.2)
7.0 (0.2)
G-R 7.2 (0.6)
(0.7) 7.5 (0.2)
0.5x0.7x0.2 = 0.07
A&S
(0.5) 7.0 (0.2)
Char. 7.2 (0.6)
(0.3) 7.5 (0.2)
Probabilistic Seismic Hazard Analysis
Logic tree methods
Attenuation Magnitude
Model Distribution
Mmax w
7.0 (0.2) 0.07
G-R 7.2 (0.6) 0.21
(0.7) 7.5 (0.2) 0.07
BJF
(0.5) 7.0 (0.2) 0.03
Char. 7.2 (0.6) 0.09
(0.3) 7.5 (0.2) 0.03
7.0 (0.2) 0.07
Final
Finalvalue
valueofofYYis
is G-R 7.2 (0.6) 0.21
obtained
obtainedas as (0.7) 7.5 (0.2) 0.07
weighted
weightedaverage
averageof of A&S
all
allvalues
valuesgiven
givenby
by (0.5) 7.0 (0.2) 0.03
terminal
terminalbranches
branches Char. 7.2 (0.6) 0.09
of
oflogic
logictree
tree (0.3) 7.5 (0.2) 0.03
Probabilistic Seismic Hazard Analysis
Logic tree methods
Cascadia interplate
2 attenuation relationships
2 updip boundaries
3 downdip boundaries
2 return periods
4 segmentation models
2 maximum magnitude approaches
192 terminal branches
Probabilistic Seismic Hazard Analysis
Logic tree methods
Cascadia intraplate
2 intraslab geometries
3 maximum magnitudes
2 a-values
2 b-values
24 terminal branches
Probabilistic Seismic Hazard Analysis
Logic tree methods
Total
TotalPSHA
PSHArequired
requiredanalysis
analysisof
of 612
612combinations
combinations