Risk Targeting in Seismic Design Codes: The State of The Art, Outstanding Issues and Possible Paths Forward
Risk Targeting in Seismic Design Codes: The State of The Art, Outstanding Issues and Possible Paths Forward
Risk Targeting in Seismic Design Codes: The State of The Art, Outstanding Issues and Possible Paths Forward
Abstract Over the past decade there have been various studies on the development
of seismic design maps using the principle of “risk-targeting”. The basis of these
studies is the calculation of the seismic risk by convolution of a seismic hazard
curve for a given location (derived using probabilistic seismic hazard analysis) with
a fragility curve for a code-designed structure (ideally derived from structural
modelling). The ground-motion level that the structure is designed for is chosen so
that the structure has a pre-defined probability of achieving a certain performance
level (e.g. non-collapse). At present, seismic design maps developed using this
approach are only widely applied in practice in the US but studies have also been
conducted on a national basis for France, Romania, Canada and Indonesia, as well
as for the whole of Europe using the European Seismic Hazard Model. This short
article presents a review of the state of the art of this technique, highlighting efforts
to constrain better some of the input parameters. In addition, we discuss the diffi-
culties of applying this method in practice as well as possible paths forward,
including an empirical method to estimate an upper bound for the acceptable col-
lapse and yield risk.
Keywords Seismic hazard Earthquake engineering Fragility curves
Risk targeting Design Acceptable risk
1 Introduction
Current seismic building codes (e.g. Eurocode 8), based on results from a proba-
bilistic seismic hazard analysis (PSHA), generally adopt a constant hazard approach
to define the ground motions used for design. In other words, the peak ground
acceleration (PGA, or other intensity measure, IM, e.g. spectral acceleration) used
for design in one location has the same probability of being exceeded in a given
year as the design PGA in another location. Often this annual probability is
1/475 = 0.0021 (equivalent to 10% in 50 years or a return period of 475 years
assuming a Poisson process). Ten years ago, Luco et al. (2007) proposed a new
approach that targets a constant risk level across a territory. This has three principal
advantages over the use of design levels defined in the traditional way: trans-
parency, a uniform risk level across a territory and the ability to compare (and
ideally control) risk for different types of hazard (e.g. earthquake and wind). It does
come, however, with the disadvantage of making more choices explicit, rather than
implicitly assumed because of convention (e.g. the choice of 475 years as the
design return period).
The procedure of Luco et al. (2007), although often using different input
parameters (see below), has been applied to France (Douglas et al. 2013), Romania
(Vacareanu et al. 2017), Indonesia (SNI 2012), Canada (Allen et al. 2015) and at a
European scale (Silva et al. 2016), as well as forming the basis of the current US
seismic design code (ASCE 2010). Despite its numerous attractions (see above) and
the fact that it is a relatively simple procedure to implement, there are a number of
outstanding issues. For example, Douglas et al. (2013) note that the collapse
probabilities targeted by Luco et al. (2007) appear to be at least an order of mag-
nitude too high when compared with observed damage in previous earthquakes
(also see Sect. 3.1 of this article).
The next section presents an overview of the risk targeting approach and dis-
cusses previous choices of the critical input parameters. Section 3 highlights the
outstanding problems and some potential solutions, which are currently being
investigated by the authors.
The risk of collapse (or other level of structural damage) of a building at a given site
from earthquake shaking can be estimated by convolving the seismic hazard curve,
expressing the probability of different levels of ground motion, with the fragility
curve, expressing the probability of collapse given these ground motions (e.g.
Kennedy 2011). This so-called “risk integral” forms the basis of the risk-targeting
approach. For this approach there needs to be a link between the design acceleration
and the fragility curve used to compute the risk of collapse. For a standard fragility
curve based on the lognormal distribution, a single point on this curve (if the
Risk Targeting in Seismic Design Codes … 213
Fig. 1 Flowchart of the method to find risk-targeted design parameters. l is the mean of the
fragility curve assuming a lognormal distribution
214 J. Douglas and A. Gkimprixis
A summary from the literature of the three key inputs to the procedure is given in
this section (we are assuming that the seismic hazard curve has been correctly
defined by a recent PSHA for a wide range of probabilities of exceedance and that
lognormal fragility curves are used).
ASCE (2005) provides design criteria for nuclear power plants and critical
facilities based on a risk-targeting framework. The targeted risk ranges from 10−5 to
10−4, while for lower design categories values up to 10−3 are given in Braverman
et al. (2007). To comply with the code, a nuclear plant must have a smaller than 1%
chance of unacceptable performance for the design IM and less than 10% for 1.5
times the design IM.
Kennedy (2011) comments on the ASCE (2005) approach and implements it for 28
US nuclear plants. With a b in the range 0.3–0.6, the Seismic Core Damage
Frequencies (SCDFs) are between 6 10−7 and 6 10−6. Based on this, the United
States Nuclear Regulatory Commission (2007) has set a target SCDF of 10−5.
Luco et al. (2007) use generic fragility curves with b = 0.8. They assume a
probability of collapse of 10% under the ground motion corresponding to a
2475-year return period. The ATC-63 project initially proposed this value after
analysing buildings designed with ASCE 7-05. Luco et al. (2007) find that the 2003
NEHRP design ground motions led to non-uniform risk across US territory. By
targeting the average collapse probability in 50 years, which they found to be
around 1%, Luco et al. (2007) calculated new design ground motions with ratios of
0.7–1.15 of the 2003 NEHRP design ground motions.
The results of that report were considered for the ASCE 7-10 and FEMA (2009)
NEHRP Provisions. The latter propose the convolution of a hazard curve with a
fragility curve with b = 0.8 (decreased to 0.6 in ASCE 7-10), X = 0.1 and a target
probability of collapse of 1% in 50 years. Liel et al. (2015) comment on ASCE 7-10
and suggest modifications to consider subduction earthquakes and near-fault effects.
With the same X and b as in the regulation, the probability of collapse in 50 years
varied from 0.21 to 0.62%, for areas affected by subduction earthquakes, while at
near-fault sites, the risk was much higher, reaching in some cases 6%.
Allen et al. (2015) comment on the implementation of risk-targeting in Canada
for future versions of its National Building Code (NBCC). For comparison, ground
motion values derived from the risk-targeting approach of Luco et al. (2007) are
divided with those proposed in 2015 NBCC (for a 2%-in-50-year exceedance
probability). The resulting risk coefficients for spectral accelerations at 1.0 s are
between 0.84 and 1.00 in southwest Canada.
In Goulet et al. (2007) the collapse probabilities for a four-storey
reinforced-concrete (RC) frame under the 2%-in-50-year ground motion were cal-
culated as between 0 and 2%; after considering structural uncertainties the values
increased to 2–7%. They also found a variation of Y between 0.1 10−4 and
0.5 10−4 and between 0.4 10−4 and 1.4 10−4, when they considered struc-
tural uncertainties.
Risk Targeting in Seismic Design Codes … 215
As can be seen by Sect. 2 there is a rapidly increasing number of studies that have
attempted to apply the risk-targeting approach to different areas or sought to con-
strain the various inputs upon which this technique relies. Despite these many
studies a number of problems remain, which we highlight in this section along with
a brief discussion of potential solutions.
Risk Targeting in Seismic Design Codes … 217
As noted above, when evaluating the risk integral, ground motions for very low
probabilities of exceedance sometimes need to be estimated. A power-law
extrapolation appears to work well for most examples we have studied to date.
Further research on the technicalities of the calculations, however, needs to be
conducted to define a stable procedure that works for all possible inputs.
Targeting a non-zero value of collapse risk accepts that some buildings will
collapse (potentially leading to human casualties) in earthquakes, even when they
are designed in compliance with the building code. From a moral point of view this
is problematic and it leads to the difficulty of trying to define what risk is “ac-
ceptable”. One potential solution to this is to estimate an upper bound on the risk
that has been “accepted” historically based on the levels of observed damage in
previous earthquakes, as attempted by Labbé (2010). This “risk” is an upper bound,
as generally after every damaging earthquake the population lament the damage
that occurred. In the following Sect. (3.1), some preliminary assessments using this
technique are presented. Another potential solution to this problem is to target a
damage state that is less severe than collapse, for example structural yielding.
Targeting this limit state is less morally problematic and the level of acceptable risk
in this context could be defined using, for example, cost-benefit analysis based on
the cost of reducing the risk further (perhaps adopting an ALARP “as low as
reasonably practicable” philosophy). This choice also has other benefits, e.g. it is
generally easier to assess in numerical modelling when a structure yields rather than
when it collapses.
Once the targeted risk is chosen, unless great changes to the accelerations cur-
rently used for design are accepted by practicing engineers, the probability of
collapse at the design acceleration is automatically implied, as shown by Fig. 3 of
Douglas et al. (2013). Here there is a trilemma: any two out of the three input
parameters, design IM, X and Y, can be chosen independently but not all three.
To make this point clearer, we can use the work of Žižmond and Dolšek (2017).
In their article, they follow a risk-targeting approach for the derivation of the q
factor for use in force-based design using Eurocode 8 (CEN 2004b). In Eq. 4 of that
paper an analytical formula is presented for the median acceleration of the
risk-targeted lognormal fragility curve, l. Adapting it to our nomenclature and
recalling that our design acceleration ag implies a probability of collapse of X, we
have the following equations, which show that the various input parameters are
connected:
k2 b2
!1k
k0 e 2
l¼ ð1Þ
Y
where ag is the risk-targeted design acceleration and k0, k are parameters that define
the idealized hazard curve. The solution to this potential problem is to check all
three values are physically reasonable. This trilemma may be the reason for the
218 J. Douglas and A. Gkimprixis
apparently high target collapse probabilities (Y) used by Luco et al. (2007), as they
were forced to adopt them once X had been defined and they did not want the
design IM to change greatly from the previous code.
There is a need to derive fragility curves for a wide variety of code-designed
structures with different geometries and materials. Previous studies have adopted
generic fragility curves that scale constantly with the design acceleration so that the
iterative procedure used to converge to the targeted risk is simple. As shown by
Fig. 3 of Ulrich et al. (2014b), however, this desired feature appears not to be true
following current design practice, as they were not created with risk targeting in
mind. It is possible, however, to generate a suite of fragility curves for all potential
design accelerations and then to use the appropriate ones when iterating to find the
actual design acceleration for a location.
Two issues concerning the used fragility functions are: what value of b to use
and, indeed, whether the lognormal distribution should continue to be used at all?
Fragility curves that are used within the risk targeting calculations can imply very
strong buildings (e.g. Table 1 of Douglas et al. 2013) as well as non-negligible
chances of collapse for very low ground accelerations. This is a consequence of the
high values of b that need to be used to account for different types of structures of
varying geometries. Two solutions to this problem are: (a) adopt a different func-
tional form for the fragility curve that equals zero for low ground accelerations and
unity for very high accelerations, or (b) move from generic fragility curves for all
types and geometries to a curve covering only a small set of structures. The second
of these changes would mean a change of philosophy of design codes from being
associated with a single map giving the design accelerations for all structures to
potentially many maps giving design accelerations for different structural types and
geometries. This additional complexity, however, appears to be necessary if the
risk-targeting approach is to imply physically realistic buildings and levels of risk.
In this section, we use field observations of building damage from recent earth-
quakes and an empirical Fermi-type approach to estimate an upper bound on the
“acceptable” probability of collapse (Y) discussed above. In particular, the number
of RC buildings that collapsed due to earthquakes during recent time periods in
Italy and Greece are used to estimate the observed annual collapse rate.
The technique relies on the ratio of the total number of observed cases of a given
damage level in a given period of time to the total number of buildings that could
have been affected by earthquakes, which is then normalized to give an annual rate.
Estimating both of these totals for a given country requires: (a) post-earthquake
field mission reports that provide the number of cases of damage and (b) building
census data indicating how many structures were present when the earthquake
occurred.
Risk Targeting in Seismic Design Codes … 219
annual collapse rates of between 2 10−6 and 1 10−5, which are much lower
than some of the probabilities assumed in previous risk-targeting exercises (e.g.
Luco et al. 2007).
For the period 1978–2003 inclusive the following RC building collapses due to
earthquakes in Greece have been reported: Thessaloniki 1978 (4 buildings)
(CEQID), Gulf of Corinth 1981 (15) (Carydis et al. 1982), Kalamata 1986 (5),
Aigion 1995 (1), Athens 1999 (69) (CEQID) and Lefkada 2003 (0) (Karababa and
Pomonis 2011), giving a total of 94 total or partial collapses. Using information
from the Greek statistical service (http://www.statistics.gr) and the TABULA pro-
ject (www.building-typology.eu) for the total number of RC buildings in Greece,
and assuming that the building total remained constant, leads to estimates of the
collapse rate between 1 10−6 and 2 10−6, again much lower than some
assumed probabilities for the targeted risk.
The accuracy of these estimates could be improved by using a more complete
database both for the numerator (the damage) and the denominator (the number of
buildings at risk) and by extending the time period covered so that more potential
earthquake scenarios are considered. Nevertheless, such an analysis is constrained
by the long recurrence intervals of damaging earthquakes, which means that the
assessed rates will always be associated with large uncertainties. In addition,
variations due to building type (e.g. masonry versus reinforced concrete), location
(e.g. developed versus less-developed countries) and damage level (e.g. collapse
versus yield) could be studied as well.
The same procedure was repeated for the same two countries considering a
different damage state. A broader class was defined to represent ‘yield’, as con-
sidered by Ulrich et al. (2014a, b), which includes the total number of RC buildings
that are neither not-damaged nor collapsed (partially or in total). The same limi-
tations discussed before are also the case here. In addition, complete data were
available for a shorter period of 18 years for Greece (but still 30 years for Italy).
With these assumptions, the annual percentage of yield is found equal to 3 10−5
for Italy and 1 10−4 for Greece. Apart from providing useful constraints on the
inputs to risk-targeting this type of analysis could have other benefits as well.
Firstly, it would provide an observational-based assessment of earthquake risk that
can be compared to estimates from computer modelling, which require a large
number of inputs and can be opaque to decision makers. If there are large differ-
ences in these estimates it may indicate that the computer models require calibra-
tion. Secondly, these estimates will enable differences in earthquake risk levels
between countries to be judged more easily and used to understand what the impact
of reducing risk could be in terms of the average number of buildings that could
collapse per year, for example.
Risk Targeting in Seismic Design Codes … 221
4 Conclusions
Crowley et al. (2013) propose that ‘Risk-targeted seismic design actions’ should be
considered for future versions of the Eurocodes. A call echoed by Formichi et al.
(2016) in a report on the background and application of Eurocodes, published on
behalf of the European Commission. This report also proposes that changes made in
other international seismic design codes should be considered when updating the
Eurocodes. It is, therefore, clear that application of the risk-targeting approach is
being seriously discussed in Europe. In consequence, additional research effort to
this end would provide valuable input for the development of risk-targeted design
maps for new buildings in Europe.
In this brief article, we have highlighted the critical issues that we believe need
to be solved before the risk-targeting approach for the development of seismic
design codes can be employed in practice. Some of these (e.g. development of
appropriate fragility functions) solely require engineering calculations but others
(e.g. choice of the acceptable level of risk) need input from other domains,
including decision makers. In the coming years we plan to tackle these issues,
particularly with respect to future Eurocodes.
References
Allen TI, Adams J, Halchuk S (2015) The seismic hazard model for Canada: past, present and
future. In: Proceedings of the tenth pacific conference on earthquake engineering building an
earthquake-resilient pacific, Sydney, Australia (Paper number 100)
ASCE (2005) Seismic design criteria for structures, systems and components in nuclear facilities.
American Society of Civil Engineers, ASCE Standard 43-05
ASCE (2010) Minimum design loads for buildings and other structures. ASCE Standard 7-10.
American Society of Civil Engineers, Reston, VA
Braverman JI, Xu J, Ellingwood BR, Costantino CJ, Morante RJ, Hofmayer CH (2007) Evaluation
of the seismic design criteria in ASCE/SEI standard 43-05 for application to nuclear power
plants. Agencywide Documents Access and Management System (ADAMS)—USNRC
Carydis PG, Tilford NR, Brandow GE, Jirsa JO (1982) The central Greece earthquakes of
February–March 1981. Report No. CETS-CND-018 of the Earthquake Engineering Research
Institute, Berkeley, CA, National Academy Press, Washington, DC
CEN (2004a) EN 1992-1-1:2004 Eurocode 2: design of concrete structures—Part 1-1: general
rules and rules for buildings. European Committee for Standardization, Brussels
CEN (2004b) EN 1998-1:2004 Eurocode 8: design of structures for earthquake resistance—Part 1:
general rules, seismic actions and rules for buildings. European Committee for Standardization,
Brussels
222 J. Douglas and A. Gkimprixis
Silva V, Crowley H, Bazzurro P (2016) Exploring risk-targeted hazard maps for Europe. Earthq
Spectra 32(2):1165–1186
SNI (2012) Tata cara perencanaan ketahanan gempa untuk struktur bangunan gedung dan non
gedung. Badan Standardisasi Nasional, ICS 91.120.25; 91.080.01. Report 1726:2012
(In Indonesian)
Tsang HH, Wenzel F (2016) Setting structural safety requirement for controlling earthquake
mortality risk. Saf Sci 86:174–183
Tsang HH, Lumantarna E, Lam NTK, Wilson JL, Gad E (2017) Annualised collapse risk of
soft-storey building with precast RC columns in Australia. Advancements and Challenges,
Mechanics of Structures and Materials, pp 1681–1686
Ulrich T, Negulescu C, Douglas J (2014a) Fragility curves for risk-targeted seismic design maps.
Bull Earthq Eng 12(4):1479–1491
Ulrich T, Douglas J, Negulescu C (2014b) Seismic risk maps for Eurocode-8 designed buildings.
In: Proceedings of the second European conference on earthquake engineering and seismology
United States Nuclear Regulatory Commission (2007) A performance-based approach to define the
site specific earthquake ground motion. Tech report 1:208
Vacareanu R, Pavel F, Craciun I, Coliba V, Arion C, Aldea A, Neagu C (2017) Risk-targeted maps
for Romania. J Seismolog. https://doi.org/10.1007/s10950-017-9713-x
Žižmond J, Dolšek M (2017) The formulation of risk-targeted behaviour factor and its application
to reinforced concrete buildings. In: Proceedings of the 16th world conference on earthquake
engineering, Santiago, Chile, 9–13 January (Paper no. 1659)