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

<p><strong>Directory &quot;waveform_data&quot;</strong>&nbsp;includes 13 tsunami time series data from different instruments: TM1, TM2,&nbsp;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&nbsp;amplitude (m).</p> <p><strong>Directory &ldquo;kin_models&rdquo;</strong> includes the following:</p> <p>1. Seafloor Mesh Geometry</p> <ul> <li>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 &ldquo;*Pt{1,2}-SM2.PointCoord.txt&rdquo; 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&nbsp;are calculated in UTM projection system, relative to an arbitrary reference point.</li> <li>The files &ldquo;*.ClipPath.txt&rdquo; includes the ID/lon/lat of mesh&nbsp;boundary nodes, which can be used for plotting.</li> <li>Visualization of the mesh parts are&nbsp;provided in PDF&nbsp;files.</li> <li>The file &ldquo;*Total-SM2.PointCoord.txt&rdquo; excludes boundary nodes and contains seafloor locations&nbsp;(504 nodes) that are directly used in&nbsp;tsunami arrival time calculations.</li> </ul> <p>2. Posterior Mean Models</p> <ul> <li>Ensemble-averaged models&nbsp;of seafloor displacements and uncertainty estimates, with&nbsp;no spatial averaging (&ldquo;0R&rdquo; in the file name) or with one-ring spatial averaging (&ldquo;1R&rdquo;). These models are shown in Figures 5 and 6 of&nbsp;<em>Jiang and Simons</em>&nbsp;(2016). The data files &ldquo;posterior_mean_{0,1}R.txt&rdquo; have three columns for (1) vertical seafloor displacement (m), (2) one-sigma standard deviation of displacement (m), and (3) corresponding resolution length (km). The&nbsp;model values&nbsp;(558 rows)&nbsp;correspond to nodes in&nbsp;files &ldquo;*Pt{1,2}-SM2.PointCoord.txt&rdquo; concatenated in sequential order.</li> <li>Tsunami arrival times&nbsp;(in sec) are calculated from the posterior mean values of propagation speeds, with zero sec&nbsp;at the earthquake epicenter. The&nbsp;coordinates (508&nbsp;nodes) are included in geometry file &ldquo;*Total-SM2.PointCoord.txt.&rdquo;</li> </ul> <p><strong>Directory &ldquo;kin_ensemble&rdquo;</strong> includes the entire posterior model ensemble (98304 samples) in HDF5 format. Using a Linux command <em>h5dump</em>&nbsp;will show the following information about the contained datasets, with their names and dimensions. The main datasets are: (1) Covariance (1008&times;1008); (2) Data Log-likelihood (98304&times;1); (3) Posterior Log-likelihood (98304&times;1); and (4) Sample Set (98304&times;1008). Each model has 1008 parameters (504 for displacement and 504 for propagation speeds). The source coordinates&nbsp;(504 nodes) are included in geometry file &quot;*Total-SM2-Parameter.PointCoord.txt.&quot;</p> <p><strong>Note:</strong>&nbsp;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).&nbsp;</p>