Skip to content

Data formats > Catalogs > Seismic catalog

Catalog - description

Catalog – General Description

The catalogue is a variable in the Matlab format file and it is kept in a file MAT. The structure is array with named fields that can contain data of various types and sizes. In the file there is only one variable, the file name and variable name are optional.

The variable describing the catalogue is a vector of structures, consisting of fields:

  • field – name of field in the catalogue (text value);
  • type – type of field in the catalogue and way of showing the field (numeric value); the numbers description is shown below
  • val – column array of For the text the column is an array type cell with text fields. For the remaining value the column is a numeric column.
  • unit – description of unit for individual data (text value).
  • description – short description of the parameter (text value).
  • fieldType – semantic meaning of the field. When some field values are similar/related then fieldType name is entered and for another case is entered.

Table The general parameters in catalogue MAT format

fieldtypevalunitdescriptionfieldTypeCommentsData format
ID3data vectorEvent IDrequired fieldtext
Time5data vectorEvent origin timerequired field, Matlab serial numerical timedouble
Lat14,15,24,25data vectordegLatitudedeg – North positivedouble
Long14,15,24,25,34,35data vectordegLongitudedeg – East positivedouble
Depth11-13data vectorkmHypocenter depth measured from the ground levelHypocenter depth counted downwards from surface (positive values below surface)double
Elevation13data vectorkmHypocenter elevation measured over the sea levelHypocenter elevation counted upwards from sea level (positive values above sea level, negative values below sea level)double
X10data vectorm, kmX coordinateOriginal coordinates if other than geographical. Description of coordinates in the metadatadouble
Y10data vectorm, kmY coordinatedouble
Z10data vectorm, kmZ coordinatedouble
EPI_err10data vectorMEpicentral errordouble
Depth_err10data vectormDepth errordouble
Nl2data vectorNo of stations used in the localisationdouble
M0222data vectorNmScalar momentdouble
Mw4data vectorMoment magnitude'Magnitude'double 0.1[1]
ML4data vectorLocal magnitude'Magnitude'double 0.1
Ns_decomp2data vectorNo of stations used in MT inversiondouble
DecompMethod3data vectorcharMethod used to decompose moment tensortext
MTrr222data vectorNmFull solution: Moment tensor rr component (r - up)double
MTss222data vectorNmFull solution: Moment tensor ss component (s - South)double
MTee222data vectorNmFull solution: Moment tensor ee component (e - East)double
MTrs222data vectorNmFull solution: Moment tensor rs componentdouble
MTre222data vectorNmFull solution: Moment tensor re componentdouble
MTse222data vectorNmFull solution: Moment tensor se componentdouble
MT_err222data vectorNmFull solution: Moment tensor errordouble
ISO120data vector%Isotropic MT component% - positive or negativedouble
CLVD120data vector%CLVD component% - positive or negativedouble
DC20data vector%Double-Couple component% - only positivedouble
StrikeA30data vectordegStrike of nodal plane AThe values range from 0 to 360double
DipA20data vectordegDip of nodal plane AThe values range from 0 to 90double
RakeA130data vectordegRake of nodal plane AThe values range from -180 to 180double
SlopeA20data vectordegInclination for nodal plane AThe values range from 0 to 90double
StrikeB30data vectordegStrike of nodal plane BThe values range from 0 to 360double
DipB20data vectordegDip of nodal plane BThe values range from 0 to 90double
RakeB130data vectordegRake of nodal plane BThe values range from ‑180 to 180double
SlopeB20data vectordegInclination for nodal plane BThe values range from 0 to 90double
Strike_err10data vectordegStrike errordouble
Dip_err10data vectordegDip errordouble
Rake_err10data vectordegRake errordouble
Slope_err10data vectordegInclination errordouble
Plunge_T10data vectordegPlunge of T-axisThe values range from 0 to 360double
PlungeT_err10data vectordegT-axis plunge errordouble
Trend_T10data vectordegTrend of T-axisThe values range from 0 to 90double
TrendT_err10data vectordegT-axis trend errordouble
Plunge_P10data vectordegPlunge of P-axisThe values range from 0 to 360double
PlungeP_err10data vectordegP-axis plunge errordouble
Trend_P10data vectordegTrend of P-axisThe values range from 0 to 90double
TrendP_err10data vectordegP-axis trend errordouble
DCrr222data vectorNmDouble-Couple solution: Moment tensor rr component (r - up)double
DCss222data vectorNmDouble-Couple solution: Moment tensor ss component (s - South)double
DCee222data vectorNmDouble-Couple solution: Moment tensor ee component (e - East)double
DCrs222data vectorNmDouble-Couple solution: Moment tensor rs componentdouble
DCre222data vectorNmDouble-Couple solution: Moment tensor re componentdouble
DCse222data vectorNmDouble-Couple solution: Moment tensor se componentdouble
DC_err222data vectorNmDouble-Couple solution: Moment tensor errordouble
DCStrikeA30data vectordegDouble-Couple solution: Strike of nodal plane AThe values range from 0 to 360double
DCDipA20data vectordegDouble-Couple solution: Dip of nodal plane AThe values range from 0 to 90double
DCRakeA130data vectordegDouble-Couple solution: Rake of nodal plane AThe values range from -180 to 180double
DCStrikeB20data vectordegDouble-Couple solution: Strike of nodal plane BThe values range from 0 to 90double
DCDipB30data vectordegDouble-Couple solution: Dip of nodal plane BThe values range from 0 to 360double
DCRakeB20data vectordegDouble-Couple solution: Rake of nodal plane BThe values range from 0 to 90double
DCStrike_err10data vectordegDouble-Couple solution: Strike errordouble
DCDip_err10data vectordegDouble-Couple solution: Dip errordouble
DCRake_err10data vectordegDouble-Couple solution: Rake errordouble
DCPlunge_T10data vectordegDouble-Couple solution: Plunge of T-axisThe values range from 0 to 90double
DCPlungeT_err10data vectordegDouble-Couple solution: T-axis plunge errordouble
DCTrend_T10data vectordegDouble-Couple solution: Trend of T-axisThe values range from 0 to 360double
DCTrendT_err10data vectordegDouble-Couple solution: T-axis trend errordouble
DCPlunge_P10data vectordegDouble-Couple solution: Plunge of P-axisThe values range from 0 to 90double
DCPlungeP_err10data vectordegDouble-Couple solution: P-axis plunge errordouble
DCTrend_P10data vectordegDouble-Couple solution: Trend of P-axisThe values range from 0 to 360double
DCTrendP_err10data vectordegDouble-Couple solution: P-axis trend errordouble
TNrr222data vectorNmTN solution: Moment tensor rr component (r - up)double
TNss222data vectorNmTN solution: Moment tensor ss component (s - South)double
TNee222data vectorNmTN solution: Moment tensor ee component (e - East)double
TNrs222data vectorNmTN solution: Moment tensor rs componentdouble
TNre222data vectorNmTN solution: Moment tensor re componentdouble
TNse222data vectorNmTN solution: Moment tensor se componentdouble
TN_err222data vectorNmTN solution: Moment tensor errordouble
TNStrikeA30data vectordegTN solution: Strike of nodal plane AThe value range from 0 to 360double
TNDipA20data vectordegTN solution: Dip of nodal plane AThe value range from 0 to 90double
TNRakeA130data vectordegTN solution: Rake of nodal plane AThe value range from -180 to 180double
TNStrikeB20data vectordegTN solution: Strike of nodal plane BThe value range from 0 to 90double
TNDipB30data vectordegTN solution: Dip of nodal plane BThe value range from 0 to 360double
TNRakeB20data vectordegTN solution: Rake of nodal plane BThe value range from 0 to 90double
TNStrike_err10data vectordegTN solution: Strike errordouble
TNDip_err10data vectordegTN solution: Dip errordouble
TNRake_err10data vectordegTN solution: Rake errordouble
TNPlunge_T20data vectordegTN solution: Plunge of T-axisThe value range from 0 to 90double
TNPlungeT_err10data vectordegTN solution: T-axis plunge errordouble
TNTrend_T30data vectordegTN solution: Trend of T-axisThe value range from 0 to 360double
TNTrendT_err10data vectordegTN solution: T-axis trend errordouble
TNPlunge_P20data vectordegTN solution: Plunge of P-axisThe value range from 0 to 90double
TNPlungeP_err10data vectordegTN solution: P-axis plunge errordouble
TNTrend_P30data vectordegTN solution: Trend of P-axisThe value range from 0 to 360double
TNTrendP_err10data vectordegTN solution: P-axis trend errordouble
NsP2data vectorNo of stations used in the P-wave spectral analysisdouble
E222data vectorJTotal seismic energydouble
E_err222data vectorJTotal seismic energy errordouble
Ep222data vectorJP-wave energydouble
Ep_err222data vectorJP-wave energy errordouble
fp12data vectorHzP-wave corner frequencydouble
fp_err12data vectorHzP-wave corner frequency errordouble
rad_eff_P12data vectorRadiation efficiency Pdouble
Qp10data vectorQuality factor Pwavesdouble
NsS2data vectorNo of stations used in the S-wave spectral analysisdouble
Es222data vectorJS-wave energydouble
Es_err222data vectorJS-wave energy error Jdouble
fs12data vectorHzS-wave corner frequency Hzdouble
fs_err12data vectorHzS-wave corner frequency error Hzdouble
Qs10data vectorQuality factor Swavesdouble
rad_eff_S12data vectorRadiation efficiency Sdouble
R10data vectormSource radiusdouble
R_err10data vectormSource radius errordouble
R_model3data vectorcharSource radius model used (Brune, Madariaga, Sato&Hirasawa)text
rad_eff12data vectorRadiation efficiencydouble
sigma_a13data vectorMPaApparent stressdouble
delta_sigma13data vectorMPaStatic stress dropdouble
sigma_d13data vectorMPaDynamic stress dropdouble
sigma_rms13data vectorMPaRMS dynamic stress dropdouble
vr10data vectorm/sRupture velocitydouble
vr_model3data vectorcharRupture velocity model (unilateral etc.)text
SW_eff12data vectorSavage-Wood efficiencydouble
u12data vectormFault slipdouble
RMS_time_residual13data vectorsecTime residual after location in a 3D modeldouble
Hypocenter_quality_index2data vectorAverage number of station, azimuthal coverage, stability of location against noiseApart from the regular parameters used to select best earthquake location, including min RMS values, max recording stations,  max phase numbers and gap  values, we further used for lacq catalogue the following criteria as derived from 3D inversion of hypocenters using REL3D code (e.g. Guyoton et al GRL1992, Roecker JGR1982,  Boyer 1996).

