The March 11, 2011 M9 Tohoku-Oki Earthquake, which is believed to be the largest event recorded in Japanese history, created a major tsunami that caused numerous deaths and enormous destruction on the nearby Honshu coast. Various tsunami sources were developed for this event, based on inverting seismic or GPS data, often using very simple underlying fault models (e.g., Okada, 1985). Tsunami simulations with such sources can predict deep water and far-field observations quite well, but coastal impact is not as well predicted, being over- or under-estimated at many locations. In this work, we developed a new tsunami source, similarly based on inverting onshore and offshore geodetic (GPS) data, but using 3D Finite Element Models (FEM) that simulate elastic dislocations along the plate boundary interface separating the stiff subducting Pacific Plate, and relatively weak forearc and volcanic arc of the overriding Eurasian plate. Due in part to the simulated weak forearc materials, such sources produce significant shallow slip along the updip portion of the rupture near the trench (several tens of meters). We assess the accuracy of the new approach by comparing numerical simulations to observations of the tsunami far- and near-field coastal impact using:
one of the standard seismic inversion sources, which we found provided the best prediction of tsunami near-field impact in our model (UCSB; Shao et al., 2011); and
the new FEM source. Specifically, we compare numerical results to DART buoy, GPS tide gage, and inundation/runup measurements. Numerical simulations are performed using the fully nonlinear and dispersive Boussinesq wave model FUNWAVE-TVD, which is parallelized and available in Cartesian or spherical coordinates. We use a series of nested model grids, with varying resolution (down to 250 m nearshore) and size, and assess effects on results of the latter and of model physics (such as when including dispersion or not).