GCMT提供目录的打包下载。下载地址是 http://www.globalcmt.org/CMTfiles.html

下载下来的ndk文件其实是一个文本文件。每一个地震有5行的内容:

PDE  2005/01/01 01:20:05.4  13.78  -88.78 193.1 5.0 0.0 EL SALVADOR             
C200501010120A   B:  4    4  40 S: 27   33  50 M:  0    0   0 CMT: 1 TRIHD:  0.6
CENTROID:     -0.3 0.9  13.76 0.06  -89.08 0.09 162.8 12.5 FREE S-20050322125201
23  0.838 0.201 -0.005 0.231 -0.833 0.270  1.050 0.121 -0.369 0.161  0.044 0.240
V10   1.581 56  12  -0.537 23 140  -1.044 24 241   1.312   9 29  142 133 72   66

首先要说明的是,本质上ndk文件是用字符数来决定表达的含义,而不是用空格分割来决定的。 不过人还是习惯用空格分隔来阅读。

第一行:

PDE      2005/01/01 01:20:05.4  13.78  -88.78 193.1 5.0    0.0    EL SALVADOR
参考目录  发震日期     发震时刻     纬度    经度    深度   mb震级 MS震级  地理位置

[1-4] 第 1 个字符到第 4 个字符表示参考目录名称,(比如 PDE 表示 USGS,ISC 表示 ISC,SWE表示 surface-wave location,[Ekstrom, BSSA, 2006])

[6-15] 发震日期

[17-26] 发震时刻

[28-33] 纬度

[35-41] 经度

[43-47] 深度

[49-55] 震级,一般是 mb 震级和 MS 震级

[57-80] 地理位置(24 个字符)

第二行:

C200501010120A   B:  4    4  40 S: 27   33  50 M:  0    0   0 CMT: 1 TRIHD:  0.6
CMT 事件编号       反演所用的数据                                 反演类型 震源速度函数

[1-16] CMT 事件编号。老格式的编号是 8 位,现在是 14 位。老格式为 XMMDDYYZ,其中 MMDDYY 是日期,Z 是用于区分同一天的事件的字符。X 是指反演用数据的类型。 新格式:XYYYYMMDDhhmmZ,第一个字符 X 只可能是:B - 仅体波反演,S - 仅面波反演,M - 仅地幔波(mantle waves)反演,C - 组合。

[18-61] 反演用的数据。三个数据类型可能是长周期体波反演(B),中等周期的面波反演(S),长周期的地幔波(mantle waves)反演。 对于每一个数据类型,会给出三个值:台站数量、分量数量和使用的最短周期。

[63-68] 反演类型。0:矩张量,1:带零迹的矩张量,2:双力偶。通常我们说矩张量就是包含了零迹的,我猜测也许这里说的矩张量是 DC+CLVD

[70-80] 震源速度函数(moment-rate function)的类型和长度。TRIHD 是三角形的,BOXHD 是矩形的。所给的值是真实的一半。 长度正比于标量矩的立方根。如果标量矩是 10^24,长度是 1.05 秒。如果标量矩是 10^27,长度是 10.5 秒。

第三行:

CENTROID:     -0.3 0.9  13.76 0.06  -89.08 0.09 162.8 12.5 FREE    S-20050322125201
反演中的质心参数                                              深度类型 时间戳

[1-58] 反演中的质心参数。和地震参考目录相比的参数。依次为时间、纬度、经度、深度和他们的估计标准差。如果标准差是 0,则说明直接用的参考目录。

[60-63] 深度类型。FREE 表示为反演的结果。FIX 表示经过过修正。BDY 表示经过了宽频带 P 波的修正。

[65-80] 时间戳. 16 个字符。Q- 表示事件后几小时内计算的。S- 表示几小时之后计算的。

第四行:

23          0.838 0.201 -0.005 0.231 -0.833 0.270  1.050 0.121 -0.369 0.161  0.044 0.240
标量矩的系数  USE 坐标系下的矩张量元素和他们的估计标准差

[1-2] 标量矩的系数。如果是 24,矩张量的元素如果单位是 dyne-cm, 那么就需要乘以 10^24。

[3-80] 矩张量的元素,USE 坐标系: Muu Mss Mee Mus Mue Mse

第五行:

V10   1.581 56  12  -0.537 23 140  -1.044 24 241   1.312   9 29  142 133 72   66
版本号 主轴坐标系的矩张量                              标量矩   断层面解

[1-3] 版本号。

[4-48] 主轴坐标系的矩张量:三个特征向量的特征值、俯角和方位角。特征值需要乘以标量矩系数(见第四行)。

[50-56] 标量矩。还需要乘以标量矩系数(见第四行)。

[58-80] Aki and Richards 描述的两个节面的 Strike、 dip 和 rake