Inversion of Teleseismic Body Waves and Geodetic Data Including Seafloor Observation Data: Source Model with Enhanced Reliability by Using Objectively Determined Inversion Settings

joint inversion of teleseismic body waves and geodetic data using multiple time-window analysis. To enhance the reliability of the inversion solution, two important inversion-setting parameters that control the spatiotemporal slip distribution, the maximum slip duration of each subfault (Tsd) and the propagation velocity of the first time window (Vftw), are determined by objective criteria. By performing a semblance analysis of small-aperture seismic array data of near-source strong motion waveforms, Tsd is determined to be 100 s. Vftw is determined to be 2:0 km=s from the matching of the teleseismic and geodetic slip models. The spatial resolution of geodetic inversion is significantly improved by using seafloor crustal deformation data in the source area in addition to the terrestrial Global Positioning System network data. The estimated seismic moment and the maximum slip are 3:4 × 1022 N·m (Mw 9:0) and 43 m, respectively. The total rupture duration is ∼150 s. Large slip is seen in the shallow part of the fault surface, including the hypocenter, and its spatial extent along the trench axis is ∼300 km. This model is consistent with the results of bathymetric surveys and the slip models from tsunami data. The relationships among the fault parameters of the characterized source model of the 2011 Tohoku earthquake suggest the possibility that saturation of fault width occurs for this huge earthquake. Online Material: Text S1 denotes the explanation of how the rupture propagation in the dip direction affects the estimated value of maximum slip duration (Tsd). Table S1 denotes the detailed information of the teleseismic stations. Table S2 denotes the detailed information of the terrestrial geodetic stations