New Zealand Tsunami on Nov. 13, 2016
We have simulated the tsunami generated from the New Zealand earthquake (42.757°S, 173.077°E, depth=23.0 km, M = 7.8 at 11:02:56 UTC according to USGS) on November 13, 2016 (Fig. 1). The assumed tsunami source is located within the aftershock area (Fig. 1). The fault length and width are 150 km × 40 km. The focal mechanism is strike:219º, dip:38º, slip:128º from the USGS's W-phase moment tensor solution. The top depth of the fault was assumed to 5 km. The average slip on the fault is 2 m. The seismic moment is 6.0 x 10**29 Nm (Mw = 7.8) assuming the rigidity of 5 x 10**10 N/m**2. As the initial condition for tsunami, static deformation of the seafloor is calculated for a rectangular fault model [Okada, 1985] using the source model. The used bathymetry data is the 30 arc-second grid data from GEBCO for coastal tide gauges. For the near-field computation we used 12 arc-second grid data resampled from the GEBCO data. To calculate tsunami propagation, the linear shallow-water, or long-wave, equations were numerically solved by using a finite-difference method [Satake, 1995]. We have downloaded the Tide gauge data from UNESCO/IOC's web site, and compared the simulated tsunami waveforms and the observed ones (Fig. 2). We can see the tsunami propagation in the animation (Fig. 3).
Fig.1 Tsunami Source Model
The red lines indicate uplift with the contour interval of 0.1 m, while the blue dotted lines indicate subsidence with the contour interval of 0.1 m. Aftershocks determined by USGS are shown by red circles. The focal mechanism determined by USGS is also shown.
Fig.2 Maximum Height of Simulated Tsunami and Tsunami Waveforms
Solid lines in red and blue indicate the observed tsunami waveforms and synthtic ones, respectively.
Fig.3 Tsunami Propagation (Click to start animation)
The red color means that the water surface is higher than normal sea level, while the blue means lower.
by Yushiro Fujii (IISEE, BRI) and Kenji Satake (ERI, Univ. of Tokyo)
References Okada, Y. (1985), Surface Deformation Due to Shear and Tensile Faults in a Half-Space, Bull. Seismol. Soc. Am., 75, 1135-1154.
Satake, K. (1995), Linear and Nonlinear Computations of the 1992 Nicaragua Earthquake Tsunami, Pure and Appl. Geophys., 144, 455-470.
Last Updated on 2016/11/18