3 quality levels are used for the Lacq seismicity hypocenters
quality level 1:
-conditioning value <100
Specifically , the conditioning value threshold aims to remove ill-conditioned matrix patterns.
- max location error (e.g. Tarantola Valette 1982) as estimated for 3D velocity model solutions <1 km
- last iteration step for convergence before location completed < 0.5 km
quality level 2:
- conditioning value <100
quality level 3:
- all selected events (in the lacq case study it correspond to Nmin=6 stations and  a residual variance <0.2 s)

Reference:
Boyer, E. (1996). Sismicité induite et production pétrolière. Mémoire de diplôme d’ingénieur, Université Joseph Fourier - Grenoble I.
Guyoton F, J.R.Grasso, and P.Volant,Interrelation between induced seismic instabilities and complex geological structure Geophys. Res. Lett., 19, 705-708, 1992.
Roecker, S., Velocity structure  of the Pamir-Hindu Kush region: Possible evidence for subjected crust, JGR, 87, 945-959, 1982
Tarantola, A., & Valette, B. (1982). Generalized nonlinear inverse problems solved using the least squares criterion. Reviews of Geophysics, 20(2), 219-232.
double
Comments3charParticularities of the respective event (main shock, aftershock etc. ) and parameter estimation issues/uncertainties descriptiontext
M0_p222data vectorNmP-wave scalar momentdouble
M0_s222data vectorNmS-wave scalar momentdouble
M0_Dev_p222data vectorNmP-wave scalar moment errordouble
M0_Dev_s222data vectorNmS-wave scalar moment errordouble
M0_Dev222data vectorNmScalar moment errordouble
P222data vectorm^3Seismic potencydouble
X-relErr4data vectormRelative error in X coordintatedouble
Y-relErr4data vectormRelative error in Y coordintatedouble
Z-relErr4data vectormRelative error in Z coordintatedouble
ML-rel4data vectorRelative magnitude calibrated with local magnitude of larger event in the sequence'Magnitude'double
NCCP2data vectorNumber of cross-correlated P-wave datadouble
NCCS2data vectorNumber of cross-correlated S-wave datadouble
NCTP2data vectorNumber of catalog P-wave datadouble
NCTS2data vectorNumber of catalog S-wave datadouble
RMSCC13data vectorRMS residual for CC-derived differential arrivalsdouble
RMSCT13data vectorRMS residual for catolog-derived differential arrivalsdouble
CID2data vectorCluster IDdouble
usedStationLocation3data vectorName of stations from which data was used for localization computationtext
unusedStationLocation3data vectorName of stations from which data was not used for localization computationtext

