Landslide Tsunami Hazard Along the Upper US East Coast: Effects of Slide Deformation, Bottom Friction, and Frequency Dispersion

Abstract

Numerical simulations of Submarine Mass Failures (SMFs) are performed along the upper US East Coast to assess the effect of slide deformation on predicted tsunami hazard. Tsunami generation is simulated using the three-dimensional non-hydrostatic model NHWAVE. For rigid slumps, the geometry and law of motion are specified as bottom boundary conditions. Deforming slide motion is modeled using a depth-integrated bottom layer of dense Newtonian fluid, fully coupled to the overlying fluid motion. Once the SMFs are no-longer tsunamigenic, tsunami propagation simulations are performed using the Boussinesq wave model FUNWAVE-TVD, using nested grids of increasingly fine resolution towards shore and employing a one-way coupling methodology. Probable maximum tsunamis are simulated for Currituck SMF proxies sited in four areas of the shelf break slope that have enough sediment accumulation to cause large failures. Deforming slides have a slightly larger initial acceleration, but still generate a smaller tsunami than rigid slumps due to their spreading and thinning out during motion, which gradually makes them less tsunamigenic. Comparing the maximum envelope of surface elevations along a 5 m isobath, consistent with earlier work, the bathymetry of the wide shelf is found to strongly control the spatial distribution of tsunami inundation. Overall, tsunamis caused by rigid slumps are worst case scenarios, providing up to 50% more inundation than for deforming slides having a moderate level of viscosity set in the upper range of debris flows. Tsunamis from both types of SMFs are able to cause water withdrawal to the 5 m isobath or deeper. Bottom friction effects are assessed by performing some of the simulations using two different Manning coefficients, one 50% larger than the other. With increased bottom friction, the largest tsunami inundations at the coast are reduced by up to 15%. Selected simulations are rerun by turning off dispersion in the model, which leads to moderate changes in maximum surface elevations nearshore (− 10 to + 5% changes), but to more significant effects in the far field (− 40 to 80% changes). Onshore, dispersion causes the appearance of short period undular bores that eventually break nearshore without significantly affecting inundation at the coast. However, these bores increase wave-induced maximum flow velocity and impulse forces, the latter by up to 40%, which may affect the design of coastal structures.