Data and Models for "Probabilistic Imaging of Tsunamigenic Seafloor Deformation During the 2011 Tohoku-oki Earthquake"

Directory "waveform_data" includes 13 tsunami time series data from different instruments: TM1, TM2, KPG1, KPG2, GB801, GB802, GB803, GB804, GB806, GB807, D21401, D21413, and D21418. Each data file (*.dat) has two columns for (1) the time since earthquake initiation (min) and (2) ocean surface or seafloor displacement amplitude (m). Directory “kin_models” includes the following: 1. Seafloor Mesh Geometry The entire seafloor mesh consists of two separate parts (422 and 136 nodes each; 558 in total) due to the need to resolve potential discontinuity at the trench. The files “*Pt{1,2}-SM2.PointCoord.txt” consists of six columns for the ID, longitude (deg), latitude (deg), East (km), North (km) and depth (km) of the nodes in each triangular mesh. The E/N coordinates are calculated in UTM projection system, relative to an arbitrary reference point. The files “*.ClipPath.txt” includes the ID/lon/lat of mesh boundary nodes, which can be used for plotting. Visualization of the mesh parts are provided in PDF files. The file “*Total-SM2.PointCoord.txt” excludes boundary nodes and contains seafloor locations (504 nodes) that are directly used in tsunami arrival time calculations. 2. Posterior Mean Models Ensemble-averaged models of seafloor displacements and uncertainty estimates, with no spatial averaging (“0R” in the file name) or with one-ring spatial averaging (“1R”). These models are shown in Figures 5 and 6 of Jiang and Simons (2016). The data files “posterior_mean_{0,1}R.txt” have three columns for (1) vertical seafloor displacement (m), (2) one-sigma standard deviation of displacement (m), and (3) corresponding resolution length (km). The model values (558 rows) correspond to nodes in files “*Pt{1,2}-SM2.PointCoord.txt” concatenated in sequential order. Tsunami arrival times (in sec) are calculated from the posterior mean values of propagation speeds, with zero sec at the earthquake epicenter. The coordinates (508 nodes) are included in geometry file “*Total-SM2.PointCoord.txt.” Directory “kin_ensemble” includes the entire posterior model ensemble (98304 samples) in HDF5 format. Using a Linux command h5dump will show the following information about the contained datasets, with their names and dimensions. The main datasets are: (1) Covariance (1008×1008); (2) Data Log-likelihood (98304×1); (3) Posterior Log-likelihood (98304×1); and (4) Sample Set (98304×1008). Each model has 1008 parameters (504 for displacement and 504 for propagation speeds). The source coordinates (504 nodes) are included in geometry file "*Total-SM2-Parameter.PointCoord.txt." Note: three different geometries files above are used for (1) posterior mean displacements (558 nodes), (2) arrival time calculation (508 nodes), and (3) source inversion models (504 nodes)."waveform_data"目录包含13组来自不同观测仪器的海啸时间序列数据，仪器标识分别为TM1、TM2、KPG1、KPG2、GB801、GB802、GB803、GB804、GB806、GB807、D21401、D21413以及D21418。每个数据文件（后缀为*.dat）包含两列数据：(1) 地震起始后的时长（单位：分钟）；(2) 海面或海底位移振幅（单位：米）。 "kin_models"目录包含以下内容： 1. 海底网格几何（Seafloor Mesh Geometry） 整个海底网格由两个独立部分组成，各自包含422和136个节点，总计558个节点，以解决海沟处可能存在的不连续性问题。文件"*Pt{1,2}-SM2.PointCoord.txt"包含六列数据，分别对应每个三角网格节点的编号、经度（单位：度）、纬度（单位：度）、东向坐标（单位：千米）、北向坐标（单位：千米）以及深度（单位：千米）。其中东/北向坐标采用通用横轴墨卡托（UTM）投影坐标系，相对于任意参考点计算得到。 文件"*.ClipPath.txt"包含网格边界节点的编号、经度与纬度，可用于绘图可视化。网格各部分的可视化结果已保存为PDF文件。文件"*Total-SM2.PointCoord.txt"不含边界节点，包含直接用于海啸到达时间计算的海底位置信息（共504个节点）。 2. 后验均值模型（Posterior Mean Models） 包含海底位移的集合平均模型与不确定性估计结果，分为未进行空间平均（文件名中标注"0R"）与采用单环空间平均（文件名中标注"1R"）两种类型。此类模型对应Jiang与Simons（2016）论文中的图5和图6。数据文件"posterior_mean_{0,1}R.txt"包含三列数据：(1) 海底垂直位移（单位：米）；(2) 位移的一倍标准差（单位：米）；(3) 对应的分辨率长度（单位：千米）。模型数据共558行，与文件"*Pt{1,2}-SM2.PointCoord.txt"中按顺序拼接的节点一一对应。 基于传播速度的后验均值计算得到海啸到达时间（单位：秒），以地震震中处的到达时间为0秒。用于计算的坐标（共508个节点）已包含在几何文件"*Total-SM2.PointCoord.txt"中。 "kin_ensemble"目录包含以HDF5格式存储的完整后验模型集合（共98304个样本）。在Linux系统中可通过`h5dump`命令查看其中包含的数据集信息，包括各数据集的名称与维度。主要数据集包括：(1) 协方差矩阵（1008×1008）；(2) 数据对数似然（98304×1）；(3) 后验对数似然（98304×1）；(4) 样本集（98304×1008）。每个模型包含1008个参数，其中504个对应海底位移，504个对应传播速度。源反演模型的坐标（共504个节点）已包含在几何文件"*Total-SM2-Parameter.PointCoord.txt"中。 注意：上述三类不同的几何文件分别用于：(1) 后验均值位移计算（558个节点）；(2) 到达时间计算（508个节点）；(3) 源反演模型构建（504个节点）。