FiSH: MATLAB Tools to Turn Fault Data into Seismic‐Hazard Models

Faults have been increasingly integrated into seismic‐hazard assessments. We have developed a package of MATLAB ( , last accessed January 2016) tools (called FiSH ), designed to help seismic‐hazard modelers analyze fault data. These tools enable the derivation of expected earthquake rates, given common fault data, and allow researchers to test the consistency between the magnitude–frequency distributions (MFDs) assigned to a fault and some available observations. The basic assumption of FiSH is that the geometric and kinematic features of a fault are the expression of its seismogenic potential. Three tools have been designed to integrate the variable levels of information available: (1) the first tool allows users to convert fault geometry and slip rates into a global budget of the seismic moment released in a given time frame, taking uncertainties into account; (2) the second tool computes the recurrence parameters and associated uncertainties from historical and/or paleoseismological data; and (3) the third tool outputs time‐independent or time‐dependent earthquake rates for different MFD models. We present a test case to illustrate the capabilities of FiSH , on the Paganica normal fault in central Italy that ruptured during the 2009 L’Aquila earthquake sequence (mainshock M w 6.3). The source codes are open, and we encourage users to handle the scripts, communicate with us regarding bugs, and/or suggest further improvements. Our intent is to distribute these tools in order to help researchers to pinpoint potential inconsistencies and obtain reliable fault‐based seismic‐hazard evaluations.

[1]  C. Allin Cornell,et al.  A hybrid recurrence model and its implication on seismic hazard results , 1995 .

[2]  Bruno Pace,et al.  Layered Seismogenic Source Model and Probabilistic Seismic-Hazard Analyses in Central Italy , 2006 .

[3]  B. E. Shaw,et al.  Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3)—The Time‐Independent Model , 2014 .

[4]  D. Wells,et al.  New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement , 1994, Bulletin of the Seismological Society of America.

[5]  B. Gutenberg,et al.  Frequency of Earthquakes in California , 1944, Nature.

[6]  David D. Jackson,et al.  A mutually consistent seismic-hazard source model for southern California , 1999 .

[7]  C. Scholz.  The Mechanics of Earthquakes and Faulting , 1990 .

[8]  A. Michetti,et al.  Fault scaling relationships, deformation rates and seismic hazards: an example from the Lazio-Abruzzo Apennines, central Italy , 2004 .

[9]  B. Pace,et al.  Sensitivity analysis for seismic source characteristics to probabilistic seismic hazard assessment in central Apennines (Abruzzo area). , 2002 .

[10]  Carlo Meletti,et al.  CRISIS2008: A Flexible Tool to Perform Probabilistic Seismic Hazard Assessment , 2013 .

[11]  D. Owen,et al.  Handbook of statistical distributions , 1978 .

[12]  David M. Perkins,et al.  SEISRISK III; a computer program for seismic hazard estimation , 1987 .

[13]  F. Visini,et al.  Fault-Based Earthquake Rupture Forecast in Central Italy: Remarks after the L'Aquila Mw 6.3 Event , 2011 .

[14]  William L. Ellsworth,et al.  A Brownian Model for Recurrent Earthquakes , 2002 .

[15]  Carlo Alberto Brunori,et al.  Evidence for surface faulting events along the Paganica fault prior to the 6 April 2009 L'Aquila earthquake (central Italy) , 2011 .

[16]  Timothy E. Dawson,et al.  Long‐Term Time‐Dependent Probabilities for the Third Uniform California Earthquake Rupture Forecast (UCERF3) , 2015 .

[17]  T. Jordan,et al.  OpenSHA: A Developing Community-modeling Environment for Seismic Hazard Analysis , 2003 .

[18]  G. Lavecchia,et al.  Defining a model of 3D seismogenic sources for Seismic Hazard Assessment applications: The case of central Apennines (Italy) , 2004 .

[19]  G. Weatherill,et al.  OpenQuake Engine: An Open Hazard (and Risk) Software for the Global Earthquake Model , 2014 .

[20]  Tom Parsons.  Monte Carlo method for determining earthquake recurrence parameters from short paleoseismic catalogs: Example calculations for California , 2008 .

[21]  Alessandro Maria Michetti,et al.  Sintesi delle conoscenze sulle faglie attive in Italia Centrale: parametrizzazione ai fini della caratterizzazione della pericolosità sismica , 2000 .

[22]  Stefan Wiemer,et al.  A Software Package to Analyze Seismicity: ZMAP , 2001 .

[23]  Pilar Villamor,et al.  National Seismic Hazard Model for New Zealand: 2010 Update , 2012 .

[24]  Y. Kagan.  Seismic moment distribution revisited: I. Statistical results , 2002 .

[25]  G. Lavecchia,et al.  Coseismic ground deformation of the 6 April 2009 L'Aquila earthquake (central Italy, Mw6.3) , 2010 .

[26]  M. Leonard.  Earthquake Fault Scaling: Self‐Consistent Relating of Rupture Length, Width, Average Displacement, and Moment Release , 2010 .

[27]  LASSCI2009.2: layered earthquake rupture forecast model for central Italy, submitted to the CSEP project , 2010 .

[28]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[29]  F. Cavallini,et al.  Error propagation in time-dependent probability of occurrence for characteristic earthquakes in Italy , 2010 .

[30]  Christopher Bronk Ramsey,et al.  Recent and Planned Developments of the Program OxCal , 2013 .