The Numbers of Data type:

1 – the real data without limits,

2 – the integer data,

3 – text value,

4 – the real number rounded to 0.1 (shown as 11),

5 – time in Matlab format serial time – the time display format; seconds with accuracy 1/10,

6 – the real data display in an engineering manner with one decimal place, e.g.: 3.5E6, (obsolete, recommended 2cd)

7 – the real data display in an engineering manner with two decimal place, (obsolete, recommended 2cd)

bc – (b and c are code digits) the real data display in fix-point manner with at minimum b places before decimal and c decimal place

e.g. For number 3.149.

10:          „3”

11:          „3.1”

12:          „3.15”

20:          „03”

23:          „03.149”

1bc– the same manner as bc, but with place for a sign (space for sign „+”, sign - for sign „-”)

2cd– (c and d are code digits), the real data is displayed in an engineering manner, with place for a sign (space for sign „+”, sign ‘-’ for sign „-”), with c decimal place and exponent expressed by d places. The sign in exponent is always displayed.

e.g. For number 0.001:

211:       „1.0E-3”

221:       „1.00E-3”

212:       „1.0E-03”

222:       „1.00E-03”

e.g. For number 1000:

211:       „1.0E+3”

221:       „1.00E+3”

212:       „1.0E+03”

222:       „1.00E+03”

