CRISIS2008: A Flexible Tool to Perform Probabilistic Seismic Hazard Assessment

In the frame of the Italian research project INGV-DPC S2 (http://nuovoprogettoesse2.stru.polimi.it/), funded by the Dipartimento della Protezione Civile (DPC; National Civil Protection Department) within the agreement 2007-2009, a tool for probabilistic seismic hazard assessment (PSHA) was developed. The main goal of the project was to provide a flexible computational tool for PSHA; the requirements considered essential for the success of the project included:
• ability to handle both stationary and non-stationary earthquake time-occurrence models;
• ability to use ground-motion prediction models that are not parametric equations but probabilistic "footprints" of the intensities generated by earthquakes of known magnitude and focal characteristics. Usually, these footprints are results of ground motion simulations.

Some commonly used programs (e.g., FRISK, by McGuire, 1978; SEISRISK III, by Bender and Perkins, 1987) and more recent and state-of-the-art tools (e.g. OpenSHA, by Field et al., 2003, http://www.opensha.org; OpenQuake, http://openquake.org) for PSHA were analyzed. It was decided to focus on CRISIS2007, which was already a mature and well known application (e.g., Kalyan Kumar and Dodagoudar, 2011; Teraphan et al., 2011; D’Amico et al., 2012; see also http://ecapra.org/CRISIS-2007), but also suitable for additional development and evolution since its source code is freely available on request. The computational tool resulted in an extensive redesign and renovation of the previous CRISIS2007 version.

CRISIS is a computer program for PSHA, originally developed in the late 1980's using Fortran as programming language (Ordaz, 1991). In this format, still without a graphical user interface (GUI), it was distributed as part of SEISAN tools (Ottemöller et al., 2011).
Ten years later, a GUI was constructed, generating what was called CRISIS99 (Ordaz, 1999). In this version, all the graphic features were written in Visual Basic, but the computation engine remained a Fortran dynamic link library. The reason for the use of mixed-language programming was that computations in Visual Basic were extremely slow.
Around 2007 the program was upgraded, in view of the advantages offered by the object-oriented technologies. An object-oriented programming language was required and the natural choice was Visual Basic.Net. In the new version (called CRISIS2007), both the GUI and the computation engine were written in the same language.
Finally, in the frame of the mentioned S2 project, starting from 2008, the program was split into two logical layers: core (CRISIS Core Library) and presentation (CRISIS2008). In addition, a new presentation layer was developed for accessing the same functionalities via Web (CRISISWeb).
It is worth noting that CRISIS has been mainly written by people that are, at the same time, PSHA practitioners. Therefore, the development loop has been relatively short, and most of the modifications and improvements have been made to satisfy the needs of the developers themselves.

Italian Presidenza del Consiglio dei Ministri, Dipartimento della Protezione Civile (DPC).

Abrahamson, N. A., and W. J Silva (1997). Empirical response spectral attenuation relations for shallow crustal earthquakes, Seismological Research Letters 68, 94–127.

Akkar, S., and J. J. Bommer (2010). Empirical equations for the prediction of PGA, PGV, and spectral accelerations in Europe, the Mediterranean region, and the Middle East, Seismological Research Letters 81, 195-206.

Arroyo, D., D. García, M. Ordaz, M. A. Mora, and S. K. Singh (2010). Strong ground-motion relations for Mexican interplate earthquakes, Journal of Seismology 14, 769-785.

Atkinson, G. M., and D. M. Boore (2003). Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions, Bulletin of the Seismological Society of America 93, 1703-1729.

Bazzurro, P., and C. A. Cornell (1999). Disaggregation of seismic hazard, Bulletin of the Seismological Society of America 89, 501-520.

Bender, B., and D. M. Perkins (1987). SEISRISK III: a computer program for seismic hazard estimation, U.S. Geological Survey Bulletin 1772, 48 pp.

Boore, D. M., and G. M. Atkinson (2008). Ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods between 0.01 s and 10.0 s, Earthquake Spectra 24, 99-138.

Campbell, K. W., and Y. Bozorgnia (2003). Updated near-source ground-motion (attenuation) relations for the horizontal and vertical components of peak ground acceleration and acceleration response spectra, Bulletin of the Seismological Society of America 93, 314–331.

Cauzzi, C., and E. Faccioli (2008). Broadband (0.05 to 20 s) prediction of displacement response spectra based on worldwide digital records, Journal of Seismology 12, 453–475.

Cornell, C.A. (1968). Engineering seismic risk analysis, Bulletin of the Seismological Society of America 58, 1583-1606.

D’Amico, V., C. Meletti, and F. Martinelli (2012). Probabilistic seismic hazard assessment in the high-risk area of south-eastern Sicily (Italy), Bollettino di Geofisica Teorica ed Applicata 53, 19-36.

Esteva, L. (1967). Criteria for the construction of spectra for seismic design, 3rd Pan-American Symposium of Structures, Caracas, Venezuela, 3–8 July (in Spanish).

Field, E.H., T. H. Jordan, and C. A. Cornell (2003). OpenSHA: A developing community-modeling environment for seismic hazard analysis, Seismological Research Letters 74, 406-419.

García, D., S. K. Singh, M. Herráiz, M. Ordaz, and J. F. Pacheco (2005). Inslab earthquakes of Central Mexico: peak ground-motion parameters and response spectra, Bulletin of the Seismological Society of America 95, 2272–2282.

Kalyan Kumar, G., and G. R. Dodagoudar (2011). Seismic input motion for Kanchipuram, South India, International Journal of Earth Sciences and Engineering 4, 6 SPL, 189-192.

Kiremidjian, A. S., and T. Anagnos (1984). Stochastic slip-predictable model for earthquake occurrences, Bulletin of the Seismological Society of America 74, 739-755.

Marzocchi, W., and A. M. Lombardi (2008). A double branching model for earthquake occurrence, Journal of Geophysical Research 113, B08317, doi:10.1029/2007JB005472.

Matthews, M. V., W. L. Ellsworth, and P. A. Reasenberg (2002). A Brownian model for recurrent earthquakes, Bulletin of the Seismological Society of America 92, 2233-2250.

McGuire, R. K. (1978). FRISK: A Computer Program for Seismic Risk Analysis Using Faults as Earthquake Sources, U.S. Geological Survey Open-File Rept. 76-67, 90 pp.

Meletti, C., V. D’Amico, and F. Martinelli (2009). Module for ER model based on Poisson applied to ZS9. Progetto INGV-DPC S2, Deliverable D2.1, http://nuovoprogettoesse2.stru.polimi.it/Deliverables.html

Ordaz, M. (1991). CRISIS. Brief description of program CRISIS. Internal report, Institute of Solid Earth Physics, University of Bergen, Norway. 16 pp.

Ordaz, M. (1999). User's manual for program CRISIS-99. Technical report, Universidad Nacional Autonoma de Mexico, Mexico City.

Ottemöller, L., P. Voss, and J. Havskov (2011). Seisan earthquake analysis software for Windows, Solaris, Linux and MacOsX. http://seis.geus.net/software/seisan/seisan.html

Pasolini, C., D. Albarello, P. Gasperini, V. D’Amico, and B. Lolli (2008). The attenuation of seismic intensity in Italy part II: modeling and validation, Bulletin of the Seismological Society of America 98, 692-708.

Sabetta, F., and A. Pugliese (1996). Estimation of response spectra and simulation of nonstationary earthquake ground motions, Bulletin of the Seismological Society of America 86, 337-352.

Spudich, P., W. B. Joyner, A. G. Lindh, D. M. Boore, B. M. Margaris, and J. B. Fletcher (1999). SEA99: A revised ground motion prediction relation for use in extensional tectonic regimes, Bulletin of the Seismological Society of America 89, 1156-1170.

Stucchi, M., C. Meletti, V. Montaldo, H. Crowley, G. M. Calvi, and E. Boschi (2011). Seismic hazard assessment (2003-2009) for the Italian building code, Bulletin of the Seismological Society of America 101, 1885-1911.

Stupazzini, M., R. Paolucci, and H. Igel (2009). Near-fault earthquake ground motion simulation in the Grenoble Valley by a high-performance spectral element code, Bulletin of the Seismological Society of America 99, 286-301.

Teraphan, O., J. Douglas, R. Sigbjörnsson, and C. G. Lai (2011). Assessment of Ground Motion Variability and Its Effects on Seismic Hazard Analysis: A Case Study for Iceland, Bulletin of Earthquake Engineering 4, 931-953.

Thomas, P., Wong, I. and N. Abrahamson (2010). Verification of Probabilistic Seismic Hazard Analysis Computer Programs, PEER Report 2010/106, Pacific Earthquake Engineering Research Center, College of Engineering, University of California, Berkeley.

Youngs, R. R., S. J. Chiou, W. J. Silva, and J. R. Humphrey (1997). Strong ground motion attenuation relationships for subduction zone earthquakes, Seismological Research Letters 68, 58-73.