Modeling coastal tsunami hazard from submarine mass failures: effect of slide rheology,experimental validation,and case studies off the US East Coast

Natural Hazards - We perform numerical simulations to assess how coastal tsunami hazard from submarine mass failures (SMFs) is affected by slide kinematics and rheology. Two types of two-layer SMF...     

Tsunami detection by high-frequency radar in British Columbia: performance assessment of the time-correlation algorithm for synthetic and real events

The authors recently proposed a new method for detecting tsunamis using high-frequency (HF) radar observations, referred to as “time-correlation algorithm” (TCA; Grilli et al. Pure Appl Geophys 173(12):3895–3934, 2016a, 174(1): 3003–3028, 2017). Unlike standard algorithms that detect surface current patterns, the TCA is based on analyzing space-time correlations of radar signal time series in pairs of radar cells, which does not require inverting radial surface currents. This was done by calculating a contrast function, which quantifies the change in pattern of the mean correlation between pairs of neighboring cells upon tsunami arrival, with respect to a reference correlation computed in the recent past. In earlier work, the TCA was successfully validated based on realistic numerical simulations of both the radar signal and tsunami wave trains. Here, this algorithm is adapted to apply to actual data from a HF radar installed in Tofino, BC, for three test cases: (1) a simulated far-field tsunami generated in the Semidi Subduction Zone in the Aleutian Arc; (2) a simulated near-field tsunami from a submarine mass failure on the continental slope off of Tofino; and (3) an event believed to be a meteotsunami, which occurred on October 14th, 2016, off of the Pacific West Coast and was measured by the radar. In the first two cases, the synthetic tsunami signal is superimposed onto the radar signal by way of a current memory term; in the third case, the tsunami signature is present within the radar data. In light of these test cases, we develop a detection methodology based on the TCA, using a correlation contrast function, and show that in all three cases the algorithm is able to trigger a timely early warning.     

Numerical Simulation of the 2011 Tohoku Tsunami Based on a New Transient FEM Co-seismic Source: Comparison to Far- and Near-Field Observations

In this work, we simulate the 2011 M9 Tohoku-Oki tsunami using new coseismic tsunami sources based on inverting onshore and offshore geodetic data, using 3D Finite Element Models (FEM). Such FEMs simulate elastic dislocations along the plate boundary interface separating the stiff subducting Pacific Plate from the 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 (several tens of meters) along the updip portion of the rupture near the trench. To assess the accuracy of the new approach, we compare observations and numerical simulations of the tsunami's far- and near-field coastal impact for: (i) one of the standard seismic inversion sources (UCSB; Shao et al. 2011); and (ii) the new FEM sources. Specifically, results of numerical simulations for both sources, performed using the fully nonlinear and dispersive Boussinesq wave model FUNWAVE-TVD, are compared to DART buoy, GPS tide gauge, and inundation/runup measurements. We use a series of nested model grids with varying resolution (down to 250 m nearshore) and size, and assess effects on model results of the latter and of model physics (such as when including dispersion or not). We also assess the effects of triggering the tsunami sources in the propagation model: (i) either at once as a hot start, or with the spatiotemporal sequence derived from seismic inversion; and (ii) as a specified surface elevation or as a more realistic time and space-varying bottom boundary condition (in the latter case, we compute the initial tsunami generation up to 300 s using the non-hydrostatic model NHWAVE). Although additional refinements are expected in the near future, results based on the current FEM sources better explain long wave near-field observations at DART and GPS buoys near Japan, and measured tsunami inundation, while they simulate observations at distant DART buoys as well or better than the UCSB source. None of the sources, however, are able to explain the largest runup and inundation measured between 39.5° and 40.25°N, which could be due to insufficient model resolution in this region (Sanriku/Ria) of complex bathymetry/topography, and/or to additional tsunami generation mechanisms not represented in the coseismic sources (e.g., splay faults, submarine mass failure). This will be the object of future work.     

Progress on Nonlinear-Wave-Forced Sediment Transport Simulation

In this paper, we report on the use of a numerical wave tank (NWT), based on fully nonlinear potential flow (FNPF) equations, in driving simulations of flow and sediment transport around partially buried obstacles. The suspended sediment transport is modeled in the near-field in a Navier-Stokes (NS) model using an immersed-boundary method and an attached sediment transport simulation module. Turbulence is represented by large eddy simulation (LES). The NWT is based on a higher order boundary element method (BEM), with an explicit second-order time stepping. Hence, only the NWT boundary is discretized. The solution for the velocity potential and its derivatives along the boundary is obtained in the BEM, which subsequently provides a solution at any required internal point within the domain. At initial time, the NS-LES model domain is initialized with the 3-D velocity field provided by the NWT and driven for later time by the pressure gradient field obtained in the NWT. Incident wave fields, as specified in the NWT to drive sediment transport, can be arbitrary. Applications are presented here for single frequency waves, such as produced by a harmonic piston wavemaker in the laboratory, and modulated frequency wave groups. The feasibility of coupling the irrotational flow and NS solutions is demonstrated.