Examples catalogs

Bobrek Catalog

fieldtypevalunitdescriptionfieldType
ID3data vectorEvent ID
Time5data vectorEvent occurrence time
Lat25data vectordegLatitude
Long25data vectordegLongitude
Depth13data vectorkmHypocenter depth measured from the ground level
Elevation13data vectorkmHypocenter elevation measured over the see level
X10data vectormX coordinate
Y10data vectormY coordinate
Z10data vectormZ coordinate
ML4data vectorLocal magnitude'Magnitude'
E222data vectorJTotal seismic energy

LGCD Catalog

fieldtypevalunitdescriptionfieldType
ID3data vectorEvent ID
Time5data vectorEvent occurrence time
Lat24data vectordegLatitude
Long24data vectordegLongitude
Depth13data vectorkmHypocenter depth measured from the ground level
Elevation13data vectorkmHypocenter elevation measured over the see level
M0222data vectorNmScalar moment
Mw4data vectorMoment magnitude'Magnitude'
MTrr222data vectorNmFull solution: Moment tensor rr component (r - up)
MTss222data vectorNmFull solution: Moment tensor ss component (s - South)
MTee222data vectorNmFull solution: Moment tensor ee component (e - East)
MTrs222data vectorNmFull solution: Moment tensor rs component
MTre222data vectorNmFull solution: Moment tensor re component
MTse222data vectorNmFull solution: Moment tensor se component
MT_err222data vectorNmFull solution: Moment tensor error
ISO120data vector%Isotropic MT component
CLVD120data vector%CLVD component
DC20data vector%Double-Couple component
StrikeA30data vectordegStrike of nodal plane A
DipA20data vectordegDip of nodal plane A
RakeA130data vectordegRake of nodal plane A
SlopeA20data vectordegInclination for nodal plane A
StrikeB30data vectordegStrike of nodal plane B
DipB20data vectordegDip of nodal plane B
RakeB130data vectordegRake of nodal plane B
fp12data vectorHzP-wave corner frequency
fs12data vectorHzS-wave corner frequency Hz

Song Tranh Catalog

fieldtypevalunitdescriptionfieldType
ID3data vectorEvent ID
Time5data vectorEvent occurrence time
Lat24data vectordegLatitude
Long24data vectordegLongitude
Depth13data vectorkmHypocenter depth measured from the ground level
Elevation13data vectorkmHypocenter elevation measured over the see level
M0222data vectorNmScalar moment
Mw4data vectorMoment magnitude'Magnitude'
ML4data vectorLocal magnitude'Magnitude'
MTrr222data vectorNmFull solution: Moment tensor rr component (r - up)
MTss222data vectorNmFull solution: Moment tensor ss component (s - South)
MTee222data vectorNmFull solution: Moment tensor ee component (e - East)
MTrs222data vectorNmFull solution: Moment tensor rs component
MTre222data vectorNmFull solution: Moment tensor re component
MTse222data vectorNmFull solution: Moment tensor se component
MT_err222data vectorNmFull solution: Moment tensor error
ISO120data vector%Isotropic MT component
CLVD120data vector%CLVD component
DC20data vector%Double-Couple component
StrikeA30data vectordegStrike of nodal plane A
DipA20data vectordegDip of nodal plane A
RakeA130data vectordegRake of nodal plane A
SlopeA20data vectordegInclination for nodal plane A
StrikeB30data vectordegStrike of nodal plane B
DipB20data vectordegDip of nodal plane B
RakeB130data vectordegRake of nodal plane B
fp12data vectorHzP-wave corner frequency
fs12data vectorHzS-wave corner frequency Hz

[1] The values rounded to 0.1.