Earthquake Source Model Database
FSP File Format
  • The *fsp (finite-source parameter) provides a comprehensive description of modeling/inversion parameters and the finite-source rupture parameters, starting with a header block, then listing source parameters on the fault plane, potentially segment by segment (in case of multi-segment events).
  • The length of the header block is variable, depending on the velocity-density structure, but can be inferred from the information given in the header.
  • The rupture model is given in a column-format, each column presenting a different quantity; each row represents one point (subfault) on the rupture plane, arranged starting from the top left corner of the fault plane, and progressing first along-strike, then down-dip, to the bottom-right of the faults.
  • The number of columns (source parameters) may be variable for different rupture models, depending on what source parameters are available. The minimum number of columns is SIX (Lat, Lon, X, Y, Z, SLIP), while the maximum number is unlimited. Additional columns may comprise RAKE, RISE time, TRUP (rupture onset time), or slip (and perhaps rake) in subsequent time windows (indicated as TW1 and rkTW1, TW2 rkTW2, etc ....
  • If the source model has been obtained using various fault segments, the data are given segment-by-segment in the same manner as just described, where each segment is preceeded by a 9-line header block with segment specific meta-data.
  • The header comprises basic earthquake source parameters and meta-data on he modeling /inversion approach; besides the parameters described above for the evTAG.slp format, some additional entries are given:
  • Ntw, Nsg: number of time-windows, and number of fault segments, respectively
  • LEN, SHF: length (in sec) of slip-function, and the overlap (in sec) of successive time-windows
  • SVF: type of slip-velocity function used in the inversion
  • In the Data-section, the number of stations are listed for the different types of data used in the inversion; PHmx is a rough estimate of the maximum angle between any two neighboring stations (ie. PHmx > 180 deg means the event was outside the observational array). Rmin estimates the distance of the closest station to the nearest point on the fault.
    We make a special case for InSAR data wherein NoS refers to "number of InSAR data points", and PHImx corresponds to "number of independent interferograms".
  • Data on the velocity-density structure at the source is preceeded by the corresponding number of layers.
  • Nsbfs: number of sub-faults, either for the single-plane rupture or given for each segment separately.
The format for 'single segment' model and 'multi-segment' model differs slightly.

An example for single segment model:

% ---------------------------------- FINITE-SOURCE RUPTURE MODEL -------------------------------- %
% Event : Yamaguchi (Japan) 06/25/1997 [Ide (1999) ]
% EventTAG: s1997YAMAGU01IDEx
%
% Loc : LAT = 34.441 LON = 131.676 DEP = 7.5
% Size : LEN = 16.0 km WID = 12.0 km Mw = 5.81 Mo = 5.86e+17 Nm
% Mech : STRK = 235.0 DIP = 86.0 RAKE = 182.0 Htop = 0.1 km
% Rupt : HypX = 6.00 km HypZ = 7.00 km avTr = 4.2 s avVr = 3.0 km/s
%
% ---------------------------------- inversion-related parameters --------------------------------
%
% Invs : Nx = 17 Nz = 13 Fmin = 0.10 Hz Fmax = 2.00 Hz
% Invs : Dx = 1.00 km Dz = 1.00 km
% Invs : Ntw = 13 Nsg = 1 (# of time-windows,# of fault segments)
% Invs : LEN = 0.6 s SHF = 0.3 s (time-window length and time-shift)
% SVF : triang (type of slip-velocity function used)
%
% Data : SGM TELE TRIL LEVEL GPS INSAR SURF OTHER
% Data : 0 10 0 0 0 0 0 0
% PHImx: 0 79 0 0 0 0 0 0
% Rmin : 0 79 0 0 0 0 0 0
%
% --------------------------------------------------------------------------------------------------
%
% VELOCITY-DENSITY STRUCTURE
% No. of layers = 5
%
% DEPTH P-VEL S-VEL DENS QP QS
% [km] [km/s] [km/s] [g/cm^3]
% 0.00 5.60 3.20 2.50 400 200
% 3.00 6.00 3.47 2.70 600 300
% 16.00 6.60 3.82 3.00 800 400
% 30.00 7.80 4.50 3.20 1000 500
% 70.00 8.00 4.62 3.25 1000 500
%
% --------------------------------------------------------------------------------------------------
% 23-Jul-2007 PMM (mai@sed.ethz.ch)
% --------------------------------------------------------------------------------------------------
%
% SOURCE MODEL PARAMETERS
% Nsbfs = 221 subfaults
% X,Y,Z coordinates in km; SLIP in m
% if applicable: RAKE in deg, RISE in s, TRUP in s, slip in each TW in m
%
% Coordinates are given for top-center of each subfault or segment: |'|
% Origin of local coordinate system at epicenter: X (EW) = 0, Y (NS) = 0
% LAT LON X==EW Y==NS Z SLIP TRUP
% --------------------------------------------------------------------------------------------------
34.4710 131.7372 5.6046 3.3283 0.1000 0.0000 3.0730
34.4658 131.7282 4.7854 2.7547 0.1000 0.0000 2.8670
34.4606 131.7193 3.9663 2.1811 0.1000 0.0000 2.6870
34.4555 131.7103 3.1471 1.6075 0.1000 0.0000 2.5390
........ many rows omitted here ....

An example for multiple segment model:

% ---------------------------- FINITE-SOURCE RUPTURE MODEL --------------------------------
%
% Event : Kobe (Japan) 01/17/1995 [Sekiguchi etal (2000)]
% EventTAG: s1995KOBEJ1seki
%
% Loc : LAT = 34.598 LON = 135.044 DEP = 16.37
% Size : LEN = 63.55 km WID = 20.50 km Mw = 6.99 Mo = 3.44e+019 Nm
% Mech : STRK = 45 DIP = 85 RAKE = 180 Htop = 1.00 km
% Rupt : avTr = 3.3 s avVr = 3.1 km/s
%
% ---------------------------- inversion-related parameters --------------------------------
%
% Invs : Nx = 31 Nz = 10 Fmin = 0.10 Hz Fmax = 1.0 Hz
% Invs : Dx = 2.05 km Dz = 2.05 km
% Invs : Ntw = 6 Nsg = 5 (# of time-windows,# of fault segments)
% Invs : LEN = 0.80 s SHF = 0.30 s (time-window length and time-shift)
% SVF : tanhyp (type of slip-velocity function used)
%
% Data : SGM TELE GPS tril level
% NoS : 25 0 0 0 0
% PHImx: 110 0 0 0 0
% Rmin : 0.2 0.0 0.0 0.0 0.0
%
% -----------------------------------------------------------------------------------------
%
% VELOCITY-DENSITY STRUCTURE
% No. of layers = 4
%
% DEPTH P-VEL S-VEL DENS QP QS
% [km] [km/s] [km/s] [g/cm^3]
% 0.00 5.50 3.20 2.60 9999 9999
% 5.10 6.00 3.46 2.70 9999 9999
% 18.00 6.70 3.87 2.80 9999 9999
% 34.50 7.50 4.33 3.00 9999 9999
%
% -----------------------------------------------------------------------------------------
% 18-Jul-2007, created by P.Martin Mai (mai@sed.ethz.ch)
% -----------------------------------------------------------------------------------------
%
% SOURCE MODEL PARAMETERS
% X,Y,Z coordinates in km; SLIP in m
% if applicable: RAKE in deg, RISE in s, TRUP in s, slip in each TW in m
%
% Coordinates are given for top-center of each subfault or segment: |'|
% Origin of local coordinate system at epicenter: X (NS) = 0, Y (EW) = 0
% ------------------------------------------------------------------------------------------
% --------------------------- MULTISEGMENT MODEL -------------------------------------------
% ------------------------------------------------------------------------------------------
% SEGMENT # 1: STRIKE = 45 deg DIP = 78.0 deg
% LEN = 20.50 km WID = 20.50 km
% depth to top: Z2top = 1.00 km
% coordinates of top-center:
% LAT = 34.5530, LON = 134.9402
% hypocenter on SEG # 1 : along-strike (X) = 20.50, down-dip (Z) = 15.37
% Nsbfs = 100 subfaults
% -------------------------------------------------------------------------------------------
% LAT LON X==EW Y==NS Z SLIP RAKE TW1 rkTW1 TW2 rkTW2 ...
% -------------------------------------------------------------------------------------------
34.4943 134.8689 -11.5234 -16.0192 1.0000 0.224 166.504 0.057 108.435 0.000 90.000 ...
34.5073 134.8847 -10.0738 -14.5696 1.0000 0.876 221.820 0.205 222.797 0.148 250.201 ...
34.5204 134.9006 -8.6242 -13.1200 1.0000 0.647 241.314 0.005 90.000 0.128 220.101 ...

........ many rows omitted here ....

% -------------------------------------------------------------------------------------------
% SEGMENT # 2: STRIKE = 50 deg DIP = 90.0 deg
% LEN = 14.35 km WID = 20.50 km
% depth to top: Z2top = 1.00 km
% coordinates of top-center:
% LAT = 34.6395, LON = 135.1041
% hypocenter on SEG # 1 : along-strike (X) = 20.50, down-dip (Z) = 15.37
% Nsbfs = 70 subfaults
% -------------------------------------------------------------------------------------------
% LAT LON X==EW Y==NS Z SLIP RAKE TW1 rkTW1 TW2 rkTW2 ...

% -------------------------------------------------------------------------------------------
34.6039 135.0526 0.6584 0.7858 1.0000 0.891 118.610 0.206 121.827 0.122 110.225 ...

........ many rows omitted here ...

©2012 Paul Martin Mai ♦ Disclaimer/Acknowlegments: see About page.